Field Survey of the March 28, 2005 Nias-Simeulue Earthquake and Tsunami

On the evening of March 28, 2005 at 11:09?p.m. local time (16:09 UTC), a large earthquake occurred offshore of West Sumatra, Indonesia. With a moment magnitude (M _w) of 8.6, the event caused substantial shaking damage and land level changes between Simeulue Island in the north and the Batu Islands in the south. The earthquake also generated a tsunami, which was observed throughout the source region as well as on distant tide gauges. While the tsunami was not as extreme as the tsunami of December 26th, 2004, it did cause significant flooding and damage at some locations. The spatial and temporal proximity of the two events led to a unique set of observational data from the earthquake and tsunami as well as insights relevant to tsunami hazard planning and education efforts.     

Comparison of theory and experiment for solute transport in weakly heterogeneous bimodal porous media

In this work, the influence of non-equilibrium effects on solute transport in a weakly heterogeneous medium is discussed. Three macro-scale models (upscaled via the volume averaging technique) are investigated: (i) the two-equation non-equilibrium model, (ii) the one-equation asymptotic model and (iii) the one-equation local equilibrium model. The relevance of each of these models to the experimental system conditions (duration of the pulse injection, dispersivity values…) is analyzed. The numerical results predicted by these macroscale models are compared directly with the experimental data (breakthrough curves). Our results suggest that the preasymptotic zone (for which a non-Fickian model is required) increases as the solute input pulse time decreases. Beyond this limit, the asymptotic regime is recovered. A comparison with the results issued from the stochastic theory for this regime is performed. Results predicted by both approaches (volume averaging method and stochastic analysis) are found to be consistent.     

The North Devon Basin: a Devonian passive margin shelf succession

The North Devon Basin, situated in a more proximal passive margin regime than the rift basins to the south, is not constrained but its succession is thought to represent in large part the sediments debouched from a northerly hinterland. Rather than that immediate source being South Wales an original location of the basin well to the south-east and west of the Ardennes massif is considered probable, with its present position being attained by Carboniferous displacement along the Bristol Channel-Bray Fault. The basin's thick (6000 m) succession comprises terrestrial and marine deposits that form two major sedimentary cycles, which are apparently closely linked to rift basin formation to the south. The GCR sites span a relatively straightforward shelf succession that extends from the late Early Devonian to the Carboniferous. The sedimentology, palaeontology, and depositional environments of terrestrial and marine facies lithostratigraphical units are detailed, some sites providing the macrofossil assemblages important in the identification and definition by Sedgwick and Murchison of the Devonian System.     

Petrogenesis of the anorthosite-chromitite association: crystal-chemical and petrological insights from the Rum Layered Suite, NW Scotland

Thirteen Cr-bearing spinels from major horizons of magma replenishment in the open-system Rum Layered Suite have been analysed by X-ray single crystal diffraction and electron microprobe analyses. On the basis of the structural parameters and the chemistry of these spinels the so-called Rum trend, in which Al-content increases at the expense of Cr and Fe³⁺, has been easily recognised. In addition, natural spinels with Fe³⁺ content similar to synthetic spinels on the MgCr₂O₄?CMgFe₂O₄ join have been analysed for the first time. Layers of chromitite, anorthosite and peridotite situated within several cm of one another have yielded different intracrystalline exchange temperatures using an intercrystalline spinel-olivine thermometer. The Rum anorthosite Cr-spinels are interpreted as having crystallised within the cumulus pile following rejuvenation of the crystal mush. Their low Al-content is a function of simultaneous plagioclase crystallisation, reducing the amount of Al³⁺ present for the Cr-spinel. By contrast, Cr-spinels in well-known Archean anorthosites (e.g. Ujaragssuit nunat and Fisken?sset, western Greenland) and Sittampundi (southern India) are very aluminous in composition, attributed to crystallisation of Cr-spinel from high-alumina basalts in lower crustal magma chambers and linked to the control exerted by plagioclase crystallisation on Al content of the melt, in the absence of clinopyroxene crystallisation. The compositional differences between the Rum anorthosite Cr-spinels and the Fisken?sset and Sittampundi Cr-spinels suggest that postcumulus reaction of Cr-spinel and melt to low (800?C900°C) temperatures, as invoked for the Rum crystals, may not have been as important a process in the Archean anorthosites.     

Seasonal variability in the input of lead, barium and indium to Law Dome, Antarctica

Lead (Pb) isotopic compositions and concentrations, and barium (Ba) and indium (In) concentrations have been determined at monthly resolution in five Law Dome (coastal Eastern Antarctica) ice core sections dated from ∼1757 AD to ∼1898 AD. ‘Natural’ background Pb concentrations in ∼1757 AD average ∼0.2 pg g⁻¹ and can be attributed to mineral dust and volcanic emissions, with ²⁰⁶Pb/²⁰⁷Pb ratios reaching up to 1.266 ± 0.002. From ∼1887 AD to ∼1898 AD, Pb concentrations reached ∼5 pg g⁻¹ and ²⁰⁶Pb/²⁰⁷Pb ratios decreased to 1.058 ± 0.001 as a result of additional inputs of Pb from anthropogenic sources. Seasonal variability in the late 1880s has been investigated by decoupling volcanic Pb from the total measured Pb concentrations, revealing spring and autumn maxima, and consistent winter minima, in anthropogenic Pb and mineral dust (Ba) concentrations. We link this variability to the annual cycle in the position and strength of the Antarctic Circumpolar Trough and, the Southern Ocean westerly winds to the north of the trough region. During the autumn and spring seasons, these systems increase in strength, transporting more impurity laden air from the Southern Hemisphere continental regions to Eastern Antarctica and Law Dome. As this Pb is isotopically identical to that emitted from south-eastern Australia (Broken Hill, Port Pirie) this implies a relatively direct air trajectory pathway from southern Australia to Law Dome (Eastern Antarctica).     

Highly depleted oceanic lithosphere in the Rheic Ocean: Implications for Paleozoic plate reconstructions

The Rheic Ocean formed at ca. 500 Ma, when several peri-Gondwanan terranes (e.g. Avalonia and Carolinia) drifted from the northern margin of Gondwana, and were consumed during the Late Carboniferous collision between Laurussia and Gondwana, a key event in the formation of Pangea. Several mafic complexes ranging in age from ca. 400–330 Ma preserve many of the lithotectonic and/or chemical characteristics of ophiolites. They are characterized by anomalously high εNd values that are typically either between or above the widely accepted model depleted mantle curves. These data indicate derivation from a highly depleted (HD) mantle and imply that (i) the mantle source of these complexes displays time-integrated depletion in Nd relative to Sm, and (ii) depletion is the result of an earlier melting event in the mantle from which basalt was extracted. The extent of mantle depletion indicates that this melting event occurred in the Neoproterozoic, possibly up to 500 million years before the Rheic Ocean formed. If so, the mantle lithosphere that gave rise to the Rheic Ocean mafic complexes must have been captured from an adjacent, older oceanic tract. The transfer of this captured lithosphere to the upper plate enabled it to become preferentially preserved. Possible Mesozoic–Cenozoic analogues include the capture of the Caribbean plate or the Scotia plate from the Pacific to the Atlantic oceanic realm. Our model implies that virtually all of the oceanic lithosphere generated during the opening phase of the Rheic Ocean was consumed by subduction during Laurentia–Gondwana convergence.     

Effect of sea-level lowering on ELA depression during the LGM

Decreases in equilibrium-line altitudes (ELAs) varied geographically during the last glacial maximum (LGM), with a mid-range value of ~ 900 m commonly deduced from altitude ratio and accumulation–area ratio calculations. Sea level, however, was 120 m lower during the LGM, so the ELA lowering relative to sea level would only be 780 m for a 900-m absolute lowering. With a lapse rate of 0.006°C m⁻¹, this implies a 4.7°C lowering of global temperature. It has been argued that this correction for sea-level change is unnecessary, but the logic on which this is based requires adiabatic compression to apply over much longer time scales than is typically invoked. We find that the correction is necessary. In addition, geometric changes in the atmosphere during the LGM, pointed out by Osmaston (2006), could lead to 0.4°C decrease in the average temperature of the troposphere. Additionally, orographic effects could significantly change the snow distribution on mountain masses near sea level.     

Hydrogeology,Chemical and Microbial Activity Measurement Through Deep Permafrost

Little is known about hydrogeochemical conditions beneath thick permafrost, particularly in fractured crystalline rock, due to difficulty in accessing this environment. The purpose of this investigation was to develop methods to obtain physical, chemical, and microbial information about the subpermafrost environment from a surface‐drilled borehole. Using a U‐tube, gas and water samples were collected, along with temperature, pressure, and hydraulic conductivity measurements, 420 m below ground surface, within a 535 m long, angled borehole at High Lake, Nunavut, Canada, in an area with 460‐m‐thick permafrost. Piezometric head was well above the base of the permafrost, near land surface. Initial water samples were contaminated with drill fluid, with later samples <40% drill fluid. The salinity of the non‐drill fluid component was <20,000 mg/L, had a Ca/Na ratio above 1, with δ¹⁸O values ～5‰ lower than the local surface water. The fluid isotopic composition was affected by the permafrost‐formation process. Nonbacteriogenic CH₄ was present and the sample location was within methane hydrate stability field. Sampling lines froze before uncontaminated samples from the subpermafrost environment could be obtained, yet the available time to obtain water samples was extended compared to previous studies. Temperature measurements collected from a distributed temperature sensor indicated that this issue can be overcome easily in the future. The lack of methanogenic CH₄ is consistent with the high sulfate concentrations observed in cores. The combined surface‐drilled borehole/U‐tube approach can provide a large amount of physical, chemical, and microbial data from the subpermafrost environment with few, controllable, sources of contamination.     

Hydrogeophysical Methods for Analyzing Aquifer Storage and Recovery Systems

Hydrogeophysical methods are presented that support the siting and monitoring of aquifer storage and recovery (ASR) systems. These methods are presented as numerical simulations in the context of a proposed ASR experiment in Kuwait, although the techniques are applicable to numerous ASR projects. Bulk geophysical properties are calculated directly from ASR flow and solute transport simulations using standard petrophysical relationships and are used to simulate the dynamic geophysical response to ASR. This strategy provides a quantitative framework for determining site‐specific geophysical methods and data acquisition geometries that can provide the most useful information about the ASR implementation. An axisymmetric, coupled fluid flow and solute transport model simulates injection, storage, and withdrawal of fresh water (salinity ～500 ppm) into the Dammam aquifer, a tertiary carbonate formation with native salinity approximately 6000 ppm. Sensitivity of the flow simulations to the correlation length of aquifer heterogeneity, aquifer dispersivity, and hydraulic permeability of the confining layer are investigated. The geophysical response using electrical resistivity, time‐domain electromagnetic (TEM), and seismic methods is computed at regular intervals during the ASR simulation to investigate the sensitivity of these different techniques to changes in subsurface properties. For the electrical and electromagnetic methods, fluid electric conductivity is derived from the modeled salinity and is combined with an assumed porosity model to compute a bulk electrical resistivity structure. The seismic response is computed from the porosity model and changes in effective stress due to fluid pressure variations during injection/recovery, while changes in fluid properties are introduced through Gassmann fluid substitution.