Assessing basin storage: Comparison of hydrometric- and tracer-based indices of dynamic and total storage

Storage is a fundamental but elusive component of drainage basin function, influencing synchronization between precipitation input and streamflow output and mediating basin sensitivity to climate and land use/land cover (LULC) change. We compare hydrometric and isotopic approaches to estimate indices of dynamic and total basin storage, respectively, and assess inter-basin differences in these indices across the Oak Ridges Moraine (ORM) region of southern Ontario, Canada. Dynamic storage indices for the 20 study basins included the ratio of baseflow to total streamflow (baseflow index BFI), Q ₉₉ flow and flow duration curve (FDC) slope. Ratios of the standard deviation of the streamflow stable isotope signal relative to that of precipitation were determined for each basin from a 1 year bi-weekly sampling program and used as indicators of total storage. Smaller ratios imply longer water travel times, smaller young water fractions (F _yw, < ~2–3 months in age) in streamflow and greater basin storage. Ratios were inversely related to BFI and Q ₉₉, and positively related to FDC slope, suggesting longer travel times and smaller F _yw for basins with stable baseflow-dominated streamflow regimes. Inter-basin differences in all indices reflected topographic, hydrogeologic and LULC controls on storage, which was greatest in steep, forest-covered headwaters underlain by permeable deposits with thick and relatively uniform unsaturated zones. Nevertheless, differential sensitivity of indices to controls on storage indicates the value of using several indices to capture more completely how basin characteristics influence storage. Regression relationships between storage indices and basin characteristics provided reasonable predictions of aspects of the streamflow regime of test basins in the ORM region. Such relationships and the underlying knowledge of controls on basin storage in this landscape provide the foundation for initial predictions of relative differences in streamflow response to regional changes in climate and LULC.     

Mercury cycling during acid rain recovery at the forested Lesní potok catchment,Czech Republic

From 2011 to 2019, mercury (Hg) stores and fluxes were studied in the small forested catchment Lesní potok (LES) in the central Czech Republic using the watershed mass balance approach together with internal measurements. Mean input fluxes of Hg via open bulk deposition, beech throughfall and spruce throughfall during the periodwere 2.9, 3.9 and 7.6 μg m⁻² year⁻¹, respectively. These values were considerably lower than corresponding deposition Hg fluxes reported in the early years of the 21st century from catchments in Germany. Current bulk precipitation inputs at unimpacted Czech mountainous sites were lower than those in Germany. The largest Hg inputs to the catchment were via litterfall, averaging 22.6 and 17.8 μg m⁻² year⁻¹ for beech and spruce stands. The average Hg input, based on the sum of mean litterfall and throughfall deposition, was 23.0 μg m⁻² year⁻¹, compared to the estimated Hg output in runoff of 0.5 μg m⁻² year⁻¹, which is low compared to other reported values. Thus, only ~2% of Hg input is exported in stream runoff. Stream water Hg was only weakly related to dissolved organic carbon (DOC) but both concentrations were positively correlated with water temperature. The estimated total soil Hg pool averaged 47.5 mg m⁻², only 4% of which was in the O-horizon. Thus Hg in the O-horizon pool represents 72 years of deposition at the current input flux and 3800 years of export at the current runoff flux. Age-dating by ¹⁴C suggested that organic soil contains Hg from recent deposition, while mineral soil at 40–80 cm depth contained 4400-year old carbon, suggesting the soil had accumulated atmospheric Hg inputs through millennia to reach the highest soil Hg pool of the soil profile. These findings suggest that industrial era intensification of the Hg cycle is superimposed on a slower-paced Hg cycle during most of the Holocene.     

Detrital zircon age spectra of middle and upper Eocene outcrop belts,U.S. Gulf Coast region

Recently reported detrital zircon (DZ) data help to associate the Paleogene strata of the Gulf of Mexico region to various provenance areas. By far, recent work has emphasised upper Paleocene‐lower Eocene and upper Oligocene strata that were deposited during the two episodes of the highest sediment supply in the Paleogene. The data reveal a dynamic drainage history, including (1) initial routing of western Cordilleran drainages towards the Gulf of Mexico in the Paleocene, (2) an eastward shift of the western continental divide, from the Jura‐Cretaceous cordilleran arc to the eastern edge of the Laramide province after the Paleocene and (3) a southward shift, along the eastern Laramide province, of the headwaters of river systems draining to the Mississippi and Houston embayments at some time between the early Eocene and Oligocene. However, DZ characterisation of most (~20 Myr) of the middle Eocene‐lower Oligocene section remains limited. We present 60 DZ age spectra, most of which are from the middle or upper Eocene outcrop belts, with 50–200‐km spacing. We define six to eight distinct groups of DZ age spectra for middle and upper Eocene strata. Data from this and other studies resolve at least six substantial temporal changes in age spectra at various positions along the continental margin. The evolving age spectra constrain the middle and upper Eocene drainage patterns of large parts of interior North America. The most well‐resolved aspects of these drainage patterns include (1) persistent rivers that flowed from erosional landscapes across the Paleozoic Appalachian orogen either into the low‐lying Mississippi embayment or directly into the eastern Gulf; (2) at least during marine regressions, a trunk channel that likely flowed southward along the axial part of Mississippi Embayment and integrated tributaries from the east and west; and (3) rivers that flowed to the Houston embayment in the middle Eocene that likely originated in the Laramide province in central Colorado and southern Wyoming, as Precambrian basement highs in those source areas were being unroofed.     

Field testing cellulose-water oxygen isotope relations in periphyton for paleohydrological reconstructions

Journal of Paleolimnology - The oxygen isotope composition (δ18O) of aquatic cellulose extracted from lake sediment cores has long been used to reconstruct past hydrological responses to...     

A Novel Approach for Evaluating Nonstationary Response of Dynamic Systems to Stochastic Excitation

The transient state of a dynamic system, such as offshore structures, to random excitation is always nonstationary. Many studies have contributed to evaluating response covariances at the transient state of a linear multi-degree-of-freedom(MDOF) system to random excitations, but a closed-form solution was not available unless the excitation was assumed to be a physically unrealizable white noise process. This study derives explicit, closed-form solutions for the response covariances at the transient state by using a pole-residue(PR) approach operated in the Laplace domain when the excitations are assumed to be stationary random processes described by physically realizable spectral density functions. By using the PR method, we can analytically solve the triple integral in evaluating the nonstationary response covariance. As this approach uses the poles and residues of system transfer functions, rather than the conventional mode superposition technique, the method is applicable to MDOF systems with non-classical damping models. Particular application of the proposed method is demonstrated for multi-story shear buildings to stochastic ground acceleration characterized by the Kanai–Tajimi spectral density function model, and a numerical example is provided to illustrate the detailed steps. No numerical integrations are required for computing the response covariances as the exact closed-form solution has been derived. The correctness of the proposed method is numerically verified by Monte Carlo simulations.     

Crossing Boundaries in a Collaborative Modeling Workspace

There is substantial literature on the importance of bridging across disciplinary and science–management boundaries. One of the ways commonly suggested to cross boundaries is for participants from both sides of the boundary to jointly produce information (i.e., knowledge co-production). But simply providing tools or bringing people together in the same room is not sufficient. Here we present a case study documenting the mechanisms by which managers and scientists collaborated to incorporate climate change projections into Colorado’s State Wildlife Action Plan. A critical component of the project was the use of a collaborative modeling and visualization workspace: the U.S. Geological Survey’s Resource for Advanced Modeling (RAM). Using video analysis and pre/post surveys from this case study, we examine how the RAM facilitated cognitive and social processes that co-produced a more salient and credible end product. This case provides practical suggestions to scientists and practitioners who want to implement actionable science.     

Planning for Sea Level Rise and Shore Protection Under Climate Uncertainty

Attention is focused here on the effect of additional sources of uncertainty derived from climate change on the cost-benefit procedures applied by coastal planners to evaluate shoreline protection projects. The largest effect would be felt if planners were trying to achieve the first best economic optimum. Given the current view that the seas will rise by significantly less than one meter through the year 2100, present procedures should work reasonably well assuming (1) informed vigilance in monitoring the pace of future greenhouse induced sea level rise, (2) careful attention to the time required for market-based adaptation to minimize the economic cost of abandonment, and (3) firm support of the credibility of an announced policy to proceed with plans to retreat from the sea when warranted. Assumptions (1) and (2) might be satisfied in reality, even cursory review of existing policy makes it clear that meeting (3) is a "long shot" at the very best. In any case, planners should periodically revisit potential protection sites, especially in the wake of catastrophic events, to assess the impact of the most recent information on sea level rise trajectories, local development patterns, and protection costs on the decision calculus.     

Evidence for Holocene environmental change from C/N ratios, and δ13C and δ15N values in Swan Lake sediments, western Sand Hills, Nebraska

Profiles of percent carbon and nitrogen, carbon/nitrogen (C/N) ratios and stable carbon (¹³C), and nitrogen (¹⁵N) isotopic ratios in organic matter from an 11.6 m core were used to reconstruct environments of deposition in the Swan Lake basin during the past 5300 YBP. The upper 6.5 m consisted of gyttja containing variable amounts of reddish brown-colored fine organic matter and calcium carbonate. It was followed by a 0.5 m sandy silt, which was followed by a 3.6 m reduced layer characterized by large quantities of black organic plant remains, sapropel, and then by another sapropel layer consisting mainly of well-sorted sapropelic sand with relatively low organic matter content. The C- and N-contents in the organic matter in the sediment profile ranged from 0.5 to 23% and from 0.02 to 2%, respectively. Carbon content were positively correlated to both N and clay content while carbon content was negatively correlated to sand content. Two major environmental phases in Swan Lake were apparent from large differences in the C and N data of the sediment organic matter. These include the sapropel (marsh) stage that stretched from approximately 5330 to 3930 YBP, and the following gyttja (open water stage). During the sapropel marsh plants identified in a previous pollen study as cattails and sedges proliferated and produced copious amounts of well-preserved organic matter. C/N ratios, ¹³C values, and ¹⁵N values in the sapropel were significantly different from those that characterized organic matter in the gyttja. During the gyttja ¹³C values indicated that deep primary producers have dominated lake biomass. By utilizing bicarbonate as their C-source, the accumulating biomass became relatively enriched ¹³C values. The presence of high sediment CaCO₃ contents indicated more alkaline and deeper water conditions prevailed during the gyttja. Further refinement of the data suggested that each major phase initially contained an identifiable transition stage. During the sapropelic (initial marsh stage) which occurred before 5330 YBP, sand content gradually decreased as organic matter increased. As reflected by high C/N ratios and slightly enriched ¹³C values, these sands appear to have contained sufficient permeability to promote partial mineralization of accumulated organic-N containing compounds. A short initial gyttja transition period from about 3930–3830 YBP occurred in which the sediment silt content was anomalously high relative that measured in the surrounding layers. The silt content suggests that this turbid transition layer can not be completely explained by sediment mixing via bioturbation. The silts appeared to have been associated with the sharp climate change that resulted in higher water-table conditions during the gyttja stage.