Search

Abstract Long-term streamflow observations contain essential information for understanding hydrological changes and managing water resources. A continental-scale dataset or analysis of temporal streamflow change is still lacking across hydrologic gauges in the Conterminous United States (CONUS). Here, we compiled 70 years of streamflow records from 1951 to 2021 at ~ 8000 hydrologic stations across the CONUS and characterized temporal trends, regime shifts, and extreme events using a Bayesian time series analysis algorithm. We found that the occurrences of sudden streamflow changes (e.g., regime shifts and extreme events) have been increasing with time across the CONUS. In addition, we derived 181 streamflow indicators that are valuable for hydrological and biological applications, such as the duration and frequency of high or low streamflow events. The Mississippi River Basin, especially the middle and lower parts, was a hot spot of high-frequency high-flow events. Overall, we anticipate the dataset generated here offers valuable information for understanding and quantifying changes in water resources across the CONUS.

Abstract Riverine nitrogen is a pivotal determinant influencing water quality in inland and coastal waters. Despite the recognized utility, no spatially-explicit data on riverine nitrogen yield is available for large parts of the world, thus hindering our ability to identify the contributors to riverine nitrogen and understand aquatic nitrogen cycling. To fill the data gap for the United States, here we (1) compiled 294,996 total nitrogen (TN), 225,827 nitrate (NO3 −), 204,015 ammonium (NH4 +), and 158,837 total organic nitrogen (TON) concentrations, with concurrent streamflow data, across the Conterminous United States (CONUS), (2) estimated riverine nitrogen loads for over 1,800 hydrological stations, (3) derived the spatial distribution of annual riverine nitrogen yield by leveraging river and catchment connectivity information contained in the National Hydrography Dataset plus (NHDPlus), and (4) characterized nonpoint-source TN loads by excluding point-source loads. This new spatial dataset quantifies spatial sources of nitrogen yield from point and non-point sources (e.g., up to 36% from point sources across the U.S.) and serves as ground-truthing to validate water quality models.

Abstract The time it takes for an ecosystem to recover from a disturbance is a key to environmental management. Conventionally, recovery is defined as a return to the pre‐disturbance state, assuming ecosystem stationarity. However, this view does not account for the impact of external forces like climate change, imposing non‐stationarity and trends. Alternatively, the counterfactual approach views recovery as the state the ecosystem would have achieved if the disturbance had not occurred, accounting for external forces. Here, we present a simple method to estimate the counterfactual recovery time. By implementing our method to the greening of the Arctic region, we showed that counterfactual greening recovery is twice as long as conventional recovery over the region. We argue that the well‐documented greening of the region acts as an external force, leading to such a large difference. We advocate for embracing the counterfactual definition of recovery, as it aligns with realistic decision‐making processes.

Harmful algal blooms (HABs) threaten lake ecosystems and public health. Early HAB detection is possible by monitoring chlorophyll-a (Chl-a) concentration. Ground-based Chl-a data have limited spatial and temporal coverage but can be geo-registered with temporally coincident satellite imagery to calibrate a remote sensing-based predictive model for regional mapping over time. When matching ground and satellite data, positional and temporal discrepancies are unavoidable due particularly to dynamic lake surfaces, thereby biasing the model calibration. This limitation has long been recognized but so far has not been addressed explicitly. To mitigate such effects of data mismatching, we proposed an Akaike Information Criterion (AIC)-like weighted regression algorithm that relies on an error-based heuristic to automatically favor “good” data points and downplay “bad” points. We evaluated the algorithm for estimating Chl-a over inland lakes in Ohio using Harmonized Landsat Sentinel-2. The AIC-like weighted regression estimates showed superior performance with an R2 of 0.91 and an error variance (σE2) of 0.29 μg/L, outperforming linear regression (R2 = 0.34, σE2 = 2.34 μg/L) and random forest (R2 = 0.82, σE2 = 0.92 μg/L). We also noticed the poorest performance occurred in the spring due to low reflectance variation in clear water and low Chl-a concentration. Our weighted regression scheme is adaptive and generically applicable. Future studies may adopt our scheme to tackle other remote sensing estimation problems (e.g., terrestrial applications) for alleviating the adverse effects of geolocation errors and temporal discrepancies.

Cyanobacterial harmful algal blooms release toxins and form thick blanket layers on the water surface causing widespread problems, including serious threats to human health, water ecosystem, economics, and recreation. To identify the potential drivers for the bloom, there is a need for extensive observations of the water sources with bloom occurrences. However, the traditional methods for monitoring water sources, such as collection of point ground samples, have proven limited due to spatial and temporal variability of water resources, and the cost associated with collecting samples that accurately represent this variability. These limitations can be addressed through the use of high-frequency satellite data. In this study, we explored the use of Random Forest (RF), which is one of the widely used machine learning architectures, to evaluate the performance of Sentinel-3 OLCI (Ocean and Land Color Imager) images in predicting bloom proxies in the western region of Lake Erie. The sixteen available bands of Sentinel-3 images were used as the predictor variables, while four proxies of the cyanobacterial masses, including Chlorophyll-a, Microcystin, Phycocyanin, and Secchi-depth, were considered as response variables in the RF models, with one RF model per proxy. Each of the proxies comes with a unique set of traits that can help with bloom detection. Among four RF models, the model for Chlorophyll-a performed the best with R2 = 0.55 and RMSE = 20.84 µg/L, while R2 performance for the rest of the other proxies was less than 0.5. This is because Chlorophyll-a is the most dominant and optically active pigment in water, while Phycocyanin, which is a strong indicator of harmful bloom, is present in low concentrations. Additionally, Microcystin, responsible for bloom toxicity, has limited spectral sensitivity, and Secchi-depth could be influenced by various factors besides blooms, such as colored dissolved organic and inorganic matter. On further examining the relationship between the proxies, Microcystin and Secchi-depth were significantly correlated with Chlorophyll-a, which enhances the usefulness of Chlorophyll-a in accurately identifying the presence of algal blooms.

Multi-angular remote sensing observation contains crucial information on forest structure parameters. Here, our goal is to examine the ability of multi-angular indices, which are constructed by the typical-angular reflectances in red and NIR bands from MODIS observations, for the retrieval of forest biomass based on the field-measured above-ground biomass (AGB) data. Specifically, we employed the updated version of the MCD43A1 BRDF parameter product as an input for BRDF models to reconstruct the MODIS typical-angular reflectances. Furthermore, we evaluated the effects of different configurations of BRDF models and solar zenith angles (SZA) on forest AGB estimation using our developed multi-angular indices. The semivariogram analysis strategy combined with Landsat ground-surface reflectance data was employed to determine the MODIS pixel heterogeneity; the survey data from field sites of homogeneous pixels was used in our analysis and validation. The results show that our developed multi-angular indices based on a hot-revised BRDF model, under a SZA of 45°, when combined with forest cover information, can account for up to 72% of the variation forest AGB, with an RMSE = 45 Mg/ha. We also found that different kernels for the BRDF models influenced the weight parameters of the biomass inversion equation but did not significantly affect the estimated AGB. In conclusion, our method can enable the better usage of MODIS multi-angular observations for forest AGB estimation.

Aerosol loadings and their spatial distribution are among the most important atmospheric information needed for a range of applications such as air quality monitoring, climate research, and public health. A key measure of aerosol quantity is aerosol optical depth (AOD) and it has been routinely observed from space by Earth observing satellites/instrument, especially the Moderate Resolution Imaging Spectroradiometer (MODIS). Despite its global coverage and daily temporal resolution, the MODIS Multi-Angle Implementation of Atmospheric Correction (MAIAC) AOD product is fraught with missing values, severely limiting its use. A gap-filling method which is suitable for large-area application with high efficiency to obtain gapless AOD with reasonable spatial pattern and complete coverage is still lacking. Here, we proposed a novel spatiotemporal fitting algorithm to gap-fill the daily MODIS AOD product. Our algorithm is a multi-stage method aimed to address the non-stationary nature of AOD time series. First, the trend of daily AOD in a year in each pixel was fitted via smoothing splines and the residual was derived based on the original data and the trend. Second, the residual was spatially interpolated, leveraging the spatial correlation between the target pixel and the neighboring pixels. Third, the actual AOD was calculated as the sum of the trend and interpolated residual. We tested the algorithm against ground-based AOD data from 2011 to 2018 in China and further evaluated it via cross-validation at the global scale based on 10 selected MODIS tiles. Compared to the ground-reference AOD, the RMSE of our gapless datasets were 0.24 and 0.27 for Terra and Aqua, respectively; and the cross-validation showed a RMSE ranging from 0.045 to 0.055 (Terra) and 0.047 to 0.057 (Aqua) under different missing ratios. The novel gap-filling method outperforms the Interpolation-based Correlation Weighting (ICW) and Inverse Distance Weighting (IDW) algorithms in accuracy. Meanwhile, the gapless AOD using the novel algorithm shows lower accuracy than original MAIAC AOD, similar accuracy with the AOD from the Long-term Gap-free High-resolution Air Pollutants (LGHAP) concentration dataset, higher accuracy than the AOD from the Modern-Era Retrospective Analysis for Research and Applications, Version 2 (MERRA-2). Overall, the accuracy of gapless AOD using this algorithm meets the need of typical applications in relevant studies. The proposed algorithm is transferable to other regions, with the potential to be used even operationally and efficiently for generating accurate gapless global daily AOD datasets with the input of only MODIS MAIAC AOD data.

Earlier snowmelt, warmer temperatures and herbivory are among the factors that influence high-latitude tundra productivity near the town of Utqiaġvik in northern Alaska. However, our understanding of the potential interactions between these factors is limited. MODIS observations provide cover fractions of vegetation, snow, standing water, and soil, and fractional absorption of photosynthetically active radiation by canopy chlorophyll (fAPAR _chl ) per pixel. Here, we evaluated a recent time-period (2001–2014) that the tundra experienced large interannual variability in vegetation productivity metrics (i.e. fAPAR _chl and APAR _chl ), which was explainable by both abiotic and biotic factors. We found earlier snowmelt to increase soil and vegetation cover, and productivity in June, while warmer temperatures significantly increased monthly productivity. However, abiotic factors failed to explain stark decreases in productivity during August of 2008, which coincided with a severe lemming outbreak. MODIS observations found this tundra ecosystem to completely recover two years later, resulting in elevated productivity. This study highlights the potential roles of both climate and herbivory in modulating the interannual variability of remotely retrieved plant productivity metrics in Arctic coastal tundra ecosystems.

Multi-angle optical reflectance measurements such as those from the NASA moderate resolution imaging spectroradiometer (MODIS) are sensitive to forest 3D structures, potentially serving as a useful proxy to estimate forest structural variables such as aboveground biomass (AGB)—a potential theoretically recognized but rarely explored. In this paper, we examined the effectiveness of the reconstructed MODIS typical-angle reflectances—reflectances observed from the hotspot, darkspot, and nadir directions—for estimating forest AGB from both theoretical and practical perspectives. To gain theoretical insights, we first tested the sensitivities of typical-angle reflectances to forest AGB through simulations using the 4-scale bidirectional reflectance distribution function (BRDF) model. We then built statistical models to fit the relationship between MODIS multi-angle observations and field-measured deciduous-broadleaf/mixed-temperate forest AGB at five sites in the eastern USA, assisted by a semivariogram analysis to determine the effect of pixel heterogeneity on the MODIS–AGB relationship. We also determined the effects of terrain and season on the predictive relationships. Our results indicated that multi-angle reflectances with fewer visible shadows yielded better AGB estimates (hotspot: R2 = 0.63, RMSE = 54.28 Mg/ha; nadir: R2 = 0.55, RMSE = 59.95 Mg/ha; darkspot: R2 = 0.46, RMSE = 65.66 Mg/ha) after filtering out the effects of complex terrain and pixel heterogeneity; the MODIS typical-angle reflectances in the NIR band were the most sensitive to forest AGB. We also found strong sensitivities of estimated accuracies to MODIS image acquisition dates or season. Overall, our results suggest that the current practice of leveraging only single-angle MODIS data can be a suboptimal strategy for AGB estimation. We advocate the use of MODIS multi-angle reflectances for optical remote sensing of forest AGB or potentially other ecological applications requiring forest structure information.

Accurate estimation of forest aboveground biomass (AGB) is vital for informing ecosystem and carbon management. The Global Ecosystem Dynamics Investigation (GEDI) instrument—a new-generation spaceborne lidar system from NASA—provides the first global coverage of high-resolution 3D altimetry data aimed specifically for mapping Earth’s forests, but its performance is yet to be tested for large parts of the world. Here, our goal is to evaluate the accuracies of GEDI in measuring terrain, forest vertical structures, and AGB in reference to independent airborne lidar data over temperate and tropical forests in North America. We compared GEDI-derived elevations and canopy heights (e.g., relative height percentiles such as RH50 and RH100) with those from the Shuttle Radar Topography Mission (SRTM) or from two airborne lidar systems: the Laser Vegetation Imaging Sensor (LVIS) and Goddard’s Lidar, Hyperspectral and Thermal system (G-LiHT). We also estimated GEDI’s geolocation errors by matching GEDI waveforms and G-LiHT pseudo-waveforms. We assessed the predictive power of GEDI metrics in estimating AGB using Random Forests regression. Results showed that GEDI-derived ground elevations correlated strongly those from LVIS, G-LiHT, and LVIS (R2 > 0.91), but with nonnegligible RMSEs of 5.7 m (G-LiHT), 3.1 m (LVIS), and 10.9 m (SRTM). GEDI canopy heights had poorer correlation with LVIS (e.g., R2 = 0.44 for RH100) than with G-LiHT (e.g., R2 = 0.60 for RH100). The estimated horizontal geolocation errors of GEDI footprints averaged 6.5 meters, comparable to the nominal accuracy of 9 m. Correction for the locational errors improved the correlation of GEDI vs G-LiHT canopy heights significantly, on average by 53% (e.g., R2 from 0.57 to 0.82 for RH50). GEDI canopy metrics were useful for predicting AGB (R2 = 0.82 and RMSE = 19.1 Mg/Ha), with the maximum canopy height RH100 being the most useful predictor. Our results highlight the importance of accommodating or correcting for GEDI geolocation errors for estimating forest characteristics and provide empirical evidence on the utility of GEDI for monitoring global biomass dynamics from space.

Predictive modeling with remotely sensed data requires an accurate representation of spatial variability by ground truth data. In this study, we assessed the reliability of the size and location of ground truth data in capturing the landscape spatial variability embedded in the Airborne Visible Infrared Imaging Spectrometer-Next Generation (AVIRIS-NG) hyperspectral image in an agricultural region in Anand, India. We derived simulated spectral vegetation and soil indices using Gaussian simulation from AVIRIS-NG image for two point-location datasets, (1) ground truth points from adaptive sampling and (2) points from conditional Latin Hypercube Sampling (cLHS). We compared values of the simulated image indices against the actual image indices (measured) through the analysis of mean absolute errors. Modeling the variogram of the measured indices with the hyperspectral image in high spatial resolution (4m), is an effective way to characterize the spatial heterogeneity at the landscape level. We used geostatistical techniques to analyze the shapes of experimental variograms in order to assess whether or not the ground truth points, when compared against the cLHS-derived points, captured the spatial structures and variability of the studied agricultural area using measured indices. In addition, we explored the capability of the variogram by running tests in different point sample sizes. The ground truth and cLHS datasets were able to derive equivalent values for field spatial variability from image indices, according to our findings. Furthermore, this research presents a methodology for selecting spectral indices and determining the best sample size for efficiently replicating spatial patterns in hyperspectral images.

Change detection (CD) is essential for accurate understanding of land surface changes with multitemporal Earth observation data. Due to the great advantages in spatial information modeling, Morphological Attribute Profiles (MAPs) are becoming increasingly popular for improving the recognition ability in CD applications. However, most of the MAPs-based CD methods are implemented by setting the scale parameters of Attribute Profiles (APs) manually and ignoring the uncertainty of change information from different sources. To address these issues, a novel method for CD in high-resolution remote sensing (HRRS) images based on morphological attribute profiles and decision fusion is proposed in this study. By establishing the objective function based on the minimum of average interscale correlation, a morphological attribute profile with adaptive scale parameters (ASP-MAPs) is presented to exploit the spatial structure information. On this basis, a multifeature decision fusion framework based on the Dempster–Shafer (D-S) theory is constructed for obtaining the CD map. Experiments of multitemporal HRRS images from different sensors have shown that the proposed method outperforms the other advanced comparison CD methods, and the overall accuracy (OA) can reach more than 83.9%.

Hydropower in the Brazilian Amazon is a prevalent form of development, but dams have widespread and long-term environmental impacts that include deforestation in the areas surrounding the dams. Small hydropower plants (SHPs) are often perceived as having reduced environmental impacts compared to the large ones. In Brazil, SHPs are licensed by state governments, which have less strict requirements than the federal environmental agency. Brazil’s definition of ‘small’ dams has grown with successive increases in the maximum installed capacity from 10 to 30 to 50 megawatts (MW). This expanding loophole has increased the attractiveness of investing in multiple small dams rather than a single large dam, with resulting proliferation of SHPs. Forest dynamics surrounding the clustered SHPs when compared to single large dams are not well documented. In this study, we capitalized on a dense time series of satellite images to quantify and compare forest loss in the regions (over 110 000 km ^2 in area) surrounding 15 SHPs and 7 large dams at multiple watershed and buffer scales in the Brazilian Amazon for nearly two decades (2000–2018). The landscapes containing SHP clusters had lower cumulative forest loss as compared to those with large dams. However, when deforestation and hydroelectric generating capacity were jointly considered (i.e. forest loss per megawatt installed), we discovered an opposite trend. The regions surrounding the SHP clusters exhibited significant impacts ranging from 1.9 to 2.5 times that of the regions surrounding large dams across 5 km to sub-basin scales. Due to the considerable consequences of SHPs on deforestation, we argue that the rapid expansion of small hydropower should be approached with caution and requires more stringent environmental assessments.

Vegetation indices are widely used to derive land surface phenology (LSP). However, due to inconsistent illumination geometries, reflectance varies with solar zenith angles (SZA), which in turn affects the vegetation indices, and thus the derived LSP. To examine the SZA effect on LSP, the MODIS bidirectional reflectance distribution function (BRDF) product and a BRDF model were employed to derive LSPs under several constant SZAs (i.e., 0°, 15°, 30°, 45°, and 60°) in the Harvard Forest, Massachusetts, USA. The LSPs derived under varying SZAs from the MODIS nadir BRDF-adjusted reflectance (NBAR) and MODIS vegetation index products were used as baselines. The results show that with increasing SZA, NDVI increases but EVI decreases. The magnitude of SZA-induced NDVI/EVI changes suggests that EVI is more sensitive to varying SZAs than NDVI. NDVI and EVI are comparable in deriving the start of season (SOS), but EVI is more accurate when deriving the end of season (EOS). Specifically, NDVI/EVI-derived SOSs are relatively close to those derived from ground measurements, with an absolute mean difference of 8.01 days for NDVI-derived SOSs and 9.07 days for EVI-derived SOSs over ten years. However, a considerable lag exists for EOSs derived from vegetation indices, especially from the NDVI time series, with an absolute mean difference of 14.67 days relative to that derived from ground measurements. The SOSs derived from NDVI time series are generally earlier, while those from EVI time series are delayed. In contrast, the EOSs derived from NDVI time series are delayed; those derived from the simulated EVI time series under a fixed illumination geometry are also delayed, but those derived from the products with varying illumination geometries (i.e., MODIS NBAR product and MODIS vegetation index product) are advanced. LSPs derived from varying illumination geometries could lead to a difference spanning from a few days to a month in this case study, which highlights the importance of normalizing the illumination geometry when deriving LSP from NDVI/EVI time series.

Cover cropping is a conservation practice that helps to alleviate soil health problems and reduce nutrient losses. Understanding the spatial variability in historic and current adoption of cover cropping practices and their impacts on soil, water, and nutrient dynamics at a landscape scale is an important step in determining and prioritizing areas in a watershed to effectively utilize this practice. However, such data are lacking. Our objective was to develop a spatial and temporal inventory of winter cover cropping practices in the Maumee River watershed using images collected by Landsat satellites (Landsat 5, 7 and 8) from 2008 to 2019 in Google Earth Engine (GEE) platform. Each year, satellite images collected during cover crop growing season (i.e., between October and April) were converted into two seasonal composites based on cover crop phenology. Using these composites, various image-based covariates were extracted for 628 ground-truth (field) data. By integrating ground-truth and image-based covariates, a cover crop classification model based on a random forest (RF) algorithm was developed, trained and validated in GEE platform. Our classification scheme differentiated four cover crop categories: Winter Hardy, Winter Kill, Spring Emergent, and No Cover. The overall classification accuracy was 75%, with a kappa coefficient of 0.63. The results showed that more than 50% of the corn-soybean areas in the Maumee River watershed were without winter crops during 2008–2019 period. It was found that 2019/2020 and 2009/2010 were the years with the largest and lowest cover crop areas, with 34% and 10% in the watershed, respectively. The total cover cropping area was then assessed in relation to fall precipitation and cumulative growing degree days (GDD). There was no apparent increasing trend in cover crop areas between 2008 and 2019, but the variability in cover crops areas was found to be related to higher accumulated GDD and fall precipitation. A detailed understanding of the spatial and temporal distribution of cover crops using GEE could help in promoting site-specific management practices to enhance their environmental benefits. This also has significance to policy makers and funding agencies as they could use the information to localize areas in need of interventions for supporting adoption of cover cropping practice.

Clumping index (CI) is a canopy structural variable important for modeling the terrestrial biosphere, but its retrieval from remote sensing data remains one of the least reliable. The majority of regional or global CI products available so far were generated from multiangle optical reflectance data. However, these reflectance-based estimates have well-known limitations, such as the mere use of a linear relationship between the normalized difference hotspot and darkspot (NDHD) and CI, uncertainties in bidirectional reflectance distribution function (BRDF) models used to calculate the NDHD, and coarse spatial resolutions (e.g., hundreds of meters to several kilometers). To remedy these limitations and develop alternative methods for large-scale CI mapping, here we explored the use of spaceborne lidar—the Geoscience Laser Altimeter System (GLAS)—and proposed a semi-physical algorithm to estimate CI at the footprint level. Our algorithm was formulated to leverage the full vertical canopy profile information of the GLAS full-waveform data; it converted raw waveforms to forest canopy gap distributions and gap fractions of random canopies, which was used to estimate CI based on the radiative transfer theory and a revised Beer–Lambert model. We tested our algorithm over two areas in China—the Saihanba National Forest Park and Heilongjiang Province—and assessed its relative accuracies against field-measured CI and MODIS CI products. We found that reliable estimation of CI was possible only for GLAS waveforms with high signal-to-noise ratios (e.g., >65) and at gentle slopes (e.g., 2 = 0.72, RMSE = 0.07, and bias = 0.02), but they showed less correlation to MODIS CI (e.g., R2 = 0.26, RMSE = 0.12, and bias = 0.04). The difference highlights the impact of the scale effect in conducting comparisons of products with huge differences resolution. Overall, our analyses represent the first attempt to use spaceborne lidar to retrieve high-resolution forest CI and our algorithm holds promise for mapping CI globally.

The vertical foliage profile (VFP) and leaf area index (LAI) are critical descriptors in terrestrial ecosystem modeling. Although light detection and ranging (lidar) observations have been proven to have potential for deriving the VFP and LAI, existing methods depend only on the received waveform information and are sensitive to additional input parameters, such as the ratio of canopy to ground reflectance. In this study, we proposed a new method for retrieving forest VFP and LAI from Ice, Cloud and land Elevation Satellite (ICESat) Geoscience Laser Altimeter System (GLAS) data over two sites similar in their biophysical parameters. Our method utilized the information from not only the interaction between the laser and the forest but also the sensor configuration, which brought the benefit that our method was free from an empirical input parameter. Specifically, we first derived the transmitted energy profile (TEP) through the lidar 1-D radiative transfer model. Then, the obtained TEP was utilized to calculate the vertical gap distribution. Finally, the vertical gap distribution was taken as input to derive the VFP based on the Beer–Lambert law, and the LAI was calculated by integrating the VFP. Extensive validations of our method were carried out based on the discrete anisotropic radiative transfer (DART) simulation data, ground-based measurements, and the Moderate Resolution Imaging Spectroradiometer (MODIS) LAI product. The validation based on the DART simulation data showed that our method could effectively characterize the VFP and LAI under various canopy architecture scenarios, including homogeneous turbid and discrete individual-tree scenes. The ground-based validation also proved the feasibility of our method: the VFP retrieved from the GLAS data showed a similar trend with the foliage distribution density in the GLAS footprints; the GLAS LAI had a high correlation with the field measurements, with a determination coefficient (R2) of 0.79, root mean square error (RMSE) of 0.49, and bias of 0.17. Once the outliers caused by low data quality and large slope were identified and removed, the accuracy was further improved, with R2 = 0.85, RMSE = 0.35, and bias = 0.10. However, the MODIS LAI product did not present a good relationship with the GLAS LAI. Relative to the GLAS LAI, the MODIS LAI showed an overestimation in the low and middle ranges of the LAI and a saturation at high values of approximately LAI = 5.5. Overall, this method has the potential to produce continental- and global-scale VFP and LAI datasets from the spaceborne lidar system.

Accurate estimation of the fraction of absorbed photosynthetically active radiation (fPAR) for maize canopies are important for maize growth monitoring and yield estimation. The goal of this study is to explore the potential of using airborne LiDAR and hyperspectral data to better estimate maize fPAR. This study focuses on estimating maize fPAR from (1) height and coverage metrics derived from airborne LiDAR point cloud data; (2) vegetation indices derived from hyperspectral imagery; and (3) a combination of these metrics. Pearson correlation analyses were conducted to evaluate the relationships among LiDAR metrics, hyperspectral metrics, and field-measured fPAR values. Then, multiple linear regression (MLR) models were developed using these metrics. Results showed that (1) LiDAR height and coverage metrics provided good explanatory power (i.e., R2 = 0.81); (2) hyperspectral vegetation indices provided moderate interpretability (i.e., R2 = 0.50); and (3) the combination of LiDAR metrics and hyperspectral metrics improved the LiDAR model (i.e., R2 = 0.88). These results indicate that LiDAR model seems to offer a reliable method for estimating maize fPAR at a high spatial resolution and it can be used for farmland management. Combining LiDAR and hyperspectral metrics led to better performance of maize fPAR estimation than LiDAR or hyperspectral metrics alone, which means that maize fPAR retrieval can benefit from the complementary nature of LiDAR-detected canopy structure characteristics and hyperspectral-captured vegetation spectral information.

A novel adaptive morphological attribute profile under object boundary constraint (AMAP–OBC) method is proposed in this study for automatic building extraction from high-resolution remote sensing (HRRS) images. By investigating the associated attributes in morphological attribute profiles (MAPs), the proposed method establishes corresponding relationships between AMAP–OBC and building characteristics in HRRS images. In the preprocessing step, the candidate object set is extracted by a group of rules for screening of non-building objects. Second, based on the proposed adaptive scale parameter extraction and object boundary constraint strategies, AMAP–OBC is conducted to obtain the initial building set. Finally, a further identification strategy with adaptive threshold combination is proposed to obtain the final building extraction results. Through experiments of multiple groups of HRRS images from different sensors, the proposed method shows outstanding performance in terms of automatic building extraction from diverse geographic objects in urban scenes.

Figuring out the effect of the built-up environment on artificial light at night is essential for better understanding nighttime luminosity in both socioeconomic and ecological perspectives. However, there are few studies linking artificial surface properties to nighttime light (NTL). This study uses a statistical method to investigate effects of construction region environments on nighttime brightness and its variation with building height and regional economic development level. First, we extracted footprint-level target heights from Geoscience Laser Altimeter System (GLAS) waveform light detection and ranging (LiDAR) data. Then, we proposed a set of built-up environment properties, including building coverage, vegetation fraction, building height, and surface-area index, and then extracted these properties from GLAS-derived height, GlobeLand30 land-cover data, and DMSP/OLS radiance-calibrated NTL data. Next, the effects of non-building areas on NTL data were removed based on a supervised method. Finally, linear regression analyses were conducted to analyze the relationships between nighttime lights and built-up environment properties. Results showed that building coverage and vegetation fraction have weak correlations with nighttime lights (R2 < 0.2), building height has a moderate correlation with nighttime lights (R2 = 0.48), and surface-area index has a significant correlation with nighttime lights (R2 = 0.64). The results suggest that surface-area index is a more reasonable measure for estimating light number and intensity of NTL because it takes into account both building coverage and height, i.e., building surface area. Meanwhile, building height contributed to nighttime lights greater than building coverage. Further analysis showed the correlation between NTL and surface-area index becomes stronger with the increase of building height, while it is the weakest when the regional economic development level is the highest. In conclusion, these results can help us better understand the determinants of nighttime lights.

Full waveform (FW) LiDAR holds great potential for retrieving vegetation structure parameters at a high level of detail, but this prospect is constrained by practical factors such as the lack of available handy processing tools and the technical intricacy of waveform processing. This study introduces a new product named the Hyper Point Cloud (HPC), derived from FW LiDAR data, and explores its potential applications, such as tree crown delineation using the HPC-based intensity and percentile height (PH) surfaces, which shows promise as a solution to the constraints of using FW LiDAR data. The results of the HPC present a new direction for handling FW LiDAR data and offer prospects for studying the mid-story and understory of vegetation with high point density (~182 points/m2). The intensity-derived digital surface model (DSM) generated from the HPC shows that the ground region has higher maximum intensity (MAXI) and mean intensity (MI) than the vegetation region, while having lower total intensity (TI) and number of intensities (NI) at a given grid cell. Our analysis of intensity distribution contours at the individual tree level exhibit similar patterns, indicating that the MAXI and MI decrease from the tree crown center to the tree boundary, while a rising trend is observed for TI and NI. These intensity variable contours provide a theoretical justification for using HPC-based intensity surfaces to segment tree crowns and exploit their potential for extracting tree attributes. The HPC-based intensity surfaces and the HPC-based PH Canopy Height Models (CHM) demonstrate promising tree segmentation results comparable to the LiDAR-derived CHM for estimating tree attributes such as tree locations, crown widths and tree heights. We envision that products such as the HPC and the HPC-based intensity and height surfaces introduced in this study can open new perspectives for the use of FW LiDAR data and alleviate the technical barrier of exploring FW LiDAR data for detailed vegetation structure characterization.

This paper reviews LiDAR ground filtering algorithms used in the process of creating Digital Elevation Models. We discuss critical issues for the development and application of LiDAR ground filtering algorithms, including filtering procedures for different feature types, and criteria for study site selection, accuracy assessment, and algorithm classification. This review highlights three feature types for which current ground filtering algorithms are suboptimal, and which can be improved upon in future studies: surfaces with rough terrain or discontinuous slope, dense forest areas that laser beams cannot penetrate, and regions with low vegetation that is often ignored by ground filters.

Transition of freeze/thaw (F/T) affects land-atmospheric interactions and other biospheric dynamics. Global F/T statuses are normally monitored using microwave remote sensing, but at coarse resolutions (e.g., 25 km). Integration of coarse microwave remote sensing data with finer satellite products represents an opportunity to further enhance our ability to map F/T statuses regionally and globally. Here, we implemented and tested an approach to generate daily F/T status maps at a 5-km spatial resolution through the fusion of passive microwave data from AMSR2 and land surface temperature products from MODIS, using China as our study area for the year 2013 and 2014. Moreover, possible influences from elevation, vegetation, seasonality, etc., were also analyzed, as such analysis provides a direction to improve the approach. Overall, our freeze/thaw maps agreed well with ground reference observations, with an accuracy of ~86.6%. The new F/T maps helped to identify regions subject to frequent F/T transitions through the year, such as the Qinghai-Tibetan Plateau, Xinjiang, Gansu, Heilongjiang, Jilin, and Liaoning Province. This study indicates that the combination of AMSR2 and MODIS observations provides an effective method to obtain finer F/T maps (5-km or lower) for extensive regions. The finer F/T maps improve our knowledge of the F/T state detected by satellite remote sensing, and have a wide range of applications in regional studies considering land surface heterogeneity and models (e.g., community land models).

Photogrammetric UAV sees a surge in use for high-resolution mapping, but its use to map terrain under dense vegetation cover remains challenging due to a lack of exposed ground surfaces. This paper presents a novel object-oriented classification ensemble algorithm to leverage height, texture and contextual information of UAV data to improve landscape classification and terrain estimation. Its implementation incorporates multiple heuristics, such as multi-input machine learning-based classification, object-oriented ensemble, and integration of UAV and GPS surveys for terrain correction. Experiments based on a densely vegetated wetland restoration site showed classification improvement from 83.98% to 96.12% in overall accuracy and from 0.7806 to 0.947 in kappa value. Use of standard and existing UAV terrain mapping algorithms and software produced reliable digital terrain model only over exposed bare grounds (mean error = −0.019 m and RMSE = 0.035 m) but severely overestimated the terrain by ~80% of mean vegetation height in vegetated areas. The terrain correction method successfully reduced the mean error from 0.302 m to −0.002 m (RMSE from 0.342 m to 0.177 m) in low vegetation and from 1.305 m to 0.057 m (RMSE from 1.399 m to 0.550 m) in tall vegetation. Overall, this research validated a feasible solution to integrate UAV and RTK GPS for terrain mapping in densely vegetated environments.

BACKGROUND: A common challenge to the study of several infectious diseases consists in combining limited cross-sectional survey data, collected with a more sensitive detection method, with a more extensive (but biased) syndromic sentinel surveillance data, collected with a less sensitive method. Our article describes a novel modeling framework that overcomes this challenge, resulting in enhanced understanding of malaria in the Western Brazilian Amazon. METHODOLOGY/PRINCIPAL FINDINGS: A cohort of 486 individuals was monitored using four cross-sectional surveys, where all participants were sampled regardless of symptoms (aggressive-active case detection), resulting in 1,383 microscopy and 1,400 polymerase chain reaction tests. Data on the same individuals were also obtained from the local surveillance facility (i.e., passive and active case detection), totaling 1,694 microscopy tests. Our model accommodates these multiple pathogen and case detection methods. This model is shown to outperform logistic regression in terms of interpretability of its parameters, ability to recover the true parameter values, and predictive performance. We reveal that the main infection determinant was the extent of forest, particularly during the rainy season and in close proximity to water bodies, and participation on forest activities. We find that time residing in Acrelandia (as a proxy for past malaria exposure) decreases infection risk but surprisingly increases the likelihood of reporting symptoms once infected, possibly because non-naïve settlers are only susceptible to more virulent Plasmodium strains. We suggest that the search for asymptomatic carriers should focus on those at greater risk of being infected but lower risk of reporting symptoms once infected. CONCLUSIONS/SIGNIFICANCE: The modeling framework presented here combines cross-sectional survey data and syndromic sentinel surveillance data to shed light on several aspects of malaria that are critical for public health policy. This framework can be adapted to enhance inference on infectious diseases whenever asymptomatic carriers are important and multiple datasets are available.

Urbanization, a major driver of global change, profoundly impacts our physical and social world, for example, altering not just water and carbon cycling, biodiversity, and climate, but also demography, public health, and economy. Understanding these consequences for better scientific insights and effective decision-making unarguably requires accurate information on urban extent and its spatial distributions. We developed a method to map the urban extent from the defense meteorological satellite program/operational linescan system nighttime stable-light data at the global level and created a new global 1 km urban extent map for the year 2000. Our map shows that globally, urban is about 0.5% of total land area but ranges widely at the regional level, from 0.1% in Oceania to 2.3% in Europe. At the country level, urbanized land varies from about 0.01 to 10%, but is lower than 1% for most (70%) countries. Urbanization follows land mass distribution, as anticipated, with the highest concentration between 30° N and 45° N latitude and the largest longitudinal peak around 80° W. Based on a sensitivity analysis and comparison with other global urban area products, we found that our global product of urban areas provides a reliable estimate of global urban areas and offers the potential for producing a time-series of urban area maps for temporal dynamics analyses.