1. Introduction
Naturally occurring chemical signatures in precipitation (e.g., Bailey et al. 2018; Bowen et al. 2019; Gibson et al. 2005; Kendall and McDonnell 2012; Moerman et al. 2013; West et al. 2010; Wiederhold 2015) are frequently used as hydrometeorological tracers, especially when inferring transport or chemical transformations through terrestrial, aquatic, and atmospheric environments (e.g., Abbott et al. 2016; Brooks et al. 2014; Good et al. 2015; Gupta et al. 2020; Kanner et al. 2014; Remondi et al. 2018). Tracer-enabled modeling allows for process-level inference based not only on the size of fluxes, but also on the spatial and temporal transport and mixing of the geochemical signatures associated with the fluxes, thereby facilitating improved understanding and multiresponse model evaluation (Bowen and Good 2015; Krause et al. 2005; McGuire and McDonnell 2006; Sprenger et al. 2019; Turnadge and Smerdon 2014). Researchers have used tracers within global climate models to evaluate processes that are challenging to observe (e.g., ageostrophic circulations, convection, and turbulence) or are modeled at subgrid scales and are therefore not explicitly simulated but parameterized (e.g., Gupta et al. 2020; Orbe et al. 2020; Rosa et al. 2012). For instance, isotope-enabled general circulation models (GCMs) have explicitly simulated water isotope ratios within the critical zone on subdaily time scales [e.g., a version of the Community Earth System Model (iCESM1); Brady et al. 2019; Nusbaumer et al. 2017; Wong et al. 2017] and provide outputs which have been evaluated against observational datasets at various scales (e.g., Hoffmann et al. 2000; Nusbaumer et al. 2017; Risi et al. 2012; Steen-Larsen et al. 2017; Wong et al. 2017).
In many modeling applications, observed and modeled temporal resolutions are different and, in these cases, a downscaling method is required in order to use observed datasets within a model to evaluate processes with dynamic fluctuations over short temporal intervals (Ebtehaj and Foufoula-Georgiou 2013). The temporal resolution at which many geochemical tracers are collected, a result of analytical or logistical cost or the need to aggregate them in time to achieve a measurable signal due to low tracer concentrations, contrasts with the time steps typical of many hydrometeorological models (Rosa et al. 2012; Gupta et al. 2020). Accordingly, a method is needed to generate higher-frequency datasets of precipitation chemistry from low-frequency collections.
Statistical downscaling leverages relationships observed in both fine- and coarse-scale measurements to predict fine-scale variations where only coarse-scale data are available (Ebtehaj and Foufoula-Georgiou 2013; Goncu and Albek 2016). Extensive work has focused on downscaling precipitation rate, including through use of temporal neural networks (e.g., Coulibaly et al. 2005), stochastic methods (e.g., Bordoy and Burlando 2013; D’Onofrio et al. 2014; Poduje and Haberlandt 2017), and conditional multivariate statistical models (e.g., Yang et al. 2010). However, past studies have not temporally downscaled precipitation chemistry data, as is warranted for tracer applications.
Precipitation stable isotope ratios (δ2H and δ18O) are an ideal test case for developing a downscaling method that can benefit tracer applications, if the downscaling can preserve multiscale statistical properties (Ebtehaj and Foufoula-Georgiou 2013). Not only is such downscaling in demand, decades of research demonstrate patterns in precipitation isotope ratios that could be leveraged in downscaling; specifically, precipitation amount often covaries with isotopic composition, attributable to the interplay of diverse climatological, physiographical, and meteorological factors in the evaporation, condensation, and transport of atmospheric moisture (e.g., Aggarwal et al. 2016; Aggarwal et al. 2012; Bowen et al. 2019; Ingraham 1998; Konecky et al. 2019; Lee and Fung 2008; Moore et al. 2016; Risi et al. 2008; West et al. 2010). This (typically inverse) covariation between precipitation rates and isotope ratios, often referred to as an “amount effect,” represents partially systematic variations at subseasonal, monthly, and event time scales (Celle-Jeanton et al. 2001; Conroy et al. 2016; Craig 1961; Craig and Gordon 1965; Gat 1996; Lee and Fung 2008; Moore et al. 2014; Tharammal et al. 2017). If these amount effects share statistical similarities across various time scales, they could support downscaling methods to predict short-term fluctuations. Hypothetically, relationships inferred from sporadic or brief datasets could be used to predict short-term variations in precipitation isotopic composition. Those patterns could be superimposed on the longer time scale seasonal patterns, which tend to follow regional patterns (Bowen et al. 2019; Dansgaard 1964; Feng et al. 2009; Allen et al. 2019), to potentially generate realistic, continuous, high-frequency time series of precipitation isotope ratios.
In this study, we developed and evaluated a downscaling method that uses the statistical structure of observed stable isotope time series to downscale and generate stable isotope time series at finer resolutions. We used daily observations of precipitation amounts and isotope ratios from 27 monitoring stations across the globe. The daily data were artificially aggregated to weekly, biweekly, and monthly scales, using amount-weighted running means to simulate coarser-scale datasets on which to apply the method. These aggregated time series were evaluated for statistical trends, specifically characterizing how the time series means, standard deviations and correlation structures changed as the temporal sampling interval increased. Then, the statistical downscaling method was applied on each of the weekly, biweekly, and monthly aggregated tracer time series to generate downscaled tracer values. An ensemble of downscaled realizations was generated at each site, the statistics of which were compared to those of the original daily observations. Our objective was to generate downscaled realizations that accurately preserved the observed daily δ2H and δ18O means and standard deviations and the correlation structure between precipitation amount, δ2H, and δ18O, so that these realizations could be suitable for various potential modeling applications.
2. Data and methods
a. Site information and tracer datasets
Daily precipitation stable water isotope time series were downloaded from the International Atomic Energy’s (IAEA) Global Network of Isotopes in Precipitation (GNIP) and Water Isotope System of Data Analysis, Visualization and Electronic Retrieval (WISER) database (IAEA 2020). Each time series was filtered to ensure precipitation values were greater than zero and had corresponding δ2H and δ18O isotope ratios. All time series with greater than one year of observations were selected, resulting in the 27 datasets used in the subsequent analysis; details pertaining to each site are included in Table S1 in the online supplemental material. A minimum time series length of one year was chosen because we wanted to account for site-specific seasonal precipitation patterns in the generated downscaled tracer time series. We acknowledge seasonality is usually characterized over time scales greater than one year, however for this analysis we decided on a minimum of one year so the downscaling method could be applied to as many datasets as possible. In section 4, the downscaling method’s performance was evaluated against the number of years represented in the time series and the frequency of collection, i.e., the number of recorded precipitation events divided by the total number of days represented in the time series. The time series lengths ranged from 1.22 to 15.94 years, with an average of 5.34 years. The total number of samples in a time series ranged from 33 to 1026, with an average of ~210. The site with 33 samples (Barasat, Kolkata; Table S1) was sampled over 1.33 years.

The large map displays the 27 GNIP site locations and their average δ2H precipitation measurements. The smaller figure is a dual isotope plot with the mean and standard deviations of all daily precipitation stable water isotope measurements (δ2H, δ18O) at the 27 GNIP sites. The Global Meteoric Water Line (GMWL) is included in the subplot. Refer to Table S1 in the supplemental material for more site-specific characteristics.
Citation: Journal of Hydrometeorology 22, 6; 10.1175/JHM-D-20-0142.1

The large map displays the 27 GNIP site locations and their average δ2H precipitation measurements. The smaller figure is a dual isotope plot with the mean and standard deviations of all daily precipitation stable water isotope measurements (δ2H, δ18O) at the 27 GNIP sites. The Global Meteoric Water Line (GMWL) is included in the subplot. Refer to Table S1 in the supplemental material for more site-specific characteristics.
Citation: Journal of Hydrometeorology 22, 6; 10.1175/JHM-D-20-0142.1
The large map displays the 27 GNIP site locations and their average δ2H precipitation measurements. The smaller figure is a dual isotope plot with the mean and standard deviations of all daily precipitation stable water isotope measurements (δ2H, δ18O) at the 27 GNIP sites. The Global Meteoric Water Line (GMWL) is included in the subplot. Refer to Table S1 in the supplemental material for more site-specific characteristics.
Citation: Journal of Hydrometeorology 22, 6; 10.1175/JHM-D-20-0142.1
Isotope ratios are often evaluated relative to the Global Meteoric Water Line (GMWL), which is defined as δ2H = 8δ18O + 10‰ (Craig 1961). Deuterium excess [d-excess (‰) = δ2H − 8δ18O] measures the deviation of a water sample’s composition from the GMWL (Dansgaard 1964) and is a useful secondary tracer in that it varies with respect to the evaporation and mixing history of air masses (e.g., Benetti et al. 2014; Fröhlich et al. 2002; Pfahl and Sodemann 2014). One can use d-excess to understand both the source of precipitation and the evolution of moisture during transport (Fröhlich et al. 2002; Good et al. 2014). We aimed to preserve a site’s d-excess in the downscaled time series because it can be informative for a variety of hydrological and meteorological applications.
b. Constructing low-resolution datasets
Time series statistics were evaluated across a range of temporal intervals. Moerman et al. (2013) investigated the correlation structure between precipitation amount and δ18O at Mulu Meteo, Sarawak (Table S1) at daily to 12-week (84 days) time scales. Following their approach, we evaluated trends in the mean μ, standard deviation σ, and Pearson correlation coefficient ρ at different temporal intervals to capture the time series response and prediction accuracy. The ρ measures the linear correlation between two variables and has a value between −1 and 1, where 1 is a total positive linear correlation, 0 is no linear correlation, and −1 is total negative linear correlation. At the daily scale for the 27 GNIP sites, the average (±standard deviation) ρ(P, δ2H) was −0.18 (±0.18), ρ(P, δ18O) was −0.20 (±0.17), and ρ(δ2H, δ18O) was 0.96 (±0.03).
c. Statistical precipitation tracer downscaling method
1) Removal of the deterministic time series component
Each aggregated weekly, biweekly, and monthly time series [Eq. (2)] was treated as an example of a low-resolution dataset on which to apply the downscaling method. We considered each tracer time series to have a deterministic component and a stochastic component. In the first step, the deterministic component was characterized by the seasonality in the precipitation signal and was removed from each set of observations. Isotope ratios in precipitation frequently have been observed to exhibit distinct seasonal signals. These can be approximated as a combination of sinusoidal functions through Fourier decomposition (Allen et al. 2018, 2019; Dutton et al. 2005; Feng et al. 2009; Halder et al. 2015; Vachon et al. 2007; Wilkinson and Ivany 2002). Sinusoidal functions effectively describe the collinear structure and fluctuations in the covariation of δ18H and δ18O relative to the GMWL (Fig. 1; Allen et al. 2018; Craig 1961, Dansgaard 1964). The sine curve parameters (amplitude, phase, and offset) are often predictable in space (Allen et al. 2018; Jasechko et al. 2016) and succinctly represent temporal dynamics because they express continuous, cyclic time series. Allen et al. (2019) used monthly isotopes in precipitation GNIP datasets from across the globe to capture patterns in the precipitation isotope seasonality using sinusoidal functions. When predicting the isotope seasonality, the values of the sine parameters can be described as functions of climate and geography. Additionally, sine curves are useful when describing the propagation of cyclic signals, this has been done to infer catchment-scale mixing processes using the dampening ratio of seasonal isotope amplitudes in streamflow versus precipitation (Kirchner 2016a,b; also see Clow et al. 2018; von Freyberg et al. 2018; Gallart et al. 2020; Jacobs et al. 2018; Jasechko et al. 2016; Lutz et al. 2018; Song et al. 2017).
2) Generation of stochastic tracer realizations
Next, the daily statistics of the stochastic hydrogen (δ2H*) and oxygen (δ18O*) isotope time series were estimated by using the relationship between the observed daily stochastic statistics and the stochastic signal’s statistics across a range of aggregation intervals (t) multiplied by each site’s specific precipitation frequency (λ; defined as the number of days with precipitation divided by the total number of days in a time series). The statistics of the stochastic signal at aggregation interval t were denoted with * as

The x axis is λ (recorded events/number of days in the time series) multiplied by daily to 12-week aggregation intervals (days) and the y axes were the deviations for each of the 27 sites in their stochastic time series (a),(b) means (
Citation: Journal of Hydrometeorology 22, 6; 10.1175/JHM-D-20-0142.1

The x axis is λ (recorded events/number of days in the time series) multiplied by daily to 12-week aggregation intervals (days) and the y axes were the deviations for each of the 27 sites in their stochastic time series (a),(b) means (
Citation: Journal of Hydrometeorology 22, 6; 10.1175/JHM-D-20-0142.1
The x axis is λ (recorded events/number of days in the time series) multiplied by daily to 12-week aggregation intervals (days) and the y axes were the deviations for each of the 27 sites in their stochastic time series (a),(b) means (
Citation: Journal of Hydrometeorology 22, 6; 10.1175/JHM-D-20-0142.1
The ratio of
Pseudo-random numbers were generated using a Gaussian copula (Sklar 1959), defined by the estimated daily statistics,
For each observed precipitation amount, values of δ2H* and δ18O* (a two-number sample representing the stochastic signal) were drawn from a multivariate Gaussian distribution using Python’s (v3.7.6) SciPy Library (v1.2.1) with parameters described by
3) Full synthetic time series generation
The deterministic component, the sinusoidal function from Eq. (3), was added to each generated stochastic time series. The result was a downscaled tracer time series, which captured site-specific daily precipitation amount effects, seasonal signals, and stochastic variability. Finally, we applied a residual correction on the downscaled synthetic series to preserve the observed aggregated weighted tracer values. For each synthetic value within each aggregation interval, an interval-specific, single correction factor was subtracted from the downscaled values so that there was no difference between that period’s downscaled synthetic values aggregated for that interval and the observed coarse-resolution interval’s value. In doing so, the precipitation-weighted values of the synthetic time series then equaled the known aggregated value. This property is particularly important as it closes the tracer mass balance. The statistical downscaling method applied to a dataset with a biweekly sampling frequency was summarized and visualized in section ii of the supplemental material.
d. Evaluation of precipitation tracer downscaling mythology
The statistical method was iterated over 100 times generating an ensemble of downscaled isotope time series at each of the 27 GNIP locations. The large number of time series generated for each ensemble allowed for us to quantify the performance of the downscaling method. Each ensemble was expected to capture the observed site-specific tracer means and standard deviations and the correlation coefficients between precipitation amount, δ2H and δ18O. The statistical downscaling method was evaluated using multiple techniques, detailed in the subsequent paragraphs.
After removing the deterministic components, the stochastic time series were expected to have means of approximately zero, a predicable decrease in standard deviation [Eq. (4)] and Pearson correlation coefficients at low temporal resolutions appropriately defining daily covariate structures. To test this,
The average of ensemble means
3. Results
a. Evaluation of estimated daily stochastic signal statistics
The estimated daily stochastic signal statistics from weekly, biweekly and monthly aggregation intervals accurately described the observed statistics (Fig. 3). The method best predicted

The estimated (a),(b) standard deviations and (c)–(e) Pearson correlation coefficients of the stochastic signal from downscaled weekly, biweekly, and monthly time series compared to the observed daily stochastic statistics. Each data point is one site location, and the black lines are the 1:1 lines. The means were not shown because they are approximately zero (refer to Figs. S2a,b).
Citation: Journal of Hydrometeorology 22, 6; 10.1175/JHM-D-20-0142.1

The estimated (a),(b) standard deviations and (c)–(e) Pearson correlation coefficients of the stochastic signal from downscaled weekly, biweekly, and monthly time series compared to the observed daily stochastic statistics. Each data point is one site location, and the black lines are the 1:1 lines. The means were not shown because they are approximately zero (refer to Figs. S2a,b).
Citation: Journal of Hydrometeorology 22, 6; 10.1175/JHM-D-20-0142.1
The estimated (a),(b) standard deviations and (c)–(e) Pearson correlation coefficients of the stochastic signal from downscaled weekly, biweekly, and monthly time series compared to the observed daily stochastic statistics. Each data point is one site location, and the black lines are the 1:1 lines. The means were not shown because they are approximately zero (refer to Figs. S2a,b).
Citation: Journal of Hydrometeorology 22, 6; 10.1175/JHM-D-20-0142.1
b. Evaluation of the downscaled tracer Realizations
The average of each ensemble’s means

The average (a),(b) means; (c),(d) standard deviations; and (e)–(g) Pearson correlation coefficients of the downscaled ensembles from the weekly, biweekly, and monthly time series compared to the observed daily site statistics. Each data point is one location, and the black lines are the 1:1 lines.
Citation: Journal of Hydrometeorology 22, 6; 10.1175/JHM-D-20-0142.1

The average (a),(b) means; (c),(d) standard deviations; and (e)–(g) Pearson correlation coefficients of the downscaled ensembles from the weekly, biweekly, and monthly time series compared to the observed daily site statistics. Each data point is one location, and the black lines are the 1:1 lines.
Citation: Journal of Hydrometeorology 22, 6; 10.1175/JHM-D-20-0142.1
The average (a),(b) means; (c),(d) standard deviations; and (e)–(g) Pearson correlation coefficients of the downscaled ensembles from the weekly, biweekly, and monthly time series compared to the observed daily site statistics. Each data point is one location, and the black lines are the 1:1 lines.
Citation: Journal of Hydrometeorology 22, 6; 10.1175/JHM-D-20-0142.1

The average (a),(b) means; (c),(d) standard deviations; and (e)–(g) Pearson correlation coefficients of the residual corrected downscaled ensembles from the weekly, biweekly, and monthly time series compared to the observed daily site statistics. Each data point is one location, and the black lines are the 1:1 lines.
Citation: Journal of Hydrometeorology 22, 6; 10.1175/JHM-D-20-0142.1

The average (a),(b) means; (c),(d) standard deviations; and (e)–(g) Pearson correlation coefficients of the residual corrected downscaled ensembles from the weekly, biweekly, and monthly time series compared to the observed daily site statistics. Each data point is one location, and the black lines are the 1:1 lines.
Citation: Journal of Hydrometeorology 22, 6; 10.1175/JHM-D-20-0142.1
The average (a),(b) means; (c),(d) standard deviations; and (e)–(g) Pearson correlation coefficients of the residual corrected downscaled ensembles from the weekly, biweekly, and monthly time series compared to the observed daily site statistics. Each data point is one location, and the black lines are the 1:1 lines.
Citation: Journal of Hydrometeorology 22, 6; 10.1175/JHM-D-20-0142.1
Average bias (predicted − observed statistic) (±standard deviation) for downscaled and residual corrected downscaled (in italics) ensembles.


c. Conserved processes with the method
An analysis of the observed time series demonstrates strong autocorrelation; when the seasonal signal is removed, the observed autocorrelation is nearly all removed (Fig. 6). In fact, the median autocorrelation of the observed time series stochastic signals falls below 5% after 3-day lags and are approximately zero at 4-day lags, supporting the assumption that the sinusoidal function adequately described the deterministic component and the residual was stationary (i.e., white noise). The Shapiro–Wilk and the D’Agostino’s K2 normality tests suggest that we could not reject the assumption of normality in the weekly, biweekly, and monthly time series (p value > 0.05). Histograms of the stochastic signals for both isotope ratios across all 27 sites are provided in the supplemental material (Fig. S3). Next, we calculated autocorrelations from 1- to 20-day lags of the residual corrected downscaled ensembles. The autocorrelations mimicked the observed temporal trends and memory of the daily time series. Autocorrelations for δ2H and δ2H* (Fig. 6) and are highly correlated with trends observed in the autocorrelations for δ18O and δ18O* (refer to the supplemental material, Fig. S4). Based on the results from the autocorrelation analysis and normality tests, we concluded the addition of the seasonal signal to the generated stochastic time series captured the large majority of the observed tracer memory in the system.

(a),(c),(e) Median autocorrelation of the observed daily δ2H datasets and the daily residual corrected ensembles (solid lines). (b),(d),(f) Median autocorrelations of the δ2H* stochastic signals for the observations and the downscaled ensembles. The 5th–95th percentiles of the observed and ensemble autocorrelations are represented as shaded regions. The horizontal red dashed lines indicate where ρ is ±5%.
Citation: Journal of Hydrometeorology 22, 6; 10.1175/JHM-D-20-0142.1

(a),(c),(e) Median autocorrelation of the observed daily δ2H datasets and the daily residual corrected ensembles (solid lines). (b),(d),(f) Median autocorrelations of the δ2H* stochastic signals for the observations and the downscaled ensembles. The 5th–95th percentiles of the observed and ensemble autocorrelations are represented as shaded regions. The horizontal red dashed lines indicate where ρ is ±5%.
Citation: Journal of Hydrometeorology 22, 6; 10.1175/JHM-D-20-0142.1
(a),(c),(e) Median autocorrelation of the observed daily δ2H datasets and the daily residual corrected ensembles (solid lines). (b),(d),(f) Median autocorrelations of the δ2H* stochastic signals for the observations and the downscaled ensembles. The 5th–95th percentiles of the observed and ensemble autocorrelations are represented as shaded regions. The horizontal red dashed lines indicate where ρ is ±5%.
Citation: Journal of Hydrometeorology 22, 6; 10.1175/JHM-D-20-0142.1
The means and standard deviations in d-excess were accurately captured in the resulting downscaled time series (Fig. 7). At each site, d-excess was calculated for the observed daily series and each ensemble from the downscaled weekly, biweekly, and monthly time scales. The downscaled d-excess was overestimated for the three sites with lowest observed d-excess, indicating potential effects to the downscaling method’s performance when precipitation is predominantly composed of evaporated waters. These d-excess estimates provide a metric for evaluating the downscaled series relative to the GMWL and increases the applicability of this method for tracing meteorological forcing variables and their constituents through modeling environments. Alternative downscaling approaches that independently model δ2H and δ18O may not preserve d-excess signals and thus would provide precipitation predictions that should not be used in simulations that leverage the information provided by dual-isotope analyses.

The average d-excess of the residual corrected downscaled ensemble at each site location compared to the average observed d-excess. Each data point is one site location, and the black line is the 1:1 line.
Citation: Journal of Hydrometeorology 22, 6; 10.1175/JHM-D-20-0142.1

The average d-excess of the residual corrected downscaled ensemble at each site location compared to the average observed d-excess. Each data point is one site location, and the black line is the 1:1 line.
Citation: Journal of Hydrometeorology 22, 6; 10.1175/JHM-D-20-0142.1
The average d-excess of the residual corrected downscaled ensemble at each site location compared to the average observed d-excess. Each data point is one site location, and the black line is the 1:1 line.
Citation: Journal of Hydrometeorology 22, 6; 10.1175/JHM-D-20-0142.1
4. Discussion
a. Method evaluation for select site characteristics
We compared the absolute error, calculated by taking the absolute value of the mean of the downscaled ensembles minus the observed mean, to the site’s latitude, calculated rainfall frequency λ, and total length of the time series in years (Fig. 8). The largest absolute errors of the mean resulted from downscaling calculations that used monthly aggregated data, yielding average (±standard deviation) absolute errors of 2.26‰ (2.54‰) for δ2H and 0.33‰ (0.35‰) for δ18O. Linear regressions between site latitude and absolute errors in the means (derived from monthly, biweekly, and weekly ensembles) showed no strong correlations, suggesting that performance may be partially climate independent (Fig. 8). Absolute errors were also not related to the strength of the seasonal isotopic variation, nor were they related to the overall variability in isotopic composition (as quantified by the standard deviation; Fig. 9). Alternatively, a weak, but significant relationship was observed between absolute error and λ (R2 = 0.25 and p value = 0.0001 for δ2H, R2 = 0.23 and p value = 0.0002 for δ18O) and average recorded precipitation amount (R2 = 0.12 and p value = 0.009 for δ2H, R2 = 0.13 and p value = 0.0003 for δ18O) for downscaled weekly ensembles, but not for downscaled biweekly or monthly ensembles (Figs. 8 and 9). Although not a site characteristic, time series length significantly influenced absolute errors of downscaled biweekly and monthly ensembles. Longer time series spanning many years support better accounting for interannual variability and removing potential biases toward certain seasons. Nonlinear effects (e.g., continentality; Dansgaard 1964; Rozanski et al. 1993) may be contributing to relatively high absolute errors, especially at the subtropics and midlatitudes (Figs. 8a,b). When applying the downscaling method to datasets from these regions, one can adapt the copula framework to account for other influential site-specific characteristics (refer to section 2b).

Absolute error E of the residual corrected ensemble means compared to various site-specific characteristics: (a),(b) latitude; (c),(d) λ; and (e),(f) total length of the time series. Each data point is one site location, dashed lines represent p values > 0.05, and solid lines represent p values < 0.05.
Citation: Journal of Hydrometeorology 22, 6; 10.1175/JHM-D-20-0142.1

Absolute error E of the residual corrected ensemble means compared to various site-specific characteristics: (a),(b) latitude; (c),(d) λ; and (e),(f) total length of the time series. Each data point is one site location, dashed lines represent p values > 0.05, and solid lines represent p values < 0.05.
Citation: Journal of Hydrometeorology 22, 6; 10.1175/JHM-D-20-0142.1
Absolute error E of the residual corrected ensemble means compared to various site-specific characteristics: (a),(b) latitude; (c),(d) λ; and (e),(f) total length of the time series. Each data point is one site location, dashed lines represent p values > 0.05, and solid lines represent p values < 0.05.
Citation: Journal of Hydrometeorology 22, 6; 10.1175/JHM-D-20-0142.1

Absolute error E of the residual corrected ensemble means compared to various site-specific characteristics: (a),(b) the sinusoidal function’s estimated amplitude; (c),(d) standard deviation of each isotope ratio; and (e),(f) average daily precipitation. Each data point is one site location, dashed lines represent p values > 0.05, and solid lines represent p values < 0.05.
Citation: Journal of Hydrometeorology 22, 6; 10.1175/JHM-D-20-0142.1

Absolute error E of the residual corrected ensemble means compared to various site-specific characteristics: (a),(b) the sinusoidal function’s estimated amplitude; (c),(d) standard deviation of each isotope ratio; and (e),(f) average daily precipitation. Each data point is one site location, dashed lines represent p values > 0.05, and solid lines represent p values < 0.05.
Citation: Journal of Hydrometeorology 22, 6; 10.1175/JHM-D-20-0142.1
Absolute error E of the residual corrected ensemble means compared to various site-specific characteristics: (a),(b) the sinusoidal function’s estimated amplitude; (c),(d) standard deviation of each isotope ratio; and (e),(f) average daily precipitation. Each data point is one site location, dashed lines represent p values > 0.05, and solid lines represent p values < 0.05.
Citation: Journal of Hydrometeorology 22, 6; 10.1175/JHM-D-20-0142.1
b. Method adaptation for broader applications
In this study, weekly, biweekly, and monthly data were used to generate daily observations, but more sophisticated applications could potentially be supported by different datasets. As a general rule, the deterministic time series component can be more accurately estimated with increased tracer sampling frequencies (Figs. 8c,d) and samples collected over longer time frames (Figs. 8e,f). Accurately representing the deterministic component increases the likelihood of a downscaled synthetic time series effectively representing the underlying seasonal patterns and interannual variability at a site. Depending on the application, one may increase or decrease the temporal downscaling intervals beyond daily or 12-week time scales. While not evaluated in this study, one could predict subdaily datasets with appropriate observation datasets or known statistical properties (i.e., mean, standard deviation, covariance structure of precipitation and its tracer composition) of a site at subdaily scales (e.g., diurnal cycle).
Theoretically, the downscaling methods used in this study can be expanded to higher dimensions and account for other tracer covariates including site conditions such as air temperature and relative humidity. At sites where the method under or overestimates the site statistics, other meteorological variables, such as air temperature, may correlate more strongly with isotope signals than precipitation amount. To do this, one needs to increase the number of covariates accounted for and the matrix dimensions within the copula framework [refer to supplemental material, section i: “Definition of a Gaussian copula”]. In these instances, adding more known dimensions to Eqs. (6) and (7) (in the supplemental material) will incorporate additive information into the generated downscaled time series. Including additional known covariates within the copula framework may improve the representation of nonlinear effects at sites in the subtropics and midlatitudes if meteorological variables (e.g., relative humidity, air temperature) are highly correlated with changes in tracer concentrations (Figs. 8a,b).
Not only would a downscaled time series facilitate running more detailed models that improve process understanding, but they also allow for better tracking of uncertainties associated with inferences drawn from those models. We compared the mean of the observed biweekly series and the mean of the downscaled biweekly ensemble aggregated to biweekly time scales using Eq. (2). The absolute error of the mean across all sites was 0.90‰ for δ2H and 0.14‰ for δ18O. This suggests models using downscaled tracers would mimic temporal trends observed at biweekly time scales, while also tracking processes and uncertainties only discernible at finer time scales. As expected, when the residual corrected downscaled biweekly ensemble was aggregated to biweekly time scales, the absolute error of the mean was approximately zero. To evaluate how the downscaling method compared to a naive downscale with no high-frequency statistical information, we created a daily time series where all precipitation events that occurred within each 14-day interval had the same isotopic composition equal to the observed biweekly values. The absolute error of the mean across all sites (±standard deviation) was 2.74‰ (2.24‰) for δ2H and 0.39‰ (0.31‰) for δ18O, which was higher than the absolute error of the mean calculated for all the downscaled biweekly ensembles [1.69‰ (1.61‰) for δ2H and 0.23‰ (0.24‰) for δ18O].
Due to limited data, all of the above analyses used the entire dataset to calculate the statistics, fit the models and apply the downscaling method. At sites with more than 5 years of data, we used the first 4 years to build a downscaling model to apply on the fifth year’s precipitation time series. We generated an ensemble of 100 downscaled δ2H and δ18O time series at each site and compared it to the observed δ2H and δ18O from the fifth year of the time series. Based on training and testing sizes, eight sites were used in this analysis and the absolute error of the mean (±standard deviation) for δ2H was 4.80‰ (3.17‰), 4.89‰ (3.26‰), and 5.52‰ (3.50‰) downscaled from weekly, biweekly, and monthly series, respectively. The absolute error of the mean (±standard deviation) for δ18O was 0.79‰ (0.59‰), 0.78‰ (0.63‰), and 0.85‰ (0.67‰) downscaled from weekly, biweekly, and monthly series, respectively. Based on these promising results, our downscaling method could be built using several years of precipitation data with a known concentration and then applied to years where only precipitation amount is available.
This method can be broadly applied to produce ensembles of downscaled datasets for various geochemical modeling applications. Ensembles decrease the risk of tying conclusions to one specific time series. The downscaled ensembles can be generated using the same statistics (like shown here) or multiple ensembles can be generated with varying statistical properties. Examples of different ensembles include time series generated from downscaled statistics estimated from different aggregation intervals (e.g., weekly and biweekly), employing a non-Gaussian copula framework (e.g., Gumbel copula, extreme-value copula) to populate a conditioned stochastic signal’s time series, and increasing dimensions of the copula framework by including additive meteorological variables (e.g., air temperature). Correspondingly, the geochemical tracer ensembles could be used for model selection and with numerous model and parameter sensitivity and uncertainty analyses. Ensembles could be useful in developing frameworks for model–data fusion by merging observational data with model outputs to improve model quality and characterize its uncertainty.
This downscaling approach could be extended across large spatial extents for use in global isotopic models or empirically based geographic simulations to represent sites with limited or no high-frequency observations available. To do this, one could generate downscaled geochemical tracers correlating with precipitation inputs at the grid scale. Last, the methodology can be applied to other geochemical tracers for understanding site-specific dynamics (e.g., chemical leaching, sediment transport, and loading) or climatological applications (e.g., nitrogen deposition, carbon sequestration).
5. Conclusions
This statistical downscaling method generates datasets that maintain informative site-specific correlation structures between covariates and the geochemical tracer and retains the statistical properties of underlying processes (e.g., d-excess, amount effects). By modeling hydrologic dynamics using downscaled tracers, researchers can enhance understanding of physical processes without collecting fine temporal in situ data. While an individual realization of this downscaling approach may generate reasonable estimates of true high-frequency values, iterating analyses using an ensemble of realizations allows for uncertainties in generated time series to be propagated through subsequent modeling and tracer-based analyses. The method is sufficiently general and can be applied for a variety of applications to generate downscaled ensembles for use in meteorological and hydrometeorological models to evaluate model performance, investigate system processes across spatial scales and is additive to model-data fusion frameworks.
Acknowledgments
The authors acknowledge the support of the U.S. National Science Foundation (DEB-1802885, DEB-1802880 and AGS-1954660) and thank the anonymous reviewers for their time and effort spent reviewing this manuscript. The authors state no conflict of interest.
Data availability statement
Python code is provided in the supplemental material to generate geochemical synthetic time series based on the user’s site-specific time series statistics. The code is intended to be easily adaptable to higher dimensions or other user specific applications. Additional materials can be made available upon request.
REFERENCES
Abbott, B. W., and Coauthors, 2016: Using multi-tracer inference to move beyond single catchment ecohydrology. Earth-Sci. Rev., 160, 19–42, https://doi.org/10.1016/j.earscirev.2016.06.014.
Aggarwal, P. K., O. A. Alduchov, K. O. Froehlich, L. J. Araguas-Araguas, N. C. Sturchio, and N. Kurita, 2012: Stable isotopes in global precipitation: A unified interpretation based on atmospheric moisture residence time. Geophys. Res. Lett., 39, L11705, https://doi.org/10.1029/2012GL051937.
Aggarwal, P. K., U. Romatschke, L. Araguas-Araguas, D. Belachew, F. J. Longstaffe, P. Berg, C. Schumacher, and A. Funk, 2016: Proportions of convective and stratiform precipitation revealed in water isotope ratios. Nat. Geosci., 9, 624–629, https://doi.org/10.1038/ngeo2739.
Allen, S. T., J. W. Kirchner, and G. R. Goldsmith, 2018: Predicting spatial patterns in precipitation isotope (δ2H and δ18O) seasonality using sinusoidal isoscapes. Geophys. Res., 45, 4859–4868, https://doi.org/10.1029/2018GL077458.
Allen, S. T., S. Jasechko, W. R. Berghuijs, J. M. Welker, G. R. Goldsmith, and J. W. Kirchner, 2019: Global sinusoidal seasonality in precipitation isotopes. Hydrol. Earth Syst. Sci., 23, 3423–3436, https://doi.org/10.5194/hess-23-3423-2019.
Bailey, A., E. Posmentier, and X. Feng, 2018: Patterns of evaporation and precipitation drive global isotopic changes in atmospheric moisture. Geophys. Res. Lett., 45, 7093–7101, https://doi.org/10.1029/2018GL078254.
Bárdossy, A., and G. G. S. Pegram, 2009: Copula based multisite model for daily precipitation simulation. Hydrol. Earth Syst. Sci., 13, 2299–2314, https://doi.org/10.5194/hess-13-2299-2009.
Benetti, M., G. Reverdin, C. Pierre, L. Merlivat, C. Risi, H. C. Steen-Larsen, and F. Vimeux, 2014: Deuterium excess in marine water vapor: Dependency on relative humidity and surface wind speed during evaporation. J. Geophys. Res. Atmos., 119, 583–593, https://doi.org/10.1002/2013JD020535.
Bordoy, R., and P. Burlando, 2013: Stochastic downscaling of climate model precipitation outputs in orographically complex regions: 2. Downscaling methodology. Water Resour. Res., 50, 562–579, https://doi.org/10.1002/wrcr.20443.
Bowen, G. J., and S. P. Good, 2015: Incorporating water isotopes in hydrological and water resource investigations. Wiley Interdiscip. Rev.: Water, 2, 107–119, https://doi.org/10.1002/wat2.1069.
Bowen, G. J., Z. Cai, R. P. Fiorella, and A. L. Putman, 2019: Isotopes in the water cycle: Regional- to global-scale patterns and applications. Annu. Rev. Earth Planet. Sci., 47, 453–479, https://doi.org/10.1146/annurev-earth-053018-060220.
Brady, E., and Coauthors, 2019: The connected isotopic water cycle in the community Earth system model version 1. J. Adv. Model. Earth Syst., 11, 2547–2566, https://doi.org/10.1029/2019MS001663.
Brooks, J. R., J. J. Gibson, S. J. Birks, M. H. Weber, K. D. Rodecap, and J. L. Stoddard, 2014: Stable isotope estimates of evaporation: Inflow and water residence time for lakes across the United States as a tool for national lake water quality assessments. Limnol. Oceanogr., 59, 2150–2165, https://doi.org/10.4319/lo.2014.59.6.2150.
Celle-Jeanton, H., Y. Travi, and B. Blavoux, 2001: Isotopic typology of the precipitation in the western Mediterranean region at three different time scales. Geophys. Res. Lett., 28, 1215–1218, https://doi.org/10.1029/2000GL012407.
Clow, D. W., M. A. Mast, and J. O. Sickman, 2018: Linking transit times to catchment sensitivity to atmospheric deposition of acidity and nitrogen in mountains of the western United States. Hydrol. Processes, 32, 2456–2470, https://doi.org/10.1002/hyp.13183.
Cong, R., and M. Brady, 2012: The interdependence between rainfall and temperature: Copula analysis. Sci. World J., 2012, 1–11, https://doi.org/10.1100/2012/405675.
Conroy, J. L., D. Noone, K. M. Cobb, J. W. Moerman, and B. L. Konecky, 2016: Paired stable isotopologues in precipitation and vapor: A case study of the amount effect within western tropical Pacific storms. J. Geophys. Res. Atmos., 121, 3290–3303, https://doi.org/10.1002/2015JD023844.
Coulibaly, P., Y. B. Dibike, and F. Anctil, 2005: Downscaling precipitation and temperature with temporal neutral networks. J. Hydrometeor., 6, 483–496, https://doi.org/10.1175/JHM409.1.
Craig, H., 1961: Isotopic variations in meteoric waters. Science, 133, 1702–1703, https://doi.org/10.1126/science.133.3465.1702.
Craig, H., and L. I. Gordon, 1965: Deuterium and oxygen-18 variations in the ocean and the marine atmosphere. Stable Isotopes in Oceanographic Studies and Paleotemperatures, E. Tongiorgi, Ed., Consiglio Nazionale della Ricerche, 9–130.
Dansgaard, W., 1964: Stable isotopes in precipitation. Tellus, 16, 436–468, https://doi.org/10.3402/tellusa.v16i4.8993.
D’Onofrio, D., E. Palazzi, J. von Hardenberg, A. Provenzale, and S. Calmanti,2014: Stochastic rainfall downscaling of climate models. J. Hydrometeor., 15, 830–843, https://doi.org/10.1175/JHM-D-13-096.1.
Dutton, A., B. H. Wilkinson, J. M. Welker, G. J. Bowen, and K. C. Lohmann, 2005: Spatial distribution and seasonal variation in 18O/16O of modern precipitation and river water across the conterminous USA. Hydrol. Processes, 19, 4121–4146, https://doi.org/10.1002/hyp.5876.
Ebtehaj, A. M., and E. Foufoula-Georgiou, 2013: On variational downscaling, fusion, and assimilation of hydrometeorological states: A unified framework via regularization. Water Resour. Res., 49, 5944–5963, https://doi.org/10.1002/wrcr.20424.
Feng, X., A. M. Faiia, and E. S. Posmentier, 2009: Seasonality of isotopes in precipitation: A global perspective. J. Geophys. Res., 114, D08116, https://doi.org/10.1029/2008JD011279.
Fröhlich, K., J. J. Gibson, and P. K. Aggarwal, 2002: Deuterium excess in precipitation and its climatological significance. Study of environmental change using isotope techniques, C&S Papers Series 13/P, International Atomic Energy Agency, 54–65.
Gallart, F., M. Valiente, P. Llorens, C. Cayuela, M. Sprenger, and J. Latron, 2020: Investigating young water fractions in a small Mediterranean mountain catchment: Both precipitation forcing and sampling frequency matter. Hydrol. Processes, 34, 3618–3634, https://doi.org/10.1002/hyp.13806.
Gao, C., Y. Xu, Q. Zhu, Z. Bai, and L. Liu, 2018: Stochastic generation of daily rainfall events: A single-site rainfall model with copula-based joint simulation of rainfall characteristics and classification and simulation of rainfall patterns. J. Hydrol., 564, 41–58, https://doi.org/10.1016/j.jhydrol.2018.06.073.
Gat, J. R., 1996: Oxygen and hydrogen isotopes in the hydrologic cycle. Annu. Rev. Earth Planet. Sci., 24, 225–262, https://doi.org/10.1146/annurev.earth.24.1.225.
Gibson, J. J., T. W. D. Edwards, S. J. Birks, N. A. St Amour, W. M. Buhay, P. McEachern, B. B. Wolfe, and D. L. Peters, 2005: Progress in isotope tracer hydrology in Canada. Hydrol. Processes, 19, 303–327, https://doi.org/10.1002/hyp.5766.
Goncu, S., and E. Albek, 2016: Statistical downscaling of meteorological time series and climatic projections in a watershed in Turkey. Theor. Appl. Climatol., 126, 191–211, https://doi.org/10.1007/s00704-015-1563-2.
Good, S. P., D. V. Mallia, J. C. Lin, and G. J. Bowen, 2014: Stable isotope analysis of precipitation samples obtained via crowdsourcing reveals the spatiotemporal evolution of superstorm sandy. PLOS ONE, 9, e91117, https://doi.org/10.1371/journal.pone.0091117.
Good, S. P., D. Noone, and G. Bowen, 2015: Hydrologic connectivity constrains partitioning of global terrestrial water fluxes. Science, 349, 175–177, https://doi.org/10.1126/science.aaa5931.
Gupta, A., E. P. Gerber, and P. H. Lauritzen, 2020: Numerical impacts on tracer transport: A proposed intercomparison test of atmospheric general circulation models. Quart. J. Roy. Meteor. Soc., 146, 3937–3964, https://doi.org/10.1002/qj.3881.
Gyasi-Agyei, Y., 2011: Copula-based daily rainfall disaggregation model. Water Resour. Res., 47, W07535, https://doi.org/10.1029/2011WR010519.
Halder, J., S. Terzer, L. I. Wassenaar, L. J. Araguás-Araguás, and P. K. Aggarwal, 2015: The Global Network of Isotopes in Rivers (GNIR): Integration of water isotopes in watershed observation and riverine research. Hydrol. Earth Syst. Sci., 19, 3419–3431, https://doi.org/10.5194/hess-19-3419-2015.
Haslauer, C. P., P. Guthke, A. Bádossy, and E. A. Sudicky, 2012: Effects of non-Gaussian copula-based hydraulic conductivity fields on macrodispersion. Water Resour. Res., 48, W07507, https://doi.org/10.1029/2011WR011425.
Hoffmann, G., J. Jouzel, and V. Masson, 2000: Stable water isotopes in atmospheric general circulation models. Hydrol. Processes, 14, 1385–1406, https://doi.org/10.1002/1099-1085(20000615)14:8<1385::AID-HYP989>3.0.CO;2-1.
Huang, J., J. Zhang, Z. Zhang, S. Sun, and J. Yao, 2012: Simulation of extreme precipitation indices in the Yangtze River basin by using statistical downscaling method (SDSM). Theor. Appl. Climatol., 108, 325–343, https://doi.org/10.1007/s00704-011-0536-3.
IAEA, 2020: Global network of isotopes in precipitation. https://www.iaea.org/services/networks/gnip.
Ingraham, N. L., 1998: Isotopic variations in precipitation. Isotope Tracers in Catchment Hydrology, C. Kendall and J. J. McDonnell, Eds., Elsevier, 87–118, https://doi.org/10.1016/B978-0-444-81546-0.50010-0.
Jacobs, S. R., E. Timbe, B. Weeser, M. C. Rufino, K. Butterbach-Bahl, and L. Breuer, 2018: Assessment of hydrological pathways in East African montane catchments under different land use. Hydrol. Earth Syst. Sci., 22, 4981–5000, https://doi.org/10.5194/hess-22-4981-2018.
Jasechko, S., and R. Taylor, 2015: Intensive rainfall recharges tropical groundwaters. Environ. Res. Lett., 10, 124015, https://doi.org/10.1088/1748-9326/10/12/124015.
Jasechko, S., J. W. Kirchner, J. M. Welker, and J. J. McDonnell, 2016: Substantial proportion of global streamflow less than three months old. Nat. Geosci., 9, 126–129, https://doi.org/10.1038/ngeo2636.
Kanner, L. C., N. H. Buenning, L. D. Stott, A. Timmermann, and D. Noone, 2014: The role of soil processes in δ18O. Global Biogeochem. Cycles, 28, 239–252, https://doi.org/10.1002/2013GB004742.
Kendall, C., and J. J. McDonnell, 2012: Isotope Tracers in Catchment Hydrology. Elsevier, 870 pp.
Khedun, C. P., A. K. Mishra, V. P. Singh, and J. R. Giardino, 2014: A copula-based precipitation forecasting model: Investigating the interdecadal modulation of ENSO’s impacts on monthly precipitation. Water Resour. Res., 50, 580–600, https://doi.org/10.1002/2013WR013763.
Kirchner, J. W., 2016a: Aggregation in environmental systems—Part 1: Seasonal tracer cycles quantify young water fractions, but not mean transit times, in spatially heterogeneous catchments. Hydrol. Earth Syst. Sci., 20, 279–297, https://doi.org/10.5194/hess-20-279-2016.
Kirchner, J. W., 2016b: Aggregation in environmental systems—Part 2: Catchment mean transit times and young water fractions under hydrologic nonstationarity. Hydrol. Earth Syst. Sci., 20, 299–328, https://doi.org/10.5194/hess-20-299-2016.
Konecky, B. L., D. C. Noone, and K. M. Cobb, 2019: The influence of competing hydroclimate processes on stable isotope ratios in tropical rainfall. Geophys. Res. Lett., 46, 1622–1633, https://doi.org/10.1029/2018GL080188.
Krause, P., D. P. Boyle, and F. Bäse, 2005: Comparison of different efficiency criteria for hydrological model assessment. Adv. Geosci., 5, 89–97, https://doi.org/10.5194/adgeo-5-89-2005.
Kuhn, G., S. Khan, A. R. Ganguly, and M. L. Branstetter, 2007: Geospatial-temporal dependence among weekly precipitation extremes with applications to observations and climate model simulations in South America. Adv. Water Resour., 30, 2401–2423, https://doi.org/10.1016/j.advwatres.2007.05.006.
Laux, P., S. Wagner, A. Wagner, J. Jacobeit, A. Bárdossy, and H. Kunstmann, 2009: Modelling daily precipitation features in the Volta basin of West Africa. Int. J. Climatol., 29, 937–954, https://doi.org/10.1002/joc.1852.
Laux, P., S. Vogl, W. Qiu, H. R. Knoche, and H. Kunstmann, 2011: Copula-based statistical refinement of precipitation in RCM simulations over complex terrain. Hydrol. Earth Syst. Sci., 15, 2401–2419, https://doi.org/10.5194/hess-15-2401-2011.
Lee, J. E., and I. Fung, 2008: “Amount effect” of water isotopes and quantitative analysis of post-condensation processes. Hydrol. Processes, 22, 1–8, https://doi.org/10.1002/hyp.6637.
Lutz, S. R., R. Krieg, C. Müller, M. Zink, K. Knöller, L. Samaniego, and R. Merz, 2018: Spatial patterns of water age: Using young water fractions to improve the characterization of transit times in contrasting catchments. Water Resour. Res., 54, 4767–4784, https://doi.org/10.1029/2017WR022216.
McGuire, K. J., and J. J. McDonnell, 2006: A review and evaluation of catchment transit time modeling. J. Hydrol., 330, 543–563, https://doi.org/10.1016/j.jhydrol.2006.04.020.
Moerman, J. W., K. M. Cobb, J. F. Adkins, H. Sodemann, B. Clark, and A. A. Tuen, 2013: Diurnal to interannual rainfall δ18O variations in northern Borneo driven by regional hydrology. Earth Planet. Sci. Lett., 369–370, 108–109, https://doi.org/10.1016/j.epsl.2013.03.014.
Moore, M., Z. Kuang, and P. N. Blossey, 2014: A moisture budget perspective of the amount effect. Geophys. Res. Lett., 41, 1329–1335, https://doi.org/10.1002/2013GL058302.
Moore, M., P. N. Blossey, A. Muhlbauer, and Z. Kuang, 2016: Microphysical controls on the isotopic composition of wintertime orographic precipitation. J. Geophys. Res., 121, 7235–7253, https://doi.org/10.1002/2015JD023763.
Nusbaumer, J., T. E. Wong, C. Bardeen, and D. Noone, 2017: Evaluating hydrological processes in the Community Atmosphere Model Version 5 (CAM5) using stable isotope ratios of water. J. Adv. Model. Earth Syst., 9, 949–977, https://doi.org/10.1002/2016MS000839.
Orbe, C., and Coauthors, 2020: GISS model E2.2: A climate model optimized for the middle atmosphere—2. Validation of large-scale transport and evaluation of climate response. J. Geophys. Res. Atmos., 125, e2020JD033151, https://doi.org/10.1029/2020JD033151.
Pfahl, S., and H. Sodemann, 2014: What controls deuterium excess in global precipitation? Climate Past, 10, 771–781, https://doi.org/10.5194/cp-10-771-2014.
Poduje, A. C. C., and U. Haberlandt, 2017: Short time step continuous rainfall modeling and simulation of extreme events. J. Hydrol., 552, 182–197, https://doi.org/10.1016/j.jhydrol.2017.06.036.
Remondi, F., J. W. Kircher, P. Burlando, and S. Fatichi, 2018: Water flux tracking with a distributed hydrologic model to quantify controls on the spatio-temporal variability of transit time distributions. Water Resour. Res., 54, 3081–3099, https://doi.org/10.1002/2017WR021689.
Risi, C., S. Bony, and F. Vimeux, 2008: Influence of convective processes on the isotopic composition (δ18O and δD) of precipitation and water vapor in the tropics: 2. Physical interpretation of the amount effect. J. Geophys. Res., 113, D19306, https://doi.org/10.1029/2008JD009943.
Risi, C., and Coauthors, 2012: Process-evaluation of tropospheric humidity simulated by general circulation models using water vapor isotopic observations: 2. Using isotopic diagnostics to understand the mid and upper tropospheric moist bias in the tropics and subtropics. J. Geophys. Res., 117, D05304, https://doi.org/10.1029/2011JD016623.
Rosa, D., J. F. Lamarque, and W. D. Collins, 2012: Global transport of passive tracers in conventional and superparameterized climate models: Evaluation of multi-scale methods. J. Adv. Model. Earth Syst., 4, M10003, https://doi.org/10.1029/2012MS000206.
Rozanski, K., L. Araguás-Araguás, and R. Gonfiantini, 1993: Isotopic patterns in modern global precipitation. Climate Change in Continental Isotopic Records, Geophys. Monogr., Vol. 78, Amer. Geophys. Union, 1–36, https://doi.org/10.1029/GM078p0001.
Schneider, M., and F. Ramos, 2014: Transductive learning for multi-task copula processes. arXiv, 8 pp., https://arxiv.org/abs/1406.0304.
Schölzel, C., and P. Friederichs, 2008: Multivariate non-normally distributed random variables in climate research - Introduction to the copula approach. Nonlinear Processes Geophys., 15, 761–772, https://doi.org/10.5194/npg-15-761-2008.
Sklar, A., 1959: Fonctions de répartition à n dimensions et leurs marges. Publ. Inst. Stat. Univ. Paris, 8, 229–231.
So, B., J. Kim, H. Kwon, and C. H. R. Lima, 2017: Stochastic extreme downscaling model for an assessment of changes in rainfall intensity-duration-frequency curves over South Korea using multiple regional climate models. J. Hydrol., 553, 321–337, https://doi.org/10.1016/j.jhydrol.2017.07.061.
Song, C., G. Wang, G. Liu, T. Mao, X. Sun, and X. Chen, 2017: Stable isotope variations of precipitation and streamflow reveal the young water fraction of a permafrost watershed. Hydrol. Processes, 31, 935–947, https://doi.org/10.1002/hyp.11077.
Sprenger, M., and Coauthors, 2019: The demographics of water: A review of water ages in the critical zone. Rev. Geophys., 57, 800–834, https://doi.org/10.1029/2018RG000633.
Steen-Larsen, H. C., C. Risi, M. Werner, K. Yoshimura, and V. Masson-Delmotte, 2017: Evaluating the skills of isotope-enabled general circulation models against in situ atmospheric water vapor isotope observations. J. Geophys. Res. Atmos., 122, 246–263, https://doi.org/10.1002/2016JD025443.
Tharammal, T., G. Bala, and D. Noone, 2017: Impact of deep convection on the isotopic amount effect in tropical precipitation. J. Geophys. Res. Atmos., 122, 1505–1523, https://doi.org/10.1002/2016JD025555.
Turnadge, C., and B. D. Smerdon, 2014: A review of methods for modelling environmental tracers in groundwater: Advantages of tracer concentration simulation. J. Hydrol., 519, 3674–3689, https://doi.org/10.1016/j.jhydrol.2014.10.056.
Vachon, R. W., J. W. C. White, E. Gutmann, and J. M. Welker, 2007: Amount-weighted annual isotopic (δ18O) values are affected by the seasonality of precipitation: A sensitivity study. Geophys. Res. Lett., 34, L21707, https://doi.org/10.1029/2007GL030547.
von Freyberg, J., S. T. Allen, S. Seeger, M. Weiler, and J. W. Kirchner, 2018: Sensitivity of young water fractions to hydroclimatic forcing and landscape properties across 22 Swiss catchments. Hydrol. Earth Syst. Sci., 22, 3841–3861, https://doi.org/10.5194/hess-22-3841-2018.
West, J. B., G. J. Bowen, T. E. Dawson, and K. P. Tu, 2010: Isoscapes: Understanding Movement, Pattern, and Process on Earth through Isotope Mapping. Springer, 487 pp.
Wiederhold, J. G., 2015: Metal stable isotope signatures as tracers in environmental geochemistry. Environ. Sci. Technol., 49, 2606–2624, https://doi.org/10.1021/es504683e.
Wilkinson, B. H., and L. C. Ivany, 2002: Paleoclimatic inference from stable isotope profiles of accretionary biogenic hardparts – A quantitative approach to the evaluation of incomplete data. Palaeogeogr. Palaeoclimatol. Palaeoecol., 185, 95–114, https://doi.org/10.1016/S0031-0182(02)00279-1.
Wong, T. E., J. Nusbaumer, and D. C. Noone, 2017: Evaluation of modeled land-atmosphere exchanges with a comprehensive water isotope fractionation scheme in version 4 of the Community Land Model. J. Adv. Model. Earth Syst., 9, 978–1001, https://doi.org/10.1002/2016MS000842.
Yang, W., A. Bárdossy, and H. J. Caspary, 2010: Downscaling daily precipitation time series using a combined circulation- and regression-based approach. Theor. Appl. Climatol., 102, 439–454, https://doi.org/10.1007/s00704-010-0272-0.