Innovative Environmental Tracer Techniques for Evaluating Sources of Spring Discharge from a Carbonate Aquifer Bisected by a River

Littlefield Springs discharge about 1.6 m³/s along a 10‐km reach of the Virgin River in northwestern Arizona. Understanding their source is important for salinity control in the Colorado River Basin. Environmental tracers suggest that Littlefield Springs are a mixture of older groundwater from the regional Great Basin carbonate aquifer and modern (post‐1950s) seepage from the Virgin River. While corrected ¹⁴C apparent ages range from 1 to 9 ka, large amounts of nucleogenic ⁴He and low ³He/⁴He ratios suggest that the carbonate aquifer component is likely even older Pleistocene recharge. Modeled infiltration of precipitation, hydrogeologic cross sections, and hydraulic gradients all indicate recharge to the carbonate aquifer likely occurs in the Clover and Bull Valley Mountains along the northern part of the watershed, rather than in the nearby Virgin Mountains. This high‐altitude recharge is supported by relatively cool noble‐gas recharge temperatures and isotopically depleted δ²H and δ¹⁸O. Excess (crustal) SF₆ and ⁴He precluded dating of the modern component of water from Littlefield Springs using SF₆ and ³H/³He methods. Assuming a lumped‐parameter model with a binary mixture of two piston‐flow components, Cl^?/Br^?, Cl^?/F^?, δ²H, and CFCs indicate the mixture is about 60% Virgin River water and 40% groundwater from the carbonate aquifer, with an approximately 30‐year groundwater travel time for Virgin River seepage to re‐emerge at Littlefield Springs. This suggests that removal of high‐salinity sources upstream of the Virgin River Gorge would reduce the salinity of water discharging from Littlefield Springs into the Virgin River within a few decades.     

Geophysical constraints on deep weathering and water storage potential in the Southern Sierra Critical Zone Observatory

The conversion of bedrock to regolith marks the inception of critical zone processes, but the factors that regulate it remain poorly understood. Although the thickness and degree of weathering of regolith are widely thought to be important regulators of the development of regolith and its water‐storage potential, the functional relationships between regolith properties and the processes that generate it remain poorly documented. This is due in part to the fact that regolith is difficult to characterize by direct observations over the broad scales needed for process‐based understanding of the critical zone. Here we use seismic refraction and resistivity imaging techniques to estimate variations in regolith thickness and porosity across a forested slope and swampy meadow in the Southern Sierra Critical Zone Observatory (SSCZO). Inferred seismic velocities and electrical resistivities image a weathering zone ranging in thickness from 10 to 35 m (average = 23 m) along one intensively studied transect. The inferred weathering zone consists of roughly equal thicknesses of saprolite (P‐velocity < 2 km s^?1) and moderately weathered bedrock (P‐velocity = 2–4 km s^?1). A minimum‐porosity model assuming dry pore space shows porosities as high as 50% near the surface, decreasing to near zero at the base of weathered rock. Physical properties of saprolite samples from hand augering and push cores are consistent with our rock physics model when variations in pore saturation are taken into account. Our results indicate that saprolite is a crucial reservoir of water, potentially storing an average of 3 m³ m^?2 of water along a forested slope in the headwaters of the SSCZO. When coupled with published erosion rates from cosmogenic nuclides, our geophysical estimates of weathering zone thickness imply regolith residence times on the order of 10⁵ years. Thus, soils at the surface today may integrate weathering over glacial–interglacial fluctuations in climate. Copyright © 2013 John Wiley & Sons, Ltd.     

Identification of soil degradation during earthquake excitations by Bayesian inference

A Bayesian inference approach is introduced to identify soil degradation behaviours at four downhole array sites. The approach of inference is based on a parametric time‐varying infinite impulse response filter model. The approach is shown to be adaptive to the changes of filter parameters and noise amplitudes. Four sites, including the Lotung (Taiwan), Chiba (Japan), Garner Valley (California), and Treasure Island (California) sites with downhole seismic arrays are analysed. Our results show two major types of soil degradation behaviour: the well‐known strain‐dependent softening, and reduction in stiffness that is not instantaneously recoverable. It is also found that both types of soil degradation are more pronounced in sandy soils than in clayey soils. The mechanism for the second type of soil degradation is not yet clear to the authors and suggested to be further studied. Copyright © 2003 John Wiley & Sons, Ltd.     

Connectivity of the Apalachicola River flow variability and the physical and bio-optical oceanic properties of the northern West Florida Shelf

Maps of satellite-derived estimates of monthly averaged chlorophyll a concentration over the northern West Florida Shelf show interannual variations concentrated near the coastline, but also extending offshore over the shelf in a tongue-like pattern from the Apalachicola River during the late winter and early spring. These anomalies are significantly correlated with interannual variability in the flow rate of the Apalachicola River, which is linked to the precipitation anomalies over the watershed, over a region extending 150–200 km offshore out to roughly the 100 m isobath. This study examines the variability of the Apalachicola River and its impacts on the variability of water properties over the northern West Florida Shelf. A series of numerical model experiments show that episodic wind-driven offshore transport of the Apalachicola River plume is a likely physical mechanism for connecting the variability of the river discharge with oceanic variability over the middle and outer shelf.     

ENCORE: the effect of nutrient enrichment on coral reefs. Synthesis of results and conclusions

Coral reef degradation resulting from nutrient enrichment of coastal waters is of increasing global concern. Although effects of nutrients on coral reef organisms have been demonstrated in the laboratory, there is little direct evidence of nutrient effects on coral reef biota in situ. The ENCORE experiment investigated responses of coral reef organisms and processes to controlled additions of dissolved inorganic nitrogen (N) and/or phosphorus (P) on an offshore reef (One Tree Island) at the southern end of the Great Barrier Reef, Australia. A multi-disciplinary team assessed a variety of factors focusing on nutrient dynamics and biotic responses. A controlled and replicated experiment was conducted over two years using twelve small patch reefs ponded at low tide by a coral rim. Treatments included three control reefs (no nutrient addition) and three + N reefs (NH4Cl added), three + P reefs (KH2PO4 added), and three + N + P reefs. Nutrients were added as pulses at each low tide (ca twice per day) by remotely operated units. There were two phases of nutrient additions. During the initial, low-loading phase of the experiment nutrient pulses (mean dose = 11.5 microM NH4+; 2.3 microM PO4(-3)) rapidly declined, reaching near-background levels (mean = 0.9 microM NH4+; 0.5 microM PO4(-3)) within 2-3 h. A variety of biotic processes, assessed over a year during this initial nutrient loading phase, were not significantly affected, with the exception of coral reproduction, which was affected in all nutrient treatments. In Acropora longicyathus and A. aspera, fewer successfully developed embryos were formed, and in A. longicyathus fertilization rates and lipid levels decreased. In the second, high-loading, phase of ENCORE an increased nutrient dosage (mean dose = 36.2 microM NH4+; 5.1 microM PO4(-3)) declining to means of 11.3 microM NH4+ and 2.4 microM PO4(-3) at the end of low tide) was used for a further year, and a variety of significant biotic responses occurred. Encrusting algae incorporated virtually none of the added nutrients. Organisms containing endosymbiotic zooxanthellae (corals and giant clams) assimilated dissolved nutrients rapidly and were responsive to added nutrients. Coral mortality, not detected during the initial low-loading phase, became evident with increased nutrient dosage, particularly in Pocillopora damicornis. Nitrogen additions stunted coral growth, and phosphorus additions had a variable effect. Coral calcification rate and linear extension increased in the presence of added phosphorus but skeletal density was reduced, making corals more susceptible to breakage. Settlement of all coral larvae was reduced in nitrogen treatments, yet settlement of larvae from brooded species was enhanced in phosphorus treatments. Recruitment of stomatopods, benthic crustaceans living in coral rubble, was reduced in nitrogen and nitrogen plus phosphorus treatments. Grazing rates and reproductive effort of various fish species were not affected by the nutrient treatments. Microbial nitrogen transformations in sediments were responsive to nutrient loading with nitrogen fixation significantly increased in phosphorus treatments and denitrification increased in all treatments to which nitrogen had been added. Rates of bioerosion and grazing showed no significant effects of added nutrients. ENCORE has shown that reef organisms and processes investigated in situ were impacted by elevated nutrients. Impacts were dependent on dose level, whether nitrogen and/or phosphorus were elevated and were often species-specific. The impacts were generally sub-lethal and subtle and the treated reefs at the end of the experiment were visually similar to control reefs. Rapid nutrient uptake indicates that nutrient concentrations alone are not adequate to assess nutrient condition of reefs. Sensitive and quantifiable biological indicators need to be developed for coral reef ecosystems. The potential bioindicators identified in ENCORE should be tested in future research on coral reef/nutrient interactions. Synergistic and cumulative effects of elevated nutrients and other environmental parameters, comparative studies of intact vs. disturbed reefs, offshore vs. inshore reefs, or the ability of a nutrient-stressed reef to respond to natural disturbances require elucidation. An expanded understanding of coral reef responses to anthropogenic impacts is necessary, particularly regarding the subtle, sub-lethal effects detected in the ENCORE studies.     

Implications of model uncertainty for the mapping of hillslope‐scale soil erosion predictions

An integrated modelling approach (MIRSED) which utilizes the process‐based soil erosion model WEPP (Water Erosion Prediction Project) is presented for the assessment of hillslope‐scale soil erosion at five sites throughout England and Wales. The methodology draws upon previous uncertainty analysis of the WEPP hillslope soil erosion model by the authors to qualify model results within an uncertainty framework. A method for incorporating model uncertainty from a range of sources is discussed as a first step towards using and learning from results produced through the GLUE (Generalized Likelihood Uncertainty Estimation) technique. Results are presented and compared to available observed data, which illustrate that levels of uncertainty are significant and must be taken into account if a meaningful understanding of output from models such as WEPP is to be achieved. Furthermore, the collection of quality, observed data is underlined for two reasons: as an essential tool in the development of soil erosion modelling and also to allow further constraint of model uncertainty. Copyright © 2001 John Wiley & Sons, Ltd.     

The mid-ocean ridge axial valley as a steady-state neck

In this paper the mid-ocean ridge axial valley is modelled as a steady-state lithospheric neck in which lithospheric stretching balances lithospheric accretion. Conversely, the axial high is a steady-state lithospheric bulge. The lithosphere is modelled as a thin plate with a Newtonian rheology. It is shown that an axial valley will occur if the rate of viscosity increase away from the ridge axis is faster than the rate at which accretion decreases. An axial high will occur if the opposite condition holds. This is consistent with the observation that axial valleys occur at low spreading rates and axial highs at high spreading rates. By fitting our model to profiles across the Mid-Atlantic Ridge and the East Pacific Rise and assuming the lithospheric thickness at the ridge axis to be 5 km, we find accretion widths of 6–8 km. We find the width over which there is a significant increase in lithospheric viscosity to be also 6–8 km.     

Fluid ascent through the solid lithosphere and its relation to earthquakes

The Earth is continuously expelling gases and liquids from great depths—juvenile volatiles from the mantle and recycled metamorphic products. Some of these fluids ascend through liquid rock in volcanic processes, but others utilize fractures and faults as conduits through the solid lithosphere. The latter process may have a major influence on earthquakes, since fluids at near lithostatic pressures appear to be required to activate deep faults that would otherwise remain locked.Fluids can be driven upward through solid rock by buoyancy, but only if present in sufficient concentration to form large-scale domains occupying interconnected fracture porosity. A growing fluid domain becomes so mobilized only when it attains the critical vertical dimension required for hydrostatic instability. This dimension, depending on the ultimate compressive yield strength of the rock, may be as much as several kilometers.Any column of fluid ascending through fractures in the solid lithosphere from a prolific deep source must become organized into a vertical sequence of discrete domains, separated by fluid-pressure discontinuities. This is required because a continuous hydrostatic-fluid-pressure profile extending from an arbitrarily deep source to the surface cannot be permitted by the finite strength of rock. A vertically stacked sequence of domains allows the internal fluid-pressure profile to approximate the external rock-stress profile in a stepwise fashion. The pressure discontinuity below the base of the uppermost hydrostatic domain may be responsible for some occurrences of so-called anomalous geopressures. An ascending stream of fluid that percolates upward from a deep source through a column of domains must encounter a sequence of abrupt pressure decreases at the transitions between successive domains. If supercritical gases act as solvents, the dissolved substances may drop out of solution at such pressure discontinuities, resulting in a local concentration of minerals and other substances.At great depths, brittle fracture would normally be prevented by high pressure and temperature, with all excessive stress discharged by ductile flow. Rock strata invaded by an ascending fluid domain are weakened, however, because cracks generated or reactivated by the high-pressure fluid can support the overburden, greatly reducing internal friction. This reduction of strength may cause a previously stressed rock to fail, resulting in hydraulic shear fracture. Thus, earthquakes may be triggered by the buoyant migration of deep-source fluids.The actual timing of the failure that leads to such an earthquake may be determined by the relatively rapid inflation of a fluid domain and not by any significant increase in the probably much slower rate of regional tectonic strain. Many earthquake precursory phenomena may be secondary symptoms of an increase in pore-fluid pressure, and certain coseismic phenomena may result from the venting of high-pressure fluids when faults break the surface. Instabilities in the migration of such fluid domains may also contribute to or cause the eruption of mud volcanoes, magma volcanoes, and kimberlite pipes.