Very high-frequency and seasonal cave atmosphere PCO2 variability: Implications for stalagmite growth and oxygen isotope-based paleoclimate records

Cave air P_CO2 at two Irish sites varied dramatically on daily to seasonal timescales, potentially affecting the timing of calcite deposition and consequently climate proxy records derived from stalagmites collected at the same sites. Temperature-dependent biochemical processes in the soil control CO₂ production, resulting in high summer P_CO2 values and low winter values at both sites. Large Large-amplitude, high-frequency variations superimposed on this seasonal cycle reflect cave air circulation. Here we model stalagmite growth rates, which are controlled partly by CO₂ degassing rates from drip water, by considering both the seasonal and high-frequency cave air P_CO2 variations. Modeled hourly growth rates for stalagmite CC-Bil from Crag Cave in SW Ireland reach maxima in late December (0.063 μm h^− 1) and minima in late June/early July (0.033 μm h^− 1). For well-mixed ‘diffuse flow’ cave drips such as those that feed CC-Bil, high summer cave air P_CO2 depresses summer calcite deposition, while low winter P_CO2 promotes degassing and enhances deposition rates. In stalagmites fed by well-mixed drips lacking seasonal variations in δ¹⁸O, integrated annual stalagmite calcite δ¹⁸O is unaffected; however, seasonality in cave air P_CO2 may influence non-conservative geochemical climate proxies (e.g., δ¹³C, Sr/Ca). Stalagmites fed by ‘seasonal’ drips whose hydrochemical properties vary in response to seasonality may have higher growth rates in summer because soil air P_CO2 may increase relative to cave air P_CO2 due to higher soil temperatures. This in turn may bias stalagmite calcite δ¹⁸O records towards isotopically heavier summer drip water δ¹⁸O values, resulting in elevated calcite δ¹⁸O values compared to the ‘equilibrium’ values predicted by calcite–water isotope fractionation equations. Interpretations of stalagmite-based paleoclimate proxies should therefore consider the consequences of cave air P_CO2 variability and the resulting intra-annual variability in calcite deposition rates.     

Quantifying bank storage of variably saturated aquifers

Numerical simulations were conducted to quantify bank storage in a variably saturated, homogenous, and anisotropic aquifer abutting a stream during rising stream stage. Seepage faces and bank slopes ranging from 1/3 to 100/3 were simulated. The initial conditions were assumed steady-state flow with water draining toward the stream. Then, the stream level rose at a constant rate to the specified elevation of the water table given by the landward boundary condition and stayed there until the system reached a new steady state. This represents a highly simplified version of a real world hydrograph. For the specific examples considered, the following conclusions can be made. The volume of surface water entering the bank increased with the rate of stream level rise, became negligible when the rate of rise was slow, and approached a positive constant when the rate was large. Also, the volume decreased with the dimensionless parameter M (the product of the anisotropy ratio and the square of the domain's aspect ratio). When M was large (>10), bank storage was small because most pore space was initially saturated with ground water due to the presence of a significant seepage face. When M was small, the seepage face became insignificant and capillarity began to play a role. The weaker the capillary effect, the easier for surface water to enter the bank. The effect of the capillary forces on the volume of surface water entering the bank was significant and could not be neglected.     

Tectonic influences on ground water quality: insight from complementary methods

A study using multiple techniques provided insight into tectonic influences on ground water systems; the results can help to understand ground water systems in the tectonically active western United States and other parts of the world. Ground water in the San Bernardino Valley (Arizona, United States and Sonora, Mexico) is the main source of water for domestic use, cattle ranching (the primary industry), and the preservation of threatened and endangered species. To improve the understanding of ground water occurrence, movement, and sustainability, an investigation was conducted using a number of complementary methods, including major ion geochemistry, isotope hydrology, analysis of gases dissolved in ground water, aquifer testing, geophysics, and an examination of surface and subsurface geology. By combining information from multiple lines of investigation, a more complete picture of the basin hydrogeology was assembled than would have been possible using fewer methods. The results show that the hydrogeology of the San Bernardino Valley is markedly different than that of its four neighboring basins in the United States. The differences include water quality, chemical evolution, storage, and residence time. The differences result from the locally unique geology of the San Bernardino Valley, which is due to the presence of a magmatically active accommodation zone (a zone separating two regions of normal faults with opposite dips). The geological differences and the resultant hydrological differences between the San Bernardino Valley and its neighboring basins may serve as a model for the distinctive nature of chemical evolution of ground water in other basins with locally distinct tectonic histories.     

An evaluation of the health status of the Lavaca Bay, Texas ecosystem using Crassostrea virginica as the sentinel species

Locational risks for compromised ecosystem health for the eastern oysters (Crassostrea virginica) harvested from Lavaca Bay, Texas were estimated. Flow cytometric evaluation of variations in DNA content and the lysosomal destabilization assay were used for evaluation of genotoxicity and stress, respectively. Bayesian geo-statistical methods were utilized to estimate and evaluate spatial effects. For models with spatial risks, continuous surface maps of predicted parameter values were created to evaluate risk location. Lysosomal destabilization assay results were spatially oriented whereas flow cytometry results were fit best with the random effects model. While not spatially oriented, the highest levels of variations in DNA content were also present near industrial facilities. Locational risks of increased biomarkers of genotoxicity and stress in the eastern oyster (Crassostrea virginica) were increased with proximity to industrial facilities     

Polycyclic aromatic hydrocarbons in coastal sediments of southwest Taiwan: An appraisal of diagnostic ratios in source recognition

Fifty-seven surface sediment samples were collected from the coast of southwest Taiwan and analyzed for polycyclic aromatic hydrocarbons (PAHs). Concentrations of total PAHs (28 PAH compounds) ranged from 15 to 907 ng g⁻¹ dry weight. Diagnostic ratios showed that PAHs in the sediments of the Gaoping estuary were predominantly of petroleum origin, whereas sediments from the Kaohsiung coast contained principally combustion-derived PAHs. Principal component analysis indicated that emissions from automobiles and coal burning were the main sources of combustion-derived PAHs. The relatively high ratios of perylene/penta-aromatic PAH isomers in sediments from the Tainan coast and some off-shore stations on the Kaohsiung coast suggest a significant diagenetic PAH contribution. The study shows that certain diagnostic ratios are useful and sensitive in delineating the distribution of PAHs from specific sources in southwest Taiwan. The phenanthrene/anthracene ratio is a better indicator than the methylphenanthrenes/phenanthrene ratio for tracing petrogenic PAHs, and the benzo(a)anthracene/chrysene and indeno(1,2,3-c,d)pyrene/benzo(g,h,i)perylene ratios are more specific than the benzo(a)pyrene/benzo(e)pyrene and benzo(b)fluoranthcene/benzo(k)fluoranthcene ratios in distinguishing PAHs from various pyrogenic sources.     

Characterization of Damage in Sandstones along the Mojave Section of the San Andreas Fault: Implications for the Shallow Extent of Damage Generation

Following theoretical calculations that suggest shallow generation of rock damage during an earthquake rupture, we measure the degree of fracture damage in young sedimentary rocks from the Juniper Hills Formation (JHF) that were displaced 21 km along the Mojave section of San Andreas Fault (SAF) and were not exhumed significantly during their displacement. In exposures adjacent to the fault, the JHF typically displays original sedimentary fabrics and little evidence of bulk shear strain at the mesoscopic scale. The formation is, however, pervasively fractured at the microscopic scale over a zone that is about a 100 m wide on the southwest side of the SAF near Little Rock. The abundance of open fractures, the poor consolidation, and the shallow inferred burial depth imply that the damage was generated close to the surface of the Earth. The spatial correlation of this damage with a seismically active trace of the SAF suggests that it was generated by SAF slip events that by assumption were of a seismic nature throughout the displacement history of the JHF. Thus the JHF provides a very shallow upper bound for the generation of brittle damage in a seismic fault zone. The fracture fabric is characterized by preferred orientations of fractures that split grains between contact points and is consistent with overall deformation under directed compression. However, the available results cannot be used to distinguish between proposed off-fault damage mechanisms. Fracture orientations are compatible with a maximum compressive stress oriented at a high angle to the fault at about 10 m, and at a lower, more variable angle farther away from the fault. The fracture distribution and fabric are consistent with observations made of the microscale damage characteristics of the Hungry Valley Formation in the northwestern section of the SAF in the Mojave, and with previous observations of exhumed, ancestral strands of the SAF.     

Nonplanar Faults: Mechanics of Slip and Off-fault Damage

Stress interactions and sliding characteristics of faults with random fractal waviness in a purely elastic medium differ both qualitatively and quantitatively from those of faults with planar surfaces. With nonplanar fault models, solutions for slip diverge as resolution of the fractal features increases, and the scaling of fault slip with fault rupture dimension becomes nonlinear. We show that the nonlinear scaling of slip and divergence of solutions arise because stresses from geometric interactions at irregularities along nonplanar faults grow with increasing slip and produce backstresses that progressively impede slip. However, in real materials with finite strength, yielding will halt the growth of the interaction stresses, which will profoundly affect slip of nonplanar faults. We infer that in the brittle seismogenic portion of the Earth’s crust, off-fault yielding occurs on pervasive secondary faults. Predicted rates of stress relaxation with distance from major faults with random fractal roughness follow a power-law relationship that is consistent with reported clustering of background seismicity up to 15 kilometers from faults.