The distribution and stable carbon isotopic composition of dissolved organic carbon in estuaries

The distributions of dissolved organic carbon (DOC) and the natural carbon isotope ratio of DOC (DO¹³C) in estuaries reflect the predominant sources and sinks of organic matter from both allochthonous and autochthonous origins. The traditional view is that DOC in land-margin ecosystems reflects mainly the mixing of land-derived and oceanic DOC. However, this view is not consistent with the bulk of our data from a survey of DOC and DO¹³C distributions in estuaries on the East and Gulf coasts of the USA. While it is accurate that the DOC in estuaries includes material derived from land and from the ocean, the distributions of DOC and DO¹³C in several estuaries reflect additional DOC inputs from estuarine phytoplankton and tidal marshes. Even when DOC concentrations were distributed conservatively, the isotopic composition of the DOC revealed the existence of a dynamic cycle of DOC input and removal in some systems.     

Hurricanes 2004: An overview of their characteristics and coastal change

Four hurricanes battered the state of Florida during 2004, the most affecting any state since Texas endured four in 1884. Each of the storms changed the coast differently. Average shoreline change within the right front quadrant of hurricane force winds varied from 1 m of shoreline advance to 20 m of retreat, whereas average sand volume change varied from 11 to 66 m³ m⁻¹ of net loss (erosion). These changes did not scale simply with hurricane intensity as described by the Saffir-Simpson Hurricane Scale. The strongest storm of the season, category 4 Hurricane Charley, had the least shoreline retreat. This was likely because of other factors like the storm's rapid forward speed and small size that generated a lower storm surge than expected. Two of the storms, Hurricanes Frances and Jeanne, affected nearly the same area on the Florida east coast just 3 wk apart. The first storm, Frances, although weaker than the second, caused greater shoreline retreat and sand volume erosion. As a consequence, Hurricane Frances may have stripped away protective beach and exposed dunes to direct wave attack during Jeanne, although there was significant dune erosion during both storms. The maximum shoreline change for all four hurricanes occurred during Ivan on the coasts of eastern Alabama and the Florida Panhandle. The net volume change across a barrier island within the Ivan impact zone approached zero because of massive overwash that approximately balanced erosion of the beach. These data from the 2004 hurricane season will prove useful in developing new ways to scale and predict coastal-change effects during hurricanes.     

Groundwater-supported evapotranspiration within glaciated watersheds under conditions of climate change

This paper analyzes the effects of geology and geomorphology on surface-water/-groundwater interactions, evapotranspiration, and recharge under conditions of long-term climatic change. Our analysis uses hydrologic data from the glaciated Crow Wing watershed in central Minnesota, USA, combined with a hydrologic model of transient coupled unsaturated/saturated flow (HYDRAT2D). Analysis of historical water-table (1970–1993) and lake-level (1924–2002) records indicates that larger amplitude and longer period fluctuations occur within the upland portions of watersheds due to the response of the aquifer system to relatively short-term climatic fluctuations. Under drought conditions, lake and water-table levels fell by as much as 2–4 m in the uplands but by 1 m in the lowlands. The same pattern can be seen on millennial time scales. Analysis of Holocene lake-core records indicates that Moody Lake, located near the outlet of the Crow Wing watershed, fell by as much as 4 m between about 4400 and 7000 yr BP. During the same time, water levels in Lake Mina, located near the upland watershed divide, fell by about 15 m. Reconstructed Holocene climate as represented by HYDRAT2D gives somewhat larger drops (6 and 24 m for Moody Lake and Lake Mina, respectively). The discrepancy is probably due to the effect of three-dimensional flow. A sensitivity analysis was also carried out to study how aquifer hydraulic conductivity and land-surface topography can influence water-table fluctuations, wetlands formation, and evapotranspiration. The models were run by recycling a wet year (1985, 87 cm annual precipitation) over a 10-year period followed by 20 years of drier and warmer climate (1976, 38 cm precipitation). Model results indicated that groundwater-supported evapotranspiration accounted for as much as 12% (10 cm) of evapotranspiration. The aquifers of highest hydraulic conductivity had the least amount of groundwater-supported evapotranspiration owing to a deep water table. Recharge was even more sensitive to aquifer hydraulic conductivity, especially in the lowland regions. These findings have important implications for paleoclimatic studies, because the hydrologic response of a surface-water body will vary across the watershed to a given climate signal.     

Fault slip in the 1997 Manyi, Tibet earthquake from linear elastic modelling of InSAR displacements

The M _w 7.6 1997 Manyi earthquake occurred in an area of central northern Tibet where sparse vegetation coverage and a lack of human habitation provide excellent conditions for Interferometric Synthetic Aperture Radar (InSAR) studies. We use coseismic pairs of radar images acquired by the ESA ERS-2 satellite to construct interferograms of the surface displacement field due to the earthquake. The location and extent of the coseismic fault rupture are mapped using a combination of optical satellite imagery, high-resolution digital topography, interferometric correlation and azimuth offset measurements; in so doing, we are able to relate prominent geomorphic features in the fault zone to bends in the fault.
Using elastic dislocation models consistent with this mapped fault trace, we then test a range of fault geometries and slip conditions to find the combination which best explains the InSAR displacements. Our favoured model contains a reversal in fault dip, approximately halfway along its length, occurring at the location of a restraining bend. Slip on this model fault is heterogeneous, with two areas of peak slip of 7 m or greater, and components of dip-slip displacement which vary significantly along-strike. The success of this model in fitting the data implies that an observed asymmetry in the coseismic interferograms can be explained in terms of the local fault geometry, rather than by using non-linear elastic rheologies as suggested by earlier authors.     

Breadcrust bombs as indicators of Vulcanian eruption dynamics at Guagua Pichincha volcano,Ecuador

Vulcanian eruptions are common at many volcanoes around the world. Vulcanian activity occurs as either isolated sequences of eruptions or as precursors to sustained explosive events and is interpreted as clearing of shallow plugs from volcanic conduits. Breadcrust bombs characteristic of Vulcanian eruptions represent samples of different parts of these plugs and preserve information that can be used to infer parameters of pre-eruption magma ascent. The morphology and preserved volatile contents of breadcrust bombs erupted in 1999 from Guagua Pichincha volcano, Ecuador, thus allow us to constrain the physical processes responsible for Vulcanian eruption sequences of this volcano. Morphologically, breadcrust bombs differ in the thickness of glassy surface rinds and in the orientation and density of crack networks. Thick rinds fracture to create deep, widely spaced cracks that form large rectangular domains of surface crust. In contrast, thin rinds form polygonal networks of closely spaced shallow cracks. Rind thickness, in turn, is inversely correlated with matrix glass water content in the rind. Assuming that all rinds cooled at the same rate, this correlation suggests increasing bubble nucleation delay times with decreasing pre-fragmentation water content of the melt. A critical bubble nucleation threshold of 0.4–0.9 wt% water exists, below which bubble nucleation does not occur and resultant bombs are dense. At pre-fragmentation melt H₂O contents of >∼0.9 wt%, only glassy rinds are dense and bomb interiors vesiculate after fragmentation. For matrix glass H₂O contents of ≥1.4 wt%, rinds are thin and vesicular instead of thick and non-vesicular. A maximum measured H₂O content of 3.1 wt% establishes the maximum pressure (63 MPa) and depth (2.5 km) of magma that may have been tapped during a single eruptive event. More common H₂O contents of ≤1.5 wt% suggest that most eruptions involved evacuation of ≤1.5 km of the conduit. As we expect that substantial overpressures existed in the conduit prior to eruption, these depth estimates based on magmastatic pressure are maxima. Moreover, the presence of measurable CO₂ (≤17 ppm) in quenched glass of highly degassed magma is inconsistent with simple models of either open- or closed-system degassing, and leads us instead to suggest re-equilibration of the melt with gas derived from a deeper magmatic source. Together, these observations suggest a model for the repeated Vulcanian eruptions that includes (1) evacuation of the shallow conduit during an individual eruption, (2) depressurization of magma remaining in the conduit accompanied by open-system degassing through permeable bubble networks, (3) rapid conduit re-filling, and (4) dome formation prior to the subsequent explosion. An important part of this process is densification of upper conduit magma to allow repressurization between explosions. At a critical overpressure, trapped pressurized gas fragments the nascent impermeable cap to repeat the process.     

The identification of CTA 21

CTA 21 has long been the outstanding example of a strong, compact radio source lacking an optical identification. We report the discovery at infrared wavelengths of the counterpart of CTA 21, and show that its spectral shape is unusual.     

Arsenic incorporation into FeS2 pyrite and its influence on dissolution: A DFT study

FeS₂ pyrite can incorporate large amounts of arsenic (up to ca. 10 wt%) and hence has a strong impact on the mobility of this toxic metalloid. Focussing on the lowest arsenic concentrations for which the incorporation occurs in solid solution, the substitution mechanisms involved have been investigated by assuming simple incorporation reactions in both oxidising and reducing conditions. The solution energies were calculated by Density Functional Theory (DFT) calculations and we predict that the formation of AsS dianion groups is the most energetically favourable mechanism. The results also suggest that the presence of arsenic will accelerate the dissolution and thus the generation of acid drainage, when the crystal dissolves in oxidising conditions.