Depositional history and evolution of the Paso del Indio site,Vega Baja,Puerto Rico

Potshards discovered during excavation of bridge pilasters for a major expressway over the Rio Indio floodplain, a stream incised within the karsts of north‐central Puerto Rico, required large‐scale archaeological excavation. Five‐meter‐deep bridge pilaster excavations in the alluvial valley provide a 4500‐year history of deposition. Stratigraphic analysis of the exposed pilaster walls in combination with textural and organic carbon analyses of sediment cores obtained over a much broader area suggest a fluvial system dominated by overbank deposition. Six sequences of alternating light and dark layers of sediment were identified. The darker layers are largely composed of silts and clays, whereas the lighter layers are rich in sand‐sized sediment. Archaeological evidence indicates the organic‐rich dark layers, believed to be buried A horizons, coincide with pre‐historic occupation by Cedrosan Saladoid, Elenan Ostionoid, and Chican Ostionoid, extending from A.D. 450 to A.D. 1500. Lighter layers below the dark soil horizons are interpreted as overbank deposits from large magnitude flood events. The floodplain aggraded discontinuously with rapid deposition of sand followed by gradual accumulation of silt, clay, and organic material. An approximately 1‐m‐thick layer of coarse sand and gravel halfway up the stratigraphic column represents an episode of more frequent and severe floods. Based on radiocarbon ages, this layer aggraded between A.D. 1000 and A.D. 1100, which is well within the Elenan Ostionoid era (A.D. 900–1200). Rates of sedimentation during this period were approximately 8 mm per year, ten times greater than the estimates of sedimentation rates before and after this flood sequence. The cause for the change in deposition is unknown. Nonetheless the Elenan Ostionoid would have had to endure frequent loss of habitation structures and crops during these events. © 2003 Wiley Periodicals, Inc.     

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.     

Extent and magnitude of copper contamination in marinas of the San Diego region, California, USA

Marinas are areas of special water quality concern because of the potential for pollutant accumulation within their protected waters. Perhaps the largest contaminant source to marinas is antifouling paints that leach copper to prevent the growth of encrusting organisms on vessel bottoms. Very little monitoring of marinas is typically conducted despite the potential environmental risk, particularly in the San Diego region of California, USA where as many as 17,000 recreational vessels are berthed. The objective of this study was twofold: (1) determine the extent and magnitude of dissolved copper concentrations in marinas throughout the San Diego region, and (2) determine if elevated copper concentrations in marinas of the San Diego region are resulting in adverse biological impacts. A probabilistic study design was used to sample water column copper concentrations and toxicity (using Mytilus galloprovincialis) at 30 stations. Results indicated that exceedence of state water quality objectives was widespread (86% of marina area), but that toxicity was much less prevalent (21% of marina area). Toxicity identification evaluations (TIEs) conducted at the most toxic sites indicated that toxicity was largely due to trace metals, most likely copper. Toxicity was reduced using TIE treatments that chelated trace metals such as cation exchange column, ethylenediaminetetraacetic acid (EDTA), and sodium thiosulfate (STS). Moreover, increasing dissolved copper concentrations correlated with increasing toxicity and these copper concentrations were high enough to account for virtually all of the observed toxicity.     

Real-time measurement of volcanic SO2 emissions: validation of a new UV correlation spectrometer (FLYSPEC)

A miniaturized, lightweight and low-cost UV correlation spectrometer, the FLYSPEC, has been developed as an alternative for the COSPEC, which has long been the mainstay for monitoring volcanic sulfur dioxide fluxes. Field experiments have been conducted with the FLYSPEC at diverse volcanic systems, including Masaya (Nicaragua), Poás (Costa Rica), Stromboli, Etna and Vulcano (Italy), Villarica (Chile) and Kilauea (USA). We present here those validation measurements that were made simultaneously with COSPEC at Kilauea between March 2002 and February 2003. These experiments, with source emission rates that ranged from 95 to 1,560 t d⁻¹, showed statistically identical results from both instruments. SO₂ path-concentrations ranged from 0 to >1,000 ppm-m with average correlation coefficients greater than r ²=0.946. The small size and low cost create the opportunity for FLYSPEC to be used in novel deployment modes that have the potential to revolutionize the manner in which volcanic and industrial monitoring is performed.     

Remote sensing documentation of historic rangeland remediation treatments in southern New Mexico

The Jornada Experimental Range and the New Mexico State University Chihuahuan Desert Rangeland Research Center are fruitful areas to study the long-term effects of rangeland remediation treatments which started in the 1930s. A number of diverse manipulations were completed under the direction of federal agency and university scientists, and abundant remote sensing imagery is available to assist in relocating the treatments and evaluating their success. This is particularly important because few of the treatments were maintained following the loss of scientific personnel coinciding with the start of World War II, and most records of Civilian Conservation Corps scientific work were lost with the disbanding of the agency in 1942. Aerial photography, which was systematically used to image the United States beginning in the 1930s, can be used to identify types of treatments, measure areal coverage, estimate longevity, and help plan locations for new experiments. No long-lasting vegetation response could be determined for contour terraces, brush water spreaders, strips grubbed free of shrubs (despite the fact that these strips have remained visible for 65 years), and mechanical rootplowing and seeding. Distinct positive, long-term vegetation responses could be seen in aerial photos for water retention dikes, certain fenced exclosures, and some boundaries where different land management practices meet. It appears from both aerial photos and existing conventional records that experimental manipulation of rangelands has often been ineffective on the landscape scale because treatments are not performed over large enough contiguous areas and hydrological and ecological processes overwhelm the treatments. In addition, treatments are not maintained over time, treatment evaluation periods are sometimes too short, multi-purpose treatments are not used to maximize effects, and treatments are often not located in appropriate sites.     

An assessment of large-scale flood modelling based on LiDAR data

Large-scale flood modelling approaches designed for regional to continental scales usually rely on relatively simple assumptions to represent the potentially highly complex river bathymetry at the watershed scale based on digital elevation models (DEMs) with a resolution in the range of 25–30 m. Here, high-resolution (1 m) LiDAR DEMs are employed to present a novel large-scale methodology using a more realistic estimation of bathymetry based on hydrogeomorphological GIS tools to extract water surface slope. The large-scale 1D/2D flood model LISFLOOD-FP is applied to validate the simulated flood levels using detailed water level data in four different watersheds in Quebec (Canada), including continuous profiles over extensive distances measured with the HydroBall technology. A GIS-automated procedure allows to obtain the average width required to run LISFLOOD-FP. The GIS-automated procedure to estimate bathymetry from LiDAR water surface data uses a hydraulic inverse problem based on discharge at the time of acquisition of LiDAR data. A tiling approach, allowing several small independent hydraulic simulations to cover an entire watershed, greatly improves processing time to simulate large watersheds with a 10-m resampled LiDAR DEM. Results show significant improvements to large-scale flood modelling at the watershed scale with standard deviation in the range of 0.30 m and an average fit of around 90%. The main advantage of the proposed approach is to avoid the need to collect expensive bathymetry data to efficiently and accurately simulate flood levels over extensive areas.     

Hydrological Outlook UK: an operational streamflow and groundwater level forecasting system at monthly to seasonal time scales

This paper describes the development of the first operational seasonal hydrological forecasting service for the UK, the Hydrological Outlook UK (HOUK). Since June 2013, this service has delivered monthly forecasts of streamflow and groundwater levels, with an emphasis on forecasting hydrological conditions over the next three months, accompanied by outlooks over longer time horizons. This system is based on three complementary approaches combined to produce the outlooks: (i) national-scale modelling of streamflow and groundwater levels based on dynamic seasonal rainfall forecasts, (ii) catchment-scale modelling where streamflow and groundwater level models are driven by historical meteorological forcings (i.e. the Ensemble Streamflow Prediction, ESP, approach), and (iii) a catchment-scale statistical method based on persistence and historical analogues. This paper provides the background to the Hydrological Outlook, describes the various component methods in detail and then considers the impact and usefulness of the product. As an example of a multi-method, operational seasonal hydrological forecasting system, it is hoped that this overview provides useful information and context for other forecasting initiatives around the world.     

Using UAV acquired photography and structure from motion techniques for studying glacier landforms: application to the glacial flutes at Isfallsglaciären

Glacier and ice sheet retreat exposes freshly deglaciated terrain which often contains small‐scale fragile geomorphological features which could provide insight into subglacial or submarginal processes. Subaerial exposure results in potentially rapid landscape modification or even disappearance of the minor‐relief landforms as wind, weather, water and vegetation impact on the newly exposed surface. Ongoing retreat of many ice masses means there is a growing opportunity to obtain high resolution geospatial data from glacier forelands to aid in the understanding of recent subglacial and submarginal processes. Here we used an unmanned aerial vehicle to capture close‐range aerial photography of the foreland of Isfallsglaciären, a small polythermal glacier situated in Swedish Lapland. An orthophoto and a digital elevation model with ~2 cm horizontal resolution were created from this photography using structure from motion software. These geospatial data was used to create a geomorphological map of the foreland, documenting moraines, fans, channels and flutes. The unprecedented resolution of the data enabled us to derive morphological metrics (length, width and relief) of the smallest flutes, which is not possible with other data products normally used for glacial landform metrics mapping. The map and flute metrics compare well with previous studies, highlighting the potential of this technique for rapidly documenting glacier foreland geomorphology at an unprecedented scale and resolution. The vast majority of flutes were found to have an associated stoss‐side boulder, with the remainder having a likely explanation for boulder absence (burial or erosion). Furthermore, the size of this boulder was found to strongly correlate with the width and relief of the lee‐side flute. This is consistent with the lee‐side cavity infill model of flute formation. Whether this model is applicable to all flutes, or multiple mechanisms are required, awaits further study. © 2016 The Authors. Earth Surface Processes and Landforms published by John Wiley & Sons Ltd.     

Carbon,nitrogen, and water stable isotopes in plant tissue and soils across a moisture gradient in Puerto Rico

Stable isotopes in the water molecule (²H or D and ¹⁸O), carbon, and nitrogen are useful tracers and integrators of processes in plant ecohydrological systems across scales. Over the last few years, there has been growing interest in regional to continental scale synthesis of stable isotope data with a view to elucidating biogeochemical and ecohydrological patterns. Published datasets from the humid tropics, however, are limited. To be able to contribute to bridging the “data gap” in the humid tropics, here, we publish a relatively novel and unique suite of δ¹³C, δ¹⁵N, δ²H, and δ¹⁸O isotope data from three sites across a moisture gradient and contrasting land use in Puerto Rico. Plant tissue (xylem and leaf) samples from two species of mahogany (Swietenia macrophylla and Swietenia mahagoni) and soil samples down to 60 cm in the soil profile were collected in relatively “wet” (July 2012) and “dry” (February 2013) periods at two sites in northeastern (Luquillo) and southwestern (Susua) Puerto Rico. The same sampling suite is also being made available from a highly urbanized site in the capital San Juan. Leaf samples taken in July 2012 and February 2013 were analyzed for δ¹³C and δ¹⁵N; all xylem and bulk soil samples were analyzed for δ²H and δ¹⁸O. Soil samples taken in July 2012 were analyzed for δ¹³C and δ¹⁵N. Leaf δ¹⁵N and δ¹³C dataset showed patterns that are possibly associated with site differences. While spatial patterns were also apparent in soil δ¹⁵N and δ¹³C dataset, the positively linear δ¹⁵N –δ¹³C relationship tends to weaken with site moisture. Soil depth and site moisture patterns were also observed in the δ²H and δ¹⁸O datasets of bulk soil and xylem samples. The purpose of these datasets is to provide baseline information on soil–plant water (δ²H and δ¹⁸O, N = 319), δ¹³C (N = 272), and δ¹⁵N (N = 269) that may be useful in a wide range of research questions from ecohydrological relations to biogeochemical patterns in soils and vegetation.