skip to main content


Title: A Statistical Method for Generating Temporally Downscaled Geochemical Tracers in Precipitation
Abstract Sampling intervals of precipitation geochemistry measurements are often coarser than those required by fine-scale hydrometeorological models. This study presents a statistical method to temporally downscale geochemical tracer signals in precipitation so that they can be used in high-resolution, tracer-enabled applications. In this method, we separated the deterministic component of the time series and the remaining daily stochastic component, which was approximated by a conditional multivariate Gaussian distribution. Specifically, statistics of the stochastic component could be explained from coarser data using a newly identified power-law decay function, which relates data aggregation intervals to changes in tracer concentration variance and correlations with precipitation amounts. These statistics were used within a copula framework to generate synthetic tracer values from the deterministic and stochastic time series components based on daily precipitation amounts. The method was evaluated at 27 sites located worldwide using daily precipitation isotope ratios, which were aggregated in time to provide low resolution testing datasets with known daily values. At each site, the downscaling method was applied on weekly, biweekly and monthly aggregated series to yield an ensemble of daily tracer realizations. Daily tracer concentrations downscaled from a biweekly series had average (+/- standard deviation) absolute errors of 1.69‰ (1.61‰) for δ 2 H and 0.23‰ (0.24‰) for δ 18 O relative to observations. The results suggest coarsely sampled precipitation tracers can be accurately downscaled to daily values. This method may be extended to other geochemical tracers in order to generate downscaled datasets needed to drive complex, fine-scale models of hydrometeorological processes.  more » « less
Award ID(s):
1802885 1954660
NSF-PAR ID:
10281865
Author(s) / Creator(s):
; ; ; ;
Date Published:
Journal Name:
Journal of Hydrometeorology
ISSN:
1525-755X
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    The National Ecological Observatory Network (NEON) provides open-access measurements of stable isotope ratios in atmospheric water vapor (δ2H, δ18O) and carbon dioxide (δ13C) at different tower heights, as well as aggregated biweekly precipitation samples (δ2H, δ18O) across the United States. These measurements were used to create the NEON Daily Isotopic Composition of Environmental Exchanges (NEON-DICEE) dataset estimating precipitation (P; δ2H, δ18O), evapotranspiration (ET; δ2H, δ18O), and net ecosystem exchange (NEE; δ13C) isotope ratios. Statistically downscaled precipitation datasets were generated to be consistent with the estimated covariance between isotope ratios and precipitation amounts at daily time scales. Isotope ratios in ET and NEE fluxes were estimated using a mixing-model approach with calibrated NEON tower measurements. NEON-DICEE is publicly available on HydroShare and can be reproduced or modified to fit user specific applications or include additional NEON data records as they become available. The NEON-DICEE dataset can facilitate understanding of terrestrial ecosystem processes through their incorporation into environmental investigations that require daily δ2H, δ18O, and δ13C flux data.

     
    more » « less
  2. Downscaling coarse global and regional climate models allows researchers to access weather and climate data at finer temporal and spatial resolution, but there remains a need to compare these models with empirical data sources to assess model accuracy. Here, we validate a widely used software for generating North American downscaled climate data, ClimateNA, with a novel empirical data source, 20th century weather journals kept by Admiralty Island, Alaska homesteader, Allen Hasselborg. Using Hasselborg’s journals, we calculated monthly precipitation and monthly mean of the maximum daily air temperature across the years 1926 to 1954 and compared these to ClimateNA data generated from the Hasselborg homestead location and adjacent areas. To demonstrate the utility and potential implications of this validation for other disciplines such as hydrology, we used an established regression equation to generate time series of 95% low duration flow estimates for the month of August using mean annual precipitation from ClimateNA predictions and Hasselborg data. Across 279 months, we found strong correlation between modeled and observed measurements of monthly precipitation ( ρ  = 0.74) and monthly mean of the maximum daily air temperature ( ρ  = 0.98). Monthly precipitation residuals (calculated as ClimateNA data - Hasselborg data) generally demonstrated heteroscedasticity around zero, but a negative trend in residual values starting during the last decade of observations may have been due to a shift to the cold-phase Pacific Decadal Oscillation. Air temperature residuals demonstrated a consistent but small positive bias, with ClimateNA tending to overestimate air temperature relative to Hasselborg’s journals. The degree of correlation between weather patterns observed at the Hasselborg homestead site and ClimateNA data extracted from spatial grid cells across the region varied by wet and dry climate years. Monthly precipitation from both data sources tended to be more similar across a larger area during wet years (mean ρ across grid cells = 0.73) compared to dry years (mean ρ across grid cells = 0.65). The time series of annual 95% low duration flow estimates for the month of August generated using ClimateNA and Hasselborg data were moderately correlated ( ρ  = 0.55). Our analysis supports previous research in other regions which also found ClimateNA to be a robust source for past climate data estimates. 
    more » « less
  3. null (Ed.)
    International Ocean Discovery Program (IODP) Expedition 357 successfully cored an east–west transect across the southern wall of Atlantis Massif on the western flank of the Mid-Atlantic Ridge to study the links between serpentinization processes and microbial activity in the shallow subsurface of highly altered ultramafic and mafic sequences that have been uplifted to the seafloor along a major detachment fault zone. The primary goals of this expedition were to (1) examine the role of serpentinization in driving hydrothermal systems, sustaining microbial communities, and sequestering carbon; (2) characterize the tectonomagmatic processes that lead to lithospheric heterogeneities and detachment faulting; and (3) assess how abiotic and biotic processes change with variations in rock type and progressive exposure on the seafloor. To accomplish these objectives, we developed a coring and sampling strategy based around the use of seabed rock drills—the first time that such systems have been used in the scientific ocean drilling programs. This technology was chosen in hopes of achieving high recovery of the carbonate cap sequences and intact contact and deformation relationships. The expedition plans also included several engineering developments to assess geochemical parameters during drilling; sample bottom water before and after drilling; supply synthetic tracers during drilling for contamination assessment; gather downhole electrical resistivity and magnetic susceptibility logs for assessing fractures, fluid flow, and extent of serpentinization; and seal boreholes to provide opportunities for future experiments. Seventeen holes were drilled at nine sites across Atlantis Massif, with two sites on the eastern end of the southern wall (Sites M0068 and M0075), three sites in the central section of the southern wall north of the Lost City hydrothermal field (Sites M0069, M0072, and M0076), two sites on the western end (Sites M0071 and M0073), and two sites north of the southern wall in the direction of the central dome of the massif and Integrated Ocean Drilling Program Site U1309 (Sites M0070 and M0074). Use of seabed rock drills enabled collection of more than 57 m of core, with borehole penetration ranging from 1.3 to 16.44 meters below seafloor and core recoveries as high as 75% of total penetration. This high level of recovery of shallow mantle sequences is unprecedented in the history of ocean drilling. The cores recovered along the southern wall of Atlantis Massif have highly heterogeneous lithologies, types of alteration, and degrees of deformation. The ultramafic rocks are dominated by harzburgites with intervals of dunite and minor pyroxenite veins, as well as gabbroic rocks occurring as melt impregnations and veins, all of which provide information about early magmatic processes and the magmatic evolution in the southernmost portion of Atlantis Massif. Dolerite dikes and basaltic rocks represent the latest stage of magmatic activity. Overall, the ultramafic rocks recovered during Expedition 357 revealed a high degree of serpentinization, as well as metasomatic talc-amphibole-chlorite overprinting and local rodingitization. Metasomatism postdates an early phase of serpentinization but predates late-stage intrusion and alteration of dolerite dikes and the extrusion of basalt. The intensity of alteration is generally lower in the gabbroic and doleritic rocks. Chilled margins in dolerite intruded into talc-amphibole-chlorite schists are observed at the most eastern Site M0075. Deformation in Expedition 357 cores is variable and dominated by brecciation and formation of localized shear zones; the degree of carbonate veining was lower than anticipated. All types of variably altered and deformed ultramafic and mafic rocks occur as components in sedimentary breccias and as fault scarp rubble. The sedimentary cap rocks include basaltic breccias with a carbonate sand matrix and/or fossiliferous carbonate. Fresh glass on basaltic components was observed in some of the breccias. The expedition also successfully applied new technologies, namely (1) extensively using an in situ sensor package and water sampling system on the seabed drills for evaluating real-time dissolved oxygen and methane, pH, oxidation-reduction potential, temperature, and conductivity during drilling; (2) deploying a borehole plug system for sealing seabed drill boreholes at four sites to allow access for future sampling; and (3) proving that tracers can be delivered into drilling fluids when using seabed drills. The rock drill sensor packages and water sampling enabled detection of elevated dissolved methane and hydrogen concentrations during and/or after drilling, with “hot spots” of hydrogen observed over Sites M0068–M0072 and methane over Sites M0070–M0072. Shipboard determination of contamination tracer delivery confirmed appropriate sample handling procedures for microbiological and geochemical analyses, which will aid all subsequent microbiological investigations that are part of the science party sampling plans, as well as verify this new tracer delivery technology for seabed drill rigs. Shipboard investigation of biomass density in select samples revealed relatively low and variable cell densities, and enrichment experiments set up shipboard reveal growth. Thus, we anticipate achieving many of the deep biosphere–related objectives of the expedition through continued scientific investigation in the coming years. Finally, although not an objective of the expedition, we were serendipitously able to generate a high-resolution (20 m per pixel) multibeam bathymetry map across the entire Atlantis Massif and the nearby fracture zone, Mid-Atlantic Ridge, and eastern conjugate, taking advantage of weather and operational downtime. This will assist science party members in evaluating and interpreting tectonic and mass-wasting processes at Atlantis Massif. 
    more » « less
  4. Tracer-aided hydrological models (TAHMs) are one of the most powerful tools to identify new (event) and old (pre-event) water fractions contributing to stormflow because they account both for streamflow and tracer mixing dynamics in model calibration. Nevertheless, their representativeness of hydrograph dynamics is often limited due to the unavailability of high-resolution conservative tracer data (e.g., water stable isotopes or chloride). Hence, there is a need to identify alternative tracers yielding similar flow partitioning results than “ideal” ones while requiring fewer financial resources for high-frequency monitoring (e.g., sub-hourly). Here, wecompare flow partitioning results of a TAHM calibrated using high-frequency electrical conductivity (EC) and water stable isotope (18O) data collected during 37 rainfall-runoff events monitored during variable hydrometeorological conditions in the Zhurucay Ecohydrological Observatory, a tropical alpine catchment located in southern Ecuador. When the model was calibrated using the sampling resolution of stables isotopes (6-hours to 1-hour), no statistically significant differences of pre-event water fractions (PEWFs) using both tracers for model calibration were found. PEWF differences between both tracers for 89% of the events were < 20% regardless of the events’ antecedent moisture and rainfall conditions. Model transfer functions were also similar suggesting that catchment internal processes inferred using both tracers are comparable. Events presenting larger differences (n = 4; up to 27% PEWF difference) had no samples collected during peak flow. Calibration of the model using EC data collected at sub-hourly intervals (every 5-minutes) showed a significant increase in model performance as compared to the frequency of collection of isotopic data. Similarity in flow partitioning results can be attributed to a quasi-conservative nature of EC due to the presence of organic-rich riparian soils (peat-type) overlying compact bedrock across the catchment. Findings also highlight the importance of capturing rapidly occurring catchment mixing processes though high-temporal frequency monitoring of tracer data. Our study encourages the value of assessing the use of alternative tracers, such as EC, to identify fast occurring rainfall runoff processes, while lowering the costs needed to implement and sustain tracer data collection for long time periods. 
    more » « less
  5. Abstract

    Tracer‐aided rainfall‐runoff modelling is a promising tool for understanding catchment hydrology, particularly when tracers provide information about coupled hydrological‐biogeochemical processes. Such models allow for predicting the quality and quantity of water under changing climatic and anthropogenic conditions. Here, we present the Spatially‐distributed Tracer‐Aided Rainfall‐Runoff model with a coupled biogeochemical reactive tracer module (STARR‐DOC) to simulate dissolved organic carbon (DOC) dynamics and sources. The STARR‐DOC model was developed and tested for a humid high Andean ecosystem (páramo) using high‐resolution hourly DOC and hydrometeorological data to simulate hourly discharge and DOC at a fine spatial (10 × 10 m) resolution. Overall, the model was able to acceptably reproduce discharge (KGE ~ 0.45) and stream DOC (KGE ~ 0.69) dynamics. Spatially distributed DOC simulations were independently compared using point DOC measurements for different soil types across the catchment, which allowed for identifying DOC production hot spots and hot moments. Results showed higher hydrological connectivity between slopes and valleys with increasing precipitation. Wetter conditions also favoured DOC production (wet month = 82 mg L−1, dry month = 5 mg L−1) and transport to the stream network (DOC concentrations: during events ~15 mg L−1, during baseflows ~4 mg L−1). Our results also suggest that minor changes in meteorological conditions directly affect páramo soil water dynamics and biogeochemistry. Knowledge of when and where DOC production in mountain catchments is greatest is important for water managers to understand when they make decisions about water security, especially considering climate change predictions for the Andean region.

     
    more » « less