222Rn as Natural Tracer for LNAPL Recovery in a Crude Oil‐Contaminated Aquifer

The objective of this study was to investigate whether ²²²Rn in groundwater can be used as a tracer for light non‐aqueous phase liquid (LNAPL) quantification at a field site treated by dual‐phase LNAPL removal. After the break of a pipeline, 5 ha of soil in the nature reserve Coussouls de Crau in southern France was contaminated by 5100 m³ of crude oil. Part of this oil seeped into the underlying gravel aquifer and formed a floating oil body of about 3.9 ha. The remediation consists of plume management by hydraulic groundwater barriers and LNAPL extraction in the source zone. ²²²Rn measurements were performed in 21 wells in and outside the source zone during 15 months. In uncontaminated groundwater, the radon activity was relatively constant and remained always >11 Bq/L. The variability of radon activity measurements in wells affected by the pump‐and‐skim system was consistent with the measurements in wells that were not impacted by the system. The mean activities in wells in the source zone were, in general, significantly lower than in wells upgradient of the source zone, owing to partitioning of ²²²Rn into the oil phase. The lowest activities were found in zones with high non‐aqueous phase liquid (NAPL) recovery. LNAPL saturations around each recovery well were furthermore calculated during a period of high groundwater level, using a laboratory‐determined crude oil–water partitioning coefficient of 38.5 ± 2.9. This yielded an estimated volume of residual crude oil of 309 ± 93 m³ below the capillary fringe. We find that ²²²Rn is a useful and cheap groundwater tracer for finding zones of good LNAPL recovery in an aquifer treated by dual‐phase LNAPL removal, but that quantification of NAPL saturation using Rn is highly uncertain.     

Simulating Stable Isotope Ratios in Plumes of Groundwater Pollutants with BIOSCREEN‐AT‐ISO

BIOSCREEN is a well‐known simple tool for evaluating the transport of dissolved contaminants in groundwater, ideal for rapid screening and teaching. This work extends the BIOSCREEN model for the calculation of stable isotope ratios in contaminants. A three‐dimensional exact solution of the reactive transport from a patch source, accounting for fractionation by first‐order decay and/or sorption, is used. The results match those from a previously published isotope model but are much simpler to obtain. Two different isotopes may be computed, and dual isotope plots can be viewed. The dual isotope assessment is a rapidly emerging new approach for identifying process mechanisms in aquifers. Furthermore, deviations of isotope ratios at specific reactive positions with respect to “bulk” ratios in the whole compound can be simulated. This model is named BIOSCREEN‐AT‐ISO and will be downloadable from the journal homepage.     

Ground Surface Temperature History in Poland in the 16th–20th Centuries Derived from the Inversion of Geothermal Profiles

Ground Surface Temperature (GST) history in Poland was derived from the inversion of temperature-depth profiles in over 20 wells. Temperature histories for the period 1500 A.D. through 1977 A.D. agree well with the instrumental record of the surface-air temperature available for the last two centuries. A statistical correlation of the reconstructed histories (from the well temperature data) with the instrumental record (air temperature) from the homogeneous Warsaw series is high (>0.8). Functional space inversion (FSI) of the temperature data with depth shows that beginning in the early 19th century, temperatures warmed by 0.9 ± 0.1°C following a long period of colder climate before. The last number could be a minimal as higher warming was calculated using a simple model based on surface temperature for the observational period (homogenized Warsaw surface temperature series, Lorenc, 2000) and POM (pre-observational mean; Harris and Chapman, 1998) of –1.53^oC below the 1951–1980 mean temperature level.     

Simulating rainfall time-series: how to account for statistical variability at multiple scales?

Daily rainfall is a complex signal exhibiting alternation of dry and wet states, seasonal fluctuations and an irregular behavior at multiple scales that cannot be preserved by stationary stochastic simulation models. In this paper, we try to investigate some of the strategies devoted to preserve these features by comparing two recent algorithms for stochastic rainfall simulation: the first one is the modified Markov model, belonging to the family of Markov-chain based techniques, which introduces non-stationarity in the chain parameters to preserve the long-term behavior of rainfall. The second technique is direct sampling, based on multiple-point statistics, which aims at simulating a complex statistical structure by reproducing the same data patterns found in a training data set. The two techniques are compared by first simulating a synthetic daily rainfall time-series showing a highly irregular alternation of two regimes and then a real rainfall data set. This comparison allows analyzing the efficiency of different elements characterizing the two techniques, such as the application of a variable time dependence, the adaptive kernel smoothing or the use of low-frequency rainfall covariates. The results suggest, under different data availability scenarios, which of these elements are more appropriate to represent the rainfall amount probability distribution at different scales, the annual seasonality, the dry-wet temporal pattern, and the persistence of the rainfall events.     

Joint inversion of teleseismic P waveforms and surface-wave group velocities from ambient seismic noise in the Bohemian Massif

Joint inversion of teleseismic P-waveforms and local group velocities of surface waves retrieved from ambient seismic noise has been performed to model velocity structure of the crust and uppermost mantle of the Bohemian Massif. We analysed P-waveforms of 381 teleseismic earthquakes recorded at 54 broadband seismic stations located on the territory of the Czech Republic and in its close surroundings. Group velocities of Rayleigh and Love surface waves were obtained by cross-correlating long-term recordings of seismic noise. The basis for waveform inversion is the well-known methodology of P-to-S receiver functions constructed from converted phases. Due to instabilities in direct inversion of receiver functions caused by the necessity of applying deconvolution, we propose an alternative formulation to fit observed and calculated radial components of P waveforms. The joint inversion is transformed into a search for the minimum of the cost function defined as a weighted sum of waveform and group velocity misfits. With the use of the robust stochastic optimizer (Differential Evolution Algorithm), neither derivatives nor a starting model are needed. The task was solved for 1D layered isotropic models of the crust and the uppermost mantle. We have performed a sequence of inversions with models containing one, two, three and four layers above a half-space. By using statistical criteria (F-test) we were able to select the simplest velocity models satisfying data and representing local geological structures. Complex crustal models are typical for stations located close to boundaries of major tectonic units. The relatively low average P to S wave-velocity ratio is in agreement with the generally accepted view that the BM crust is predominantly felsic.     

Why does the Antarctic Peninsula Warm in climate simulations?

The Antarctic Peninsula has warmed significantly since the 1950s. This pronounced and isolated warming trend is collectively captured by 29 twentieth-century climate hindcasts participating in the version 3 Coupled Model Intercomparison Project. To understand the factors driving warming trends in the hindcasts, we examine trends in Peninsula region’s atmospheric heat budget in every simulation. We find that atmospheric latent heat release increases in nearly all hindcasts. These increases are generally anthropogenic in origin, and account for about 60% of the ensemble-mean warming trend in the Peninsula. They are driven primarily by well-understood features of the anthropogenic intensification of global hydrological cycle. As sea surface temperature increases, moisture contained in atmospheric flows increases. When such flows are forced to ascend the Peninsula’s topography, enhanced local latent heat release results. The mechanism driving the warming of the Antarctic Peninsula is therefore clear in the models. Evidence for a similar mechanism operating in the real world is seen in the increasing snow accumulation rates inferred from ice cores drilled in the Peninsula. However, the relative importance of this mechanism and other processes previously identified as potentially causing the observed warming, such as the recent sea ice retreat in the Bellingshausen Sea, is difficult to assess. Thus the relevance of the simulated warming mechanism to the observed warming is unclear, in spite of its robustness in the models.     

A new airborne polar Nephelometer for the measurements of optical and microphysical cloud properties. Part I: Theoretical design

A new optical sensor, the airborne Polar Nephelometer, is described. The sensor is designed to measure the optical and microphysical parameters of clouds containing either water droplets or ice crystals or a mixture of these particles ranging in size from a few micrometers to about 500 m diameter. The probe measures the scattering phase function of an ensemble of cloud particles intersecting a collimated laser beam near the focal point of a paraboloïdal mirror. The light scattered from polar angles from 3.49° to 169° is reflected onto a circular array of 33 photodiodes. The signal processing electronics and computer storage can provide one measurement of the scattering phase function every 100 ms or every 0.2 ms. The first part of the paper describes the theoretical design of a prototype version of the probe.