Estimating sagebrush cover in semi-arid environments using Landsat Thematic Mapper data

Sagebrush ecosystems of the western US provide important habitat for several ungulate and vertebrate species. As a consequence of energy development, these ecosystems in Wyoming have been subjected to a variety of anthropogenic disturbances. Land managers require methodology that will allow them to consistently catalog sagebrush ecosystems and evaluate potential impact of proposed anthropogenic activities. This study addresses the utility of remotely sensed and ancillary geospatial data to estimate sagebrush cover using ordinal logistic regression. We demonstrate statistically significant prediction of ordinal sagebrush cover categories using spectral (χ² = 113; p < 0.0001) and transformed indices (χ² = 117; p < 0.0001). Both Landsat spectral bands (c-value = 0.88) and transformed indices (c-value = 0.89) can distinguish sites with closed, moderate and open cover sagebrush cover categories from no cover. The techniques described in this study can be used for estimating categories of sagebrush cover in arid ecosystems.     

Toxicity assessment of sediments associated with various land-uses in coastal South Carolina, USA, using a meiobenthic copepod bioassay

Coastal urbanization supplies surrounding estuarine environments with urban-related contaminants such as polycyclic aromatic hydrocarbons (PAHs), metals and pesticide mixtures. In our study, adult female and male copepods (Amphiascus tenuiremis) were chronically exposed to 18 sediment samples collected from sites influenced by different land-uses. Sediment samples were collected from three major suburban areas (Hilton Head, Kiawah Island and the Okatee River watershed) and a pristine site (North Inlet) in coastal South Carolina. Three-sediment bioassays (six sites per bioassay) were conducted by culturing copepods for 14 days in quadruplicate test sediments under flow-through conditions at 20 degrees C and 12:12 LD cycles. Adult survival and copepod reproductive output were quantified. Sediment samples were also analyzed by GC-MS for low and high molecular weight PAHs. Minimal adult mortality was observed in most sediment samples. However, sediments from Hilton Head Island and the Okatee River showed significant effects on copepod reproductive output (i.e., nauplii, copepodites and clutch size). Thus, we determined that reproductive endpoints rather than adult copepod survivorship were more sensitive to effects of contaminated sediments on A. tenuiremis. Furthermore, six (33%) of the 18 sites had a >25% reduction in copepod bioassay endpoints relative to controls, suggesting a high risk to long term A. tenuiremis population maintenance.     

Resonant Periodic Orbits of Trans-Neptunian Objects

We study two and three-dimensional resonant periodic orbits, usingthe model of the restricted three-body problem with the Sun andNeptune as primaries. The position and the stability character ofthe periodic orbits determine the structure of the phase space andthis will provide useful information on the stability and longterm evolution of trans-Neptunian objects. The circular planarmodel is used as the starting point. Families of periodic orbitsare computed at the exterior resonances 1/2, 2/3 and 3/4 withNeptune and these are used as a guide to select the energy levelsfor the computation of the Poincaré maps, so that all basicresonances are included in the study. Using the circular planarmodel as the basic model, we extend our study to more realisticmodels by considering an elliptic orbit of Neptune and introducingthe inclination of the orbit. Families of symmetric periodicorbits of the planar elliptic restricted three-body problem andthe three-dimensional problem are found. All these orbitsbifurcate from the families of periodic orbits of the planarcircular problem. The stability of all orbits is studied. Althoughthe resonant structure in the circular problem is similar for allresonances, the situation changes if the eccentricity of Neptuneor the inclination of the orbit is taken into account. All theseresults are combined to explain why in some resonances there aremany bodies and other resonances are empty.     

Numerical modeling of toxic nonaqueous phase liquid removal from contaminated groundwater systems: mesh effect and discretization error estimation

Numerical modeling has now become an indispensable tool for investigating the fundamental mechanisms of toxic nonaqueous phase liquid (NAPL) removal from contaminated groundwater systems. Because the domain of a contaminated groundwater system may involve irregular shapes in geometry, it is necessary to use general quadrilateral elements, in which two neighbor sides are no longer perpendicular to each other. This can cause numerical errors on the computational simulation results due to mesh discretization effect. After the dimensionless governing equations of NAPL dissolution problems are briefly described, the propagation theory of the mesh discretization error associated with a NAPL dissolution system is first presented for a rectangular domain and then extended to a trapezoidal domain. This leads to the establishment of the finger‐amplitude growing theory that is associated with both the corner effect that takes place just at the entrance of the flow in a trapezoidal domain and the mesh discretization effect that occurs in the whole NAPL dissolution system of the trapezoidal domain. This theory can be used to make the approximate error estimation of the corresponding computational simulation results. The related theoretical analysis and numerical results have demonstrated the following: (1) both the corner effect and the mesh discretization effect can be quantitatively viewed as a kind of small perturbation, which can grow in unstable NAPL dissolution systems, so that they can have some considerable effects on the computational results of such systems; (2) the proposed finger‐amplitude growing theory associated with the corner effect at the entrance of a trapezoidal domain is useful for correctly explaining why the finger at either the top or bottom boundary grows much faster than that within the interior of the trapezoidal domain; (3) the proposed finger‐amplitude growing theory associated with the mesh discretization error in the NAPL dissolution system of a trapezoidal domain can be used for quantitatively assessing the correctness of computational simulations of NAPL dissolution front instability problems in trapezoidal domains, so that we can ensure that the computational simulation results are controlled by the physics of the NAPL dissolution system, rather than by the numerical artifacts. Copyright © 2014 John Wiley & Sons, Ltd.     

High-resolution climate simulation of the last glacial maximum

The climate of the last glacial maximum (LGM) is simulated with a high-resolution atmospheric general circulation model, the NCAR CCM3 at spectral truncation of T170, corresponding to a grid cell size of roughly 75 km. The purpose of the study is to assess whether there are significant benefits from the higher resolution simulation compared to the lower resolution simulation associated with the role of topography. The LGM simulations were forced with modified CLIMAP sea ice distribution and sea surface temperatures (SST) reduced by 1°C, ice sheet topography, reduced CO₂, and 21,000 BP orbital parameters. The high-resolution model captures modern climate reasonably well, in particular the distribution of heavy precipitation in the tropical Pacific. For the ice age case, surface temperature simulated by the high-resolution model agrees better with those of proxy estimates than does the low-resolution model. Despite the fact that tropical SSTs were only 2.1°C less than the control run, there are many lowland tropical land areas 4–6°C colder than present. Comparison of T170 model results with the best constrained proxy temperature estimates (noble gas concentrations in groundwater) now yield no significant differences between model and observations. There are also significant upland temperature changes in the best resolved tropical mountain belt (the Andes). We provisionally attribute this result in part as resulting from decreased lateral mixing between ocean and land in a model with more model grid cells. A longstanding model-data discrepancy therefore appears to be resolved without invoking any unusual model physics. The response of the Asian summer monsoon can also be more clearly linked to local geography in the high-resolution model than in the low-resolution model; this distinction should enable more confident validation of climate proxy data with the high-resolution model. Elsewhere, an inferred salinity increase in the subtropical North Atlantic may have significant implications for ocean circulation changes during the LGM. A large part of the Amazon and Congo Basins are simulated to be substantially drier in the ice age—consistent with many (but not all) paleo data. These results suggest that there are considerable benefits derived from high-resolution model regarding regional climate responses, and that observationalists can now compare their results with models that resolve geography at a resolution comparable to that which the proxy data represent.     

Transport of water in frozen soil: I: Experimental determination of soil-water diffusivity under isothermal conditions

A new experimental method for measuring the soil-water diffusivity of frozen soil under isothermal conditions is introduced. The theoretical justification of the method is presented and the feasibility of the method is demonstrated by experiments conducted using marine deposited clay. The measured values of the soil-water diffusivity are found comparable to reported experimental data.     

Calculating surface water PCO2 from foraminiferal organic δC

The δ¹³C of organic matter bound within the crystal lattice of foraminiferal calcite tests may provide a potential tracer of the isotopic composition of the surface water primary photosynthate. Using δ¹³C of the organic matter extracted from the crystal lattice and the calcite test, it is theoretically possible to estimate the paleo-surface water pCO₂. We have tailored this technique initially for the subpolar planktonic foraminifera species Globigerina bulloides. Initial surface water pCO₂ estimates from deep-sea core BOFS 5K (50°41.3′N, 21°51.9′W, water depth 3547 m) indicate that the northeast Atlantic Ocean may have been a greater sink for CO₂ during the last glacial than during the Holocene. Greatly reduced benthic foraminifera abundances, especially phytodetritus feeders, in BOFS 5K during the last glacial indicates low surface productivity. This rules out a productivity-driven CO₂ sink. The enhanced glacial CO₂ sink must, therefore, have results from a southwards shift of the centre of deep water formation.     

Oxidation of the Kaapvaal lithospheric mantle driven by metasomatism

The oxidation state, reflected in the oxygen fugacity (fO₂), of the subcratonic lithospheric mantle is laterally and vertically heterogeneous. In the garnet stability field, the Kaapvaal lithospheric mantle becomes progressively more reducing with increasing depth from Δlog fO₂ FMQ-2 at 110 km to FMQ-4 at 210 km. Oxidation accompanying metasomatism has obscured this crystal-chemical controlled depth-fO₂ trend in the mantle beneath Kimberley, South Africa. Chondrite normalized REE patterns for garnets, preserve evidence of a range in metasomatic enrichment from mild metasomatism in harzburgites to extensive metasomatism by LREE-enriched fluids and melts with fairly unfractionated LREE/HREE ratios in phlogopite-bearing lherzolites. The metasomatized xenoliths record redox conditions extending up to Δlog fO₂ = FMQ, sufficiently oxidized that magnesite would be the stable host of carbon in the most metasomatized samples. The most oxidized lherzolites, those in or near the carbonate stability field, have the greatest modal abundance of phlogopite and clinopyroxene. Clinopyroxene is modally less abundant or absent in the most reduced peridotite samples. The infiltration of metasomatic fluids/melts into diamondiferous lithospheric mantle beneath the Kaapvaal craton converted reduced, anhydrous harzburgite into variably oxidized phlogopite-bearing lherzolite. Locally, portions of the lithospheric mantle were metasomatized and oxidized to an extent that conversion of diamond into carbonate should have occurred. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Distinct modes of bidecadal and multidecadal variability in a climate reconstruction of the last centuries from a South Pacific coral

The Rarotonga coral Sr/Ca time series (Linsley et al. in Science 290:1145–1148, 2000) provides a near-monthly resolved proxy record of South Pacific climate variability over the last ~300 years. Here we show that two distinct interdecadal, quasi-periodic time components with periods of ~80 and ~25 years can be identified in this time series by Singular Spectrum Analysis. Their associated spatial patterns in the global sea surface temperature (SST) field show notable differences. Whereas the multidecadal component is associated with a global SST pattern that was recently associated with solar forcing on multidecadal timescales, the bidecadal component is associated with a well known pattern of Pacific decadal to interdecadal SST variability.