Spatial representativeness and uncertainty of eddy covariance carbon flux measurements for upscaling net ecosystem productivity to the grid scale

Eddy covariance (EC) measurements are often used to validate net ecosystem productivity (NEP) estimated from satellite remote sensing data and biogeochemical models. However, EC measurements represent an integrated flux over their footprint area, which usually differs from respective model grids or remote sensing pixels. Quantifying the uncertainties of scale mismatch associated with gridded flux estimates by upscaling single EC tower NEP measurements to the grid scale is an important but not yet fully investigated issue due to limited data availability as well as knowledge of flux variability at the grid scale. The Heihe Watershed Allied Telemetry Experimental Research (HiWATER) Multi-Scale Observation Experiment on Evapotranspiration (MUSOEXE) built a flux observation matrix that includes 17 EC towers within a 5 km × 5 km area in a heterogeneous agricultural landscape in northwestern China, providing an unprecedented opportunity to evaluate the uncertainty of upscaling due to spatial representative differences at the grid scale. Based on the HiWATER-MUSOEXE data, this study evaluated the spatial representativeness and uncertainty of EC CO2 flux measurements for upscaling to the grid scale using a scheme that combines a footprint model and a model-data fusion method. The results revealed the large spatial variability of gross primary productivity (GPP), ecosystem respiration (Re), and NEP within the study site during the growing season from 10 June to 14 September 2012. The variability of fluxes led to high variability in the representativeness of single EC towers for grid-scale NEP. The systematic underestimations of a single EC tower may reach 92(±11)%, 30(±11)%, and 165(±150)% and the overestimations may reach 25(±14)%, 20(±13)%, and 40(±33)% for GPP, Re, and NEP, respectively. This finding suggests that remotely sensed NEP at the global scale (e.g., MODIS products) should not be validated against single EC tower data in the case of heterogeneous surfaces. Any systematic bias should be addressed before upscaling EC data to grid scale. Otherwise, most of the systematic bias may be propagated to grid scale due to the scale dependence of model parameters. A systematic bias greater than 20% of the EC measurements can be corrected effectively using four indicators proposed in this study. These results will contribute to the understanding of spatial representativeness of EC towers within a heterogeneous landscape, to upscaling carbon fluxes from the footprint to the grid scale, to the selection of the location of EC towers, and to the reduction in the bias of NEP products by using an improved parameterization scheme of remote-sensing driven models, such as VPRM

Land Surface Temperature Product Validation Best Practice Protocol Version 1.0 - October, 2017

The Global Climate Observing System (GCOS) has specified the need to systematically generate andvalidate Land Surface Temperature (LST) products. This document provides recommendations on goodpractices for the validation of LST products. Internationally accepted definitions of LST, emissivity andassociated quantities are provided to ensure the compatibility across products and reference data sets. Asurvey of current validation capabilities indicates that progress is being made in terms of up-scaling and insitu measurement methods, but there is insufficient standardization with respect to performing andreporting statistically robust comparisons.Four LST validation approaches are identified: (1) Ground-based validation, which involvescomparisons with LST obtained from ground-based radiance measurements; (2) Scene-based intercomparisonof current satellite LST products with a heritage LST products; (3) Radiance-based validation,which is based on radiative transfer calculations for known atmospheric profiles and land surface emissivity;(4) Time series comparisons, which are particularly useful for detecting problems that can occur during aninstrument's life, e.g. calibration drift or unrealistic outliers due to undetected clouds. Finally, the need foran open access facility for performing LST product validation as well as accessing reference LST datasets isidentified

Evaluation of Sentinel-1, SMAP and SMOS surface soil moisture products for distributed eco-hydrological modelling in Mediterranean forest basins

[EN] Reliable distributed hydrological modeling, especially in semi-arid areas, must consider the inclusion of surface soil moisture (SSM) spatial information during the calibration process. This variable plays a key role in the evapotranspiration processes that determine the hydrological cycle. The coarse resolution of the SSM estimates by satellite remote sensing has restricted the application of this approach to only large basins, focusing most of the studies in the consideration of simply the temporal dynamics of this variable. The growing efforts in providing higher spatial resolution through disaggregating methodologies or new sensor estimates facilitates the application of this spatial approach to small basins. This paper explores the applicability of the currently available satellite surface soil moisture estimates for distributed eco-hydrological modelling in Mediterranean forest basins. On one hand, this study contributes to fill the existing research gap on the use of remote sensing SSM spatial patterns within the distributed hydrological modelling framework in small basins. On the other hand, it serves as an indirect validation method for the spatial performance of satellite SSM products. To achieve this goal, we implemented the eco-hydrological model TETIS in three case studies named: Hozgarganta (southern Spain), Ceira (western Portugal) and Carraixet (eastern Spain). The SSM estimates selected for comparison were Sentinel-1 SSM provided by the Copernicus Global Land Services (CGLS), SMAP SSM disaggregated using Sentinel-1 (SPL2SMAP_S) provided by the National Aeronautics and Space Administration (NASA), SMOS SSM provided by the Barcelona Expert Center (BEC), and SMOS and SMAP SSM disaggregated using the DISPATCH algorithm provided by Lobelia Earth. The methodology employed involved a multi-objective and multi-variable calibration in terms of remote sensing SSM spatial patterns and in-situ streamflow, using the Spatial Efficiency Metric (SPAEF) and the Nash-Sutcliffe efficiency index (NSE) respectively. Before model calibration a sensitivity analysis of the most influent variables was performed. The temporal and spatial comparison of the reference SSM products revealed inconsistencies amongst products. The disaggregating methodology determined the spatial agreement to a greater degree than the sensor itself (i.e. SMAP, SMOS). In spite of the differences amongst products, the multi-objective calibration approach proposed increased the robustness of the hydrological modelling.This study was founded by the Spanish AEI within the program WaterJPI through the project iAqueduct (PCI2019-103729) , by the EC Life project ResilientForests (LIFE17 CCA/ES/000063) , and by the project Water4Cast funded by Generalitat Valenciana (PROMETEO/2021/074) . We also acknowledge the following hydrometeorological data providers institutions: SiAR, SAIH-HIDROSUR, SAIH Jucar and SNIRH.Gomis-Cebolla, J.; Garcia-Arias, A.; Perpinyà-Vallès, M.; Francés, F. (2022). Evaluation of Sentinel-1, SMAP and SMOS surface soil moisture products for distributed eco-hydrological modelling in Mediterranean forest basins. Journal of Hydrology. 608:1-19. https://doi.org/10.1016/j.jhydrol.2022.12756911960

Spatiotemporally representative and cost-efficient sampling design for validation activities in wanglang experimental site

Altres ajuts: EC Copernicus Global Land Service (CGLOPS-1, 199494-JRC.Spatiotemporally representative Elementary Sampling Units (ESUs) are required for capturing the temporal variations in surface spatial heterogeneity through field measurements. Since inaccessibility often coexists with heterogeneity, a cost-efficient sampling design is mandatory. We proposed a sampling strategy to generate spatiotemporally representative and cost-efficient ESUs based on the conditioned Latin hypercube sampling scheme. The proposed strategy was constrained by multi-temporal Normalized Difference Vegetation Index (NDVI) imagery, and the ESUs were limited within a sampling feasible region established based on accessibility criteria. A novel criterion based on the Overlapping Area (OA) between the NDVI frequency distribution histogram from the sampled ESUs and that from the entire study area was used to assess the sampling efficiency. A case study inWanglang National Nature Reserve in China showed that the proposed strategy improves the spatiotemporally representativeness of sampling (mean annual OA = 74.7%) compared to the single-temporally constrained (OA = 68.7%) and the random sampling (OA = 63.1%) strategies. The introduction of the feasible region constraint significantly reduces in-situ labour-intensive characterization necessities at expenses of about 9% loss in the spatiotemporal representativeness of the sampling. Our study will support the validation activities in Wanglang experimental site providing a benchmark for locating the nodes of automatic observation systems (e.g., LAINet) which need a spatially distributed and temporally fixed sampling design

Monitoring Land Surface Albedo and Vegetation Dynamics Using High Spatial and Temporal Resolution Synthetic Time Series from Landsat and the MODIS BRDF/NBAR/Albedo Product

Seasonal vegetation phenology can significantly alter surface albedo which in turn affects the global energy balance and the albedo warmingcooling feedbacks that impact climate change. To monitor and quantify the surface dynamics of heterogeneous landscapes, high temporal and spatial resolution synthetic time series of albedo and the enhanced vegetation index (EVI) were generated from the 500-meter Moderate Resolution Imaging Spectroradiometer (MODIS) operational Collection V006 daily BRDF (Bidirectional Reflectance Distribution Function) / NBAR (Nadir BRDF-Adjusted Reflectance) / albedo products and 30-meter Landsat 5 albedo and near-nadir reflectance data through the use of the Spatial and Temporal Adaptive Reflectance Fusion Model (STARFM). The traditional Landsat Albedo (Shuai et al., 2011) makes use of the MODIS BRDFAlbedo products (MCD43) by assigning appropriate BRDFs from coincident MODIS products to each Landsat image to generate a 30-meter Landsat albedo product for that acquisition date. The available cloud free Landsat 5 albedos (due to clouds, generated every 16 days at best) were used in conjunction with the daily MODIS albedos to determine the appropriate 30-meter albedos for the intervening daily time steps in this study. These enhanced daily 30-meter spatial resolution synthetic time series were then used to track albedo and vegetation phenology dynamics over three Ameriflux tower sites (Harvard Forest in 2007, Santa Rita in 2011 and Walker Branch in 2005). These Ameriflux sites were chosen as they are all quite nearby new towers coming on line for the National Ecological Observatory Network (NEON), and thus represent locations which will be served by spatially paired albedo measures in the near future. The availability of data from the NEON towers will greatly expand the sources of tower albedometer data available for evaluation of satellite products. At these three Ameriflux tower sites the synthetic time series of broadband shortwave albedos were evaluated using the tower albedo measurements with a Root Mean Square Error (RMSE) less than 0.013 and a bias within the range of 0.006. These synthetic time series provide much greater spatial detail than the 500 meter gridded MODIS data, especially over more heterogeneous surfaces, which improves the efforts to characterize and monitor the spatial variation across species and communities. The mean of the difference between maximum and minimum synthetic time series of albedo within the MODIS pixels over a subset of satellite data of Harvard Forest (16 kilometers by 14 kilometers) was as high as 0.2 during the snow-covered period and reduced to around 0.1 during the snow-free period. Similarly, we have used STARFM to also couple MODIS Nadir BRDF-Adjusted Reflectances (NBAR) values with Landsat 5 reflectances to generate daily synthetic times series of NBAR and thus Enhanced Vegetation Index (NBAR-EVI) at a 30-meter resolution. While normally STARFM is used with directional reflectances, the use of the view angle corrected daily MODIS NBAR values will provide more consistent time series. These synthetic times series of EVI are shown to capture seasonal vegetation dynamics with finer spatial and temporal details, especially over heterogeneous land surfaces

Earth Observations and Integrative Models in Support of Food and Water Security

Global food production depends upon many factors that Earth observing satellites routinely measure about water, energy, weather, and ecosystems. Increasingly sophisticated, publicly-available satellite data products can improve efficiencies in resource management and provide earlier indication of environmental disruption. Satellite remote sensing provides a consistent, long-term record that can be used effectively to detect large-scale features over time, such as a developing drought. Accuracy and capabilities have increased along with the range of Earth observations and derived products that can support food security decisions with actionable information. This paper highlights major capabilities facilitated by satellite observations and physical models that have been developed and validated using remotely-sensed observations. Although we primarily focus on variables relevant to agriculture, we also include a brief description of the growing use of Earth observations in support of aquaculture and fisheries

Evaluating the two-source energy balance model using local thermal and surface flux observations in a strongly advective irrigated agricultural area

Application and validation of many thermal remote sensing-based energy balance models involve the use of local meteorological inputs of incoming solar radiation, wind speed and air temperature as well as accurate land surface temperature (LST), vegetation cover and surface flux measurements. For operational applications at large scales, such local information is not routinely available. In addition, the uncertainty in LST estimates can be several degrees due to sensor calibration issues, atmospheric effects and spatial variations in surface emissivity. Time differencing techniques using multi-temporal thermal remote sensing observations have been developed to reduce errors associated with deriving the surface- air temperature gradient, particularly in complex landscapes. The Dual-Temperature-Difference (DTD) method addresses these issues by utilizing the Two-Source Energy Balance (TSEB) model of Norman et al. (1995) [1], and is a relatively simple scheme requiring meteorological input from standard synoptic weather station networks or mesoscale modeling. A comparison of the TSEB and DTD schemes is performed using LST and flux observations from eddy covariance (EC) flux towers and large weighing lysimeters (LYs) in irrigated cotton fields collected during BEAREX08, a large-scale field experiment conducted in the semi-arid climate of the Texas High Plains as described by Evett et al. (2012) [2]. Model output of the energy fluxes (i.e., net radiation, soil heat flux, sensible and latent heat flux) generated with DTD and TSEB using local and remote meteorological observations are compared with EC and LY observations. The DTD method is found to be significantly more robust in flux estimation compared to the TSEB using the remote meteorological observations. However, discrepancies between model and measured fluxes are also found to be significantly affected by the local inputs of LST and vegetation cover and the representativeness of the remote sensing observations with the local flux measurement footprint

Intercomparison of phenological transition dates derived from the PhenoCam Dataset V1.0 and MODIS satellite remote sensing

Phenology is a valuable diagnostic of ecosystem health, and has applications to environmental monitoring and management. Here, we conduct an intercomparison analysis using phenological transition dates derived from near-surface PhenoCam imagery and MODIS satellite remote sensing. We used approximately 600 site-years of data, from 128 camera sites covering a wide range of vegetation types and climate zones. During both “greenness rising” and “greenness falling” transition phases, we found generally good agreement between PhenoCam and MODIS transition dates for agricultural, deciduous forest, and grassland sites, provided that the vegetation in the camera field of view was representative of the broader landscape. The correlation between PhenoCam and MODIS transition dates was poor for evergreen forest sites. We discuss potential reasons (including sub-pixel spatial heterogeneity, flexibility of the transition date extraction method, vegetation index sensitivity in evergreen systems, and PhenoCam geolocation uncertainty) for varying agreement between time series of vegetation indices derived from PhenoCam and MODIS imagery. This analysis increases our confidence in the ability of satellite remote sensing to accurately characterize seasonal dynamics in a range of ecosystems, and provides a basis for interpreting those dynamics in the context of tangible phenological changes occurring on the ground