Land Surface Parameters Research Articles

Simulation of ecosystem fluxes with the SCOPE model: Sensitivity to parametrization and evaluation with flux tower observations

Accurate estimates of carbon, water and energy fluxes between the Earth surface and the atmosphere are crucial for enhancing our understanding of ecosystem–climate interactions. Such estimates can be made by combining remote sensing derived land surface parameters with climate reanalysis data. We analysed to what degree generic (plant functional type (PFT)-independent) satellite-derived vegetation properties and climate reanalysis data can explain land surface fluxes and to what extent the PFT-specific information extends the flux simulations. For this purpose, we used the Soil Canopy Observation, Photochemistry and Energy fluxes (SCOPE) model, which combines radiative transfer in plant leaves and vegetation canopies with photosynthesis and energy balance in a single model representation of the vegetation. We evaluated the performance of SCOPE in simulating fluxes by comparison to 63 eddy covariance sites representing 10 PFTs. We varied the sources of maximum carboxylation capacity (Vcmax25) and BallBerrySlope values (default vs literature), the seasonality of Vcmax25 and the meteorological forcing (locally measured and climate reanalysis). The average performance of daily flux in terms of root-mean-square error (RMSE) was 2.3 ± 0.8 μmol CO2 m−2 s−1 (R2= 0.74 ± 0.12) for gross primary productivity (GPP), 24 ± 8 W m−2 (R2= 0.68 ± 0.16) for latent heat flux (λE) and 50 ± 15 W m−2 (R20.47±0.17) for sensible heat flux (H). The inter-site variability of the annual accumulated GPP flux was captured well with seasonally varying PFT-specific Vcmax25 (R2= 0.74, RMSE = 308 g C m−2 yr−1 and bias = −68 g C m−2 yr−1). The annual accumulated evapotranspiration (ET) was overestimated (R2= 0.31, RMSE = 101 mm yr−1 and bias = 37 mm yr−1), mainly in the ecosystems with subtropical Mediterranean climate, for which the soil resistance to evaporation from porous space (rss) had to be constrained from soil moisture content (SMC) or land surface temperature (LST). Overall, the study demonstrates that SCOPE model can simulate ecosystem flux with high accuracy without site-specific calibration of its parameters.

Characterization of Hydrologic Sand and Dust Storm Sources in the Middle East

Due to diverse hydroclimatic conditions and human interventions, the Middle East hosts a variety of active sources of sand and dust storms (SDS). Discrimination of different types of SDS sources is the most important factor for adopting optimal mitigation measures to combat SDS. This study employed a binary mask-based modeling framework to identify Middle East SDS sources. Accordingly, using time series of remotely sensed data of land surface and atmospheric aerosol parameters, SDS sources covering an area of 1 million Km2 were identified with an overall accuracy of 82.6%. Considering the type of land use and spatial-temporal changes in water bodies, SDS sources were categorized into seven types in terms of origin. Desert sources have the largest share (>79%), whereas hydrologic sources accounted for about 8.4%. The results showed that water bodies had a declining trend after 2000. The occurrence of two severe drought periods in 2000–2001 and 2007–2012 led to a 52% decrease in water bodies and a 14–37% increase in SDS emission compared to the pre-2000 period. The latter drought period also led to a sharp decrease in groundwater resources across the region. Our results revealed that natural circumstances and drought actively contribute to the depletion of water resources that led to the formation of SDS sources in the Middle East, while the role of anthropogenic factors is predominant in the case of hydrologic SDS sources.

Quantification of model uncertainty in sub-daily extreme precipitation projections

Sub-daily precipitation is projected to intensify more rapidly with global warming, that in turn may lead to frequent disasters such as flash floods, landslides, etc. However, the precipitation projections from Global Climate Model (GCMs) are known to be highly uncertain, with an increasing magnitude at finer spatio-temporal scales. Hence, realistic quantification of the inherent uncertainty is vital to improve the confidence in these projections, for the mitigation/adaptation applications. In the present study, the model uncertainty in sub-daily (3-hourly) precipitation extremes is estimated, by considering 20 Coupled Model Intercomparison Project Phase 5 (CMIP5) models, during four different seasons for the near (2026–45) and far-future (2081–99) periods. An uncertainty quantification approach using “mutual information based independence weight” is proposed in this study, to account the linear and nonlinear dependence between GCMs along with the equitability property. The CMCC-CM model is found to be the most independent model having highest independence weight (0.064) in the ensemble. The uncertainty in northern hemisphere is found to be more during June to August period, whereas in southern hemispheres during December to February period. Overall, there is an increase in model uncertainty in far-future, especially over the equatorial regions. However, the behaviour of different indices is not in agreement with each other. For example, extreme characteristics such as monthly maximum 3hourly-precipitation has increasing uncertainty; whereas the mean characteristic like, simple hourly-precipitation intensity index possesses unchanging uncertainty. Model uncertainty is most sensitive to the precipitation characteristics of the regions, followed by elevation, a significant geophysical factor. The study also highlighted the need for model improvement over the grids having elevation <2 km and above 4 km. The study identifies the possible regions, seasons, factors etc. where model improvement can be made. Moreover, it is anticipated that improvement of resolution, updation of convective parameterization schemes and land surface parameterization schemes in GCMs may enhance the skill of GCMs and can improve the confidence in the projection of precipitation on sub-daily scale. • Seasonal model uncertainty is estimated to assess the reliability in the projections of sub-daily precipitation extremes. • Mutual information-based weighting technique is used to establish the independence weight of models. • Different aspects of sub-daily precipitation exhibit differing uncertainty, in terms of magnitude and also temporal evolution. • Precipitation intensity remains the dominant factor behind model uncertainty, but geophysical influences cannot be ignored. • Reliability on the sub-daily mean characteristics improves with time; whereas that on extreme characteristics declines.

Sea temperature variation associated with the 2021 Haiti Mw 7.2 earthquake and possible mechanism

A catastrophic earthquake (magnitude Mw 7.2) occurred on 14 August 2021 in Haiti. Sea surface temperature (SST) and sea potential temperature (SPT) observed from satellites show pronounced changes prior to the main earthquake event. The changes are observed near the epicenter and the boundary of the Caribbean Plate. Here we have carried out a detailed analysis of ocean, atmosphere and land surface parameters (sea water salinity, chlorophyll-a, latent energy flux, sensible heat flux, specific humidity, air temperature, methane and land surface temperature) using multiple satellite data to understand the possible mechanism of sea temperature changes. The strong coupling among these geophysical and geochemical parameters has been found within 2-week prior to the 2021 Haiti Mw 7.2 earthquake due to the interaction between ocean/land and atmosphere during the last phase of earthquake preparation period caused by the tectonic movement. The spatiotemporal distribution pattern of anomalies is discussed in view of the focal mechanism of the earthquake event.

Point-to-Surface Upscaling Algorithms for Snow Depth Ground Observations

To validate the accuracy of snow depth products retrieved from passive microwave remote sensing data with a high confidence level, the verification method based on points of ground observation is subject to great uncertainty, due to the scale effect. Thus, it is necessary to use a point-to-surface scale transformation method to obtain the relative ground truth at the remote sensing pixel scale. In this study, by using the snow depth ground observations at different observation scales, the upscaling methods are conducted based on simple average (SA), geostatistical, Bayes maximum entropy (BME), and random forest (RF) algorithms. In addition, the cross-validation of the leave-one-out method is employed to validate the upscaling results. The results show that the SA algorithm is seriously inadequate for estimating snow depth variation in space, and is only suitable for regions with relatively flat terrain and small variation of snow depth. The BME algorithm can introduce prior knowledge and perform kernel smoothing on observed data, and the upscaling result is superior to geostatistical and RF algorithms, especially when the observed data is insufficient, and outliers appear. The results of the study are expected to provide a reference for developing a point-to-surface upscaling method based on snow depth ground observations, and to further solve the uncertainties caused by scale effects in snow depth and other land surface parameter inversion and validation, by using remote sensing data.

The relationship among summer atmospheric boundary layer height over the Taklimakan Desert, its land surface parameters and Eurasian circulation

AbstractTo investigate the relationship among summer atmospheric boundary layer (ABL) over the Taklimakan Desert (TD), its land surface parameters and large scale meteorological field, the Urumqi Institute of Desert Meteorology of China Meteorological Administration conducted a comprehensive observation of boundary layer in Tazhong, the hinterland of TD from July 1 to 31, 2016. The correlations between the ABL height and surface sensible heat, latent heat, air temperature, turbulent kinetic energy (TKE) and buoyancy term are analyzed by using the observed GPS sounding data, the surface layer micro meteorological data in the hinterland of the TD and ERA5 reanalysis data so as to study the correlation of the ABL height and Eurasian circulation field. The research results indicate: (1) The height of summer ABL over the TD does not fully correspond with the variation tendency of surface sensible heat flux. In addition to surface sensible heat, large scale advection and residual mixed layer (RML) are also crucial factor influencing ABL. (2) The surface TKE over the TD is consistent with ABL height variation, with the correlation coefficient of 0.62. Hence the equation of ABL height and TKE is expressed as: Y = 1609.2X4 – 12,828X3 + 35,617X2 –39,514X + 15,858. (3) When the ABL height over the TD is quite high, 500 hPa altitude in Xinjiang is mostly controlled by anticyclone circulation, and the easterly low level jet (LLJ) tends to be formed at 850 hPa the next night. (4) The high correlation areas of summer ABL height over the TD and Eurasian temperature field are Xinjiang and Balkash Lake region. Significant positive correlation is identified between the ABL height and the temperature field at 850, 700, and 500 hPa of the above regions. The maximum correlation coefficient reaches 0.7. The study deepens the understanding of the relationship among the ABL process over the TD, land surface parameters and Eurasian circulation.

Sensitivity of METRIC-based tree crop evapotranspiration estimation to meteorology, land surface parameters and domain size

Evapotranspiration (ET), a key component of crop water balance, is essential for growers to manage agricultural irrigation. METRIC model is increasingly used to estimate and map ET with satellite observations, especially at daily and field scales. However, the analysis of its uncertainties on tree crop ET estimation is very limited. We used METRIC model to estimate daily ET over the almond and pistachio orchards throughout the entire growing season in the San Joaquin Valley, California, and evaluated the impacts of variations in meteorological variables, land surface parameters and different study sizes on ET estimation. The results showed that the METRIC-based daily ET estimates of the almond and pistachio trees were most sensitive to net radiation (Rn), followed by instantaneous reference ET (ETo), land surface temperature (LST), and near-surface temperature difference (dT). When a 5% error was added to Rn, the uncertainty in daily ET estimates was up to 8.44 ± 1.24%, corresponding to a mean increase of 0.45 ± 0.02 mm/day throughout the entire growing season. A 5% or 0.5 K error in instantaneous ETo, LST, and dT also caused changes in daily ET estimates by 4.15 ± 3.75%, 2.97 ± 0.14%, and 2.31 ± 1.20%, respectively. ET estimates with METRIC were relatively much less sensitive to other variables. The results from different study sizes indicated that daily ET estimates in Size-2 were 0.33 ± 0.08 mm/day overestimation compared with those from Size-1 in the almond site. Similarly, a mean overestimation of 0.30 ± 0.14 mm/day was found between Size-1 and Size-3 in the pistachio sites throughout the entire growing season. This study demonstrates the potential to further improve the accuracy of remote sensing based ET estimation for complex tree crops with METRIC model, by optimizing key parameters and choosing a proper domain size for study area.

Applying a Wavelet Transform Technique to Optimize General Fitting Models for SM Analysis: A Case Study in Downscaling over the Qinghai–Tibet Plateau

Soil moisture (SM) is an important land-surface parameter. Although microwave remote sensing is recognized as one of the most appropriate methods for retrieving SM, such retrievals often cannot meet the requirements of specific applications because of their coarse spatial resolution and spatiotemporal data gaps. A range of general models (GMs) for SM analysis topics (e.g., gap-filling, forecasting, and downscaling) have been introduced to address these shortcomings. This work presents a novel strategy (i.e., optimized wavelet-coupled fitting method (OWCM)) to enhance the fitting accuracy of GMs by introducing a wavelet transform (WT) technique. Four separate GMs are selected, i.e., elastic network regression, area-to-area regression kriging, random forest regression, and neural network regression. The fitting procedures are then tested within a downscaling analysis implemented between aggregated Global Land Surface Satellite products (i.e., LAI, FVC, albedo), Thermal and Reanalysis Integrating Medium-resolution Spatial-seamless LST, and Random Forest Soil Moisture (RFSM) datasets in both the WT space and the regular space. Then, eight fine-resolution SM datasets mapped from the trained GMs and OWCMs are analyzed using direct comparisons with in situ SM measurements and indirect intercomparisons between the aggregated OWCM-/GM-derived SM and RFSM. The results demonstrate that OWCM-derived SM products are generally closer to the in situ SM observations, and better capture in situ SM dynamics during the unfrozen season, compared to the corresponding GM-derived SM product, which shows fewer time changes and more stable trends. Moreover, OWCM-derived SM products represent a significant improvement over corresponding GM-derived SM products in terms of their ability to spatially and temporally match RFSM. Although spatial heterogeneity still substantially impacts the fitting accuracies of both GM and OWCM SM products, the improvements of OWCMs over GMs are significant. This improvement can likely be attributed to the fitting procedure of OWCMs implemented in the WT space, which better captures high- and low-frequency image features than in the regular space.

Estimating Maize Yield in the Black Soil Region of Northeast China Using Land Surface Data Assimilation: Integrating a Crop Model and Remote Sensing.

Accurate yield estimation at the regional scale has always been a persistent challenge in the agricultural sector. With the vigorous emergence of remote sensing land surface observations in recent decades, data assimilation methodology has become an effective means to promote the accuracy and efficiency of yield estimation by integrating regional data and point-scale crop models. This paper focuses on the black soil area of Northeast China, a national strategic grain production base, applying the AquaCrop crop growth model to simulate the fractional vegetation cover (FVC) and maize yield from 2000 to 2020 and then forming a reliable FVC optimization dataset based on an ensemble Kalman filter (EnKF) assimilation algorithm with remote sensing products. Using the random forest model, the regression relationship between FVC and yield was established from the long-term time series data, which is crucial to achieve better yield estimation through the optimized FVC. The major findings include the following: (1) The R2 of the assimilated FVC and maize yield can reach 0.557. (2) When compared with the local statistical yield, our method reduced the mean absolute error (MAE) from 1.164 ton/ha (based on GLASS FVC products) to 1.004 ton/ha (based on the calibrated AquaCrop model) and then to 0.888 ton/ha (the result after assimilation). The above results show that we have proposed a yield estimation method to provide accurate yield estimations by combining data assimilation and machine learning. This study provided deep insights into understanding the variations in FVC and revealed the spatially explicit yield prediction ability from the time series land surface parameters, which has significant potential for optimizing water and soil resource management.

Fine scale variability in Green Vegetation Fraction Over the Southern Great Plains using Sentinel-2 satellite: A case study

We estimate the green vegetation fraction (GVF) at high spatial resolution (10 m) by applying a new method based on image histogram equalization and utilizing images with multiple spectral bands from Sentinel-2 satellite observations. True color and multi-spectral images at 10 m resolution are used over 50 × 50 km2 domain centered on the US Department of Energy (DOE) Atmospheric Radiation Measurement (ARM) site over the Southern Great Plains to recognize green pixels and estimate the vegetation fraction. The GVF retrieved from several spectral bands is compared with Sentinel-2 true color images, Sentinel-3, MODIS vegetation products, and observations from the US DOE research aircraft on selected days from the Holistic Interactions of Shallow Clouds, Aerosols, and Land-Ecosystems (HI-SCALE) campaign. Unlike MODIS or Sentinel-3, our approach to retrieve the GVF is independent of absolute reflectance or photosynthetic response. Spatial patterns in vegetation distribution are similar in all three satellite approaches, however MODIS underestimates the vegetation cover by about 10% compared to the Sentinel-2, -3 retrievals over the study domain. The largest contrast in MODIS and Sentinel-2, -3 is found when the GVF is less than 0.2 or greater than 0.9, likely as a result of different retrieval approaches and spectral bands used to estimate the vegetation products. The NDVI retrieved from MODIS is also underestimated compared to the aircraft observations suggesting that MODIS underestimates both NDVI and GVF. In the past, various satellite observations have been used to retrieve the vegetation fraction at medium (0.25–0.50 km) and coarser resolutions (15 km) but our approach is extended to a much higher (10 m) resolution. Our analysis from Sentinel-2 shows large spatial variability in GVF and a clear disparity with estimates from MODIS and Sentinel-3 satellites. The new approach to estimate GVF has the potential to increase spatial accuracy and precision of surface properties relevant for crop and forest management, surface and energy balance processes, and land surface parameterizations in atmospheric models.

Global estimation of above-ground biomass from spaceborne C-band scatterometer observations aided by LiDAR metrics of vegetation structure

The backscattered power recorded by a spaceborne scatterometer operating at C-band is sensitive to land surface parameters and is operationally used by some global remote sensing services, e.g., to estimate soil moisture. The estimation of forest variables, in particular above-ground biomass (AGB), from scatterometer data instead was seldom explored. Given the availability of multi-decadal sets of scatterometer observations from space, it is of interest to address the contribution of C-band scatterometer data to the quantification of carbon stocks stored in forests even if the spatial resolution of spaceborne scatterometers is very coarse. In this paper, we investigated the prospects of AGB estimation using backscatter observations by the MetOp Advanced SCATterometer (ASCAT) with a spatial resolution of 0.25°. For this study, ASCAT observations acquired in 2010 were used to be contemporary with AGB datasets selected to benchmark the performance of the estimation. A Water Cloud Model that integrates two allometric equations derived from spaceborne LiDAR data reproduced the relationship between observations of radar backscatter as a function of AGB. Estimates of AGB from individual observations were then combined with a weighted average to reduce uncertainties. Finally, a correction was introduced to compensate for the offset introduced by sloped terrain and surfaces not covered by woody vegetation on the AGB estimate of a pixel. Uncertainties associated with the scatterometer observations, and the modelling framework were propagated to obtain per-pixel values of the standard deviation of an AGB estimate. The proposed method explains much of the variance in AGB estimates when compared to measurements from inventory data (R2 = 0.72) and generated unbiased estimates globally (bias: −3.3 Mg⋅ha−1). Nonetheless, the discrepancy between estimated and plot-based AGB values tended to increase for decreasing biomass level from 20% to 60% of the reference AGB level. A further assessment related to global stocks indicated that the value estimated from the scatterometer dataset (596 Pg, 95% of which 563 Pg stored in forest land) was in line with two published estimates based on forest inventory data only (571 Pg and 600 Pg, respectively). Despite the coarse spatial resolution, our results indicate that C-band scatterometer observations from space can contribute to the characterization of terrestrial biomass pools. The record of observations starting in the early 1990s may provide an unprecedented way to look at long-term forest dynamics as well as to constrain the strength of carbon-climate cycle feedback simulated by Earth System models.

The role of land surface schemes in non‐hydrostatic RegCM on the simulation of Indian summer monsoon

AbstractThis study evaluates the performance of a non‐hydrostatic Regional Climate Model (RegCM) with two land surface parameterizations, namely Biosphere‐Atmosphere Transfer (BAT) scheme and Community Land Model (CLM) scheme on the simulation of summer monsoon over India. The initial and boundary conditions are taken from ERA‐Interim (ERA‐I) reanalysis available at a 0.75° spatial grid. The RegCM is designed with 25 km horizontal resolution and 23 vertical levels, and it integrated for the period 1st January 1982 to 30th September 2016. We used the NOAA Optimum Interpolation Weekly sea surface temperatures, the topography and land‐use data from the United States Geological Survey, and Global Land Cover Characterization. Our results reveal that the summer monsoon precipitation anomaly is reasonably well simulated with BAT and CLM schemes compared to those of India Meteorological Department (IMD) observations. The spatial distributions of temperature anomaly over India with BAT compared to CRU shows cold bias over most of the Indian sub‐continent, whereas the CLM shows large warm bias over central India and monsoon trough region. The low‐level and upper‐level winds are in good agreement with ERA‐I. Further, the analysis of extreme monsoon events (i.e., excess and deficit seasons) reveals that CLM depicts less bias in simulating the excess and deficit precipitation regions as compared to those of BAT. Conversely, BAT simulates temperature well over India in both excess and deficit monsoon seasons. The Taylor's diagram analysis, Added Value Index based on Mean Square Error, and Modified Brier Skill Score are also applied for the extreme rainfall seasons. Our quantitative analysis demonstrates that the CLM have resolved better rainfall features with more skills with respect to that of BAT. The results suggest that CLM have added value in Indian Summer Monsoon (ISM) regional climate simulation than the BAT, particularly on simulation of the rainfall characteristics during extreme rainfall seasons.

Combination of Sentinel-1 and Sentinel-2 Data for Tree Species Classification in a Central European Biosphere Reserve

Microwave and optical imaging methods react differently to different land surface parameters and, thus, provide highly complementary information. However, the contribution of individual features from these two domains of the electromagnetic spectrum for tree species classification is still unclear. For large-scale forest assessments, it is moreover important to better understand the domain-specific limitations of the two sensor families, such as the impact of cloudiness and low signal-to-noise-ratio, respectively. In this study, seven deciduous and five coniferous tree species of the Austrian Biosphere Reserve Wienerwald (105,000 ha) were classified using Breiman’s random forest classifier, labeled with help of forest enterprise data. In nine test cases, variations of Sentinel-1 and Sentinel-2 imagery were passed to the classifier to evaluate their respective contributions. By solely using a high number of Sentinel-2 scenes well spread over the growing season, an overall accuracy of 83.2% was achieved. With ample Sentinel-2 scenes available, the additional use of Sentinel-1 data improved the results by 0.5 percentage points. This changed when only a single Sentinel-2 scene was supposedly available. In this case, the full set of Sentinel-1-derived features increased the overall accuracy on average by 4.7 percentage points. The same level of accuracy could be obtained using three Sentinel-2 scenes spread over the vegetation period. On the other hand, the sole use of Sentinel-1 including phenological indicators and additional features derived from the time series did not yield satisfactory overall classification accuracies (55.7%), as only coniferous species were well separated.

Impacts of Land Surface Parameterizations on Simulations over the Arid and Semiarid Regions: The Case of the Loess Plateau in China

Abstract The arid and semiarid areas of the Loess Plateau are extremely sensitive to climate change. Land–atmosphere interactions of these regions play an important role in the regional climate. However, most present land surface models (LSMs) are not reasonable and accurate enough to describe the surface characteristics in these regions. In this study, we investigate the effects of three key land surface parameters including surface albedo, soil thermal conductivity, and additional damping on the Noah LSM in simulating the land surface characteristics. The observational data from June to September from 2007 to 2009 collected at the Semi-Arid Climate and Environment Observatory of Lanzhou University (SACOL) station in northwest China are used to validate the Noah LSM simulations. The results suggest that the retrieved values of surface albedo, soil thermal conductivity, and additional damping based on observations are in closer agreement with those of the MULT scheme for surface albedo, the J75_NOAH scheme for soil thermal conductivity, and the Y08 scheme for additional damping, respectively. Furthermore, the model performance is not obviously affected by surface albedo parameterization schemes, while the scheme of soil thermal conductivity is vital to simulations of latent heat flux and soil temperature and the scheme of additional damping is crucial for simulating net radiation flux, sensible heat flux, and surface soil temperature. A set of optimal parameterizations is proposed for the offline Noah LSM at the SACOL station when the MULT scheme for surface albedo, the J75_NOAH scheme for soil thermal conductivity, and the Y08 scheme for additional damping are combined simultaneously, especially in the case of sensible heat flux and surface soil temperature simulations.

Quasi-global machine learning-based soil moisture estimates at high spatio-temporal scales using CYGNSS and SMAP observations

Global soil moisture mapping at high spatial and temporal resolution is important for various meteorological, hydrological, and agricultural applications. Recent research shows that the land surface reflection in the forward direction of Global Navigation Satellite System (GNSS) signals at L-band can convey high-resolution land surface information, including surface soil moisture. However, these signals are often affected by complex land surface characteristics and the bistatic nature of the GNSS-Reflectometry (GNSS-R) technique, resulting in a nonlinear relationship between the signals and surface soil moisture. In this work, a machine learning (ML) approach is used to map quasi-global soil moisture using bistatic reflectance observations acquired from the recently launched Cyclone GNSS (CYGNSS) mission. Specifically, several land surface parameters are obtained from remote sensing products and integrated with Soil Moisture Active Passive (SMAP) enhanced soil moisture retrievals to facilitate daily quasi-global CYGNSS soil moisture mapping at 9 km. Based on cross-validation against SMAP data, the ML algorithm is shown to be suitable for retrieving soil moisture from CYGNSS. Median values of unbiased root-mean-square-difference for the quasi-global coverage or regions with vegetation water content less than 5 kg/m2 are 0.0395 cm3/cm3 and 0.0320 cm3/cm3, respectively. Likewise, via independent evaluation against more than 100 in-situ sites, the algorithm is shown to have an unbiased root-mean-square-error of 0.0543 cm3/cm3. CYGNSS-based retrievals contain similar spatial variability as SMAP across different seasons. Moreover, through a robust triple collocation technique, the accuracy of CYGNSS soil moisture is relatively high over moderately vegetated regions with correlations ranging from 0.4 to 0.8. Based on these validation results, we argue that derived CYGNSS soil moisture estimates can supplement current global soil moisture databases and provide more frequent retrievals at 9 km.

A global dataset on subgrid land surface climate (2015–2100) from the Community Earth System Model

AbstractSubgrid data from earth system models are a powerful, yet underutilized, data resource for investigating the climatic impacts of land use and land cover change (LULCC). In this paper, we describe a global dataset on subgrid land surface climate variables produced by the Community Earth System Model in a fully coupled mode. The simulation was conducted at a 0.9° × 1.25° resolution under the Representative Concentration Pathway (RCP) 8.5 scenario from 2015 to 2100. Data are archived for eight subgrid tiles (urban, rural, tree, grass, shrub, bare soil, crop and lake) and include variables on the physical state, surface energy fluxes, runoff and atmospheric forcing conditions. Archival intervals are monthly, daily and hourly. Meta data on land surface parameters are also available. The data files are stored in NetCDF‐4 (Network Common Data Form, version 4) format and the meta data follow the latest Coupled Model Intercomparison Project phase 6 standards. We anticipate that this dataset will become a useful resource for characterizing local climate changes due to LULCC. This dataset can be downloaded from the Harvard Dataverse (https://doi.org/10.7910/DVN/HUXAH6).

Shrub willow canopy chlorophyll content estimation from unmanned aerial systems (UAS) data: Estimation and uncertainty analysis across time, space, and scales

Remote sensing is an important technology to map land-surface parameters, with many studies demonstrating land-surface parameter estimation, but fewer testing model transferability and quantifying model parameter uncertainty. In this study, we evaluated the uncertainty across time, space, and spatial scales of shrub willow health characterization based on canopy chlorophyll content (CCC) derived from unmanned aerial systems (UAS) data. Shrub willow is a short-rotation woody crop used to produce biomass for bioproducts and CCC is a popular biophysical parameter used to indicate shrub willow health status. We analyzed time-series UAS data at three spatial scales (5 m, 10 m, and 20 m) and stratified field-observed CCC levels. Since scale effects are often related to spatial heterogeneity, we implemented a nested analysis of variance (ANOVA) to evaluate the spatial heterogeneity within different pixel sizes and found that 5 m pixels were the most homogeneous, followed by 10 m and 20 m. Results from regression modeling of shrub willow CCC as a function of red-edge normalized difference vegetation index (NDVIre) at 5 m, 10 m, and 20 m scales showed that the models built at 5 m, 10 m, and 20 m could be applied across time, space, and scales. We also quantified the uncertainty for model parameters using two different inferential frameworks, confidence interval (CI) widths from the frequentist framework and credible interval (CrI) widths from the Bayesian framework. Unlike predictive root mean square error (RMSE), which showed model output uncertainty decreased as pixel size increased, CI and CrI widths indicated that the related model parameter uncertainty increased as pixel size increased. We calculated CI and CrI widths based on time-series model building and different stratified model building to analyze model parameter uncertainty and quantified predictive RMSE for all sampling date combinations. The results showed that, in terms of model complexity and the range of ground observations, both frequentist and Bayesian regression have advantages and disadvantages and demonstrated that the uncertainty quantified by CrI widths is able to guide future experimental design to save resources.

Data-Driven Selection of Land Product Validation Station Based on Machine Learning

Validation is a crucial technique used to strengthen the application capabilities of earthobservation satellite data and solve the quality problems of remote-sensing products. Observing land-surface parameters in the field is one of the key steps of validation. Therefore, the demand for long-term stable validation stations has gradually increased. However, the current location-selection procedure of validation stations lacks a systematic and objective evaluation system. In this research, a data-driven selection of a land product validation station (DSS-LPV) based on Machine Learning is proposed. Firstly, we construct an evaluation indicator system in which all factors affecting the location of validation stations are divided into surface characteristics, atmospheric conditions and the social environment. Then, multi-scale evaluation grids are constructed and indicators are allocated for spatial evaluation. Finally, four Machine Learning (ML) methods are used to learn the established reliable stations, and different data-driven scoring models are constructed to explore the intrinsic relationship between evaluation indicators and station locations. In this article, the reliability of DSS-LPV is effectively validated by the example of China using the national-level land product validation station that has been established. After a comparison between the four ML models, the random forest (RF) with the highest accuracy was selected as the modeling method of DSS-LPV. The correlation between the regression value of test stations and the target value is 0.9133. The average score of test stations is 0.8304. The test stations are generally located within the calculated hot-spot area of the score density map, which means that it is highly consistent with the location of the built stations. Research results indicate that DSS-LPV is an effective method that can provide a reasonable geographical distribution of the stations. The location-selection results can provide scientific decision-making support for the construction of land product validation stations.

Land Surface Parameterization at Exposed Playa and Desert Region to Support Dust Emissions Estimates in Southern California, United States

Remote sensing technologies provide a unique opportunity to identify ground surfaces that are more susceptible to dust emissions at a large scale. As part of the Salton Sea Air Quality Mitigation Program (SSAQMP) of the Imperial Irrigation District (IID), efforts have been made to improve our understanding of fugitive, wind-blown dust emissions around the Salton Sea region in Southern California, United States. Field campaigns were conducted for multiple years to evaluate surface conditions and measure the dust emissions potential in the area. Data collected during the field work were coupled with remote sensing imagery and data mining techniques to map surface characteristics that are important in identifying dust emissions potential. Around the playa domain, surface crust type, sand presence, and soil moisture were estimated. Geomorphic surface types were mapped in the desert domain. Overall accuracy ranged from 91.7% to 99.4% for the crust type mapping. Sand presence mapping showed consistent and slightly better accuracy, ranging from 96.2% to 99.7%. Soil moisture assessment agreed with precipitation records. Geomorphic mapping in the desert domain achieved accuracy above 93.5%, and the spatial pattern was consistent with previous studies. These land surface condition assessments provide important information to support dust emissions estimates in the region.

Estimation of Soil Salinization by Machine Learning Algorithms in Different Arid Regions of Northwest China

Hyperspectral data has attracted considerable attention in recent years due to its high accuracy in monitoring soil salinization. At present, most existing research focuses on the saline soil in a single area without comparative analysis between regions. The regional differences in the hyperspectral characteristics of saline soil are still unclear. Thus, we chose Golmud in the cold–dry Qaidam Basin (QB–G) and Gaotai–Minghua in the relatively warm–dry Hexi Corridor (HC–GM) as the study areas, and used the deep extreme learning machine (DELM) and sine cosine algorithm–Elman (SCA–Elman) to predict soil salinity, and then selected the most suitable algorithm in these two regions. A total of 79 (QB–G) and 86 (HC–GM) soil samples were collected and tested to obtain their electrical conductivity (EC) and corresponding hyperspectral reflectance (R). We utilized the land surface parameters that affect the soil based on Landsat 8 and digital elevation model (DEM) data, selected the variables using the light gradient boosting machine (LightGBM), and built SCA–Elman and DELM from the hyperspectral reflectance data combined with land surface parameters. The results revealed the following: (1) The soil hyperspectral reflectance in QB–G was higher than that in HC–GM. The soils of QB–G are mainly the chloride type and those of HC–GM mainly belong to the sulfate type, having lower reflectance. (2) The accuracies of some of the SCA–Elman and DELM models in QB–G (the highest MAEv, RMSEv, and Rv2 were 0.09, 0.12 and 0.75, respectively) were higher than those in HC–GM (the highest MAEv, RMSEv, and Rv2 were 0.10, 0.14 and 0.73, respectively), which has flatter terrain and less obvious surface changes. The surface parameters in QB–G had higher correlation coefficients with EC due to the regular altitude change and cold–dry climate. (3) Most of the SCA–Elman results (the mean Rv2 in HC-GM and QB-G were 0.62 and 0.60, respectively) in all areas performed better than the DELM results (the mean Rv2 in HC–GM and QB–G were 0.51 and 0.49, respectively). Therefore, SCA–Elman was more suitable for the soil salinity prediction in HC–GM and QB–G. This can provide a reference for soil salinization monitoring and model selection in the future.