Shallow seismicity, triggered seismicity, and ambient noise tomography at the long-dormant Uturuncu Volcano, Bolivia

Using a network of 15 seismometers around the inflating Uturuncu Volcano from April 2009 to 2010, we find an average rate of about three local volcano-tectonic earthquakes per day, and swarms of 5–60 events a few times per month with local magnitudes ranging from −1.2 to 3.7. The earthquake depths are near sea level, more than 10 km above the geodetically inferred inflation source and the Altiplano Puna Magma Body. The Mw 8.8 Maule earthquake on 27 February 2010 triggered hundreds of earthquakes at Uturuncu with the onset of the Love and Rayleigh waves and again with the passage of the X2/X3 overtone phases of Rayleigh waves. This is one of the first incidences in which triggering has been observed from multiple surface wave trains. The earthquakes are oriented NW–SE similar to the regional faults and lineaments. The b value of the catalog is 0.49, consistent with a tectonic origin of the earthquakes. We perform ambient noise tomography using Love wave cross-correlations to image a low-velocity zone at 1.9 to 3.9 km depth below the surface centered slightly north of the summit. The low velocities are perhaps related to the hydrothermal system and the low-velocity zone is spatially correlated with earthquake locations. The earthquake rate appears to vary with time—a seismic deployment from 1996 to 1997 reveals 1–5 earthquakes per day, whereas 60 events/day were seen during 5 days using one seismometer in 2003. However, differences in analysis methods and magnitudes of completeness do not allow direct comparison of these seismicity rates. The rate of seismic activity at Uturuncu is higher than at other well-monitored inflating volcanoes during periods of repose. The frequent swarms and triggered earthquakes suggest the hydrothermal system is metastable.     

Evolution of Mid-Atlantic Coastal and Back-Barrier Estuary Environments in Response to a Hurricane: Implications for Barrier-Estuary Connectivity

Assessments of coupled barrier island-estuary storm response are rare. Hurricane Sandy made landfall during an investigation in Barnegat Bay-Little Egg Harbor estuary that included water quality monitoring, geomorphologic characterization, and numerical modeling; this provided an opportunity to characterize the storm response of the barrier island-estuary system. Barrier island morphologic response was characterized by significant changes in shoreline position, dune elevation, and beach volume; morphologic changes within the estuary were less dramatic with a net gain of only 200,000 m³ of sediment. When observed, estuarine deposition was adjacent to the back-barrier shoreline or collocated with maximum estuary depths. Estuarine sedimentologic changes correlated well with bed shear stresses derived from numerically simulated storm conditions, suggesting that change is linked to winnowing from elevated storm-related wave-current interactions rather than deposition. Rapid storm-related changes in estuarine water level, turbidity, and salinity were coincident with minima in island and estuarine widths, which may have influenced the location of two barrier island breaches. Barrier-estuary connectivity, or the transport of sediment from barrier island to estuary, was influenced by barrier island land use and width. Coupled assessments like this one provide critical information about storm-related coastal and estuarine sediment transport that may not be evident from investigations that consider only one component of the coastal system.     

Effects of instream processes,discharge, and land cover on nitrogen export from southern Appalachian Mountain catchments

Catchments with minimal disturbance usually have low dissolved inorganic nitrogen (DIN) export, but disturbances and anthropogenic inputs result in elevated DIN concentration and export and eutrophication of downstream ecosystems. We studied streams in the southern Appalachian Mountains, USA, an area dominated by hardwood deciduous forest but with areas of valley agriculture and increasing residential development. We collected weekly grab samples and storm samples from nine small catchments and three river sites. Most discharge occurred at baseflow, with baseflow indices ranging from 69% to 95%. We identified three seasonal patterns of baseflow DIN concentration. Streams in mostly forested catchments had low DIN with bimodal peaks, and summer peaks were greater than winter peaks. Streams with more agriculture and development also had bimodal peaks; however, winter peaks were the highest. In streams draining catchments with more residential development, DIN concentration had a single peak, greatest in winter and lowest in summer. Three methods for estimating DIN export produced consistent results. Annual DIN export ranged from less than 200 g ha^?1 year^?1 for the less disturbed catchments to over 2,000 g ha^?1 year^?1 in the catchments with the least forest area. Land cover was a strong predictor of DIN concentration but less significant for predicting DIN export. The two forested reference catchments appeared supply limited, the most residential catchment appeared transport limited, and export for the other catchments was significantly related to discharge. In all streams, baseflow DIN export exceeded stormflow export. Morphological and climatological variation among watersheds created complexities unexplainable by land cover. Nevertheless, regression models developed using land cover data from the small catchments reasonably predicted concentration and export for receiving rivers. Our results illustrate the complexity of mechanisms involved in DIN export in a region with a mosaic of climate, geology, topography, soils, vegetation, and past and present land use.     

The CQU Regional Seismic Network and applications to underground mining in Central Queensland,Australia

Digital seismograms continuously recorded from 1988 to 1992 by two stations of the RESNOM seismic network in northern Baja California, Mexico, were used to search for probable shear-wave anisotropic characteristics in the region of the Cerro Prieto fault. Shear-wave splitting was identified in many of the three-component records analyzed. We measured the polarization direction of the leadingS wave inside theS-wave window as well as the delay times between fast and slow phases on those records displaying shear-wave splitting. For station CPX, which is nearest the Imperial Valley region to the north, the preferred polarization direction found in this study (azimuth 180°±10°) coincides with the direction of the regional maximum compressive stress determined for the region. This polarization direction can be interpreted in terms of the “Extensive Dilatancy Anisotropy” model as the effect of vertical parallel aligned cracks. The preferred polarization direction measured at LMX, however, gives an azimuth of 45°±5°. Thus, it appears that faults and fractures aligned oblique to the main tectonic trend have a greater influence on the anisotropic characteristics of the crust south of Cerro Prieto volcano than that of the regional stress field. Time delays between slow and fastS waves observed at CPX appear constant from 1988 to 1992 while delays measured at LMX for the same interval indicate a small increase with time which cannot be attributed to azimuthal variations of paths.     

Geostatistical Mapping of Salinity Conditioned on Borehole Logs,Montebello Oil Field,California

We present a geostatistics-based stochastic salinity estimation framework for the Montebello Oil Field that capitalizes on available total dissolved solids (TDS) data from groundwater samples as well as electrical resistivity (ER) data from borehole logging. Data from TDS samples (n = 4924) was coded into an indicator framework based on falling below four selected thresholds (500, 1000, 3000, and 10,000 mg/L). Collocated TDS-ER data from the surrounding groundwater basin were then employed to produce a kernel density estimator to establish conditional probabilities for ER data (n = 8 boreholes) falling below the selected TDS thresholds within the Montebello Oil Field area. Directional variograms were estimated from these indicator coded data, and 500 TDS realizations from conditional indicator simulation were generated for the subsurface region above the Montebello Oil Field reservoir. Simulations were summarized as 3D maps of median TDS, most likely salinity class, and probability for exceeding each of the specified TDS thresholds. Results suggested TDS was below 500 mg/L in most of the study area, with a trend toward higher values (500 to 1000 mg/L) to the southwest; consistent with the average regional groundwater flow direction. Discrete localized zones of TDS greater than 1000 mg/L were observed, with one of these zones in the greater than 10,000 mg/L range; however, these areas were not prevalent. The probabilistic approach used here is adaptable and is readily modified to include additional data and types and can be employed in time-lapse salinity modeling through Bayesian updating.     

Simulating the diurnal cycle of rainfall in global climate models: resolution versus parameterization

The effects of horizontal resolution and the treatment of convection on simulation of the diurnal cycle of precipitation during boreal summer are analyzed in several innovative weather and climate model integrations. The simulations include: season-long integrations of the Non-hydrostatic Icosahedral Atmospheric Model (NICAM) with explicit clouds and convection; year-long integrations of the operational Integrated Forecast System (IFS) from the European Centre for Medium-range Weather Forecasts at three resolutions (125, 39 and 16 km); seasonal simulations of the same model at 10 km resolution; and seasonal simulations of the National Center for Atmospheric Research (NCAR) low-resolution climate model with and without an embedded two-dimensional cloud-resolving model in each grid box. NICAM with explicit convection simulates best the phase of the diurnal cycle, as well as many regional features such as rainfall triggered by advancing sea breezes or high topography. However, NICAM greatly overestimates mean rainfall and the magnitude of the diurnal cycle. Introduction of an embedded cloud model within the NCAR model significantly improves global statistics of the seasonal mean and diurnal cycle of rainfall, as well as many regional features. However, errors often remain larger than for the other higher-resolution models. Increasing resolution alone has little impact on the timing of daily rainfall in IFS with parameterized convection, yet the amplitude of the diurnal cycle does improve along with the representation of mean rainfall. Variations during the day in atmospheric prognostic fields appear quite similar among models, suggesting that the distinctive treatments of model physics account for the differences in representing the diurnal cycle of precipitation.     

Identification of genes expressed in juvenile Haliotis rufescens in response to different copper concentrations in the north of Chile under controlled conditions

This study reports molecular markers potentially associated with resistance or sensitivity to the impact of copper in juvenile red abalone, Haliotis rufescens, in the north of Chile under experimental conditions. Genomic analysis was made applying subtractive hybridization libraries (SSH) to identify genes up-and down regulated during cooper exposure in abalone over periods of 12 and 168 h exposed to 2.5 and 10 μg/L of Cu⁺². Results obtained from the SSH library revealed 368 different sequences regulated by copper, that correspond to eight major physiological functions. The validation of these sequences obtained by SSH as well as their expression kinetics were made by PCR in real time on 14 potential genes regulated by metal stress. This study provides information for the characterization of potential genomic markers that may be used in future environmental monitoring and to investigate new mechanisms of stress to copper in this commercially important marine species.