Impact of day/night time land surface temperature in soil moisture disaggregation algorithms

Since its launch in 2009, the ESA’s SMOS mission is providing global soil moisture (SM) maps at ~40 km, using the first L-band microwave radiometer on space. Its spatial resolution meets the needs of global applications, but prevents the use of the data in regional or local applications, which require higher spatial resolutions (~1-10 km). SM disaggregation algorithms based generally on the land surface temperature (LST) and vegetation indices have been developed to bridge this gap. This study analyzes the SM-LST relationship at a variety of LST acquisition times and its influence on SM disaggregation algorithms. Two years of in situ and satellite data over the central part of the river Duero basin and the Iberian Peninsula are used. In situ results show a strong anticorrelation of SM to daily maximum LST (R˜-0.5 to -0.8). This is confirmed with SMOS SM and MODIS LST Terra/Aqua at day time-overpasses (R˜-0.4 to -0.7). Better statistics are obtained when using MODIS LST day (R˜0.55 to 0.85; ubRMSD˜0.04 to 0.06 m3 /m3 ) than LST night (R˜0.45 to 0.80; ubRMSD˜0.04 to 0.07 m3 /m3 ) in the SM disaggregation. An averaged ensemble of day and night MODIS LST Terra/Aqua disaggregated SM estimates also leads to robust statistics (R˜0.55 to 0.85; ubRMSD˜0.04 to 0.07 m3 /m3 ) with a coverage improvement of ~10-20 %.Peer ReviewedPostprint (published version

Revisiting the contribution of transpiration to global terrestrial evapotranspiration

Even though knowing the contributions of transpiration (T), soil and open water evaporation (E), and interception (I) to terrestrial evapotranspiration (ET=T+E+I) is crucial for understanding the hydrological cycle and its connection to ecological processes, the fraction of T is unattainable by traditional measurement techniques over large scales. Previously reported global mean T/(E+T+I) from multiple independent sources, including satellite-based estimations, reanalysis, land surface models, and isotopic measurements, varies substantially from 24% to 90%. Here we develop a new ET partitioning algorithm, which combines global evapotranspiration estimates and relationships between leaf area index (LAI) and T/(E+T) for different vegetation types, to upscale a wide range of published site-scale measurements. We show that transpiration accounts for about 57.2% (with standard deviation6.8%) of global terrestrial ET. Our approach bridges the scale gap between site measurements and global model simulations,and can be simply implemented into current global climate models to improve biological CO2 flux simulations

Bayesian multimodel estimation of global terrestrial latent heat flux from eddy covariance, meteorological, and satellite observations

Accurate estimation of the satellite-based global terrestrial latent heat flux (LE) at high spatial and temporal scales remains a major challenge. In this study, we introduce a Bayesian model averaging (BMA) method to improve satellite-based global terrestrial LE estimation by merging five process-based algorithms. These are the Moderate Resolution Imaging Spectroradiometer (MODIS) LE product algorithm, the revised remote-sensing-based Penman-Monteith LE algorithm, the Priestley-Taylor-based LE algorithm, the modified satellite-based Priestley-Taylor LE algorithm, and the semi-empirical Penman LE algorithm. We validated the BMA method using data for 2000–2009 and by comparison with a simple model averaging (SA) method and five process-based algorithms. Validation data were collected for 240 globally distributed eddy covariance tower sites provided by FLUXNET projects. The validation results demonstrate that the five process-based algorithms used have variable uncertainty and the BMA method enhances the daily LE estimates, with smaller root mean square errors (RMSEs) than the SA method and the individual algorithms driven by tower-specific meteorology and Modern Era Retrospective Analysis for Research and Applications (MERRA) meteorological data provided by the NASA Global Modeling and Assimilation Office (GMAO), respectively. The average RMSE for the BMA method driven by daily tower-specific meteorology decreased by more than 5 W/m2 for crop and grass sites, and by more than 6 W/m2 for forest, shrub, and savanna sites. The average coefficients of determination (R2) increased by approximately 0.05 for most sites. To test the BMA method for regional mapping, we applied it for MODIS data and GMAO-MERRA meteorology to map annual global terrestrial LE averaged over 2001–2004 for spatial resolution of 0.05°. The BMA method provides a basis for generating a long-term global terrestrial LE product for characterizing global energy, hydrological, and carbon cycles

Impact of day/night time land surface temperature in soil moisture disaggregation algorithms

18 pages, 5 figures, 1 tableSince its launch in 2009, the ESA’s SMOS mission is providing global soil moisture (SM) maps at ~40 km, using the first L-band microwave radiometer on space. Its spatial resolution meets the needs of global applications, but prevents the use of the data in regional or local applications, which require higher spatial resolutions (~1-10 km). SM disaggregation algorithms based generally on the land surface temperature (LST) and vegetation indices have been developed to bridge this gap. This study analyzes the SM-LST relationship at a variety of LST acquisition times and its influence on SM disaggregation algorithms. Two years of in situ and satellite data over the central part of the river Duero basin and the Iberian Peninsula are used. In situ results show a strong anticorrelation of SM to daily maximum LST (R≈0.5 to -0.8). This is confirmed with SMOS SM and MODIS LST Terra/Aqua at day time-overpasses (R≈-0.4 to -0.7). Better statistics are obtained when using MODIS LST day (R≈0.55 to 0.85; ubRMSD≈0.04 to 0.06 m/m) than LST night (R≈0.45 to 0.80; ubRMSD≈0.04 to 0.07 m/m) in the SM disaggregation. An averaged ensemble of day and night MODIS LST Terra/Aqua disaggregated SM estimates also leads to robust statistics (R≈0.55 to 0.85; ubRMSD≈0.04 to 0.07 m/m) with a coverage improvement of~10-20 %This work was supported by the Spanish Ministry of Economy and Competitiveness, through a Formación Personal Investigador (FPI) grant BES-2011-043322, the project PROMISES: ESP2015-67549-C3, ERDF (European Regional Development Fund) and the BBVA foundationPeer Reviewe

Multiscale Assimilation of Sentinel and Landsat Data for Soil Moisture and Leaf Area Index Predictions Using an Ensemble-Kalman-Filter-Based Assimilation Approach in a Heterogeneous Ecosystem

Data assimilation techniques allow researchers to optimally merge remote sensing observations in ecohydrological models, guiding them for improving land surface fluxes predictions. Presently, freely available remote sensing products, such as those of Sentinel 1 radar, Landsat 8 sensors, and Sentinel 2 sensors, allow the monitoring of land surface variables (e.g., radar backscatter for soil moisture and the normalized difference vegetation index (NDVI) and for leaf area index (LAI)) at unprecedentedly high spatial and time resolutions, appropriate for heterogeneous ecosystems, typical of semiarid ecosystems characterized by contrasting vegetation components (grass and trees) competing for water use. A multiscale assimilation approach that assimilates radar backscatter and grass and tree NDVI in a coupled vegetation dynamic-land surface model is proposed. It is based on the ensemble Kalman filter (EnKF), and it is not limited to assimilating remote sensing data for model predictions, but it uses assimilated data for dynamically updating key model parameters (the ENKFdc approach), including saturated hydraulic conductivity and grass and tree maintenance respiration coefficients, which are highly sensitive parameters of soil-water balance and biomass budget models, respectively. The proposed EnKFdc assimilation approach facilitated good predictions of soil moisture, grass, and tree LAI in a heterogeneous ecosystem in Sardinia for a 3-year period with contrasting hydrometeorological (dry vs. wet) conditions. Contrary to the EnKF-based approach, the proposed EnKFdc approach performed well for the full range of hydrometeorological conditions and parameters, even assuming extremely biased model conditions with very high or low parameter values compared with the calibrated ("true") values. The EnKFdc approach is crucial for soil moisture and LAI predictions in winter and spring, key seasons for water resources management in Mediterranean water-limited ecosystems. The use of ENKFdc also enabled us to predict evapotranspiration and carbon flux well, with errors of less than 4% and 15%, respectively; such results were obtained even with extremely biased initial model conditions

Evaluation of global terrestrial evapotranspiration using state-of-the-art approaches in remote sensing, machine learning and land surface modeling

Evapotranspiration (ET) is critical in linking global water, carbon and energy cycles. However, direct measurement of global terrestrial ET is not feasible. Here, we first reviewed the basic theory and state-of-the-art approaches for estimating global terrestrial ET, including remote-sensing-based physical models, machine-learning algorithms and land surface models (LSMs). We then utilized 4 remote-sensing-based physical models, 2 machine-learning algorithms and 14 LSMs to analyze the spatial and temporal variations in global terrestrial ET. The results showed that the ensemble means of annual global terrestrial ET estimated by these three categories of approaches agreed well, with values ranging from 589.6 mm yr−1 (6.56×104 km3 yr−1) to 617.1 mm yr−1 (6.87×104 km3 yr−1). For the period from 1982 to 2011, both the ensembles of remote-sensing-based physical models and machine-learning algorithms suggested increasing trends in global terrestrial ET (0.62 mm yr−2 with a significance level of p0.05), although many of the individual LSMs reproduced an increasing trend. Nevertheless, all 20 models used in this study showed that anthropogenic Earth greening had a positive role in increasing terrestrial ET. The concurrent small interannual variability, i.e., relative stability, found in all estimates of global terrestrial ET, suggests that a potential planetary boundary exists in regulating global terrestrial ET, with the value of this boundary being around 600 mm yr−1. Uncertainties among approaches were identified in specific regions, particularly in the Amazon Basin and arid/semiarid regions. Improvements in parameterizing water stress and canopy dynamics, the utilization of new available satellite retrievals and deep-learning methods, and model–data fusion will advance our predictive understanding of global terrestrial ET

Incorporating canopy structure from simulated GEDI lidar into bird species distribution models

The Global Ecosystem Dynamics Investigation (GEDI) lidar began data acquisition from the International Space Station in March 2019 and is expected to make over 10 billion measurements of canopy structure and topography over two years. Previously, airborne lidar data with limited spatial coverage have been used to examine relationships between forest canopy structure and faunal diversity, most commonly bird species. GEDI’s latitudinal coverage will permit these types of analyses at larger spatial extents, over the majority of the Earth’s forests, and most importantly in areas where canopy structure is complex and/or poorly understood. In this regional study, we examined the impact that GEDI-derived Canopy Structure variables have on the performance of bird species distribution models (SDMs) in Sonoma County, California. We simulated GEDI waveforms for a two-year period and then interpolated derived Canopy Structure variables to three grid sizes of analysis. In addition to these variables, we also included Phenology, Climate, and other Auxiliary variables to predict the probability of occurrence of 25 common bird species. We used a weighted average ensemble of seven individual machine learning models to make predictions for each species and calculated variable importance. We found that Canopy Structure variables were, on average at our finest resolution of 250 m, the second most important group (32.5%) of predictor variables after Climate variables (35.3%). Canopy Structure variables were most important for predicting probability of occurrence of birds associated with Conifer forest habitat. Regarding spatial analysis scale, we found that finer-scale models more frequently performed better than coarser-scale models, and the importance of Canopy Structure variables was greater at finer spatial resolutions. Overall, GEDI Canopy Structure variables improved SDM performance for at least one spatial resolution for 19 of 25 species and thus show promise for improving models of bird species occurrence and mapping potential habitat

Assimilation of NDVI data in a land surface – Vegetation model for leaf area index predictions in a tree-grass ecosystem

Periodic observations of vegetation index, such as the normalized difference vegetation index (NDVI), can be used for data assimilation in heterogenous ecosystems. Indeed, the new Sentinel 2 Multispectral instrument and Landsat 8 Operational Land Imager sensor data are available at such high temporal and spatial resolutions that can be used to detect the patches of the main vegetation components (grass and trees) of heterogenous ecosystems, and capture their dynamics. We demonstrate the possibility to merge grass and tree NDVI observations and models, to optimally provide robust predictions of grass and tree leaf area index. The proposed assimilation approach assimilates NDVI data through the Ensemble Kalman filter (EnKF) and dynamically calibrates a key vegetation dynamic model parameter, the maintenance respiration coefficient (ma ). In the presence of large bias of the grass and tree ma base values, only the use of the proposed assimilation approach removes the large bias in the biomass balance, dynamically calibrating maintenance respiration coefficients, and corrects the model. The performance of a land surface - vegetation model was improved by assimilating observations of NDVI. The effective impact of the proposed assimilation approach on the evapotranspiration and CO2 uptake predictions in the heterogenous ecosystem is also demonstrated

Assessing the utility of remote sensing data to accurately estimate changes in groundwater storage

Accurate and timely estimates of groundwater storage changes are critical to the sustainable management of aquifers worldwide, but are hindered by the lack of in-situ groundwater measurements in most regions. Hydrologic remote sensing measurements provide a potential pathway to quantify groundwater storage changes by closing the water balance, but the degree to which remote sensing data can accurately estimate groundwater storage changes is unclear. In this study, we quantified groundwater storage changes in California\u27s Central Valley at two spatial scales for the period 2002 through 2020 using remote sensing data and an ensemble water balance method. To evaluate performance, we compared estimates of groundwater storage changes to three independent estimates: GRACE satellite data, groundwater wells and a groundwater flow model. Results suggest evapotranspiration has the highest uncertainty among water balance components, while precipitation has the lowest. We found that remote sensing-based groundwater storage estimates correlated well with independent estimates; annual trends during droughts fall within 15% of trends calculated using wells and groundwater models within the Central Valley. Remote sensing-based estimates also reliably estimated the long-term trend, seasonality, and rate of groundwater depletion during major drought events. Additionally, our study suggests that the proposed method estimate changes in groundwater at sub-annual latencies, which is not currently possible using other methods. The findings have implications for improving the understanding of aquifer dynamics and can inform regional water managers about the status of groundwater systems during droughts

Inconsistencies of interannual variability and trends in long-term satellite leaf area index products

Understanding the long-term performance of global satellite leaf area index (LAI) products is important for global change research. However, few effort has been devoted to evaluating the long-term time-series consistencies of LAI products. This study compared four long-term LAI products (GLASS, GLOBMAP, LAI3g, and TCDR) in terms of trends, interannual variabilities, and uncertainty variations from 1982 through 2011. This study also used four ancillary LAI products (GEOV1, MERIS, MODIS C5, and MODIS C6) from 2003 through 2011 to help clarify the performances of the four long-term LAI products. In general, there were marked discrepancies between the four long-term LAI products. During the pre-MODIS period (1982-1999), both linear trends and interannual variabilities of global mean LAI followed the order GLASS>LAI3g>TCDR>GLOBMAP. The GLASS linear trend and interannual variability were almost 4.5 times those of GLOBMAP. During the overlap period (2003-2011), GLASS and GLOBMAP exhibited a decreasing trend, TCDR no trend, and LAI3g an increasing trend. GEOV1, MERIS, and MODIS C6 also exhibited an increasing trend, but to a much smaller extent than that from LAI3g. During both periods, the R2 of detrended anomalies between the four long-term LAI products was smaller than 0.4 for most regions. Interannual variabilities of the four long-term LAI products were considerably different over the two periods, and the differences followed the order GLASS>LAI3g>TCDR>GLOBMAP. Uncertainty variations quantified by a collocation error model followed the same order. Our results indicate that the four long-term LAI products were neither intraconsistent over time nor interconsistent with each other. These inconsistencies may be due to NOAA satellite orbit changes and MODIS sensor degradation. Caution should be used in the interpretation of global changes derived from the four long-term LAI products