Beaver dam analogues drive heterogeneous groundwater–surface water interactions

Beaver dam analogues (BDAs) are a cost-effective stream restoration approach that leverages the recognized environmental benefits of natural beaver dams on channel stability and local hydrology. Although natural beaver dams are known to exert considerable influence on the hydrologic conditions of a stream system by mediating geomorphic processes, nutrient cycling, and groundwater–surface water interactions, the impacts of beaver-derived restoration methods on groundwater–surface water exchange are poorly characterized. To address this deficit, we monitored hyporheic exchange fluxes and streambed porewater biogeochemistry across a sequence of BDAs installed along a central Wyoming stream during the summer of 2019. Streambed fluxes were quantified by heat tracing methods and vertical hydraulic gradients. Biogeochemical activity was evaluated using major ion porewater chemistry and principal component analysis. Vertical fluxes of approximately 1.0 m/day were observed around the BDAs, as was the development of spatially heterogeneous zones of nitrate production, groundwater upwelling, and anaerobic reduction. Strong contrasts in hyporheic zone processes were observed across BDAs of differing sizes. This suggests that structures may function with size-dependent behaviour, only altering groundwater–surface water interactions after a threshold hydraulic step height is exceeded. Patterns of hyporheic exchange and biogeochemical cycling around the studied BDAs resemble those around natural beaver dams, suggesting that BDAs may provide comparable benefits to channel complexity and near-stream function over a 1-year period.     

Local ice caps in Finderup Land,North Greenland,survived the Holocene Thermal Maximum

Local glaciers and ice caps (GICs) comprise only ~5.4% of the total ice volume, but account for ~14–20% of the current ice loss in Greenland. The glacial history of GICs is not well constrained, however, and little is known about how they reacted to Holocene climate changes. Specifically, in North Greenland, there is limited knowledge about past GIC fluctuations and whether they survived the Holocene Thermal Maximum (HTM, ~8 to 5 ka). In this study, we use proglacial lake records to constrain the ice‐marginal fluctuations of three local ice caps in North Greenland including Flade Isblink, the largest ice cap in Greenland. Additionally, we have radiocarbon dated reworked marine molluscs in Little Ice Age (LIA) moraines adjacent to the Flade Isblink, which reveal when the ice cap was smaller than present. We found that outlet glaciers from Flade Isblink retreated inland of their present extent from ~9.4 to 0.2 cal. ka BP. The proglacial lake records, however, demonstrate that the lakes continued to receive glacial meltwater throughout the entire Holocene. This implies that GICs in Finderup Land survived the HTM. Our results are consistent with other observations from North Greenland but differ from locations in southern Greenland where all records show that the local ice caps at low and intermediate elevations disappeared completely during the HTM. We explain the north–south gradient in glacier response as a result of sensitivity to increased temperature and precipitation. While the increased temperatures during the HTM led to a complete melting of GICs in southern Greenland, GICs remained in North Greenland probably because the melting was counterbalanced by increased precipitation due to a reduction in Arctic sea‐ice extent and/or increased poleward moisture transport.     

Dynamic modelling of retrogressive landslides with emphasis on the role of clay sensitivity

This paper presents a detailed numerical study of the retrogressive failure of landslides in sensitive clays. The dynamic modelling of the landslides is carried out using a novel continuum approach, the particle finite element method, complemented with an elastoviscoplastic constitutive model. The multiwedge failure mode in the collapse is captured successfully, and the multiple retrogressive failures that have been widely observed in landslides in sensitive clays are reproduced with the failure mechanism, the kinematics, and the deposition being discussed in detail. Special attention has been paid to the role of the clay sensitivity on each retrogressive failure as well as on the final retrogression distance and the final run‐out distance via parametric studies. Moreover, the effects of the viscosity of sensitive clays on the failure are also investigated for different clay sensitivities.     

On the Nonlinear Thermal Evolution of Solar Corona Structures

The thermal evolution of structures is investigated for different values of the size of structure. A simplified cooling function and a constant heating mechanism are assumed. Analytical criteria for thermal instability are obtained. It is found that the response of the thermal structure not only depends on the amplitude of the disturbances, but also on whether the disturbance increases or decreases the initial steady temperature. Additionally, the evolution of the structure is examined numerically by using a time-dependent code under the assumption that the inertia terms are small. In particular, the analytical criteria obtained for thermal instability are verified.     

On the self-consistency of evolutionary synthesis models

Evolutionary synthesis models have been used to study the physical properties of unresolved populations in a wide range of scenarios. Unfortunately, their self-consistency is difficult to test and there are some theoretical open questions without an answer: (1) The change of the homology relations assumed in the computation of isochrones due to the effect of stellar winds (or rotation) and the discontinuities in the stellar evolution are not considered. (2) There is no consensus about how the isochrones must be integrated. (3) The discreteness of the stellar populations (that produce an intrinsic statistical dispersion) usually are not taken into account, and model results are interpreted in a deterministic way instead of a statistical one. The objective of this contribution is to present some inconsistencies in the computation and some cautions in the application of the results of such codes. This revised version was published online in August 2006 with corrections to the Cover Date.     

Heavy metals in sediments of the Tecate River, Mexico

Ten sites along the Tecate River, Mexico were sampled to evaluate the cadmium, lead, nickel and chromium concentrations in sediments. The result shows contamination for cadmium in most of the sites, where two sites were class 4 (polluted to strongly polluted) according to geoaccumulation index proposed by Muller. Two sites were found polluted for all the heavy metals analyzed (Cr, Cd, Pb and Ni), indicating the effect of anthropogenic activities. A correlation between Ni and Cd concentration had been found indicating a common source. These metals are usually used in electroplating industry. The results of this study can be used for decision makers to prioritize measures to control the pollution for these metals.     

Quick-look assessments to identify optimal CO2 EOR storage sites

A newly developed, multistage quick-look methodology allows for the efficient screening of an unmanageably large number of reservoirs to generate a workable set of sites that closely match the requirements for optimal CO₂ enhanced oil recovery (EOR) storage. The objective of the study is to quickly identify miscible CO₂ EOR candidates in areas that contain thousands of reservoirs and to estimate additional oil recovery and sequestration capacities of selected top options through dimensionless modeling and reservoir characterization. Quick-look assessments indicate that the CO₂ EOR resource potential along the US Gulf Coast is 4.7 billion barrels, and CO₂ sequestration capacity is 2.6 billion metric tons. In the first stage, oil reservoirs are screened and ranked in terms of technical and practical feasibility for miscible CO₂ EOR. The second stage provides quick estimates of CO₂ EOR potential and sequestration capacities. In the third stage, a dimensionless group model is applied to a selected set of sites to improve the estimates of oil recovery and storage potential using appropriate inputs for rock and fluid properties, disregarding reservoir architecture and sweep design. The fourth stage validates and refines the results by simulating flow in a model that describes the internal architecture and fluid distribution in the reservoir. The stated approach both saves time and allows more resources to be applied to the best candidate sites.     

Seal integrity and feasibility of CO<Subscript>2</Subscript> sequestration in the Teapot Dome EOR pilot: geomechanical site characterization

This paper reports a preliminary investigation of CO₂ sequestration and seal integrity at Teapot Dome oil field, Wyoming, USA, with the objective of predicting the potential risk of CO₂ leakage along reservoir-bounding faults. CO₂ injection into reservoirs creates anomalously high pore pressure at the top of the reservoir that could potentially hydraulically fracture the caprock or trigger slip on reservoir-bounding faults. The Tensleep Formation, a Pennsylvanian age eolian sandstone is evaluated as the target horizon for a pilot CO₂ EOR-carbon storage experiment, in a three-way closure trap against a bounding fault, termed the S1 fault. A preliminary geomechanical model of the Tensleep Formation has been developed to evaluate the potential for CO₂ injection inducing slip on the S1 fault and thus threatening seal integrity. Uncertainties in the stress tensor and fault geometry have been incorporated into the analysis using Monte Carlo simulation. The authors find that even the most pessimistic risk scenario would require ∼10 MPa of excess pressure to cause the S1 fault to reactivate and provide a potential leakage pathway. This would correspond to a CO₂ column height of ∼1,500 m, whereas the structural closure of the Tensleep Formation in the pilot injection area does not exceed 100 m. It is therefore apparent that CO₂ injection is not likely to compromise the S1 fault stability. Better constraint of the least principal stress is needed to establish a more reliable estimate of the maximum reservoir pressure required to hydrofracture the caprock.     

Mineralogy and geochemistry of the Texeo Cu–Co mine site (NW Spain): screening tools for environmental assessment

The Cu–Co–Ni Texeo mine has been the most important source of Cu in NW Spain since Roman times and now, approximately 40,000 m³ of wastes from mine and metallurgical operations, containing average concentrations of 9,263 mg kg⁻¹ Cu, 1,100 mg kg⁻¹ As, 549 mg kg⁻¹ Co, and 840 mg kg⁻¹ Ni, remain on-site. Since the cessation of the activity, the abandoned works, facilities and waste piles have been posing a threat to the environment, derived from the release of toxic elements. In order to assess the potential environmental pollution caused by the mining operations, a sequential sampling strategy was undertaken in wastes, soil, surface and groundwater, and sediments. First, screening field tools were used to identify hotspots, before defining formal sampling strategies; so, in the areas where anomalies were detected in a first sampling stage, a second detailed sampling campaign was undertaken. Metal concentrations in the soils are highly above the local background, reaching up to 9,921 mg kg⁻¹ Cu, 1,373 mg kg⁻¹ As, 685 mg kg⁻¹ Co, and 1,040 mg kg⁻¹ Ni, among others. Copper concentrations downstream of the mine works reach values up to 1,869 μg l⁻¹ and 240 mg kg⁻¹ in surface water and stream sediments, respectively. Computer-based risk assessment for the site gives a carcinogenic risk associated with the presence of As in surface waters and soils, and a health risk for long exposures; so, trigger levels of these elements are high enough to warrant further investigation.     

Experimental investigation of the effects of temperature and ionic strength on Mo isotope fractionation during adsorption to manganese oxides

The Mo stable isotope system is being applied to study changes in ocean redox. Such applications implicitly assume that Mo isotope fractionation in aqueous systems is relatively insensitive to frequently changing environmental variables such as temperature (T) and ionic strength (I). A major driver of fractionation is the adsorption of Mo to Mn oxyhydroxide surfaces [Barling J. and Anbar A. D. (2004) Molybdenum isotope fractionation during adsorption by manganese oxides. Earth Planet. Sci. Lett.217(3-4), 315-329]. Here, we report the results of experiments that determine the extent to which Mo isotope fractionation during adsorption of Mo to the Mn oxyhydroxide mineral birnessite is sensitive to T and I. The results are compared to new predictions from quantum chemical computations. We measured fractionation from 1 to 50 °C at I = 0.1 m and found that Δ^97/95Mo_{dissolved-adsorbed} varies from 1.9‰ to 1.6‰ over this temperature range. Experiments were also performed at 25 °C in synthetic seawater (I = 0.7); fractionation at this condition was the same within analytical error as in low ionic strength experiments. These findings confirm that the Mo isotope fractionation during adsorption to Mn oxyhydroxides is relatively insensitive to variations and T and I over environmentally relevant ranges. To relate these findings to potential mechanisms of Mo isotope fractionation, we also report results for density functional theory computations of the fractionation between and various possible structures of molybdic acid as a function of temperature. Because no plausible species fractionates from with a magnitude matching the experiments, we are left with three possibilities to explain the fractionation: (1) solvation effects on the vibrational frequencies of aqueous species considered thus far are significant, such that our calculations in vacuo yield inaccurate fractionations; (2) a trace aqueous species not yet considered fractionates from and then adsorbs to birnessite; or (3) a surface complex not present in solution forms on birnessite in which Mo is not tetrahedrally coordinated. Our findings help validate assumptions underlying paleoceanographic applications of the Mo isotope system and also lead us closer to understanding the mechanism of isotope fractionation during adsorption of Mo to Mn oxyhydroxides.