Simplification and Validation of a Spectral-Tensor Model for Turbulence Including Atmospheric Stability

A spectral-tensor model of non-neutral, atmospheric-boundary-layer turbulence is evaluated using Eulerian statistics from single-point measurements of the wind speed and temperature at heights up to 100 m, assuming constant vertical gradients of mean wind speed and temperature. The model has been previously described in terms of the dissipation rate \(\epsilon \), the length scale of energy-containing eddies \(\mathcal {L}\), a turbulence anisotropy parameter \(\varGamma \), the Richardson number Ri, and the normalized rate of destruction of temperature variance \(\eta _\theta \equiv \epsilon _\theta /\epsilon \). Here, the latter two parameters are collapsed into a single atmospheric stability parameter z / L using Monin–Obukhov similarity theory, where z is the height above the Earth’s surface, and L is the Obukhov length corresponding to \(\{Ri,\eta _\theta \}\). Model outputs of the one-dimensional velocity spectra, as well as cospectra of the streamwise and/or vertical velocity components, and/or temperature, and cross-spectra for the spatial separation of all three velocity components and temperature, are compared with measurements. As a function of the four model parameters, spectra and cospectra are reproduced quite well, but horizontal temperature fluxes are slightly underestimated in stable conditions. In moderately unstable stratification, our model reproduces spectra only up to a scale \(\sim \) 1 km. The model also overestimates coherences for vertical separations, but is less severe in unstable than in stable cases.     

Use of Rhodamine water tracer in the marshland upwelling system

Rhodamine water tracer (RWT) was used to characterize the migration of waste water within the saline subsurface of a marshland upwelling system (MUS), which is an alternative on-site waste water treatment system for coastal areas. Field tracer studies were performed to investigate the fresh waste water plume movement within the saline ground water. Pore velocities were calculated using first detection times and ranged from 0.68 to 10.7 x 10(-4) cm/sec for the loamy sandy soil matrix present at the site. Use of RWT in the field also allowed determination of main and preferential flowpaths. One- and two-dimensional laboratory experiments were performed using silica sand to investigate the interactions of the organically rich waste water with RWT within the zone surrounding the point of injection (one-dimensional) and the impact of background salinity on plume movement (two-dimensional). The results from these studies were used to help explain the field data. One-dimensional breakthrough curves revealed retardation factors for the RWT in the waste water mixture of 1.73 to 1.90. These results were similar to other researchers, indicating little interaction between the waste water and RWT. Variations in pore water salinity (5, 15, 25, and 35 ppt) were found to have a significant effect on pore water velocities of the fresh water plume (two-dimensional), indicating the need to incorporate background salinities into the design process for MUS.     

Climate change and coastal flooding in Metro Boston: impacts and adaptation strategies

Sea level rise (SLR) due to climate change will increase storm surge height along the 825 km long coastline of Metro Boston, USA. Land at risk consists of urban waterfront with piers and armoring, residential areas with and without seawalls and revetments, and undeveloped land with either rock coasts or gently sloping beachfront and low-lying coastal marshes. Risk-based analysis shows that the cumulative 100 year economic impacts on developed areas from increased storm surge flooding depend heavily upon the adaptation response, location, and estimated sea level rise. Generally it is found that it is advantageous to use expensive structural protection in areas that are highly developed and less structural approaches such as floodproofing and limiting or removing development in less developed or environmentally sensitive areas.     

A population cartogram of New Zealand

The authors use 1981 census data to create a cartogram showing the counties of New Zealand in proportion to their population.     

Characterization of uraninite nanoparticles produced by Shewanella oneidensis MR-1

The reduction of uranium(VI) by Shewanella oneidensis MR-1 was studied to examine the effects of bioreduction kinetics and background electrolyte on the physical properties and reactivity to re-oxidation of the biogenic uraninite, UO_2(s). Bioreduction experiments were conducted with uranyl acetate as the electron acceptor and sodium lactate as the electron donor under resting cell conditions in a 30 mM NaHCO₃ buffer, and in a PIPES-buffered artificial groundwater (PBAGW). MR-1 was cultured in batch mode in a defined minimal medium with a specified air-to-medium volume ratio such that electron acceptor (O₂) limiting conditions were reached just when cells were harvested for subsequent experiments. The rate of U(VI) bioreduction was manipulated by varying the cell density and the incubation temperature (1.0 × 10⁸ cell ml⁻¹ at 20 °C or 2.0 × 10⁸ cell ml⁻¹ at 37 °C) to generate U(IV) solids at “fast” and “slow” rates in the two different buffers. The presence of Ca in PBAGW buffer altered U(VI) speciation and solubility, and significantly decreased U(VI) bioreduction kinetics. High resolution transmission electron microscopy was used to measure uraninite particle size distributions produced under the four different conditions. The most common primary particle size was 2.9-3.0 nm regardless of U(VI) bioreduction rate or background electrolyte. Extended X-ray absorption fine-structure spectroscopy was also used to estimate uraninite particle size and was consistent with TEM results. The reactivity of the biogenic uraninite products with dissolved oxygen was tested, and neither U(VI) bioreduction rate nor background electrolyte had any statistical effect on oxidation rates. With MR-1, uraninite particle size was not controlled by the bioreduction rate of U(VI) or the background electrolyte. These results for MR-1, where U(VI) bioreduction rate had no discernible effect on uraninite particle size or oxidation rate, contrast with our recent research with Shewanella putrefaciens CN32, where U(VI) bioreduction rate strongly influenced both uraninite particle size and oxidation rate. These two studies with Shewanella species can be viewed as consistent if one assumes that particle size controls oxidation rates, so the similar uraninite particle sizes produced by MR-1 regardless of U(VI) bioreduction rate would result in similar oxidation rates. Factors that might explain why U(VI) bioreduction rate was an important control on uraninite particle size for CN32 but not for MR-1 are discussed.     

The provenance of Australian uranium ore concentrates by elemental and isotopic analysis

Elemental and isotopic ratio analyses of U ore concentrate samples, from the 3 operating U mining facilities in Australia, were carried out to determine if significant variations exist between their products, thereby allowing the U ore concentrate’s origin to be identified. Elemental analyses were conducted using inductively coupled plasma mass spectrometry (ICP-MS) and X-ray fluorescence spectrometry (XRF). Lead isotope ratios were measured using ICP-MS and U isotope analyses were conducted using thermal ionisation mass spectrometry (TIMS). Minute quantities of sample, such as that obtained from a swipe, were also examined for elemental concentrations using secondary ion mass spectrometry (SIMS). The results of multivariate statistical analysis show clear patterns in the trace elemental composition of the processed U ores, indicating that it is possible to use this feature as a unique identifier of an Australian U ore concentrate’s source. Secondary ion mass spectrometry analyses also allow individual particles to be differentiated using this ‘fingerprinting’ technique. Isotope ratios determined using TIMS reveal that there is a significant difference in the n(²³⁴U)/n(²³⁸U) isotope ratio between the U ore concentrate from each mine.     

Holocene climate and environmental history of East Greenland inferred from lake sediments

Prediction of future Arctic climate and environmental changes, as well as associated ice-sheet behavior, requires placing present-day warming and reduced ice extent into a long-term context. Here we present a record of Holocene climate and glacier fluctuations inferred from the paleolimnology of small lakes near Istorvet ice cap in East Greenland. Calibrated radiocarbon dates of organic remains indicate deglaciation of the region before ~10,500 years BP, after which time the ice cap receded rapidly to a position similar to or less extensive than present, and lake sediments shifted from glacio-lacustrine clay to relatively organic-rich gyttja. The lack of glacio-lacustrine sediments throughout most of the record suggests that the ice cap was similar to or smaller than present throughout most of the Holocene. This restricted ice extent suggests that climate was similar to or warmer than present, in keeping with other records from Greenland that indicate a warm early and middle Holocene. Middle Holocene magnetic susceptibility oscillations, with a ~200-year frequency in one of the lakes, may relate to solar influence on local catchment processes. Following thousands of years of restricted extent, Istorvet ice cap advanced to within 365 m of its late Holocene limit at ~AD 1150. Variability in the timing of glacial and climate fluctuations, as well as of sediment organic content changes among East Greenland lacustrine records, may be a consequence of local factors, such as elevation, continentality, water depth, turbidity, and seabirds, and highlights the need for a detailed spatial array of datasets to address questions about Holocene climate change.     

The geomagnetic field over the past 5 million years

The modern geomagnetic field is usually expressed as a spherical harmonic expansion. Although the palaeomagnetic record is very incomplete in both space and time, sufficient data are available from a span of ages to generate time-averaged spherical harmonic field models with many degrees of freedom. Here three data sets are considered: directional measurements from lavas, inclination measurements from ocean sediments, and intensity measurements from lavas. Individual data are analysed, as well as site-averages, using the same methods that have been developed for the modern field, to give models for the past 5 Myr. The normal-polarity field model has an axial-dipole intensity similar to that of the modern-day field, whilst the equatorial-dipole component is very much smaller. The field is not axisymmetric, but shows flux concentrations at the core's surface under Canada and Siberia similar to those observed in the field over historical timescales. Tests on synthetic data show that it is unlikely that these similarities result from the overprinting of the palaeomagnetic field due to inadequate cleaning of the samples. The reverse-polarity field model does not show such obvious features, but this may be due to the sparsity of the data.
The patterns observed in the normal-polarity field, with persistent features in the northern hemisphere and a smooth southern hemisphere, could be explained if the present pattern of secular variation is typical of the past several million years. This would reveal itself as large variations over time in the direction of the magnetic vector in regions of high secular variation, with relatively little change over quieter regions. However, we have been unable to find any evidence for a geographical pattern of secular variation in the data.     

Kriging with imprecise (fuzzy) variograms. II: Application

The geostatistical analysis of soil liner permeability is based on 20 measurements and imprecise prior information on nugget effect, sill, and range of the unknown variogram. Using this information, membership functions for variogram parameters are assessed and the fuzzy variogram is constructed. Both kriging estimates and estimation variances are calculated as fuzzy numbers from the fuzzy variogram and data points. Contour maps are presented, indicating values of the kriged permeability and the estimation variance corresponding to selected membership values called levels.