Geochemical tracers to evaluate hydrogeologic controls on river salinization

The salinization of rivers, as indicated by salinity increases in the downstream direction, is characteristic of arid and semiarid regions throughout the world. Historically, salinity increases have been attributed to various mechanisms, including (1) evaporation and concentration during reservoir storage, irrigation, and subsequent reuse; (2) displacement of shallow saline ground water during irrigation; (3) erosion and dissolution of natural deposits; and/or (4) inflow of deep saline and/or geothermal ground water (ground water with elevated water temperature). In this study, investigation of salinity issues focused on identification of relative salinity contributions from anthropogenic and natural sources in the Lower Rio Grande in the New Mexico-Texas border region. Based on the conceptual model of the system, the various sources of water and, therefore, salinity to the Lower Rio Grande were identified, and a sampling plan was designed to characterize these sources. Analysis results for boron (delta(11)B), sulfur (delta(34)S), oxygen (delta(18)O), hydrogen (delta(2)H), and strontium ((87)Sr/(86)Sr) isotopes, as well as basic chemical data, confirmed the hypothesis that the dominant salinity contributions are from deep ground water inflow to the Rio Grande. The stable isotopic ratios identified the deep ground water inflow as distinctive, with characteristic isotopic signatures. These analyses indicate that it is not possible to reproduce the observed salinization by evapotranspiration and agricultural processes alone. This investigation further confirms that proper application of multiple isotopic and geochemical tracers can be used to identify and constrain multiple sources of solutes in complex river systems.     

Inverse upscaling of hydraulic parameters during constant flux infiltration using borehole radar

Modeling unsaturated flow in porous media requires constitutive relations that describe the soil water retention and soil hydraulic conductivity as a function of either potential or water content. Often, the hydraulic parameters that describe these relations are directly measured on small soil cores, and many cores are needed to upscale to the entire heterogeneous flow field. An alternative to the forward upscaling method using small samples are inverse upscaling methods that incorporate soft data from geophysical measurements observed directly on the larger flow field. In this paper, we demonstrate that the hydraulic parameters can be obtained from cross borehole ground penetrating radar by measuring the first arrival travel time of electromagnetic waves (represented by raypaths) from stationary antennae during a constant flux infiltration experiment. The formulation and coupling of the hydrological and geophysical models rely on a constant velocity wetting front that causes critical refraction at the edge of the front as it passes by the antennae. During this critical refraction period, the slope of the first arrival data can be used to calculate (1) the wetting velocity and (2) the hydraulic conductivity of the wet (or saturated) soil. If the soil is undersaturated during infiltration, then an estimate of the saturated water content is needed before calculating the saturated hydraulic conductivity. The hydraulic conductivity value is then used in a nonlinear global optimization scheme to estimate the remaining two parameters of a Broadbridge and White soil.     

Dynamic adaptive pathways in downscaled climate change scenarios

The parallel scenario process enables characterization of climate-related risks and response options to climate change under different socio-economic futures and development prospects. The process is based on representative concentration pathways, shared socio-economic pathways, and shared policy assumptions. Although this scenario architecture is a powerful tool for evaluating the intersection of climate and society at the regional and global level, more specific context is needed to explore and understand risks, drivers, and enablers of change at the national and local level. We discuss the need for a stronger recognition of such national-scale characteristics to make climate change scenarios more relevant at the national and local scale, and propose ways to enrich the scenario architecture with locally relevant details that enhance salience, legitimacy, and credibility for stakeholders. Dynamic adaptive pathways are introduced as useful tools to draw out which elements of a potentially infinite scenario space connect with decision-relevant aspects of particular climate-related and non-climate-related risks and response options. Reviewing adaptation pathways for New Zealand case studies, we demonstrate how this approach could bring the global-scale scenario architecture within reach of local-scale decision-making. Such a process would enhance the utility of scenarios for mapping climate-related risks and adaptation options at the local scale, involving appropriate stakeholder involvement.     

Using chemical and microbiological indicators to track the impacts from the land application of treated municipal wastewater and other sources on groundwater quality in a karstic springs basin

Multiple chemical constituents (nutrients; N, O, H, C stable isotopes; 64 organic wastewater compounds, 16 pharmaceutical compounds) and microbiological indicators were used to assess the impact on groundwater quality from the land application of approximately 9.5 million liters per day of treated municipal sewage effluent to a sprayfield in the 960-km² Ichetucknee Springs basin, northern Florida. Enriched stable isotope signatures (δ¹⁸O and δ²H) were found in water from the effluent reservoir and a sprayfield monitoring well (MW-7) due to evaporation; however, groundwater samples downgradient from the sprayfield have δ¹⁸O and δ²H concentrations that represented recharge of meteoric water. Boron and chloride concentrations also were elevated in water from the sprayfield effluent reservoir and MW-7, but concentrations in groundwater decreased substantially with distance downgradient to background levels in the springs (about 12 km) and indicated at least a tenfold dilution factor. Nitrate-nitrogen isotope (δ¹⁵N–NO₃) values above 10 ‰ in most water samples were indicative of organic nitrogen sources except Blue Hole Spring (δ¹⁵N–NO₃ = 4.6–4.9 ‰), which indicated an inorganic source of nitrogen (fertilizers). The detection of low concentrations the insect repellent N,N-diethyl-metatoluamide (DEET), and other organic compounds associated with domestic wastewater in Devil’s Eye Spring indicated that leakage from a nearby septic tank drainfield likely has occurred. Elevated levels of fecal coliforms and enterococci were found in Blue Hole Spring during higher flow conditions, which likely resulted from hydraulic connections to upgradient sinkholes and are consistent with previoius dye-trace studies. Enteroviruses were not detected in the sprayfield effluent reservoir, but were found in low concentrations in water samples from a downgradient well and Blue Hole Spring during high-flow conditions indicating a human wastewater source. The Upper Floridan aquifer in the Ichetucknee Springs basin is highly vulnerable to contamination from multiple anthropogenic sources throughout the springs basin.     

Usage of an early warning and information system Web-site for real-time seismicity in Iceland

Iceland has been subjected to destructive earthquakes and volcanic eruptions throughout history. Such events are often preceded by changes in earthquake activity over varying timescales. Although most seismicity is confined to micro-earthquakes, large earthquakes have occurred within populated regions. Following the most recent hazardous earthquakes in 2000, the Icelandic Meteorological Office (IMO) developed an early warning and information system (EWIS) Web-site for viewing near-real-time seismicity in Iceland. Here we assess Web-site usage data in relation to earthquake activity, as recorded by the South Iceland Lowland (SIL) seismic network. Between March 2005 and May 2006 the SIL seismic network recorded 12,583 earthquakes. During this period, the EWIS Web-site logged a daily median of 91 visits. The largest onshore event (M _L 4.2) struck 20 km from Reykjavík on 06 March 2006 and was followed by an immediate, upsurge in usage resulting in a total of 1,173 unique visits to the Web-site. The greatest cluster of large (≥M _L 3) events occurred 300 km offshore from Reykjavík in May 2005. Within this swarm, 9 earthquakes ≥M _L 3 were detected on 11 May 2005, resulting in the release of a media bulletin by IMO. During the swarm, and following the media bulletin, the EWIS Web-site logged 1,234 unique visits gradually throughout the day. In summary, the data reveal a spatial and temporal relationship between Web-site usage and earthquake activity. The EWIS Web-site is accessed immediately after the occurrence of a local earthquake, whereas distant, unfelt earthquakes generate gradual interest prompted by media bulletins and, possibly, other contributing factors. We conclude that the Internet is a useful tool for displaying seismic information in near-real-time, which has the capacity to help increase public awareness of natural hazards.     

Discussion of “Geology and diamond distribution of the 140/141 kimberlite, Fort à la Corne, central Saskatchewan, Canada”, by A. Berryman, B.H. Scott-Smith and B.C. Jellicoe (Lithos v. 76, p. 99–114)

A wide variety of geological data and geological observations by numerous geoscientists do not support a two-stage crater excavation and in-fill model, or a champagne glass-shaped geometry for the 169 or 140/141 kimberlite bodies in the Fort à la Corne kimberlite field, Saskatchewan as described by Berryman, A., Scott Smith, B.H., Jellicoe, B., (2004). Rather, these kimberlite bodies are best described as polygenetic kimberlite tephra cones and tuff rings with associated feeder vents of variable geometry as shown by previous workers for the 169 kimberlite, the 140/141 kimberlite and the Star kimberlite. The domal tephra cone geometry is preserved due to burial by conformable Cretaceous marine mudstones and siltstones and is not an artifact of Quaternary glacial processes.     

Exploring LiDAR data for mapping the micro-topography and tidal hydro-dynamics of mangrove systems: An example from southeast Queensland,Australia

The aim was to explore the use of Light Detection and Ranging (LiDAR) data to map the micro-topography of an intertidal wetland in southeast Queensland Australia. The driver for this was the need to identify and map the habitats of the immature stages of an aedine disease vector mosquito (Aedes vigilax (Skuse)). We derived a high resolution digital elevation model (DEM) data set at a vertical resolution of 0.05 m from LiDAR data. The relative accuracy of the DEM across the site was tested by comparing water depth predictions derived from the DEM against in-situ water depth readings from pressure sensors over a 10-day tidal cycle, which included high spring tides. We found that the field observations of micro-topographic units important for mosquito management matched those delineated from the DEM. The micro-topography included a low berm or central ridge that was more or less continuous across the site, a shallow back basin and fringing mangroves. The fringing mangroves had unimpeded connection to the tidal source, however the central ridge blocked tidal water from the back basin for all but the highest tides. Eggshell survey indicated that the back basin was the area suitable for immature mosquitoes. We conclude that LiDAR data has application for understanding and mapping the structure of mangrove wetlands. We have also demonstrated (in a small area) that LiDAR is useful for modelling the effect of sea level changes on the coastal fringe. LiDAR may be the only method to inform research on changes to land use and ecosystems caused by sea level change.     

Effects of dredging and open-water disposal on benthic macroinvertebrates in a South Carolina estuary

Approximately 28,475 m³ of muddy sediments were dredged from a shoal in a South Carolina estuarine system and released near the surface at a nearby site having high tidal current velocities. Effects at the dredged sites included decreased macrofaunal abundance and changes in species composition. These effects appeared to be short term, with substantial recovery occurring within 3 months. Rapid recovery was primarily attributed to immigration through slumping of channel wall sediments similar to those dredged. Detrimental effects on benthic macrofauna in the area of open water disposal were minimal. Most differences noted in community structure between collection dates were attributed to sampling and seasonal variability. The absence of a major longterm disruption to the benthos in the disposal area was probably due to (1) strong tidal currents, which rapidly dispersed the moderate amount of mud sediments released; (2) surface disposal, permitting wider dispersal; and (3) disposal during late autumn, a period of low faunal recruitment.     

A 14-year monthly climatology and trend in the 35–65 km altitude range from Rayleigh Lidar temperature measurements at a low latitude station

A sodium resonance lidar at 589 nm has been operated in São José dos Campos, Brazil (23°S, 46°W) since 1972 mainly for studies related to the origin, chemistry and dynamics of the mesospheric sodium layer. Beginning in 1993, the improved laser capability has also enabled the processing of the Rayleigh signal from which the temperatures from 35 to 65 km are retrieved on a nightly mean basis. We used these nightly profiles to determine the monthly temperature profiles from 1993 to 2006. The mean temperature characteristics for each year and for the whole period are obtained. Seasonal thermal amplitude is small (6 K peak to peak at 40 and 60 km). Compared with the MSISE-90 model, a large difference is noted, with temperature lower than the model below the stratopause and higher above. Also the seasonal variation has a large difference with better agreement occurring around local winter, but with temperatures higher by 8–10 K at the equinoxes. The semiannual component is dominant over the annual at all altitudes. Linear trends with decreasing temperature of 1.09, 2.29 and 1.42 K/decade are observed at 40, 50 and 60 km, respectively.