Excel modelling of hydrological systems

This short paper describes the use of Excel spreadsheet Solver facility for deriving estimates in a solute mixing model equation with more than one unknown.     

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.     

Phase relationship between sea level and abrupt climate change

Direct traces of past sea levels are based on the elevation of old coral reefs at times of sea level highstands. However, these measurements are discontinuous and cannot be easily correlated with climate records from ice cores. In this study we show a new approach to recognizing the imprint of sea level changes in continuous sediment records taken from the continental slope at locations that were continuously submerged, even during periods of sea level lowstand. By using a sediment core precisely synchronized with Greenland ice cores, we were able to recognize major floods of the Mediterranean continental shelf over the past 270 kyr. During the last glacial period five flooding events were observed at the onset of the warmest Greenland interstadials. Consistent correspondence between warm climate episodes and eustatic sea level rises shows that these global flooding events were generated by pronounced melting of the Northern Hemisphere ice sheets, due to rapid intensification of Atlantic Meridional Overturning Circulation.The method described in this study opens a new perspective for inter-hemispheric synchronization of marine climate records if applied in other continental margins from the Southern Hemisphere or the equatorial regions.     

Predicting contaminant fate and transport in sediment caps: Mathematical modelling approaches

Sediment capping is a remedial option for managing contaminated sediments that involves the artificial placement of a layer of material over a contaminated area. Sorbent materials such as activated C and coke can be used to amend sand caps to improve cap performance. In this study, analytical and numerical modelling approaches were compared for predicting contaminant fate and transport in sediment caps using several diffusion-controlled and advection-dominated contaminant transport scenarios. An analytical tool was used to predict cap performance at steady-state. These results were compared with the results from the numerical CoReTranS model in which the effective diffusivity and degradation rates were modelled as discontinuous functions at a prescribed bioturbation depth. The numerical approach was also applied to modelling a sorptive cap. It was shown that, while the analytical approach can be used to predict steady-state contaminant transport, the numerical approach is needed to evaluate multiple sediment layers with different transport and sorption characteristics and to examine the transient performance between the time that the single layer transient is applicable (i.e., before penetration of the cap containment layer) and until steady-state in the upper layer. For the 30 cm thick sand cap that was considered in this study, the predicted time to reach steady-state conditions for a diffusion-controlled scenario is 1 ka. For an advection-dominated transport, the time to reach steady-state conditions is reduced to 100 a. The activated C-amended sand cap was more effective in isolating the contaminant within the sorbent layer for a sustained period of time (∼100 a). Results from both modelling approaches showed that capping can effectively reduce contaminant flux to the overlying water with critical variables being cap thickness, groundwater velocity, and sediment sorptivity.     

Climate of the past millennium inferred from varved proglacial lake sediments on northeast Baffin Island,Arctic Canada

This study uses ^239+240Pu-dated varved sediments from Big Round Lake, a proglacial lake on northeast Baffin Island, Arctic Canada to generate a 1000-year-long, annual-resolution record of past climate. Varve thickness is positively correlated with July–August–September temperature measured at Clyde River, 70 km to the north of the lake (r = 0.46, p < 0.001). We therefore interpret the variability and trends in varve thickness to partially represent summer temperature. The coolest Little Ice Age temperatures occurred in this record from 1575 to 1760 AD and were approximately 1.5°C cooler than today (average from 1995 to 2005 AD) and 0.2°C cooler than the last millennium (average from 1000 to 2000 AD). Pre-twentieth-century warmth occurred during two intervals, 970–1150 AD and 1375–1575 AD; temperatures were approximately 1.2°C cooler than today, but 0.1°C warmer than the last millennium. The Big Round Lake varve-thickness record contains features similar to that reconstructed elsewhere in the eastern Canadian Arctic. This high-resolution quantitative record expands our understanding of arctic climate during the past millennium.

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.