SEARCH · Search NASA

Results for “spatial variability”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 883 records · Page 49

All sources NASA NTRS DOE OSTI

DOE OSTI2024

Development of the tangent linear and adjoint models of the global online chemical transport model MPAS-CO 2 v7.3

We describe the development of the tangent linear (TL) and adjoint models of the Model for Prediction Across Scales (MPAS)-CO 2 transport model, which is a global online chemical transport model developed upon the non-hydrostatic Model for Prediction Across Scales – Atmosphere (MPAS-A). The primary goal is to make the model system a valuable research tool for investigating atmospheric carbon transport and inverse modeling. First, we develop the TL code, encompassing all CO 2 transport processes within the MPAS-CO 2 forward model. Then, we construct the adjoint model using a combined strategy involving re-calculation and storage of the essential meteorological variables needed for CO 2 transport. This strategy allows the adjoint model to undertake a long-period integration with moderate memory demands. To ensure accuracy, the TL and adjoint models undergo vigorous verifications through a series of standard tests. The adjoint model, through backward-in-time integration, calculates the sensitivity of atmospheric CO 2 observations to surface CO 2 fluxes and the initial atmospheric CO 2 mixing ratio. To demonstrate the utility of the newly developed adjoint model, we conduct simulations for two types of atmospheric CO 2 observations, namely the tower-based in situ CO 2 mixing ratio and satellite-derived column-averaged CO 2 mixing ratio (X CO 2 ). A comparison between the sensitivity to surface flux calculated by the MPAS-CO 2 adjoint model with its counterpart from CarbonTracker–Lagrange (CT-L) reveals a spatial agreement but notable magnitude differences. These differences, particularly evident for X CO 2 , might be attributed to the two model systems' differences in the simulation configuration, spatial resolution, and treatment of vertical mixing processes. Moreover, this comparison highlights the substantial loss of information in the atmospheric CO 2 observations due to CT-L's spatial domain limitation. Furthermore, the adjoint sensitivity analysis demonstrates that the sensitivities to both surface flux and initial CO 2 conditions spread out throughout the entire Northern Hemisphere within a month. MPAS-CO 2 forward, TL, and adjoint models stand out for their calculation efficiency and variable-resolution capability, making them competitive in computational cost. In conclusion, the successful development of the MPAS-CO 2 TL and adjoint models, and their integration into the MPAS-CO 2 system, establish the possibility of using MPAS's unique features in atmospheric CO 2 transport sensitivity studies and in inverse modeling with advanced methods such as variational data assimilation.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2020

How MJO Teleconnections and ENSO Interference Impacts U.S. Precipitation

A composite analysis reveals how the Madden–Julian oscillation (MJO) impacts North American rainfall through perturbations in both the upper-tropospheric flow and regional low-level moisture availability. Upper-level divergence associated with the MJO tropical convection drives a quasi-stationary Rossby wave response to the midlatitudes. This forces a midlatitude upper-level dipolar geopotential height anomaly that is accompanied by a westward retraction of the jet stream and reduced rainfall over the central-eastern North Pacific. A reverse effect is found as the MJO propagates eastward across the Maritime Continent. These large differences in the extratropical upper-level flow, combined with anomalies in the regional supply of water vapor, have a profound impact on southeastern U.S. rainfall. The low-frequency variability, including that associated with ENSO, can modify the seasonal background flow (e.g., El Niño and La Niña basic states) affecting the distribution, strength, and propagation of the intraseasonal oscillation and the extratropical teleconnection patterns. The combined effects of the ENSO and the MJO signals result in both spatial and temporal patterns of interference and modulation of North American rainfall. The results from this study show that during a particular phase of an active MJO, the extratropical response can considerably enhance or mask the interannual ENSO signal in the United States, potentially resulting in anomalies of the opposite sign than that expected during a specific ENSO phase. Analyses of specific MJO events during an El Niño or La Niña episode reveal significant contributions to extreme events via constructive and destructive interference of the MJO and ENSO signals.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Improving the accuracy of discretisations of the vector transport equation on the lowest-order quadrilateral Raviart-Thomas finite elements

Within finite element models of fluids, vector-valued fields such as velocity or momentum variables are commonly discretised using the Raviart-Thomas elements. However, when using the lowest-order quadrilateral Raviart-Thomas elements, standard finite element discretisations of the vector transport equation typically have a low order of spatial accuracy. This paper describes two schemes that improve the accuracy of transporting such vector-valued fields on two-dimensional curved manifolds. The first scheme that is presented reconstructs the transported field in a higher-order function space, where the transport equation is then solved. The second scheme applies a mixed finite element formulation to the vector transport equation, simultaneously solving for the transported field and its vorticity. In this work, an approach to stabilising this mixed vector-vorticity formulation is presented that uses a Streamline Upwind Petrov-Galerkin (SUPG) method. These schemes are then demonstrated, along with their accuracy properties, through some numerical tests. Two new test cases are used to assess the transport of vector-valued fields on curved manifolds, solving the vector transport equation in isolation. The improvement of the schemes is also shown through two standard test cases for rotating shallow-water models.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2022

Resolve the Clear-Sky Continuous Diurnal Cycle of High-Resolution ECOSTRESS Evapotranspiration and Land Surface Temperature

Although numerous evapotranspiration (ET) estimates have been developed from remote sensing measurements, a diurnally resolved product with high spatial resolution (<100 m) is still lacking, which is critically needed for agricultural and ecological monitoring. The ECOsystem Spaceborne Thermal Radiometer Experiment on Space Station (ECOSTRESS) provides, for the first time, ET at 70 m resolution with diurnal samplings, enabling a step toward filling this data gap. This study developed the first framework to resolve the full diurnal cycle of ET and Land Surface Temperature (LST) from sporadic ECOSTRESS measurements. We first constructed 70 m diurnal LST utilizing a diurnal temperature cycle (DTC) model that fused ECOSTRESS and Geostationary Operational Environmental Satellite (GOES) LST. Next, we derived 70 m diurnal ET from the diurnally resolved LST, along with ancillary meteorological and reflectance data sets, using the Priestley-Taylor Jet Propulsion Laboratory (PT-JPL) algorithm. This framework was tested in two regions: A relatively homogeneous semiarid grassland landscape (centered at the US-Seg Ameriflux site) and a heterogeneous humid cropland landscape (centered at US-ARM). Our diurnally resolved LST and ET successfully reproduced the spatial patterns in the native ECOSTRESS measurements during overpasses for both landscapes (LST: r > 0.97, RMSE ~2°C; ET: r > 0.99, and RMSE = 1–7 W m -2 ). The constructed time series also well captured the temporal variability of the in situ measurements (especially under clear-sky conditions) at both US-Seg and US-ARM, albeit performance varying with availability and quality of ECOSTRESS samples. This study sets the stage for testing and applying our framework to broader climates, biomes, and landscapes toward eventually generating diurnally resolved 70 m global operational LST/ET products to enhance ecological/agricultural applications.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

A combined first principles study of the structural, magnetic, and phonon properties of monolayer CrI3

The first magnetic 2D material discovered, monolayer (ML) CrI3, is particularly fascinating due to its ground state ferromagnetism. However, because ML materials are difficult to probe experimentally, much remains unresolved about ML CrI3’s structural, electronic, and magnetic properties. Here, we leverage Density Functional Theory (DFT) and high-accuracy Diffusion Monte Carlo (DMC) simulations to predict lattice parameters, magnetic moments, and spin–phonon and spin–lattice coupling of ML CrI3. We exploit a recently developed surrogate Hessian DMC line search technique to determine CrI3’s ML geometry with DMC accuracy, yielding lattice parameters in good agreement with recently published STM measurements—an accomplishment given the ~10% variability in previous DFT-derived estimates depending upon the functional. Strikingly, we find that previous DFT predictions of ML CrI3’s magnetic spin moments are correct on average across a unit cell but miss critical local spatial fluctuations in the spin density revealed by more accurate DMC. DMC predicts that magnetic moments in ML CrI3 are 3.62 μB per chromium and -0.145 μB per iodine, both larger than previous DFT predictions. The large disparate moments together with the large spin–orbit coupling of CrI3’s I-p orbital suggest a ligand superexchange-dominated magnetic anisotropy in ML CrI3, corroborating recent observations of magnons in its 2D limit. We also find that ML CrI3 exhibits a substantial spin–phonon coupling of ~3.32 cm-1. Our work, thus, establishes many of ML CrI3’s key properties, while also continuing to demonstrate the pivotal role that DMC can assume in the study of magnetic and other 2D materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

VA Determinants of Health Data Curation Documentation FY25-Q2

The U.S. Department of Veterans Affairs (VA) places the health and well-being of our nation’s veterans as its top priority. VA is dedicated to offering timely access to high-quality, evidence-based mental health care that meets the needs of veterans and supports their reintegration into society. One of our core missions is to prevent suicide among veterans through innovative approaches and resources. With funding from the VA Office of Mental Health and Suicide Prevention (OMHSP), the Determinants of Health (EDH) project has developed innovative datasets associated with specific health outcomes, a methodology for transforming spatiotemporal data from one spatial reference (e.g., a 1km grid) to another (e.g., US Census Tracts), and capabilities for modeling health outcomes. These datasets represent an enhancement of the Agency for Healthcare Research and Quality (AHRQ), addressing key gaps by introducing finer spatial resolution (Census Tract) and additional geographical covariates into existing data. The curation and standardization of these datasets is a complex task since they often originate from various sources and are measured at different spatial and temporal resolutions. For example, US Census data products typically use census blocks, block groups, or counties, while data like weather data are available on 1km grids. Some economic data may only be available at the zip code level. In this context, ‘standardized’ means that all datasets share the same spatial extent (e.g., US Census Tract and/or County), and ‘curated’ implies a repeatable process with data provenance and the use of appropriate methodologies for covariate conversion. The Determinants of Health datasets draw from multiple sources, resulting in variables with varying degrees of availability, patterns of missing data, and methodological considerations across different sources, geographies, and years.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2025

VA Determinants of Health Data Curation Documentation FY25-Q3

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2025

VA Community Determinants of Health Data Curation Documentation FY25-Q4

The U.S. Department of Veterans Affairs (VA) places the health and well-being of our nation’s veterans as its top priority. VA is dedicated to offering timely access to high-quality, evidence-based mental health care that meets the needs of veterans and supports their reintegration into society. One of our core missions is to prevent suicide among veterans through innovative approaches and resources. With funding from the VA Office of Mental Health and Suicide Prevention (OMHSP), the Determinants of Health (EDH) project has developed innovative datasets associated with specific health outcomes, a methodology for transforming spatiotemporal data from one spatial reference (e.g., a 1km grid) to another (e.g., US Census Tracts), and capabilities for modeling health outcomes. These datasets represent an enhancement of the Agency for Healthcare Research and Quality (AHRQ), addressing key gaps by introducing finer spatial resolution (Census Tract) and additional geographical covariates into existing data. The curation and standardization of these datasets is a complex task since they often originate from various sources and are measured at different spatial and temporal resolutions. For example, US Census data products typically use census blocks, block groups, or counties, while data like weather data are available on 1km grids. Some economic data may only be available at the zip code level. In this context, standardized’ means that all datasets share the same spatial extent (e.g., US Census Tract and/or County), and ‘curated’ implies a repeatable process with data provenance and the use of appropriate methodologies for covariate conversion. The Determinants of Health datasets draw from multiple sources, resulting in variables with varying degrees of availability, patterns of missing data, and methodological considerations across different sources, geographies, and years.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2022

Resistivity of Surface Steps in Bulk-Insulating Topological Insulators

Electron transport in topological insulators usually involves both topologically protected surface states and trivial electronic states in the bulk material. The surface transport is particularly interesting; however, it is also susceptible to atomic defects on the surfaces, such as vacancies, impurities, and step edges. Experimental determination of scattering effects of these surface defects requires both nanoscale spatial resolution and the ability to decipher surface transport from bulk transport. Here we directly measure the resistivity of individual surface steps in the surface dominating transport process of topological insulator Bi 2 Te 2 Se. A variable probe-spacing transport spectroscopy with a multiprobe scanning tunneling microscope is used to differentiate the surface conductance from bulk conductance, allowing the identification of a surface dominant transport regime. The technique also reveals a deviation from ideal 2D transport at atomic steps. Then, a multi-probe scanning tunneling potentiometry is employed to visualize the electrochemical potentials across individual step edges. A quantitative analysis of the potential distributions enables us to acquire a resistivity of 0.530 mΩ · cm for the one quintuple-layer atomic step. The result indicates that atomic defects, despite preserving the time-reversal symmetry, can still significantly affect the transport in topological insulators.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Evaluating Entrainment–Mixing Characteristics through Direct Comparisons of Drop Size Distributions Using In Situ Observations from ACE-ENA

Abstract Constraining the impacts of entrainment and associated mixing (i.e., entrainment–mixing) on cloud properties continues to be difficult, partly due to observational uncertainties as well as a lacking consensus of which methodologies for diagnosing entrainment–mixing are most appropriate. This study introduces a novel method to evaluate the presence and degree of inhomogeneous and homogeneous mixing using ∼100 h of in situ observations from a research aircraft over the northeastern Atlantic. Specifically, drop size distributions are compared between regions containing negligible and significant entrainment for select flight legs, making a direct characterization of the degree of homogeneous and inhomogeneous mixing possible. A measure of drop concentration variance is used as a proxy variable to diagnose entrainment–mixing. Results correspond well with entrainment–mixing metrics, showing lower Damköhler numbers where drop size distributions shift toward smaller drop sizes (i.e., inhomogeneous mixing) and greater transition length scales where drop size distributions do not (i.e., homogeneous mixing). Inhomogeneous mixing occurs in most samples from Aerosol and Cloud Experiment in the Eastern North Atlantic (ACE-ENA) (regardless of homogeneous mixing frequencies increasing with increasing spatial resolution from ∼100 to ∼10 m) and is associated with decreased drop size relative dispersion and both greater aerosol and drop concentrations compared with homogeneous mixing. Precipitating clouds have a greater frequency of homogeneous mixing compared with nonprecipitating clouds. The proposed methodology is similarly applied to in situ observations of southeast Pacific stratocumulus, shallow convective clouds over central Oklahoma and low-level clouds over the Southern Ocean. All four locations are primarily dominated by inhomogeneous mixing with minimal variability among each region.

Clouds↗

DOE OSTI2025

Challenges in quantifying unparameterized spatial uncertainties in deep geologic repositories for nuclear waste

Spatially heterogeneous uncertainties are prevalent in geophysical modeling applications, such as probabilistic post-closure performance assessment (PA) of deep geologic repositories for nuclear waste. Such uncertainties are often highly influential to model outputs, so it is desirable to identify the most important mechanisms by which they influence model predictions. However, these uncertainties are often unparameterized in the sense that there is no set of parameters that can be specified to yield a particular realization of the uncertainty. Additionally, the uncertainty is not intrinsically endowed with a parameterization that captures a realization’s mechanistic influence on model outputs. Therefore, in this work we present a novel methodology to develop and assess a set of proxy variables that aim to represent this influence. We show how they can be computed, downselected, and incorporated into the construction of statistical surrogate models mapping model inputs to outputs. We present our methodology in the context of a motivating application problem in deep geologic repository PA and discuss the challenges in capturing the effects of these spatial heterogeneities in uncertainty analyses.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2023

The drivers and predictability of wildfire re-burns in the western United States (US)

Evidence is mounting that the effectiveness of using prescribed burns as a management tactic may be diminishing due to the higher incidence of wildfire re-burns. The development of predictive models of re-burns is thus essential to better understand their primary drivers so that forest management practices can be updated to account for these events. First, we assess the potential for human activity as a driver of re-burns by evaluating re-burn trends both within and outside of the wildland–urban interface (WUI) of the western US. Next, we investigate the predictability of re-burns through the application of both random forest and the explanatory machine learning non-negative matrix factorization using k-means clustering (NMFk) algorithms to predict re-burn occurrence over California based on a number of climate factors. Our findings indicate that while most states showed increasing trends within the WUI when trends were conducted over longer moving windows (e.g. 20 years), California was the only state where the rate of increase was consistently higher in the WUI, indicating a stronger potential for human activity as a driver in that location. Furthermore, we find model performance was found to be robust over most of California (Testing F1 scores = 0.688), although results were highly variable based on EPA level III Ecoregion (F1 scores = 0.0–0.778). Insights provided from this study will lead to a better understanding of climate and human activity drivers of re-burns and how these vary at broad spatial scales so that improvements in forest management practices can be tuned according to the level of change that is expected for a given region.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Determining spatially varying profit‐maximizing management practices for miscanthus and switchgrass production in the rainfed United States

Abstract Determining optimal management practices for the profitable production of perennial energy crops is critical for scaling up production beyond experimental levels. Although many experimental field studies have examined the effects of management practices on the performance of miscanthus and switchgrass, there are no recommendations for economically optimal nitrogen (N) application rates and how they should vary spatially and with the age of the energy crop as well as on optimal rotation age of the energy crop to maximize profits. We develop a modeling framework to determine economically optimal crop management decisions and simulate the variability under various scenarios for miscanthus and switchgrass production across 2287 counties in the rainfed United States. We find that profit‐maximizing N recommendations for these crops vary across maturity stages and regions and can increase the landowner's profits compared with a uniform N rate across ages and regions. We also find that the optimal rotation for these crops is shorter than the productive physical lifespan (15–20 and 10 years for miscanthus and switchgrass, respectively). Specifically, the N rate that maximizes the economic returns is negligible for miscanthus and 111 kg ha −1 for switchgrass production at age 2. The mean profit‐maximizing N rate increases with age for miscanthus, peaking at 151 kg ha −1 at age 11 before declining to 114 kg ha −1 at the optimal rotation age of 13 years while that for switchgrass is 150 kg ha −1 for middle‐aged stands and declines to 114 kg ha −1 at the optimal rotation of 8–9 years. We find that miscanthus is the most profitable energy crop in the northern region of the rainfed United States while switchgrass is most profitable in the south of the rainfed United States. Our findings are useful for improving assessments of the profitability of energy crops and guiding future management decisions by landowners.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

VA Community Determinants of Health Data Curation Documentation FY26-Q1

The U.S. Department of Veterans Affairs (VA) places the health and well-being of our nation’s veterans as its top priority. VA is dedicated to offering timely access to high-quality, evidence-based mental health care that meets the needs of veterans and supports their reintegration into society. One of our core missions is to prevent suicide among veterans through innovative approaches and resources. With funding from the VA Office of Mental Health and Suicide Prevention (OMHSP), the Community Determinants of Health (EDH) Data project has developed innovative datasets associated with specific health outcomes, a methodology for transforming spatiotemporal data from one spatial reference (e.g., a 1km grid) to another (e.g., US Census Tracts), and capabilities for modeling health outcomes. These datasets represent an enhancement of the Agency for Healthcare Research and Quality (AHRQ), addressing key gaps by introducing finer spatial resolution (Census Tract) and additional geographical covariates into existing data. The curation and standardization of these datasets is a complex task since they often originate from various sources and are measured at different spatial and temporal resolutions. For example, US Census data products typically use census blocks, block groups, or counties, while data like weather data are available on 1km grids. Some economic data may only be available at the zip code level. In this context, ‘standardized’ means that all datasets share the same spatial extent (e.g., US Census Tract and/or County), and ‘curated’ implies a repeatable process with data provenance and the use of appropriate methodologies for covariate conversion. The Community Determinants of Health datasets draw from multiple sources, resulting in variables with varying degrees of availability, patterns of missing data, and methodological considerations across different sources, geographies, and years.

99 GENERAL AND MISCELLANEOUS↗

DOE OSTI2026

VA Community Determinants of Health Data Curation Documentation FY26-Q2

The U.S. Department of Veterans Affairs (VA) places the health and well-being of our nation’s veterans as its top priority. VA is dedicated to offering timely access to high-quality, evidence-based mental health care that meets the needs of veterans and supports their reintegration into society. One of our core missions is to prevent suicide among veterans through innovative approaches and resources. With funding from the VA Office of Mental Health and Suicide Prevention (OMHSP), the Community Determinants of Health (EDH) Data project has developed innovative datasets associated with specific health outcomes, a methodology for transforming spatiotemporal data from one spatial reference (e.g., a 1 km grid) to another (e.g., U.S. Census Tracts), and capabilities for modeling health outcomes. These datasets represent an enhancement of the Agency for Healthcare Research and Quality (AHRQ), addressing key gaps by introducing finer spatial resolution (Census Tract) and additional geographical covariates into existing data. The curation and standardization of these datasets is a complex task since they often originate from various sources and are measured at different spatial and temporal resolutions. For example, U.S. Census data products typically use census blocks, block groups, or counties, while data like weather data are available on 1 km grids. Some economic data may only be available at the ZIP code level. In this context, ‘standardized’ means that all datasets share the same spatial extent (e.g., U.S. Census Tract and/or County), and ‘curated’ implies a repeatable process with data provenance and the use of appropriate methodologies for covariate conversion. The Community Determinants of Health datasets draw from multiple sources, resulting in variables with varying degrees of availability, patterns of missing data, and methodological considerations across different sources, geographies, and years.

99 GENERAL AND MISCELLANEOUS↗

DOE OSTI2025

Complexity Reduction Methods for Large-Scale Spatially Explicit Biofuels Network Design

The size and complexity of energy system optimization models have increased significantly in recent years, driven by the availability of high-resolution spatial data. We present complexity reduction and solution methods that enable us to efficiently represent high-resolution spatial data in the network design of large-scale energy systems. We aim to reduce the size and enhance the computational efficiency of network design models without sacrificing solution accuracy. Specifically, we first present how to aggregate highly granular data into larger resolutions without averaging out their specific properties through a composite-curve-based approach and then develop a method to linearly represent these curves. Second, we utilize a general clustering method to determine groups of geographically proximate biomass fields and establish a single transportation arc for all of them, reducing the number of transportation-related variables while maintaining an accurate representation of the system. Finally, we introduce a two-step algorithm that decomposes large-scale network design problems into two smaller, more manageable subproblems. We demonstrate the application of our methods using a case study of switchgrass-to-biofuels network design in the eight states of the U.S. Midwest, using realistic and highly explicit spatial data.

09 BIOMASS FUELS↗

DOE OSTI2021

Investigating the Causes and Impacts of Convective Aggregation in a High Resolution Atmospheric GCM

Abstract A ∼50 km resolution atmospheric general circulation model (GCM) is used to investigate the impact of radiative interactions on spatial organization of convection, the model's mean state, and extreme precipitation events in the presence of realistic boundary conditions. Mechanism‐denial experiments are performed in which synoptic‐scale feedbacks between radiation and dynamics are suppressed by overwriting the model‐generated atmospheric radiative cooling rates with its monthly varying climatological values. When synoptic‐scale radiative interactions are disabled, the annual mean circulation and precipitation remain almost unchanged, however tropical convection becomes less aggregated, with an increase in cloud fraction and relative humidity in the free troposphere but a decrease in both variables in the boundary layer. Changes in cloud fraction and relative humidity in the boundary layer exhibit more sensitivity to the presence of radiative interactions than variations in the degree of aggregation. The less aggregated state is associated with a decrease in the frequency of extreme precipitation events, coincident with a decrease in the dynamical contribution to the magnitude of extreme precipitation. At regional scales, the spatial contrast in radiative cooling between dry and moist regions diminishes when radiative interactions are suppressed, reducing the upgradient transport of energy, degree of aggregation, and frequency of extreme precipitation events. However, the mean width of the tropical rain belt remains almost unaffected when radiative interactions are disabled. These results offer insights into how radiation‐circulation coupling affects the spatial organization of convection, distributions of clouds and humidity, and weather extremes.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

Development of observation-based global multilayer soil moisture products for 1970 to 2016

Abstract. Soil moisture (SM) datasets are critical to understanding the global water, energy, and biogeochemical cycles and benefit extensive societal applications. However, individual sources of SM data (e.g., in situ and satellite observations, reanalysis, offline land surface model simulations, Earth system model – ESM – simulations) have source-specific limitations and biases related to the spatiotemporal continuity, resolutions, and modeling and retrieval assumptions. Here, we developed seven global, gap-free, long-term (1970–2016), multilayer (0–10, 10–30, 30–50, and 50–100 cm) SM products at monthly 0.5∘ resolution (available at https://doi.org/10.6084/m9.figshare.13661312.v1; Wang and Mao, 2021) by synthesizing a wide range of SM datasets using three statistical methods (unweighted averaging, optimal linear combination, and emergent constraint). The merged products outperformed their source datasets when evaluated with in situ observations (mean bias from −0.044 to 0.033 m3 m−3, root mean square errors from 0.076 to 0.104 m3 m−3, Pearson correlations from 0.35 to 0.67) and multiple gridded datasets that did not enter merging because of insufficient spatial, temporal, or soil layer coverage. Three of the new SM products, which were produced by applying any of the three merging methods to the source datasets excluding the ESMs, had lower bias and root mean square errors and higher correlations than the ESM-dependent merged products. The ESM-independent products also showed a better ability to capture historical large-scale drought events than the ESM-dependent products. The merged products generally showed reasonable temporal homogeneity and physically plausible global sensitivities to observed meteorological factors, except that the ESM-dependent products underestimated the low-frequency temporal variability in SM and overestimated the high-frequency variability for the 50–100 cm depth. Based on these evaluation results, the three ESM-independent products were finally recommended for future applications because of their better performances than the ESM-dependent ones. Despite uncertainties in the raw SM datasets and fusion methods, these hybrid products create added value over existing SM datasets because of the performance improvement and harmonized spatial, temporal, and vertical coverages, and they provide a new foundation for scientific investigation and resource management.

54 ENVIRONMENTAL SCIENCES↗