Simulations of a coupled hydro-chemo-mechanical system in rocks

A coupled hydro-chemo-mechanical numerical model is developed for these coupled phenomena in many engineering fields. The model has been applied to predicting the response of a stressed rockmass column to an injected reactive fluid (reagent) flow. The response includes evolutions of porosity, permeability, reagent and mineral concentrations during dissolution. In the model, the progress of dissolution is defined by the change in porosity ratio and the porosity increases with dissolution assuming there is no precipitation. The numerical evolutions of porosity, permeability, reagent and mineral concentrations during dissolution are validated against steady state solutions. The model results show that these evolutions are regulated to a certain extent by the applied external loadings: an applied extensional stress enhances the progress of the dissolution process while an applied compression stress slows the progress of the dissolution process. This revised version was published online in July 2006 with corrections to the Cover Date.     

A modelling framework for evaluation of the hydrological impacts of nature‐based approaches to flood risk management,with application to in‐channel interventions across a 29‐km2 scale catchment in the United Kingdom

Nature‐based approaches to flood risk management are increasing in popularity. Evidence for the effectiveness at the catchment scale of such spatially distributed upstream measures is inconclusive. However, it also remains an open question whether, under certain conditions, the individual impacts of a collection of flood mitigation interventions could combine to produce a detrimental effect on runoff response. A modelling framework is presented for evaluation of the impacts of hillslope and in‐channel natural flood management interventions. It couples an existing semidistributed hydrological model with a new, spatially explicit, hydraulic channel network routing model. The model is applied to assess a potential flood mitigation scheme in an agricultural catchment in North Yorkshire, United Kingdom, comprising various configurations of a single variety of in‐channel feature. The hydrological model is used to generate subsurface and surface fluxes for a flood event in 2012. The network routing model is then applied to evaluate the response to the addition of up to 59 features. Additional channel and floodplain storage of approximately 70,000 m³ is seen with a reduction of around 11% in peak discharge. Although this might be sufficient to reduce flooding in moderate events, it is inadequate to prevent flooding in the double‐peaked storm of the magnitude that caused damage within the catchment in 2012. Some strategies using features specific to this catchment are suggested in order to improve the attenuation that could be achieved by applying a nature‐based approach.     

Influence of a coarse interlayer on seawater intrusion and contaminant migration in coastal aquifers

Vertical 2D slice laboratory experiments were carried out in homogenous and layered sand tanks to elucidate the effects of a highly permeable (coarse‐grained sand) interlayer on seawater intrusion and transport of contaminants to a coastal sea. Tidal fluctuations produced oscillations in the seawater–freshwater transition zone, fluctuations of the contaminant infiltration rate and a zigzag contaminant plume outline. The seawater wedge became discontinuous at the (vertical) edges of the interlayer because of increased lateral movement of the seawater–freshwater interface within the interlayer. The contaminant plume formed a tail within the interlayer depending on the tidal stage, and similar to the wedge, its movement was accentuated. A simple analytical model that neglected vertical flow reliably predicted steady‐state seawater intrusion into the coastal aquifer. Numerical modeling was used to gain insight into the groundwater hydrodynamics and contaminant migration. The numerical results confirmed the experimental findings, i.e. that a highly permeable interlayer can provide a rapid transit path for contaminants to reach the seaward boundary and that the interlayer amplifies the effects of tidal fluctuations, resulting in wider transition zones for the seawater wedge and contaminant plume. Numerical simulations further showed that, with increasing interlayer hydraulic conductivity, the maximum seawater intrusion distance inside the interlayer increases approximately linearly. For the fixed‐head contaminant injection condition used, the model showed that contaminant infiltration increases approximately logarithmically with increasing interlayer hydraulic conductivity (other factors held fixed). Copyright © 2013 John Wiley & Sons, Ltd.     

Pattern,style and timing of British–Irish Ice Sheet advance and retreat over the last 45 000 years: evidence from NW Scotland and the adjacent continental shelf

Predicting the future response of ice sheets to climate warming and rising global sea level is important but difficult. This is especially so when fast-flowing glaciers or ice streams, buffered by ice shelves, are grounded on beds below sea level. What happens when these ice shelves are removed? And how do the ice stream and the surrounding ice sheet respond to the abruptly altered boundary conditions? To address these questions and others we present new geological, geomorphological, geophysical and geochronological data from the ice-stream-dominated NW sector of the last British–Irish Ice Sheet (BIIS). The study area covers around 45 000 km² of NW Scotland and the surrounding continental shelf. Alongside seabed geomorphological mapping and Quaternary sediment analysis, we use a suite of over 100 new absolute ages (including cosmogenic-nuclide exposure ages, optically stimulated luminescence ages and radiocarbon dates) collected from onshore and offshore, to build a sector-wide ice-sheet reconstruction combining all available evidence with Bayesian chronosequence modelling. Using this information we present a detailed assessment of ice-sheet advance/retreat history, and the glaciological connections between different areas of the NW BIIS sector, at different times during the last glacial cycle. The results show a highly dynamic, partly marine, partly terrestrial, ice-sheet sector undergoing large size variations in response to sub-millennial-scale climatic (Dansgaard–Oeschger) cycles over the last 45 000 years. Superimposed on these trends we identify internally driven instabilities, operating at higher frequency, conditioned by local topographic factors, tidewater dynamics and glaciological feedbacks during deglaciation. Specifically, our new evidence indicates extensive marine-terminating ice-sheet glaciation of the NW BIIS sector during Greenland Stadials 12 to 9 – prior to the main ‘Late Weichselian’ ice-sheet glaciation. After a period of restricted glaciation, in Greenland Interstadials 8 to 6, we find good evidence for rapid renewed ice-sheet build-up in NW Scotland, with the Minch ice-stream terminus reaching the continental shelf edge in Greenland Stadial 5, perhaps only briefly. Deglaciation of the NW sector took place in numerous stages. Several grounding-zone wedges and moraines on the mid- and inner continental shelf attest to significant stabilizations of the ice-sheet grounding line, or ice margin, during overall retreat in Greenland Stadials 3 and 2, and to the development of ice shelves. NW Lewis was the first substantial present-day land area to deglaciate, in the first half of Greenland Stadial 3 at a time of globally reduced sea-level c. 26 kabp , followed by Cape Wrath at c. 24 kabp. The topographic confinement of the Minch straits probably promoted ice-shelf development in early Greenland Stadial 2, providing the ice stream with additional support and buffering it somewhat from external drivers. However, c. 20–19 kabp , as the grounding-line migrated into shoreward deepening water, coinciding with a marked change in marine geology and bed strength, the ice stream became unstable. We find that, once underway, grounding-line retreat proceeded in an uninterrupted fashion with the rapid loss of fronting ice shelves – first in the west, then the east troughs – before eventual glacier stabilization at fjord mouths in NW Scotland by ~17 kabp. Around the same time, ~19–17 kabp , ice-sheet lobes readvanced into the East Minch – possibly a glaciological response to the marine-instability-triggered loss of adjacent ice stream (and/or ice shelf) support in the Minch trough. An independent ice cap on Lewis also experienced margin oscillations during mid-Greenland Stadial 2, with an ice-accumulation centre in West Lewis existing into the latter part of Heinrich Stadial 1. Final ice-sheet deglaciation of NW mainland Scotland was punctuated by at least one other coherent readvance at c. 15.5 kabp , before significant ice-mass losses thereafter. At the glacial termination, c. 14.5 kabp , glaciers fed outwash sediment to now-abandoned coastal deltas in NW mainland Scotland around the time of global Meltwater Pulse 1A. Overall, this work on the BIIS NW sector reconstructs a highly dynamic ice-sheet oscillating in extent and volume for much of the last 45 000 years. Periods of expansive ice-sheet glaciation dominated by ice-streaming were interspersed with periods of much more restricted ice-cap or tidewater/fjordic glaciation. Finally, this work indicates that the role of ice streams in ice-sheet evolution is complex but mechanistically important throughout the lifetime of an ice sheet – with ice streams contributing to the regulation of ice-sheet health but also to the acceleration of ice-sheet demise via marine ice-sheet instabilities.     

The evolution of the terrestrial-terminating Irish Sea glacier during the last glaciation

Here we reconstruct the last advance to maximum limits and retreat of the Irish Sea Glacier (ISG), the only land-terminating ice lobe of the western British Irish Ice Sheet. A series of reverse bedrock slopes rendered proglacial lakes endemic, forming time-transgressive moraine- and bedrock-dammed basins that evolved with ice marginal retreat. Combining, for the first time on glacial sediments, optically stimulated luminescence (OSL) bleaching profiles for cobbles with single grain and small aliquot OSL measurements on sands, has produced a coherent chronology from these heterogeneously bleached samples. This chronology constrains what is globally an early build-up of ice during late Marine Isotope Stage 3 and Greenland Stadial (GS) 5, with ice margins reaching south Lancashire by 30 ± 1.2 ka, followed by a 120-km advance at 28.3 ± 1.4 ka reaching its 26.5 ± 1.1 ka maximum extent during GS-3. Early retreat during GS-3 reflects piracy of ice sources shared with the Irish-Sea Ice Stream (ISIS), starving the ISG. With ISG retreat, an opportunistic readvance of Welsh ice during GS-2 rode over the ISG moraines occupying the space vacated, with ice margins oscillating within a substantial glacial over-deepening. Our geomorphological chronosequence shows a glacial system forced by climate but mediated by piracy of ice sources shared with the ISIS, changing flow regimes and fronting environments.     

Made in Japan: The Japanese market and herring roe production and management in Canada's southern Gulf of St. Lawrence

The herring roe fishery in Canada's southern Gulf of St. Lawrence is a productive fishery that lies outside the confines of privatized arrangements. Yet, this fishery is compromised by its dependence on the Japanese seafood market. Consideration of the efficiency of this fishery needs to account for the role assumed by this market in setting fishing practices. A more ecologically and economically efficient use of the herring resource requires greater attention to supply chain integration by local harvesters, processors and buyers prior to engaging with distant markets. This is necessary regardless of the property rights regime that is in place.     

In situ U–Pb analysis of shocked zircon from the Charlevoix impact structure,Québec,Canada

Laser ablation inductively coupled plasma‐mass spectrometry (LA‐ICP‐MS) U–Pb geochronology of shocked zircon grains in a vesicular‐fluidal impact melt rock from the ≥54 km Charlevoix impact structure, Québec, Canada, suggests an Ordovician to Silurian age of 450 ± 20 Ma for the impact. This age is anchored by concordant U–Pb results of ~450 Ma for a U‐rich, cryptocrystalline zircon grain in the melt rock, interpreted as a recrystallized metamict zircon crystal; the U–Th–Pb system of the metamict grain was seemingly chronometrically reset by the Charlevoix impact, but withstood later tectonometamorphic events. The new zircon age for Charlevoix is in agreement with a stratigraphically constrained Late Ordovician maximum age of ~453 Ma and corroborates earlier suggestions that the impact occurred most likely in the Ordovician, and not ~100 Myr later, as indicated by previous K/Ar and ⁴⁰Ar/³⁹Ar geochronologic results. The latter may reflect postimpact thermal overprint of impactites during the Salinian (Late Silurian to Early Devonian) and/or Acadian (Late Devonian) orogenies. U–Pb geochronology of zircon crystals in anorthosite exposed in the central uplift of the impact structure yielded a Grenvillian crystallization age of 1062 ± 11 Ma. The preferred Ordovician age for the Charlevoix impact structure, which is partially overthrusted by the Appalachian front, suggests the impact occurred during a phase of Taconian tectonism and an episode of enhanced asteroid bombardment of the Earth. Our results, moreover, demonstrate that (recrystallized) metamict zircon grains may be of particular interest in impact geochronology.     

Differential sensitivity of flounder (Platichthys flesus) in response to oestrogenic chemical exposure: an issue for design and interpretation of monitoring and research programmes

This study was conducted as an initial investigation of 'differential response' in one of the main sentinel organisms used for monitoring programmes in United Kingdom estuaries, the flounder Platichthys flesus. It has been hypothesised that monitoring using species with a wide geographical spread and limited migration, such as flounder, might result in the comparison of different genetic stocks and certainly of populations with differing early life stage contaminant exposure histories. Furthermore, it is probable that these pre-exposure and genetic differences could manifest themselves in an ability to respond differently to contaminant exposure, so-called 'differential response'. It is important that the extent and nature of this response is understood, if we want to be able to fully interpret the monitoring data from such programmes. During this study, flounder were collected from four separate sources; wild caught fish from the estuaries of the Rivers Alde, Mersey and Tyne, and farmed flounder from Port Erin Farm, Isle of Man. Under controlled laboratory conditions, groups of fish from each source were exposed to water-borne concentrations of the synthetic oestrogen ethynylestradiol (EE2) at a nominal concentration of 50 ng/l. Plasma was taken from each male fish after 6 and 10 days exposure and analysed for the presence of vitellogenin (VTG) using an ELISA technique. Significant levels of VTG induction were evident in fish from all sources after both 6 and 10 days exposure. Flounder from the Mersey were the only fish with significantly elevated initial background levels of VTG (day 0) and this appeared to be reflected in that these specimens showed the highest induction response after day 6. However, after day 10, fish from all other sites had a slightly higher mean VTG than those from the Mersey which showed significantly (p < 0.05) lower mean plasma VTG. It is suggested that other differential responses may have been masked by the use of a high dose of EE2 which produced maximum induction in nearly all fish. The findings of the study are discussed in terms of implications for further research into the differential response issue and how the initial plasma VTG figures contribute to a time-series from the Mersey, Tyne and Alde estuaries.