Late Pleistocene stratigraphy and sedimentary environments of the Severnaya Dvina‐Vychegda region in northwestern Russia

The Late Pleistocene stratigraphy from the Severnaya Dvina‐Vychegda region of northwestern Russia is revised based on investigations of new localities, revisiting earlier localities, introduction of about 110 new OSL dates and burial depth corrections of earlier published OSL dates, in addition to six new radiocarbon dates. Most of the OSL samples studied here are from fluvial and subaquaeous sediments, which we found to be well bleached. Six chronostratigraphical units and their sedimentary environment are described, with the oldest unit consisting of pre‐Eemian glacial beds. For the first time, Early Weichselian sediments are documented from the region and a fluvial environment with some vegetation and permafrost conditions is suggested to have persisted from the end of the Eemian until at least about 92 ka ago. The period in which a Middle Weichselian White Sea Lake could have existed is constrained to 67?62 ka, but as the lake level never reached the thresholds of the drainage basin, the lake probably existed only for a short interval within this time‐span. Blocking and reversal of fluvial drainage started again around 21?20 ka ago when the Fennoscandian Ice Sheet advanced into the area, reaching its maximum 17?15 ka ago. At that time, an ice‐dammed lake reached its maximum water level, which was around 135 m above present sea level. Drainage of the lake started shortly after 15 ka ago, and the lake was emptied within 700 years. Severe periglacial conditions, with permafrost and aeolian activity, prevailed in the area until about 10.7 ka.     

Evolution of the Taupo Volcanic Center, New Zealand: petrological and thermal constraints from the Omega dacite

The 20 ka ~0.1 km³ Omega dacite, which erupted shortly after the 26.5 ka Oruanui super-eruption, compositionally stands out among Taupo Volcanic Zone (TVZ) magmas, which are overwhelmingly characterized by rhyolites (>90 % by volume). The previously reported presence of inherited zircons in this zircon-undersaturated magma has provided unequivocal evidence for the involvement of upper-crustal material in a 1–10 year timescale prior to the Omega eruption. However, whether this crustal involvement is characterized by wholesale, melting of preexisting crust or subordinate bulk assimilation into an already differentiated magma body remains unclear. To disentangle these processes, we describe the mineral chemistry of the major phases present in the Omega dacite and determine intensive parameters describing magma chamber conditions. Dominantly unimodal populations of plagioclase (An_50–60), orthopyroxene (Mg# from 58 to 68), and clinopyroxene (Mg# from 65 to 73), along with coexisting equilibrium pairs of Fe–Ti oxides, constrain pre-eruptive temperatures to 850–950 °C, a pressure between ~3 and 7 kbars, and an oxygen fugacity of ~NNO. MELTS thermodynamic modeling suggests that this phase assemblage is in equilibrium with the bulk rock and glass compositions of the Omega dacite at these estimated P–T–fO₂ pre-eruptive conditions. Combining these petrological observations with insights into conductive thermal models of magma–crust interactions, we argue that the Omega dacite more likely formed in the mid-to-lower crust via protracted processing through fractional crystallization coupled with some assimilation (AFC). Incorporation of crustal material is likely to have occurred at various stages, with the inherited zircons (and potentially parts of glomerocrysts) representing late and subordinate upper-crustal assimilants. This petrogenetic model is consistent with the presence of a differentiating crustal column, consisting of a polybaric fractional crystallization and assimilation history. On the basis of petrological, thermal, and geophysical considerations, upper-crustal reservoirs, which feed large-scale rhyolitic volcanism in the TVZ, most likely take the form of large, long-lived crystal mush zones. Following large eruptions, such as the Oruanui event, this mush is expected to crystallize significantly (up to 70–80 vol% crystals) due to syn-eruptive decompression. Hence, the Omega dacite, immediately post-dating the Oruanui event, potentially represents incoming deeper recharge of less-evolved magma that was able to penetrate the nearly solidified upper-crustal mush. Over the past 20,000 years, similar intermediate recharge magmas have incrementally reheated, reconstructed, and reactivated the upper-crustal mush zone, allowing a gradual return to rhyolitic volcanism at the Taupo Volcanic Center.     

Pembrokeshire,South Wales

Unravelling the geology of Pembrokeshire (southwest Wales) has been a rite of passage for British undergraduate students for generations. A favourite destination for an Easter fieldtrip, this complex region contains diverse geology of late Precambrian to Carboniferous age. In the course of a five‐day excursion, a student can uncover evidence for two plate tectonic cycles: the birth of the early Palaeozoic Iapetus Ocean and its death in the Caledonian Orogen, and the reworking of the late Palaeozoic Rheic back‐arc basins on the northern margins of the Variscides. In doing so, there is an opportunity to examine highly variable palaeoenvironments including the deposits of alluvial fans, flood plains, deltas, estuaries, carbonate shelves and deeper marine settings, as well as developing fundamental field skills in geological mapping, structural analysis and sedimentary logging. In this article, I introduce the fascinating story of the geological evolution of Pembrokeshire and describe some of the classic localities useful for basic geological training.     

Anthropogenic subsidence in the Mexicali Valley,Baja California,Mexico, and slip on the Saltillo fault

Deep fluid extraction in the Cerro Prieto geothermal field (CPGF) has caused subsidence and induced slip on tectonic faults in the Mexicali Valley (Baja California, Mexico). The Mexicali Valley is located in the southern part of the Salton Trough, at the boundary between the Pacific and North American plates. The Valley is characterized by being a zone of continuous tectonic deformation, geothermal activity, and seismicity. Within the Cerro Prieto pull-apart basin, seismicity is concentrated mainly in swarms, while strong earthquakes have occurred in the Imperial and Cerro Prieto transform faults, that are the eastern and western bound of the basin. Since 1973, fluid extraction at the CPGF has influenced deformation in the area, accelerating the subsidence and causing rupture (frequently as vertical slip or creep) on the surface traces of tectonic faults. Both subsidence and fault slip are causing damage to infrastructure like roads, railroad tracks, irrigation channels, and agricultural fields. Currently, accelerated extraction in the eastern part of CPGF has shifted eastwards the area of most pronounced subsidence rate; this accelerated subsidence can be observed at the Saltillo fault, a southern branch of the Imperial fault in the Mexicali Valley. Published leveling data, together with field data from geological surveys, geotechnical instruments, and new InSAR images were used to model the observed deformation in the area in terms of fluid extraction. Since the electricity production in the CPGF is an indispensable part of Baja California economy, extraction is sure to continue and may probably increase, so that the problem of damages caused by subsidence will likely increase in the future.     

Estimation of surface roughness (z0) over a stabilizing coastal dune field based on vegetation and topography

Accurate knowledge of the surface roughness and the resultant wind speed are important for many applications, such as climatic models, wind power meteorology, agriculture and erosion hazards, especially on sand dunes in arid and semi‐arid environments, where vegetation cover is scarce. In this study we aimed at quantifying the effects of vegetation cover and topography on surface roughness over a stabilizing dune field on the southern coast of Israel. Forty‐six wind measurements were made at various distances from the coastline, ranging from 10 to 2800 m, and z₀ values were calculated from the wind measurements based on the ratio between the wind gust and the average wind speed. We estimated vegetation cover using the soil adjusted vegetation index (SAVI) from Landsat satellite images for the upwind sector at various lengths, ranging from 15 to 400 m, and based on digital elevation models and differential GPS field measurements we calculated the topographic variable of the relative heights of the stations. z₀ values were positively correlated with the winter SAVI values (r = 0·87 at an upwind length of 200 m) and negatively correlated with the relative height (r = ?0·68 at an upwind length of 200–400 m for the inland dune stations). Using these variables we were able to create a map of estimated z₀ values having an accuracy of over 64%. Such maps provide a better understanding of the spatial variability in both wind speed and sand movement over coastal dune areas. Copyright © 2007 John Wiley & Sons, Ltd.     

Numerical modelling of the potential effects of a dam on a coastal aquifer in S. Spain

This study presents the results of a three‐dimensional variable‐density numerical modelling of the Motril‐Salobreña coastal aquifer and the possible effects of the entry into service in May 2005 of the Rules Dam, located just 17 km from the coast. Present parameters of the Motril‐Salobreña aquifer show that the system's conditions are very similar to a natural regime. The dam will substantially alter aquifer recharge, as the entry flow through the alluvial sediments of the Guadalfeo River will be entirely cut off or drastically reduced. Different scenarios reproducing the possible evolution of the aquifer under operation of the Rules Dam have been modelled. In most cases, results indicate that the conditions of the aquifer would worsen, with a general advance of the freshwater–saltwater interface. The area with most risk of saltwater intrusion is the old mouth of the Guadalfeo River, where the mixing zone could advance 1200 m inland. It is proposed that maintaining a 5–6 Mm³ year^?1 ‘ecological flow’ in the Guadalfeo River could prevent this saline advance. This application demonstrates that variable‐density models are potentially useful tools for estimating the effects of dams on the hydrodynamic and hydrochemical conditions of a coastal aquifer. Copyright © 2009 John Wiley & Sons, Ltd.     

Experimental study of extreme shear stress for shallow flow under simulated rainfall

The raindrop impact and overland flow are two major factors causing soil detachment and particle transportation. In this study, the turbulent characteristics of the shallow rain‐impacted water flow were investigated using a 2‐D fibre‐optic laser Doppler velocimetry (FLDV) and an artificial rainfall simulator. The fluctuating turbulent shear stress was computed using digital data processing techniques. The experimental data showed that the Reynolds shear stress follows a probability distribution with heavy tails. The tail probability increases with an increase of rainfall intensity or raindrop diameter, and it decreases with an increase of Reynolds number. A modified empirical equation was derived using both the raindrop diameter and rainfall intensity as independent variables to provide a better prediction of the Darcy‐Weisbach friction coefficient f under rainfall conditions. Copyright © 2009 John Wiley & Sons, Ltd.     

Simultaneous identification of parameter,initial condition,and boundary condition in groundwater modelling

In this paper, we perform an inverse method to simultaneously estimate aquifer parameters, initial condition, and boundary conditions in groundwater modelling. The parameter estimation is extended to a complete inverse problem that makes the calibrated groundwater flow model more realistic. The adjoint state method, the gradient search method, and the least square error algorithm are combined to build the optimization procedure. Horizontal two‐dimensional groundwater flow in a confined aquifer is exemplified to demonstrate the correlation between unknowns, the contribution of observation, as well as the suitability of applying the inverse method. The correlation analysis shows the connection between storage coefficient and initial condition. Besides, transmissivity and boundary conditions are also highly correlated. More observations at different location and time are necessary to provide sufficient information. A time series of unsteady head is requested for estimation of storage coefficient and initial condition. Observation near boundary is very effective for boundary condition estimation. The observation at pumping well mostly contributes to the estimation of transmissivity. According to all observations, it is possible to identify parameters, initial condition, and boundary condition simultaneously. Furthermore, the results not only illustrate the traditional assumption of known boundary condition but also initial condition, which may cause an incorrect estimation. Copyright © 2009 John Wiley & Sons, Ltd.     

Bayesian analysis of stage–fall–discharge models for gauging stations affected by variable backwater

Parsimonious stage–fall–discharge rating curve models for gauging stations subject to backwater complications are developed from simple hydraulic theory. The rating curve models are compounded in order to allow for possible shifts in the hydraulics when variable backwater becomes effective. The models provide a prior scientific understanding through the relationship between the rating curve parameters and the hydraulic properties of the channel section under study. This characteristic enables prior distributions for the rating curve parameters to be easily elicited according to site‐specific information and the magnitude of well‐known hydraulic quantities. Posterior results from three Norwegian and one American twin‐gauge stations affected by variable backwater are obtained using Markov chain Monte Carlo simulation techniques. The case studies demonstrate that the proposed Bayesian rating curve assessment is appropriate for developing rating procedures for gauging stations that are subject to variable backwater. Copyright © 2009 John Wiley & Sons, Ltd.