Quantitative phosphorus measurement: A field test procedure for archaeological site analysis at Piedras Negras,Guatemala

Currently there is a wide interest in the use of chemical analyses for the evaluation of anthropogenically altered soils and other archaeological deposits. Because soil phosphorus levels increase in areas of human habitation, and leave a permanent signature that can only be removed by erosion of the soil itself, phosphorus mapping has become a popular field procedure to indicate areas of habitation where overt evidence of ancient occupance is absent. We have developed a methodology to obtain accurate acid‐extractable phosphorus concentrations (mg/kg) in calcareous soils under the primitive field conditions of Piedras Negras, Guatemala. Predicated on Mehlich‐II acid extractant and colorimetric methods, this procedure processed 36 samples per hour at very low cost per sample. Based on eight replicate measurements of a group of samples, the coefficient of variation of the procedure was 8.3%. Subsequent analysis of 35 soil samples in a controlled laboratory revealed a moderate correlation of 0.44 between the Mehlich‐extractable phosphorus and total phosphorus. The correlation was 0.91 between the Mehlich procedure and Olsen bicarbonate extractable phosphorus, indicating that Mehlich‐based results are similar to those obtainable using a traditional extractable phosphorus method on soils of neutral to alkaline pH. There was a moderate correlation between Mehlich P and ring‐test rating (r = 0.42). The wider dynamic range of the Mehlich extraction, coupled with the use of a battery‐operated colorimeter, facilitated the finding of a refuse midden within an area of phosphate enriched soils. Further tests indicated that phosphorus concentrations measured in the field deviated by only 7% from those made under controlled laboratory conditions. © 2000 John Wiley & Sons, Inc.     

Simulating the impact of the large-scale circulation on the 2-m temperature and precipitation climatology

The impact of the simulated large-scale atmospheric circulation on the regional climate is examined using the Weather Research and Forecasting (WRF) model as a regional climate model. The purpose is to understand the potential need for interior grid nudging for dynamical downscaling of global climate model (GCM) output for air quality applications under a changing climate. In this study we downscale the NCEP-Department of Energy Atmospheric Model Intercomparison Project (AMIP-II) Reanalysis using three continuous 20-year WRF simulations: one simulation without interior grid nudging and two using different interior grid nudging methods. The biases in 2-m temperature and precipitation for the simulation without interior grid nudging are unreasonably large with respect to the North American Regional Reanalysis (NARR) over the eastern half of the contiguous United States (CONUS) during the summer when air quality concerns are most relevant. This study examines how these differences arise from errors in predicting the large-scale atmospheric circulation. It is demonstrated that the Bermuda high, which strongly influences the regional climate for much of the eastern half of the CONUS during the summer, is poorly simulated without interior grid nudging. In particular, two summers when the Bermuda high was west (1993) and east (2003) of its climatological position are chosen to illustrate problems in the large-scale atmospheric circulation anomalies. For both summers, WRF without interior grid nudging fails to simulate the placement of the upper-level anticyclonic (1993) and cyclonic (2003) circulation anomalies. The displacement of the large-scale circulation impacts the lower atmosphere moisture transport and precipitable water, affecting the convective environment and precipitation. Using interior grid nudging improves the large-scale circulation aloft and moisture transport/precipitable water anomalies, thereby improving the simulated 2-m temperature and precipitation. The results demonstrate that constraining the RCM to the large-scale features in the driving fields improves the overall accuracy of the simulated regional climate, and suggest that in the absence of such a constraint, the RCM will likely misrepresent important large-scale shifts in the atmospheric circulation under a future climate.     

Dust motions in quasi-statically charged binary asteroid systems

In this paper, we discuss dust motion and investigate possible mass transfer of charged particles in a binary asteroid system, in which the asteroids are electrically charged due to solar radiation. The surface potential of the asteroids is assumed to be a piecewise function, with positive potential on the sunlit half and negative potential on the shadow half. We derive the nonautonomous equations of motion for charged particles and an analytic representation for their lofting conditions. Particle trajectories and temporary relative equilibria are examined in relation to their moving forbidden regions, a concept we define and discuss. Finally, we use a Monte Carlo simulation for a case study on mass transfer and loss rates between the asteroids.     

Radio-continuum study of the nearby Sculptor group Galaxies. Part 2: NGC 55 at λ=20, 13, 6 and 3 cm

A series of new radio-continuum (λ=20, 13, 6 and 3 cm) mosaic images focused on the NGC 55 galactic system were produced using archived observational data from the Australia Telescope Compact Array. These new images are both very sensitive (down to rms=33 μJy) and feature high angular resolution (down to <4″). Using these newly created images, 66 previously unidentified discrete sources are identified. Of these sources, 46 were classified as background sources, 11 as H?ii regions and 6 as supernova remnant candidates. This relatively low number of SNR candidates detected coupled with the low number of large H?ii regions is consistent with the estimated low star formation rate of the galaxy at 0.06 M _⊙?year^?1. Our spectral index map shows that the core of galaxy appears to have a shallow spectral index between α=?0.2 and ?0.4. This indicates that the core of the galaxy is a region of high thermal radiation output.     

Challenges and lessons from a wetland LiDAR project: a case study of the Okefenokee Swamp,Georgia, USA

Airborne LiDAR (Light Detection and Ranging) provides opportunities to generate high-quality digital elevation models (DEMs) even in wetland environments. Our project, performed over the Okefenokee Swamp in Georgia during the spring of 2010, shows that several, distinctive factors must be considered to ensure successful wetland LiDAR projects. Some of the challenges include selecting optimal flight times in accordance with weather variability and water levels, having effective and quality control protocols, applying and developing filtering and interpolation algorithms, breaklines in swamps and accounting for data striping and noise. While some of these issues are faced in any airborne LiDAR acquisition, many of these require special consideration in a low-slope wetland environment with water saturated soils, widespread shallow water, and sediments and extensive vegetation. An examination of these issues and how they were handled will help in ensuring the success of future LiDAR acquisitions and, in particular, will advance knowledge of producing quality DEMs in wetland environments.     

Developing a parallel computational implementation of AMOEBA

As geospatial researchers' access to high-performance computing clusters continues to increase alongside the availability of high-resolution spatial data, it is imperative that techniques are devised to exploit these clusters' ability to quickly process and analyze large amounts of information. This research concentrates on the parallel computation of A Multidirectional Optimal Ecotope-Based Algorithm (AMOEBA). AMOEBA is used to derive spatial weight matrices for spatial autoregressive models and as a method for identifying irregularly shaped spatial clusters. While improvements have been made to the original ‘exhaustive’ algorithm, the resulting ‘constructive’ algorithm can still take a significant amount of time to complete with large datasets. This article outlines a parallel implementation of AMOEBA (the P-AMOEBA) written in Java utilizing the message passing library MPJ Express. In order to account for differing types of spatial grid data, two decomposition methods are developed and tested. The benefits of using the new parallel algorithm are demonstrated on an example dataset. Results show that different decompositions of spatial data affect the computational load balance across multiple processors and that the parallel version of AMOEBA achieves substantially faster runtimes than those reported in related publications.     

Snowmelt runoff and soil moisture dynamics on steep subalpine hillslopes

Snowmelt water supplies streamflow and growing season soil moisture in mountain regions, yet pathways of snowmelt water and their effects on moisture patterns are still largely unknown. This study examined how flow processes during snowmelt runoff affected spatial patterns of soil moisture on two steep sub‐alpine hillslope transects in Rocky Mountain National Park, CO, USA. The transects have northeast‐facing and east‐facing aspects, and both extend from high‐elevation bedrock outcrops down to streams in valley bottoms. Spatial patterns of both snow depth and near‐surface soil moisture were surveyed along these transects in the snowmelt and summer seasons of 2008–2010. To link these patterns to flow processes, soil moisture was measured continuously on both transects and compared with the timing of discharge in nearby streams. Results indicate that both slopes generated shallow lateral subsurface flow during snowmelt through near‐surface soil, colluvium and bedrock fractures. On the northeast‐facing transect, this shallow subsurface flow emerged through mid‐slope seepage zones, in some cases producing saturation overland flow, whereas the east‐facing slope had no seepage zones or overland flow. At the hillslope scale, earlier snowmelt timing on the east‐facing slope led to drier average soil moisture conditions than on the northeast‐facing slope, but within hillslopes, snow patterns had little relation to soil moisture patterns except in areas with persistent snow drifts. Results suggest that lateral flow and exfiltration processes are key controls on soil moisture spatial patterns in this steep sub‐alpine location. Copyright © 2014 John Wiley & Sons, Ltd.     

Distance in the space of energetically bounded Keplerian orbits

In this paper we derive an explicit, analytic formula for the geodesic distance between two points in the space of bounded Keplerian orbits in a particular topology. The specific topology we use is that of a cone passing through the direct product of two spheres. The two spheres constitute the configuration manifold for the space of bounded orbits of constant energy. We scale these spheres by a factor equal to the semi-major axis of the orbit, forming a linear cone. This five-dimensional manifold inherits a Riemannian metric, which is induced from the Euclidean metric on \mathbbR⁶{\mathbb{R}^6}, the space in which it is embedded. We derive an explicit formula for the geodesic distance between any two points in this space, each point representing a physical, gravitationally bound Keplerian orbit. Finally we derive an expression for the Riemannian metric that we used in terms of classical orbital elements, which may be thought of as local coordinates on our configuration manifold.     

Characterisation of glacial sediments using geophysical methods for groundwater source protection

A sequence of glacial and alluvial deposits overlying the Cretaceous Chalk in Eastern England was characterised using two geophysical techniques: electrical resistivity imaging and electromagnetic (EM) induction. Extensive geological data were available from trenching and boreholes. Synthetic modelling of the electrical resistivity imaging technique was undertaken to identify its limitations and to optimise survey design. The EM induction method provided a quick and cost-effective reconnaissance technique for identifying large-scale lateral variation in lithology, and for siting resistivity profiles and further boreholes. The resistivity imaging technique provided detailed information on the vertical continuity of permeable units, and was able to identify permeable pathways through the sequence. Certain limitations in detecting thin sand or gravel layers underlying electrically conductive clay were seen in both the synthetic and field data. Nevertheless, the study shows that knowledge of these limitations allows interpretation for the purpose of groundwater vulnerability assessment, given that an appropriate amount of invasive investigation has been conducted.