Crustal and upper mantle seismic structure of the Australian Plate, South Island, New Zealand

Seismic reflection and refraction data were collected west of New Zealand's South Island parallel to the Pacific–Australian Plate boundary. The obliquely convergent plate boundary is marked at the surface by the Alpine Fault, which juxtaposes continental crust of each plate. The data are used to study the crustal and uppermost mantle structure and provide a link between other seismic transects which cross the plate boundary. Arrival times of wide-angle reflected and refracted events from 13 recording stations are used to construct a 380-km long crustal velocity model. The model shows that, beneath a 2–4-km thick sedimentary veneer, the crust consists of two layers. The upper layer velocities increase from 5.4–5.9 km/s at the top of the layer to 6.3 km/s at the base of the layer. The base of the layer is mainly about 20 km deep but deepens to 25 km at its southern end. The lower layer velocities range from 6.3 to 7.1 km/s, and are commonly around 6.5 km/s at the top of the layer and 6.7 km/s at the base. Beneath the lower layer, the model has velocities of 8.2–8.5 km/s, typical of mantle material. The Mohorovicic discontinuity (Moho) therefore lies at the base of the second layer. It is at a depth of around 30 km but shallows over the south–central third of the profile to about 26 km, possibly associated with a southwest dipping detachment fault. The high, variable sub-Moho velocities of 8.2 km/s to 8.5 km/s are inferred to result from strong upper mantle anisotropy. Multichannel seismic reflection data cover about 220 km of the southern part of the modelled section. Beneath the well-layered Oligocene to recent sedimentary section, the crustal section is broadly divided into two zones, which correspond to the two layers of the velocity model. The upper layer (down to about 7–9 s two-way travel time) has few reflections. The lower layer (down to about 11 s two-way time) contains many strong, subparallel reflections. The base of this reflective zone is the Moho. Bi-vergent dipping reflective zones within this lower crustal layer are interpreted as interwedging structures common in areas of crustal shortening. These structures and the strong northeast dipping reflections beneath the Moho towards the north end of the (MCS) line are interpreted to be caused by Paleozoic north-dipping subduction and terrane collision at the margin of Gondwana. Deeper mantle reflections with variable dip are observed on the wide-angle gathers. Travel-time modelling of these events by ray-tracing through the established velocity model indicates depths of 50–110 km for these events. They show little coherence in dip and may be caused side-swipe from the adjacent crustal root under the Southern Alps or from the upper mantle density anomalies inferred from teleseismic data under the crustal root.     

The Okavango giant mafic dyke swarm (NE Botswana): its structural significance within the Karoo Large Igneous Province

The structural organization of a giant mafic dyke swarm, the Okavango complex, in the northern Karoo Large Igneous Province (LIP) of NE Botswana is detailed. This N110°E-oriented dyke swarm extends for 1500 km with a maximum width of 100 km through Archaean basement terranes and Permo-Jurassic sedimentary sequences. The cornerstone of the study is the quantitative analysis of N>170 (exposed) and N>420 (detected by ground magnetics) dykes evidenced on a ca. 80-km-long section lying in crystalline host-rocks, at high-angle to the densest zone of the swarm (Shashe area). Individual dykes are generally sub-vertical and parallel to the entire swarm. Statistical analysis of width data indicates anomalous dyke frequency (few data <5.0 m) and mean dyke thickness (high value of 17 m) with respect to values classically obtained from other giant swarms. Variations of mean dyke thicknesses from 17 (N110°E swarm) to 27 m (adjoining and coeval N70°E giant swarm) are assigned to the conditions hosting fracture networks dilated as either shear or pure extensional structures, respectively, in response to an inferred NNW–SSE extension. Both fracture patterns are regarded as inherited brittle basement fabrics associated with a previous (Proterozoic) dyking event. The Okavango N110°E dyke swarm is thus a polyphase intrusive system in which total dilation caused by Karoo dykes (estimated frequency of 87%) is 12.2% (6315 m of cumulative dyke width) throughout the 52-km-long projected Shashe section. Assuming that Karoo mafic dyke swarms in NE Botswana follow inherited Proterozoic fractures, as similarly applied for most of the nearly synchronous giant dyke complexes converging towards the Nuanetsi area, leads us to consider that the resulting triple junction-like dyke/fracture pattern is not a definitive proof for a deep mantle plume in the Karoo LIP.     

Denitrification capacity in a subterranean estuary below a Rhode Island fringing salt marsh

Coastal waters are severely threatened by nitrogen (N) loading from direct groundwater discharge. The subterranean estuary, the mixing zone of fresh groundwater and sea water in a coastal aquifer, has a high potential to remove substantial N. A network of piezometers was used to characterize the denitrification capacity and groundwater flow paths in the subterranean estuary below a Rhode Island fringing salt marsh.¹⁵N-enriched nitrate was injected into the subterranean estuary (in situ push-pull method) to evaluate the denitrification capacity of the saturated zone at multiple depths (125–300 cm) below different zones (upland-marsh transition zone, high marsh, and low marsh). From the upland to low marsh, the water table became shallower, groundwater dissolved oxygen decreased, and groundwater pH, soil organic carbon, and total root biomass increased. As groundwater approached the high and low marsh, the hydraulic gradient increased and deep groundwater upwelled. In the warm season (groundwater temperature >12 °C), elevated groundwater denitrification capacity within each zone was observed. The warm season low marsh groundwater denitrification capacity was significantly higher than all other zones and depths. In the cool season (groundwater temperature <10.5 °C), elevated groundwater denitrification capacity was only found in the low marsh. Additions of dissolved organic carbon did not alter groundwater denitrification capacity suggesting that an alternative electron donor, possibly transported by tidal inundation from the root zone, may be limiting. Combining flow paths with denitrification capacity and saturated porewater residence time, we estimated that as much as 29–60 mg N could be removed from 11 of water flowing through the subterranean estuary below the low marsh, arguing for the significance of subterranean estuaries in annual watershed scale N budgets.     

Using stable isotope analyses to assess carbon dynamics in a shallow subtropical estuary

The sources of carbon, which fuel water column respiration, remain unresolved for most estuaries; our objective was to examine carbon dynamics in a shallow subtropical estuary. We sampled the Sabine-Neches estuary, Texas, during low (November 1999) and high (May 2000) freshwater inflow and measured stable carbon isotope ratios of the dissolved inorganic and orgnaic carbon (δ¹³C-DIC, δ¹³C-DOC), as well as quantifying accessory parameters (salinity, nutrients, total suspended solids, and photosynthetic pigments). Pigment analysis indicated that diatoms were the predominant phytoplankton. Data from the May 2000 sampling event exhibited conservative mixing, indicating that the system was acting as a conduit between the watershed and the Gulf of Mexico. During November, mixing was generally nonconservative indicating extensive recycling of allochthonous and autochthonous carbon sources. Our data imply that both carbon sources had similar isotope, ratios that made it impossible to unambiguously determine the dominant source supporting respiration. The nonconservative DIC concentration data indicating an autotrophic sink as well as the strong relationship between δ¹³C-DOC and chlorophylla, suggest that in situ production was an important component of the DOC pool. We hypothesize that uncharacteristically calm wind conditions during sampling may have promoted phytoplankton settling, removing autotrophs, from the water column, but leaving behind a dissolved biogeochemical signature. Interpretation of carbon dynamics may be confounded by spatial and temporal decoupling of producers and consumers from biogeochemical indicators.     

Continental rift systems and anorogenic magmatism

Precambrian Laurentia and Mesozoic Gondwana both rifted along geometric patterns that closely approximate truncated-icosahedral tessellations of the lithosphere. These large-scale, quasi-hexagonal rift patterns manifest a least-work configuration. For both Laurentia and Gondwana, continental rifting coincided with drift stagnation, and may have been driven by lithospheric extension above an insulated and thermally expanded mantle. Anorogenic magmatism, including flood basalts, dike swarms, anorthosite massifs and granite-rhyolite provinces, originated along the Laurentian and Gondwanan rift tessellations. Long-lived volcanic regions of the Atlantic and Indian Oceans, sometimes called hotspots, originated near triple junctions of the Gondwanan tessellation as the supercontinent broke apart. We suggest that some anorogenic magmatism results from decompression melting of asthenosphere beneath opening fractures, rather than from random impingement of hypothetical deep-mantle plumes.     

Mercury in Eastern Kentucky coals: Geologic aspects and possible reduction strategies

Mercury emissions from US coal-fired power plants will be regulated by the US Environmental Protection Agency (USEPA) before the end of the decade. Because of this, the control of Hg in coal is important. Control is fundamentally based on the knowledge of the amounts of Hg in mined, beneficiated, and as-fired coal. Eastern Kentucky coals, on a reserve district level, have Hg contents similar to the USA average for coal at mines. Individual coals show greater variation at the bench scale, with Hg enrichment common in the top bench, often associated with enhanced levels of pyritic sulfur. Some of the variation between parts of eastern Kentucky is also based on the position relative to major faults. The Pine Mountain thrust fault appears to be responsible for elemental enrichment, including Hg, in coals on the footwall side of the thrust.Eastern Kentucky coals shipped to power plants in 1999, the year the USEPA requested coal quality information on coal deliveries, indicate that coals shipped from the region have 0.09 ppm Hg, compared to 0.10 ppm for all delivered coals in the USA. On an equal energy basis, and given equal concentrations of Hg, the high volatile bituminous coals from eastern Kentucky would emit less Hg than lower rank coals from other USA regions.     

Representing and communicating deep uncertainty in climate-change assessments

IPCC reports provide a synthesis of the state of the science in order to inform the international policy process. This task is made difficult by the presence of deep uncertainty in the climate problem that results from long time scales and complexity. This paper focuses on how deep uncertainty can be effectively communicated. We argue that existing schemes do an inadequate job of communicating deep uncertainty and propose a simple approach that distinguishes between various levels of subjective understanding in a systematic manner. We illustrate our approach with two examples. To cite this article: M. Kandlikar et al., C. R. Geoscience 337 (2005).     

Groundwater control of mangrove surface elevation: Shrink and swell varies with soil depth

We measured monthly soil surface elevation change and determined its relationship to groundwater changes at a mangrove forest site along Shark River, Everglades National Park, Florida. We combined the use of an original design, surface elevation table with new rod-surface elevation tables to separately track changes in the mid zone (0–4 m), the shallow root zone (0–0.35 m), and the full sediment profile (0–6 m) in response to site hydrology (daily river stage and daily groundwater piezometric pressure). We calculated expansion and contraction for each of the four constituent soil zones (surface [accretion and erosion; above 0 m], shallow zone [0–0.35 m], middle zone [0.35–4 m], and bottom zone [4–6]) that comprise the entire soil column. Changes in groundwater pressure correlated strongly, with changes in soil elevation for the entire profile (Adjusted R² = 0.90); this relationship was not proportional to the depth of the soil profile sampled. The change in thickness of the bottom soil zone accounted for the majority (R² = 0.63) of the entire soil profile expansion and contraction. The influence of hydrology on specific soil zones and absolute elevation change must be considered when evaluating the effect of disturbances, sea level rise, and water management decisions on coastal wetland systems.     

The influence of edge sites on the development of surface charge on goethite nanoparticles: A molecular dynamics investigation

Large-scale molecular simulation of proton accumulations were carried out on (i) (110) and (021) slabs immersed in aqueous solution and (ii) a series of model goethite nanoparticles of dimension 2 to 8 nm with systematically varying acicularity and (110)/(021) surface areas. In the slab systems, the (021) surface exhibits 15% more proton charge per unit area than the (110) surface. In the particulate systems, the acicular particles having the highest (110)/(021) ratio accumulate the most charge, opposite to the trend expected from the slab simulations, indicating that, at length scales on the order of 10 nm, the slab results are not a good indicator of the overall charging behavior of the particles. The primary reason for the discrepancy between the particulate systems and slab systems is the preferential accumulation of protons at acute (110)-(110) intersections. Charge accumulates preferentially in this region because excess proton charge at an asperity is more effectively solvated than at a flat interface.     

Environmental significance of C/C and O/O ratios of modern land-snail shells from the southern great plains of North America

¹³C/¹²C and ¹⁸O/¹⁶O ratios of aragonite shells of modern land snails from the southern Great Plains of North America were measured for samples from twelve localities in a narrow east-west corridor that extended from the Flint Hills in North Central Oklahoma to the foothills of the Sangre de Cristo Mountains in Northern New Mexico, USA. Across the study area, shell δ¹⁸O values (PDB scale) ranged from −4.1‰ to 1.2‰, while δ¹³C values ranged from −13.2‰ to 0.0‰. δ¹⁸O values of the shell aragonite were predicted with a published, steady state, evaporative flux balance model. The predicted values differed (with one exception) by less than 1‰ from locality averages of measured δ¹⁸O values. This similarity suggests that relative humidity at the time of snail activity is an important control on the δ¹⁸O values of the aragonite and emphasizes the seasonal nature of the climatic information preserved in the shells. Correlated δ¹³C values of coexisting Vallonia and Gastrocopta suggest similar feeding habits and imply that these genera can provide information on variations in southern Great Plains plant ecology. Although there is considerable scatter, multispecies, transect average δ¹³C values of the modern aragonite shells are related to variations in the type of photosynthesis (i.e., C₃, C₄) in the local plant communities. The results of this study emphasize the desirability of obtaining isotope ratios representing averages of many shells in a locale to reduce possible biases associated with local variations among individuals, species, etc., and thus better represent the “neighborhood” scale temporal and/or spatial environmental variations of interest in studies of modern and ancient systems.