The spaceborne Global Navigation Satellite System-Reflectometry (GNSS-R) delay-Doppler map (DDM) data collected over ocean carry typical feature information about the ocean surface, which may be covered by open water, mixed water/ice, complete ice, etc. A new method based on Doppler spread analysis is proposed to remotely sense sea ice using the spaceborne GNSS-R data collected over the Northern and Southern Hemispheres. In order to extract useful information from DDM, three delay waveforms are defined and utilized. The delay waveform without Doppler shift is defined as central delay waveform (CDW), while the integration of delay waveforms of 20 different Doppler shift values is defined as integrated delay waveform (IDW). The differential waveform between normalized CDW (NCDW) and normalized IDW (NIDW) is defined as differential delay waveform (DDW), which is a new observable used to describe the difference between NCDW and NIDW, which have different Doppler spread characteristics. The difference is mainly caused by the roughness of reflected surface. First, a new data quality control method is proposed based on the standard deviation and root-mean-square error (RMSE) of the first 48 bins of DDW. Then, several different observables derived from NCDW, NIDW, and DDW are applied to distinguish sea ice from water based on their probability density function. Through validating against sea ice edge data from the Ocean and Sea Ice Satellite Application Facility, the trailing edge waveform summation of DDW achieves the best results, and its probabilities of successful detection are 98.22% and 96.65%, respectively, in the Northern and Southern Hemispheres.

The distribution of forest tree species provides crucial data for regional forest management and ecological research. Although medium-high spatial resolution remote sensing images are widely used for dynamic monitoring of forest vegetation phenology and species identification, the use of multiresolution images for similar applications remains highly uncertain. Moreover, it is necessary to explore to what extent spectral variation is responsible for the discrepancies in the estimation of forest phenology and classification of various tree species when using up-scaled images. To clarify this situation, we studied the forest area in Harqin Banner in northeast China by using year-round multiple-resolution time-series images (at four spatial resolutions: 4, 10, 16, and 30 m) and eight phenological metrics of four deciduous forest tree species in 2018, to explore potential impacts of relevant results caused by various resolutions. We also investigated the effect of using up-scaled time-series images by comparing the corresponding results that use pixel-aggregation algorithms with the four spatial resolutions. The results indicate that both phenology and classification accuracy of the dominant forest tree species are markedly affected by the spatial resolution of time-series remote sensing data (p < 0.05): the spring phenology of four deciduous forest tree species first rises and then falls as the image resolution varies from 4 to 30 m; similarly, the accuracy of tree species classification increases as the image resolution varies from 4 to 10 m, and then decreases as the image resolution gradually falls to 30 m (p < 0.05). Therefore, there remains a profound discrepancy between the results obtained by up-scaled and actual remote sensing data at the given spatial resolutions (p < 0.05). The results also suggest that combining phenological metrics and time-series NDVI data can be applied to identify the regional dominant tree species across different spatial resolutions, which would help advance the use of multiscale time-series satellite data for forest resource management.

Lambertian surfaces represent an important assumption when constructing thermal radiance transfer equations for remote sensing observations of ground surface temperatures. We identify the properties of ground surfaces in thermal infrared regions as Lambertian surfaces via field experiments. Because Lambertian surfaces present homogeneous thermal emissions levels in hemispheric directions for a specific ground surface under specific kinetic temperatures and emissions, we conducted a series of field experiments to illustrate the properties of such ground surfaces. Four typical ground surfaces were selected for the experiments to observe thermal emissions: bare soil, grass, water, and concrete. Radiance thermometers were used to observe ground emissions from seven directions: 30°, 45°, 60°, 90°, 120°, 135°, and 150°. Solar zenith angles were considered for the observation of ground emissions. Experiments were conducted in five different regions of China (Beijing, Nanjing, Xilinguole, Yongzhou, and Jiangmen) during both daytime and nighttime. To determine whether different observation angles have significantly different effects on radiance, statistical analyses (ANOVA and Friedman test) were conducted. Post hoc multiple comparison tests and pairwise multiple comparisons were also conducted to examine the various pairings of observation angles and to measure the radiance differences. Roughly half of the radiance groups of all observed sites were tested via an ANOVA, and the remaining groups with unequal variances were subjected to the Friedman test. The results indicate that statistically significant differences in the radiance levels occurred among the seven angles for almost all of the sites (39 of the 40 groups). The results of our experiments indicate that the selected ground surfaces, especially the grass and the bare soil, may not behave with Lambertian properties in the thermal infrared region. This is probably attributed to the roughness of the selected surface, because we found that roughness is an important factor affecting the observed magnitude of thermal emission from different directions of the ground surface under study. Therefore, whether or not a terrestrial surface can be assumed to be a Lambertian surface should be based on their geometric structure. When the surface is relatively smooth, we can say that it is close to the Lambertian property in thermal emission.

The influence of tree distribution and slope on canopy gap fraction (GF) and bidirectional reflectance factor (BRF) is shown here to be non-negligible. Trees are often assumed to be randomly distributed in natural forests due to random distribution of natural resources, but this assumption is not valid for forest plantations. A geometric optical model for forest plantations (GOFP) is a geometric optical model for forest plantations on horizontal surfaces based on the theory of exclusion distance among crowns. Sloping terrains change the exclusion distance among crowns, and inevitably affect the canopy GF and BRF. In this article, GOFP with a hypergeometric model (distances among trees are considered) on horizontal surfaces is modified as GOFP-T to simulate BRF for forest plantations on sloping terrains under two scenarios: (horizontal distances among crowns remain unchanged with slope) and (sloping distances among crowns remain unchanged with slope). Two three-dimensional (3-D) radiative transfer models (DART and LESS) and field measurements are used to evaluate and validate GOFP-T simulations. The results show that 1) the canopy GF, four component area ratios, and canopy BRF simulated by GOFP-T show high consistency with results from the two 3-D model: root-mean-square errors in GF, sunlit foliage, and sunlit ground are less than 0.02, 0.06, and 0.03, respectively; 2) forest coverage and canopy reflectance in GOFP-T are compared well with point cloud results from an airboard LiDAR system and Landsat 8 OLI surface reflectance products, respectively, indicating that GOFP-T has ability in simulating canopy reflectance for forest plantations on sloping terrains. GOFP-T with the hypergeometric model in this article is the first simple model to simulate canopy GF and BRF for forest plantations on sloping terrains.

As a widely used ground-based optical instrument, the LAI-2000 or LAI-2200 plant canopy analyzer (PCA) (Li-Cor, Inc., Lincoln, NE) is designed to measure the plant effective leaf area index (L_e) by measuring the canopy gap fraction at several limited or discrete view zenith angles (VZAs) (usually five VZAs: 7, 23, 38, 53, and 68°) based on Miller’s equation. Miller’s equation requires the probability of radiative transmission through the canopy to be measured over the hemisphere, i.e., VZAs in the range from 0 to 90°. However, the PCA view angle ranges are confined to several limited ranges or discrete sectors. The magnitude of the error produced by the discretization of VZAs in the leaf area index measurements remains difficult to determine. In this study, a theoretical deduction was first presented to definitely prove why the limited or discrete VZAs or ranges can affect the L_e measured with the PCA, and the specific error caused by the limited or discrete VZAs was described quantitatively. The results show that: (1) the weight coefficient of the last PCA ring is the main cause of the error; (2) the error is closely related to the leaf inclination angles (IAs)—the L_e measured with the PCA can be significantly overestimated for canopies with planophile IAs, whereas it can be underestimated for erectophile IAs; and (3) the error can be enhanced with the increment of the discrete degree of PCA rings or VZAs, such as using four or three PCA rings. Two corrections for the error are presented and validated in three crop canopies. Interestingly, although the leaf IA type cannot influence the L_e calculated by Miller’s equation in the hemispheric space, it affects the L_e measured with the PCA using the discrete form of Miller’s equation for several discrete VZAs.

Sloping terrain distorts the sun-target–sensor geometry, resulting in biases of the optical reflectance measured by remote sensors relative to flat situations. Performing topographic correction (TC) is, therefore, deemed mandatory to foster the full exploitation of satellite images worldwide to support various applications in mountainous regions. Various TC models have already been proposed and developed, while most of them were previously evaluated at local or regional scales using a few images with various evaluation criteria. Therefore, a systematic and comprehensive assessment has yet to be done on these TC models in the global mountainous regions. In the present study, 10 523 Landsat 8 Operational Land Imager (OLI) images, filtered by land cover types and seasons and sampled across global mountain regions, were corrected using ten popular TC models—statistical empirical (SE), b correction, variable empirical coefficient algorithm (VECA), C-correction (CC), sun-canopy–sensor (SCS), DymondShepherd (DS), SCS + C, path length correction (PLC), Minnaert, and Minnaert + SCS—with a unified evaluation criterion on the Google Earth Engine (GEE) platform. The outcomes are that: 1) global TC effects on Landsat 8 OLI images generally increase with sun zenith angles (SZAs) and latitudes; 2) six models (SE, b correction, CC, VECA, Minnaert, and Minnaert + SCS) show good adaptability among the ten models for the global mountainous placing a disregard to land cover types and seasons; and 3) considering permanent snow and ice, needle-leaved forests in winter, and null values might appear in b correction, SE is deemed to be with the most global adaptability. This study pioneers an evaluation of fashionable TC models concerning mountainous regions worldwide and will be useful for applying TC to Landsat images for the benefit of making global TC products in the future and a fair intercomparison of OLI surface reflectance measured in various mountainous areas of the globe.