skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


This content will become publicly available on February 10, 2026

Title: Filling the Gaps: A Bayesian Mixture Model for Imputing Missing Soil Water Content Data
ABSTRACT Soil water content (SWC) data are central to evaluating how soil moisture varies over time and space and influences critical plant and ecosystem functions, especially in water‐limited drylands. However, sensors that record SWC at high frequencies often malfunction, leading to incomplete timeseries and limiting our understanding of dryland ecosystem dynamics. We developed an analytical approach to impute missing SWC data, which we tested at six eddy flux tower sites along an elevation gradient in the southwestern United States. We impute missing data as a mixture of linearly interpolated SWC between the observed endpoints of a missing data gap and SWC simulated by an ecosystem water balance model (SOILWAT2). Within a Bayesian framework, we allowed the relative utility (mixture weight) of each component (linearly interpolated vs. SOILWAT2) to vary by depth, site and gap characteristics. We explored “fixed” weights versus “dynamic” weights that vary as a function of cumulative precipitation, average temperature, and time since the start of the gap. Both models estimated missing SWC data well (R2 = 0.70–0.88 vs. 0.75–0.91 for fixed vs. dynamic weights, respectively), but the utility of linearly interpolated versus SOILWAT2 values depended on site and depth. SOILWAT2 was more useful for more arid sites, shallower depths, longer and warmer gaps and gaps that received greater precipitation. Overall, the mixture model reliably gap‐fills SWC, while lending insight into processes governing SWC dynamics. This approach to impute missing data could be adapted to accommodate more than two mixture components and other types of environmental timeseries.  more » « less
Award ID(s):
2425290
PAR ID:
10575198
Author(s) / Creator(s):
 ;  ;  ;  ;  ;  ;  
Publisher / Repository:
Wiley Blackwell (John Wiley & Sons)
Date Published:
Journal Name:
Ecohydrology
Volume:
18
Issue:
1
ISSN:
1936-0584
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. We investigated the climatic and ecohydrological controls of the monthly methane emission fluxes from freshwater wetlands across the globe. Fluxes of methane, photosynthetically active radiation (PAR), soil temperature (TS), atmospheric pressure, latent heat flux (LE), wind speed (WS), friction velocity, vapor pressure deficit (VPD), soil water content (SWC), water table depth, and precipitation were obtained from 32 FLUXNET wetland sites. Multivariate pattern recognition techniques of principal component and factor analyses were utilized to classify and group climatic and ecological variables based on their similarity as drivers, examining their interrelation patterns across the different sites. Partial least squares regression models were developed to estimate the relative linkages of methane emission fluxes with the climatic and ecohydrological drivers. When the wetlands were flooded (i.e., positive water table depth relative to the ground), PAR, LE, VPD, and TS had the strongest controls on the methane emission fluxes. However, in the absence of flooding (i.e., negative water table depth), the methane emission fluxes were mainly controlled by SWC and WS. For the wetland sites with unavailable water table depth data, PAR, TS, and WS had the strongest controls on the methane emissions and subsequent transport. Our findings provided important knowledge and insights for predicting and managing methane emissions in freshwater wetlands at a global scale. 
    more » « less
  2. Understanding dynamics of soil water content (SWC) and pore air relative humidity (RHpa), as influenced by wetting-drying cycles, is crucial for sustaining fragile ecosystems of desert lands across the world. However, to date, such an understanding is still incomplete. The objective of this study was to examine such dynamics at a typical desert site within the Horqin Sandy Land, located in Mongolian Plateau of north China. The results indicated that vaporization primarily occurred at a depth of around 10 cm below the ground surface. The diurnal variations of the SWC and RHpa in the top 10 cm soils were much larger than those in the soils at a deeper depth. For a non-rainy day, the SWC and RHpa were mainly determined by the relative magnitude of atmospheric temperature over soil temperature, whereas, for a rainy day, the SWC and RHpa were primarily controlled by the rainfall pattern and amount. The retardation role of the top dry soil layer, which is about 10 cm thick and exists most time at the study site, can effectively prevent the beneath moist soils from being further dried up, and thus is beneficial for sustaining the desert ecosystem. 
    more » « less
  3. This package contains gap-filled daily precipitation values for the 15 NPP sites at Jornada Basin LTER in southern New Mexico, USA. Sites were selected to represent the 5 major ecosystem types in the Chihuahuan Desert (upland grasslands, playa grasslands, mesquite-dominated shrublands, creosotebush-dominated shrublands, tarbush-dominated shrublands). For each ecosystem type, three sites were selected to represent the range in variability in production and plant diversity; thus the locations are not replicates. Gap-filled daily precipitation was calculated for the period from 1980 to 2020 at each site using the closest rain gauges that provided a minimum resolution of daily precipitation data. The Methods section and attached documents describe this process in detail. The rain gauges used are described, with respect to their relationship to NPP sites, in the attached "daily_gapfill_ppt_gauge_usage.csv" file. Although automated weather stations became operational at all NPP sites in 2013 (except P-SMAL, in 2017), updates to this data package are ongoing and are intended to gap-fill any missing or invalid data from the weather stations. 
    more » « less
  4. Eddy covariance serves as one the most effective techniques for long-term monitoring of ecosystem fluxes, however long-term data integrations rely on complete timeseries, meaning that any gaps due to missing data must be reliably filled. To date, many gap-filling approaches have been proposed and extensively evaluated for mature and/or less actively managed ecosystems. Random forest regression (RFR) has been shown to be stable and perform better in these systems than alternative approaches, particularly when filling longer gaps. However, the performance of RFR gap filling remains less certain in more challenging ecosystems, e.g., actively managed agri-ecosystems and following recent land-use change due to management disturbances, ecosystems with relatively low fluxes due to low signal to noise ratios, or for trace gases other than carbon dioxide (e.g., methane). In an extension to earlier work on gap filling global carbon dioxide, water, and energy fluxes, we assess the RFR approach for gap filling methane fluxes globally. We then investigate a range of gap-filling methodologies for carbon dioxide, water, energy, and methane fluxes in challenging ecosystems, including European managed pastures, Southeast Asian converted peatlands, and North American drylands. Our findings indicate that RFR is a competent alternative to existing research standard gap-filling algorithms. The marginal distribution sampling (MDS) is still suggested for filling short (< 12 days) gaps in carbon dioxide fluxes, but RFR is better for filling longer (> 30 days) gaps in carbon dioxide fluxes and also for gap filling other fluxes (e.g. sensible heat, latent energy and methane). In addition, using RFR with globally available reanalysis environmental drivers is effective when measured drivers are unavailable. Crucially, RFR was able to reliably fill cumulative fluxes for gaps > 3 moths and, unlike other common approaches, key environment-flux responses were preserved in the gap-filled data. 
    more » « less
  5. The Re-Greening of the West African Sahel has attracted great interdisciplinary interest since it was originally detected in the mid-2000s. Studies have investigated vegetation patterns at regional scales using a time series of coarse resolution remote sensing analyses. Fewer have attempted to explain the processes behind these patterns at local scales. This research investigates bottom-up processes driving Sahelian greening in the northern Central Plateau of Burkina Faso—a region recognized as a greening hot spot. The objective was to understand the relationship between soil and water conservation (SWC) measures and the presence of trees through a comparative case study of three village terroirs, which have been the site of long-term human ecology fieldwork. Research specifically tests the hypothesis that there is a positive relationship between SWC and tree cover. Methods include remote sensing of high-resolution satellite imagery and aerial photos; GIS procedures; and chi-square statistical tests. Results indicate that, across all sites, there is a significant association between SWC and trees (chi-square = 20.144, p ≤ 0.01). Decomposing this by site, however, points out that this is not uniform. Tree cover is strongly associated with SWC investments in only one village—the one with the most tree cover (chi-square = 39.098, p ≤ 0.01). This pilot study concludes that SWC promotes tree cover but this is heavily modified by local contexts. 
    more » « less