The volume and paleobathymetry of glacial Lake Agassiz

During the last retreat of the Laurentide Ice Sheet in North America, many proglacial lakes formed as continental drainage was impounded against the southern and western ice margin. Lake Agassiz was the largest of these lakes. The bathymetry of Lake Agassiz at the Herman and Upper Campbell beach levels – formed at about 11.5–11.0 ka and 9.9–9.5 ka, respectively – was computer modelled in this study by first collecting data for the isostatically-deformed paleowater planes of the two lake levels (derived from isobase lines constructed from beach elevations), and then subtracting these from the modern topography of the former lake floor. Pixels with dimensions of 1/30 × 1/30 of a degree were used in the model. Using these data, the area and volume of the lake were also calculated: at the Herman level these were 152 500 km2 and 13 100 km3 respectively; at the Upper Campbell level these were 350 400 km2 and 38 700 km3. Contour maps showing the paleobathymetry of both periods in the lake's history were also constructed. Determining the paleobathymetry and volume of Lake Agassiz is an important step in understanding the impact that the lake had on its surrounding environment and on the rivers, lakes, and oceans into which it flowed.     

Late Pleistocene and Holocene glaciation of the Fish Lake valley,northeastern Alaska Range,Alaska

We reconstructed a chronology of glaciation spanning from the Late Pleistocene through the late Holocene for Fish Lake valley in the north‐eastern Alaska Range using ¹⁰Be surface exposure dating and lichenometry. After it attained its maximum late Wisconsin extent, the Fish Lake valley glacier began to retreat ca. 16.5 ka, and then experienced a readvance or standstill at 11.6 ± 0.3 ka. Evidence of the earliest Holocene glacial activity in the valley is a moraine immediately in front of Little Ice Age (LIA) moraines and is dated to 3.3–3.0 ka. A subsequent advance culminated at ca. AD 610–900 and several LIA moraine crests date to AD 1290, 1640, 1860 and 1910. Our results indicate that ¹⁰Be dating from high‐elevation sites can be used to help constrain late Holocene glacial histories in Alaska, even when other dating techniques are unavailable. Close agreement between ¹⁰Be and lichenometric ages reveal that ¹⁰Be ages on late Holocene moraines may be as accurate as other dating methods. Copyright © 2009 John Wiley & Sons, Ltd.     

Using the precision of the mean to estimate suitable sample sizes for monitoring total phosphorus in Australian catchments

The effect of the sample size on prediction quality is well understood. Recently, studies have assessed this relationship using near‐continuous water quality samples. However, this is rarely possible because of financial constraints, and therefore, many studies have relied on simulation‐based methods utilizing more affordable surrogates. A limitation of simulation‐based methods is the requirement of a good relationship, which is often not present. Therefore, catchment managers require a direct method to estimate the effect of sample size on the mean using historical water quality data. One measure of prediction quality is the precision with which a mean is estimated; this is the focus of this work. By characterizing the effect of sample size on the precision of the mean, it is possible for catchment managers to adjust the sample size in relation to both the cost and the precision. Historical data are often sparse and generally collected using several different sampling schemes, all without inclusion probabilities. This means that an approach is needed to obtain unbiased estimates of the variance of the mean using a model‐based approach. With the use of total phosphorus data from 17 sub‐catchments in southeastern Australia, the ability of a model‐based approach to estimate the effect of sample size on the precision of event and base‐flow mean concentrations. The results showed that for estimating annual base‐flow mean concentration, little gain in precision was achieved above 12 observations per year. Sample sizes greater than 12 samples per event improved event‐based estimates; however, the inclusion of more than 12 samples per event did not greatly reduce the event mean concentration uncertainties. The precision of the base‐flow estimates was most correlated to percentage urban cover, whereas the precision of the event mean estimates was most correlated with catchment size. The method proposed in this work could be readily applied to other water quality variables and other monitoring sites. Copyright © 2014 John Wiley & Sons, Ltd.     

Field study on influence of root characteristics on soil suction distribution in slopes vegetated with Cynodon dactylon and Schefflera heptaphylla

Vegetation evapotranspiration (ET) induced soil water suction reduces hydraulic conductivity and increases shear strength of slopes. Several field studies have been conducted to investigate suction distribution in vegetated slopes. However, these studies were conducted on natural slopes, which are prone to heterogeneity in vegetation and soil conditions. Moreover, studies quantifying the effect of different vegetation species, root characteristics (root depth and root area index) and transpiration reduction function (T_rf) on suction in slopes under natural variation are rare. This study investigated the suction distribution and root characteristics in recompacted slopes vegetated with two different species, i.e. Cynodon dactylon (Bermuda grass) and Schefflera heptaphylla (ivy tree). Bare slope served as a control. Suction distributions during different seasons and rainfall events were monitored. It is found that during the dry season, slope vegetated with young Schefflera heptaphylla seedlings have substantially higher suction within the root zone compared with bare slope and slope vegetated with Cynodon dactylon. This is because Schefflera heptaphylla has a higher root biomass, T_rf and ET than Cynodon dactylon. It was also found that suctions within root zones of vegetated slopes and bare slope were completely destroyed under rainfall events corresponding to 2 years and 20 years return period. Copyright © 2015 John Wiley & Sons, Ltd.     

Vadose CO2 gas drives dissolution at water tables in eogenetic karst aquifers more than mixing dissolution

Most models of cave formation in limestone that remains near its depositional environment and has not been deeply buried (i.e. eogenetic limestone) invoke dissolution from mixing of waters that have different ionic strengths or have equilibrated with calcite at different _pCO₂ values. In eogenetic karst aquifers lacking saline water, mixing of vadose and phreatic waters is thought to form caves. We show here calcite dissolution in a cave in eogenetic limestone occurred due to increases in vadose CO₂ gas concentrations and subsequent dissolution of CO₂ into groundwater, not by mixing dissolution. We collected high‐resolution time series measurements (1 year) of specific conductivity (SpC), temperature, meteorological data, and synoptic water chemical composition from a water table cave in central Florida (Briar Cave). We found SpC, _pCO₂ and calcite undersaturation increased through late summer, when Briar Cave experienced little ventilation by outside air, and decreased through winter, when increased ventilation lowered cave CO_2(g) concentrations. We hypothesize dissolution occurred when water flowed from aquifer regions with low _pCO₂ into the cave, which had elevated _pCO₂. Elevated _pCO₂ would be promoted by fractures connecting the soil to the water table. Simple geochemical models demonstrate that changes in _pCO₂ of less than 1% along flow paths are an order of magnitude more efficient at dissolving limestone than mixing of vadose and phreatic water. We conclude that spatially or temporally variable vadose CO_2(g) concentrations are responsible for cave formation because mixing is too slow to generate observed cave sizes in the time available for formation. While this study emphasized dissolution, gas exchange between the atmosphere and karst aquifer vadose zones that is facilitated by conduits likely exerts important controls on other geochemical processes in limestone critical zones by transporting oxygen deep into vadose zones, creating redox boundaries that would not exist in the absence of caves. Copyright © 2014 John Wiley & Sons, Ltd.     

Fluvial erosion and post-erosional processes on Titan

The surface of Titan has been revealed by Cassini observations in the infrared and radar wavelength ranges as well as locally by the Huygens lander instruments. Sand seas, recently discovered lakes, distinct landscapes and dendritic erosion patterns indicate dynamic surface processes. This study focus on erosional and depositional features that can be used to constrain the amount of liquids involved in the erosional process as well as on the compositional characteristics of depositional areas. Fluvial erosion channels on Titan as identified at the Huygens landing site and in RADAR and Visible and Infrared Mapping Spectrometer (VIMS) observations have been compared to analogous channel widths on Earth yielding average discharges of up to 1600 m³/s for short recurrence intervals that are sufficient to move centimeter-sized sediment and significantly higher discharges for long intervals. With respect to the associated drainage areas, this roughly translates to 1-150 cm/day runoff production rates with 10 years recurrence intervals and by assuming precipitation this implies 0.6-60 mm/h rainfall rates. Thus the observed surface erosion fits with the methane convective storm models as well as with the rates needed to transport sediment. During Cassini's T20 fly-by, the VIMS observed an extremely eroded area at 30° W, 7° S with resolutions of up to 500 m/pixel that extends over thousands of square kilometers. The spectral characteristics of this area change systematically, reflecting continuous compositional and/or particle size variations indicative of transported sediment settling out while flow capacities cease. To account for the estimated runoff production and widespread alluvial deposits of fine-grained material, release of area-dependent large fluid volumes are required. Only frequent storms with heavy rainfall or cryovolcanic induced melting can explain these erosional features.     

Community-level response of coastal microbial biofilms to ocean acidification in a natural carbon dioxide vent ecosystem

The impacts of ocean acidification on coastal biofilms are poorly understood. Carbon dioxide vent areas provide an opportunity to make predictions about the impacts of ocean acidification. We compared biofilms that colonised glass slides in areas exposed to ambient and elevated levels of pCO(2) along a coastal pH gradient, with biofilms grown at ambient and reduced light levels. Biofilm production was highest under ambient light levels, but under both light regimes biofilm production was enhanced in seawater with high pCO(2). Uronic acids are a component of biofilms and increased significantly with high pCO(2). Bacteria and Eukarya denaturing gradient gel electrophoresis profile analysis showed clear differences in the structures of ambient and reduced light biofilm communities, and biofilms grown at high pCO(2) compared with ambient conditions. This study characterises biofilm response to natural seabed CO(2) seeps and provides a baseline understanding of how coastal ecosystems may respond to increased pCO(2) levels.     

Modelling the tidal mixing fronts and seasonal stratification of the Northwest European Continental shelf

We investigate mixing processes under stratified conditions on the Northwest European Continental shelf using a numerical model (POLCOMS). Our results indicate that convection induced by vertical shearing of horizontal density gradients (‘shear-induced convection’) is a regularly occurring feature in the bottom and surface boundary layers in this open shelf-sea situation. Two types of turbulence models are investigated to study their capability for reproducing the observed location of tidal mixing fronts, and the physical processes occurring in seasonally stratified waters. The first model is a one-equation variant of the Mellor–Yamada model, whereas the second model combines a more recent second-momentum closure with a two-equation model. It is found that generally mean frontal positions (as estimated from ICES data) are predicted more accurately by the two-equation model. The one-equation model reproduces the mean frontal locations to 18.1 km (<3 grid spacings) and the two-equation model to 17.1 km; although in the Celtic Sea the accuracy is ∼33 and ∼12 km, respectively. Comparison with historical tide gauges, current metres, CTD stations, and thermistor chain data from the North Sea Project all show an improvement in accuracy when the two-equation model is used. This is particularly apparent in the model's ability to reproduce the spring–neap variability during stratification. We find that in the presence of shear-induced convection the routinely applied clipping of the turbulent length-scale, previously thought to be a minor ingredient in a turbulence model, has a dramatic effect on the results: if the length-scale clipping is not applied, substantial over-mixing is observed to occur. The causes and possible remedies of this effect are investigated. Overall our results demonstrate a sensitivity to the details of the turbulence model that is significantly greater than previously thought.     

The influence of precipitate formation on the chemical oxidation of TCE DNAPL with potassium permanganate

A three-dimensional two-phase flow model is coupled to a non-linear reactive transport model to study the efficacy of potassium permanganate treatment on dense, non-aqueous phase liquid (DNAPL) source removal in porous media. A linear relationship between the soil permeability (k) and concentration of manganese dioxide precipitate ([MnO_2(s)]), k = k_o + S_rind [MnO_2(s)], is utilized to simulate nodal permeability reductions due to precipitate formation. Using published experimental column studies, an S_rind = −5.5 × 10⁻¹⁶ m² L/mg was determined for trichloroethylene (TCE) DNAPL. This S_rind was then applied to treatment simulations on three-dimensional TCE DNAPL source zones comprising either DNAPL at residual saturations, or DNAPL at pooled saturations.