A lunar rock of deep crustal origin: sample 76535

Lunar sample 76535 is a coarse-grained troctolitic granulite exhibiting a texture indicative of long annealing times. It is composed of homogeneous crystals of plagioclase (58 per cent, An₉₆), olivine (37 per cent, Fo₈₈) and bronzite (4 per cent, En₈₆).Chromian spinel-bronzite-diopside (Wo₄₆En₅₀Fs₄) symplectic intergrowths commonly occur along olivine-plagioclase boundaries and as tiny inclusions within olivine grains. These symplectites apparently formed by a reaction of the type:

$\text{OI + An + Chromite → Opx + Cpx + Al-Mg-chromite}$

. The reaction is related to the experimentally determined reaction

$\text{OI + An = Opx + Cpx + Sp}$

of Kushiro and Yoder (1966). The enstatite content of the diopside coexisting with the bronzite indicates equilibration at about 1000°C. Thermodynamic calculations for 1000°C indicate that the symplectites formed at a minimum pressure of about 0.6 kb. Low alumina contents of the pyroxenes indicate equilibration near this minimum pressure.Clusters of the same assemblage found in the symplectic intergrowths, but containing accessory metal, troilite, Ca-phosphates, baddeleyite, plagioclase and/or K-feldspar occur sporadically throughout the rock. These apparent late stage products crystallized in the low temperature-high pressure region discussed above.Phase relations of co-existing metal phases indicate that the rock cooled at a few tens of degrees/my, corresponding to depths of 10–20 km below the lunar surface, in agreement with the above pressure estimate.We infer that 76535 represents an original cumulate deposited at a depth between about 10 and 30 km. The last liquid crystallized in the relatively high pressure-low temperature field opx + cpx + Al-Mg-chromite. Cooling was extremely slow and accompanied by extensive chemical and textural re-equilibration. Reaction to form the symplectites occurred during the late stages of re-equilibration.     

LNAPL Volume Calculation: Parameter Estimation by Nonlinear Regression of Saturation Profiles

An estimation of the volume of light nonaqueous phase liquids (LNAPL) is often required during site assessment, remedial design, or litigation. LNAPL volume can be estimated by a strictly empirical approach whereby core samples, distributed throughout the vertical and lateral extent of LNAPL, are analyzed for LNAPL content, and these data are then integrated to compute a volume. Alternatively, if the LNAPL has obtained vertical equilibrium, the thickness of LNAPL in monitoring wells can be used to calculate of LNAPL in monitoring wells can be used to calculate LNAPL volume at the well locations if appropriate soil and LNAPL properties can be estimated.
A method is described for estimating key soil and LNAPL properties by nonlinear regression of vertical profiles of LNAPL saturation. The methods is relatively fast, cost effective, and amenable to quantitative analysis of uncertainty. Optionally, the method allows statistical determination of best-fit values for the Van Genuchten capillary parameters (n, α_oil-water and α_oil-air), residual water saturation and ANAPL density. The sensitivity of the method was investigated by fitting field LNAPL saturation profiles and then determining the variation in misfit (mean square residual) as a function of parameter value for each parameter. Using field data from a sandy aquifer, the fitting statistics were found to be highly sensitive to LNAPL density, α_oil-water and α_oil-air moderately sensitive to the Van Genuchten n value, and weakly sensitive to residual water saturation. The regression analysis also provides information that can be used to estimate uncertainty in the estimated parameters, which can then be used to estimate uncertainty in calculated values of specific volume.     

Biological and Physical Attenuation of Endocrine Disruptors and Pharmaceuticals: Implications for Water Reuse

A select group of endocrine disrupters, pharmaceuticals, and personal care products was studied to determine the degree of biological attenuation in water reuse applications. Laboratory investigations involved both batch reactors using biologically active sand and continuous flow simulated aquifer storage and recovery experiments. All laboratory experiments were conducted using Colorado River water spiked with various target compounds at concentrations between 10 and 100 ng/L. Field studies were also conducted to determine the occurrence and attenuation of target compounds in water reuse applications. Two golf courses irrigated with reuse water were studied to determine if turf applications led to contamination of nearby ground water. A waste water treatment facility that uses rapid infiltration basins seasonally was also tested to determine the degree of attenuation of detectable target compounds along a subsurface flowpath. A qualitative structural activity relationship model was applied to the target compounds to predict the general rate of aerobic biological degradation. Significant attenuation of many target compounds was observed in both laboratory and field experiments. Conversely, several compounds displayed limited removal during these studies. Field experiments were limited to detectable compounds and various nonbiological removal or concentration effects that may impact data interpretations, which are discussed in this paper. The predictive model was found to be moderately accurate within the confines of the project scope.     

Symplectic integrators for long-term integrations in celestial mechanics

We discuss the use of symplectic integration algorithms in long-term integrations in the field of celestial mechanics. The methods' advantages and disadvantages (with respect to more common integration methods) are discussed. The numerical performance of the algorithms is evaluated using the 2-body and circular restricted 3-body problems. Symplectic integration methods have the advantages of linear phase error growth in the 2-body problem (unlike most other methods), good conservation of the integrals of the motion, good performance for moderately eccentric orbits, and ease of use. Its disadvantages include a relatively large number of force evaluations and an inability to continuously vary the step size.     

Scaling and regionalization of flood flows in British Columbia,Canada

A regionalization of flood data in British Columbia reveals a common scaling with drainage area over the range 0·5×10²<A_d<10⁴ km². This scaling is not a function of flood return period, which implies that simple scaling—consistent with a snowmelt‐dominated flow regime—applies to the province. The observed scale relation takes the form , similar to values reported in previous studies. The scaling relation identified was used to define the regional pattern of hydroclimatic variability for flood flows in British Columbia after discounting the effect of drainage area. The pattern was determined by kriging a scale‐independent runoff factor k for the mean annual flood, 5 year flood and 20 year flood. The analysis permits quantification of uncertainty of the estimates, which can be used in conjunction with the mapped k‐fields to calculate a mean and range for floods with the identified return period for ungauged basins. Owing to the sparsity of data, the precision is relatively poor. The standard error is generally less than 75% of the estimate in the southern half of the province, whereas in the northern half it is often between 75 and 100%. Examination of the relative increase in flood magnitude with increasing return period reveals spatially consistent but statistically insignificant differences. Flood magnitude tends to increase more rapidly in the western regions, where rain events may contribute to flood generation. The relative increase in flood magnitude with return period is consistently lower in the eastern mountain ranges, where snowmelt dominates the flood flow regime. Copyright © 2002 John Wiley & Sons, Ltd.     

Spatial Analysis of Environmental Change Impacts on Wheat Production in Mid-Lower North, South Australia

Three environmental change scenarios (the best scenario, the most likely scenario and the worst scenario) were used by the APSIM (Agricultural Production System sIMulator) Wheat module to study the possible impacts of future environmental change (climate change plus pCO₂ change) on wheat production in the Mid-Lower North of South Australia. GIS software was used to manage spatial-climate data and spatial-soil data and to present the results. Study results show that grain yield (kg ha⁻¹) was adversely affected under the worst environmental change scenario (−100% ∼ −42%) and the most likely environmental change scenario (−58% ∼ −3%). Grain nitrogen content (% N) either increased or decreased depending on the environmental change scenarios used and climate divisions (−25% ∼ +42%). Spatial variability was found for projected impact outcomes within climate divisions indicating the necessity of including the spatial distribution of soil properties in impact assessment.     

Climate-change impacts on sagebrush habitat and West Nile virus transmission risk and conservation implications for greater sage-grouse

Greater sage-grouse (Centrocercus urophasianus) are threatened by loss of sagebrush habitat and the spread of West Nile virus throughout much of their range in North America; yet, future impacts of climate change on these potential stressors have not been addressed. Here, we aim to quantify the potential impacts of climate change on the distribution of climatically suitable habitat for sagebrush and on transmission risk for West Nile virus in the eastern portion of the species’ range. We used Maxent to model the current and future climatically suitable habitat for two dominant sagebrush species in the study area, and we used a degree-day model to predict future West Nile virus transmission risk under likely climate-change scenarios. Our models suggest that areas with the highest future suitability for sagebrush habitat will be found in southwestern Wyoming and north-central Montana. The degree-day model suggests that greater sage-grouse in western portions of the study area, which are generally higher in elevation than where West Nile virus currently occurs, will see increasing risk of transmission in the future. We developed a spatially explicit map of suggested management actions based on our predictions that will aid in conservation of the species into the coming decades.     

The Phoenix Evening Transition Flow Experiment (TRANSFLEX)

Motivated by air quality and numerical modelling applications as well as recent theoretical advancements in the topic, a field experiment, dubbed transition flow experiment, was conducted in Phoenix, Arizona to study the evening transition in complex terrain (shift of winds from upslope to downslope). Two scenarios were considered: (i) the flow reversal due to a change of buoyancy of a cooled slab of air near the ground, and (ii) the formation of a transition front. A suite of in-situ flow, turbulence and particulate matter (PM) concentration sensors, vertically profiling tethered balloons and remote sensors were deployed, and a mesoscale numerical model provided guidance for interpreting observations. The results were consistent with the front formation mechanism, where it was also found that enhanced turbulence associated with the front increases the local PM concentration. During the transition period the flow adjustment was complex, involving the arrival of multiple fronts from different slopes, directional shear between fronts and episodic turbulent mixing events. The upward momentum diffusion from the incipient downslope flow was small because of stable stratification near the ground, and full establishment of downslope flow occurred over several hours following sunset. Episodic frontal events pose challenges to the modelling of the evening transition in complex terrain, requiring conditional parametrizations for subgrid scales. The observed increase of PM concentration during the evening transition has significant implications for the regulatory enforcement of PM standards for the area.     

Irregular satellite capture during planetary resonance passage

The passage of Jupiter and Saturn through mutual 1:2 mean-motion resonance has recently been put forward as explanation for their relatively high eccentricities [Tsiganis, K., Gomes, R., Morbidelli, A., Levison, H.F., 2005. Nature 435, 459-461] and the origin of Jupiter's Trojans [Morbidelli, A., Levison, H.F., Tsiganis, K., Gomes, R., 2005. Nature 435, 462-465]. Additional constraints on this event based on other small-body populations would be highly desirable. Since some outer satellite orbits are known to be strongly affected by the near-resonance of Jupiter and Saturn (“the Great Inequality”; ?uk, M., Burns, J.A., 2004b. Astron. J. 128, 2518-2541), the irregular satellites are natural candidates for such a connection. In order to explore this scenario, we have integrated 9200 test particles around both Jupiter and Saturn while they went through a resonance-crossing event similar to that described by Tsiganis et al. [Tsiganis, K., Gomes, R., Morbidelli, A., Levison, H.F., 2005. Nature 435, 459-461]. The test particles were positioned on a grid in semimajor axes and inclinations, while their initial pericenters were put at just 0.01 AU from their parent planets. The goal of the experiment was to find out if short-lived bodies, spiraling into the planet due to gas drag (or alternatively on orbits crossing those of the regular satellites), could have their pericenters raised by the resonant perturbations. We found that about 3% of the particles had their pericenters raised above 0.03 AU (i.e. beyond Iapetus) at Saturn, but the same happened for only 0.1% of the particles at Jupiter. The distribution of surviving particles at Saturn has strong similarities to that of the known irregular satellites. If saturnian irregular satellites had their origin during the 1:2 resonance crossing, they present an excellent probe into the early Solar System's evolution. We also explore the applicability of this mechanism for Uranus, and find that only some of the uranian irregular satellites have orbits consistent with resonant pericenter lifting. In particular, the more distant and eccentric satellites like Sycorax could be stabilized by this process, while closer-in moons with lower eccentricity orbits like Caliban probably did not evolve by this process alone.     

Thallium Mass Fraction and Stable Isotope Ratios of Sixteen Geological Reference Materials

Thallium stable isotope ratio and mass fraction measurements were performed on sixteen geological reference materials spanning three orders of magnitude in thallium mass fraction, including both whole rock and partially separated mineral powders. For stable isotope ratio measurements, a minimum of three independent digestions of each reference material was obtained. High‐precision trace element measurements (including Tl) were also performed for the majority of these RMs. The range of Tl mass fractions represented is 10 ng g^?1 to 16 μg g^?1, and Tl stable isotope ratios (reported for historical reasons as ε²⁰⁵Tl relative to NIST SRM 997) span the range ?4 to +2. With the exception – attributed to between‐bottle heterogeneity – of G‐2, the majority of data are in good agreement with published or certified values, where available. The precision of mean of independent measurement results between independent dissolutions suggests that, for the majority of materials analysed, a minimum digested mass of 100 mg is recommended to mitigate the impact of small‐scale powder heterogeneity. Of the sixteen materials analysed, we therefore recommend for use as Tl reference materials the USGS materials BCR‐2, COQ‐1, GSP‐2 and STM‐1; CRPG materials AL‐I, AN‐G, FK‐N, ISH‐G, MDO‐G, Mica‐Fe, Mica‐Mg and UB‐N; NIST SRM 607 and OREAS14P.