Turbulent velocity distributions and implied trajectory models

Well-mixed, first-order Lagrangian stochastic (LS) particle trajectory models are derived from several idealized (“toy”) turbulent velocity distributions, and their performance is compared against the observations of Project Prairie Grass, i.e., the case of a continuous point source of tracer near the ground, in the horizontally homogeneous and neutrally stratified surface layer. Although in a context of limited information a Gaussian distribution is the preferred choice, and although the Gaussian corresponds to the simplest of this set of LS models (namely, the Langevin equation), models stemming from other velocity distributions give similar, albeit distinguishable, predictions.     

Predicting plume meandering and averaging time effects on mean and fluctuating concentrations in atmospheric dispersion simulated in a water channel

Plume meandering and averaging time effects were measured directly using a high spatial resolution, high frequency, linescan laser-induced fluorescence (LIF) technique for measuring scalar concentrations in a plume dispersing in a water channel. Post-processing of the collected data removed time dependent background dye levels and corrected for attenuation across the laser beam to produce accurate measurements over long sample times in both a rough surface boundary-layer shear flow and shear free grid-generated turbulent flow. The data were used to verify the applicability of a meandering plume model for predicting the properties of mean and fluctuating concentrations. The centroid position of the crosswind concentration profile was found to have a Gaussian probability density function and the instantaneous plume spread about the centroid fluctuated log-normally. A modified travel-time power law model for averaging time adjustment was developed and compared to the widely used, but much less accurate, 0.2 power-law model.     

The settling of consecutive spheres in viscoplastic fluids

One of the factors contributing to the uncertainties involved in the estimation of particle settling velocity in viscoplastic fluids is the time-dependent effect where the viscous parameters of the fluid change as a particle flows through and shears the medium. These changes, particularly at low shear Reynolds numbers, are reflected in the settling velocity of a following sphere that is released some time after an initial one, with the following sphere having a significantly greater velocity. This study found that changes in both fall velocity and equivalent viscosity can be correlated satisfactorily by a power law equation to the dimensionless form of the time interval between releases, and the rheogram shape factor for the fluid. A collision of particles occurs in cases where the time interval between releases is small, after which the particles combine and travel at a terminal velocity. A new variable, β, which takes into account the different surficial stress of the combined spheres, was introduced to the correlation of Wilson et al. [Wilson, K.C., Horsley, R.R., Kealy, T., Reizes, J.A., Horsley, M.R., 2003. Direct prediction of fall velocities in non-Newtonian materials. Int. J. Miner. Process. 71, 17–30] β was found to depend on the rheogram shape factor for the fluid and the shear Reynolds number for the particle. The validity of this approach was supported by experimental data.     

The use of tracers to assess leakage from the sequestration of CO2 in a depleted oil reservoir,New Mexico,USA

Geological sequestration of CO₂ in depleted oil reservoirs is a potentially useful strategy for greenhouse gas management and can be combined with enhanced oil recovery. Development of methods to estimate CO₂ leakage rates is essential to assure that storage objectives are being met at sequestration facilities. Perfluorocarbon tracers (PFTs) were added as three 12 h slugs at about one week intervals during the injection of 2090 tons of CO₂ into the West Pearl Queen (WPQ) depleted oil formation, sequestration pilot study site located in SE New Mexico. The CO₂ was injected into the Permian Queen Formation. Leakage was monitored in soil–gas using a matrix of 40 capillary adsorbent tubes (CATs) left in the soil for periods ranging from days to months. The tracers, perfluoro-1,2-dimethylcyclohexane (PDCH), perfluorotrimethylcyclohexane (PTCH) and perfluorodimethylcyclobutane (PDCB), were analyzed using thermal desorption, and gas chromatography with electron capture detection. Monitoring was designed to look for immediate leakage, such as at the injection well bore and at nearby wells, and to develop the technology to estimate overall CO₂ leak rates based on the use of PFTs. Tracers were detected in soil–gas at the monitoring sites 50 m from the injection well within days of injection. Tracers continued to escape over the following years. Leakage appears to have emanated from the vicinity of the injection well in a radial pattern to about 100 m and in directional patterns to 300 m. Leakage rates were estimated for the 3 tracers from each of the 4 sets of CATs in place following the start of CO₂ injection. Leakage was fairly uniform during this period. As a first approximation, the CO₂ leak rate was estimated at about 0.0085% of the total CO₂ sequestered per annum.     

Seasonal global mean sea level change from satellite altimeter, GRACE, and geophysical models

We estimate seasonal global mean sea level changes using different data resources, including sea level anomalies from satellite radar altimetry, ocean temperature and salinity from the World Ocean Atlas 2001, time-variable gravity observations from the Gravity Recovery and Climate Experiment (GRACE) mission, and terrestrial water storage and atmospheric water vapor changes from the NASA global land data assimilation system and National Centers for Environmental Prediction reanalysis atmospheric model. The results from all estimates are consistent in amplitude and phase at the annual period, in some cases with remarkably good agreement. The results provide a good measure of average annual variation of water stored within atmospheric, land, and ocean reservoirs. We examine how varied treatments of degree-2 and degree-1 spherical harmonics from GRACE, laser ranging, and Earth rotation variations affect GRACE mean sea level change estimates. We also show that correcting the standard equilibrium ocean pole tide correction for mass conservation is needed when using satellite altimeter data in global mean sea level studies. These encouraging results indicate that is reasonable to consider estimating longer-term time series of water storage in these reservoirs, as a way of tracking climate change.     

Mapping and analysing virtual outcrops

Laser scanning is a very efficient way to generate realistic, high-resolution digital models of 3-D geological outcrops. This paper discusses the methodologies involved in the creation and analysis of virtual outcrops, based on laser scanner data. The visualisation of the laser scanner data as a photorealistic 3-D object is described. Geological features picked out on the virtual outcrop (e.g. fractures, faults or bedding planes) can be extrapolated outward, into space, and inward, into the subsurface, using tension surfaces.Electronic Supplementary Material Supplementary material is available in the online version of this article at Reviewed by: J.D. Clemens, D. Yuen     

Crustal assimilation in basalt and jotunite: Constraints from layered intrusions

To constrain the amount and rate of crustal contamination that is possible in basaltic and jotunitic magma, and to gain an insight into the physical and thermal processes of assimilation in crustal magma chambers, we have modelled published Sr and Nd isotopic data from three layered intrusions. Well-exposed sequences of cumulates with no evidence of magma recharge provide direct records of concurrent assimilation and fractional crystallization (AFC). The key to the modelling is that F, the mass fraction of magma remaining in the chamber, can be estimated from the thicknesses of the studied cumulate sequences. This allows AFC model curves to be fitted to the isotopic data by varying r, the ratio of the rate of mass assimilated to the rate of mass crystallized. The results of modelling show that r is nearly constant in 800 to 2000 m thick sequences of cumulates displaying up-section decreases in anorthite content of plagioclase, increases in whole-rock Sr₀ (initial ⁸⁷Sr/⁸⁶Sr) and decreases in whole-rock εNd₀ (initial εNd). The r-values of the layered sequences range from 0.12 in the Fongen–Hyllingen Intrusion, over 0.20 in the Bjerkreim–Sokndal Intrusion, to 0.27 in the Hasvik Intrusion. The total amount of assimilation, the bulk crust/magma ratio, reaches values of 0.08, 0.19 and 0.28 at the level of the most contaminated samples after 60% to 80% crystallisation, whereas the instantaneous crust/magma ratio of the most contaminated magmas were respectively 0.14, 0.46, and 0.70, for the three intrusions.Innumerable country rock xenoliths occur in the three layered intrusions and played a crucial role in the assimilation process. The xenoliths spalled off the roofs of the magma chambers during magma emplacement and their initial temperature and composition relate to r in the intrusions. In the Hasvik Intrusion (r = 0.27), the initial temperature of the country rocks was 450 °C and the xenoliths were fusible metasediments and therefore produced a high fraction of partial melt that could be assimilated. In the Bjerkreim–Sokndal Intrusion (r = 0.20), the country rocks were initially at temperatures of 640–880 °C but included both refractory massif-type anorthosite and fusible gneisses. In the Fongen–Hyllingen Intrusion (r = 0.12), the country rocks were cooler (300 °C) and the xenoliths include refractory metabasalt (dominant) and fusible metapelite. We argue that the refractory metabasalt and anorthosite xenoliths acted mainly as heat sinks, resulting in reduced r-values in Fongen–Hyllingen and Bjerkreim–Sokndal Intrusions.Heating of refractory and fusible xenoliths, and melting of fusible xenoliths absorbed sensible and latent heat of the magma. Energy-balanced modelling shows that up to 75% of the heat available was absorbed by xenoliths within the magma chambers, promoting higher rates of cooling and crystallisation than would have resulted from loss of heat to the envelope of country rocks alone. The high r-values reflect the amount of heat absorbed by heating and melting country rock within the magma chambers themselves, and their constancy reflects the ready availability of fusible xenoliths.