Late Cretaceous ophiolite obduction and Paleocene India-Asia collision in the westernmost Himalaya

Abstract

Collision of the Kohistan island arc with Asia at ~100 Ma resulted in N-S compression within the Neo-Tethys at a spreading center north of the Indo-Pakistani craton. Subsequent India-Asia convergence converted the Neo-Tethyan spreading center into a short-lived subduction zone. The hanging wall of the subduction zone became the Waziristan, Khost and Jalalabad igneous complexes. During the Santonian- Campanian (late Cretaceous), thrusting of the NW IndoPakistani craton beneath Albian oceanic crust and a Cenomanian volcano-sedimentary complex, generated an ophiolite-radiolarite belt. Ophiolite obduction resulted in tectonic loading and flexural subsidence of the NW Indian margin and sub-CCD deposition of shelf-derived olistostromes and turbidites in the foredeep. Campanian-Maastriehtian calci- clastic and siliciclastic sediment gravity flows derived from both margins filled the foredeep as a huge allochthon of Triassic-Jurassic rise and slope strata was thrust ahead of the ophiolites onto the Indo-Pakistani craton. Shallow to intermediate marine strata covered the foredeep during the late Maastrichtian. As ophiolite obduction neared completion during the Maastrichtian, the majority of India-Asia convergence was accommodated along the southern margin of Asia. During the Paleocene, India was thrust beneath a second allochthon that included open marine middle Maastrichtian colored mélange which represents the Asian Makran-Indus-Tsangpo accretionary prism. Latérites that formed on the eroded ophiolites and structurally higher colored mélange during the Paleocene wei’e unconformably overlapped by upper Paleocene and Middle Eocene shallow marine limestone and shale that delineate distinct episodes of Paleocene collisional and Early Eocene post-collisional deformation.     

Discretizing the fracture-matrix interface to simulate solute transport

This article examines the required spatial discretization perpendicular to the fracture-matrix interface (FMI) for numerical simulation of solute transport in discretely fractured porous media. The discrete-fracture, finite-element model HydroGeoSphere ( Therrien et al. 2005 ) and a discrete-fracture implementation of MT3DMS ( Zheng 1990 ) were used to model solute transport in a single fracture, and the results were compared to the analytical solution of Tang et al. (1981) . To match analytical results on the relatively short timescales simulated in this study, very fine grid spacing perpendicular to the FMI of the scale of the fracture aperture is necessary if advection and/or dispersion in the fracture is high compared to diffusion in the matrix. The requirement of such extremely fine spatial discretization has not been previously reported in the literature. In cases of high matrix diffusion, matching the analytical results is achieved with larger grid spacing at the FMI. Cases where matrix diffusion is lower can employ a larger grid multiplier moving away from the FMI. The very fine spatial discretization identified in this study for cases of low matrix diffusion may limit the applicability of numerical discrete-fracture models in such cases.     

Elevation adjustments of paired natural levees during flooding of the Saskatchewan River

Natural levees control the exchange of water between an alluvial channel and its floodplain, but little is known about the spatial distribution and evolution of levee heights. The summer 2005 flood of the Saskatchewan River (Cumberland Marshes, east‐central Saskatchewan) inundated large areas of floodplain for up to seven weeks, forming prominent new deposits on natural levees along main‐stem channels. Measurements of flood‐deposit thickness and crest heights of 61 levee pairs show that the thickest deposits occur on the lower pre‐flood levee in 80% of the sites, though no clear relationship exists between deposit thickness and magnitude of height difference. Only 16% of the pairs displayed thicker deposits on the higher levee, half of which occurred at sites where relatively clear floodbasin waters re‐entered turbid channels during general flooding. Difference in crest elevation (ΔE) between paired levees is approximately log‐normally distributed, both before and after the flood, though with different mean values. Supplemental observations from tank experiments indicate that during near‐bankfull flows, temporally and spatially variable deposition and erosion occur on levees due to backwater effects associated with nearby channel bars and irregular rises of the channel bed forced by channel extension. During floods, preferential deposition in lows tends to even out crest heights. Copyright © 2009 John Wiley & Sons, Ltd.     

A Survey of Techniques for Predicting Earthquake Ground Motions for Engineering Purposes

Over the past four or five decades many advances have been made in earthquake ground-motion prediction and a variety of procedures have been proposed. Some of these procedures are based on explicit physical models of the earthquake source, travel-path and recording site while others lack a strong physical basis and seek only to replicate observations. In addition, there are a number of hybrid methods that seek to combine benefits of different approaches. The various techniques proposed have their adherents and some of them are extensively used to estimate ground motions for engineering design purposes and in seismic hazard research. These methods all have their own advantages and limitations that are not often discussed by their proponents. The purposes of this article are to: summarise existing methods and the most important references, provide a family tree showing the connections between different methods and, most importantly, to discuss the advantages and disadvantages of each method.

Magnetic perturbations to the acoustic modes of roAp stars

We study the effect of a large-scale surface magnetic field on the non-radial acoustic modes of roAp stars. Special attention is given to the use of a variational principle which is used for determining the shifts in the frequencies with relative ease, enabling us to avoid having to calculate the perturbed eigenfunctions. With knowledge of the frequency shifts we then estimate the eigenfunctions in a simpler, albeit approximate way. The results indicate frequency shifts of the order of few μHz, which depend on the order, degree and azimuthal order of the mode. The loss of energy through Alfvén waves is also estimated from the imaginary parts of the frequency shifts. The results indicate that the loss is particularly high near specific frequencies. This might indicate the presence of a selection effect, which could make some modes more likely to be excited than others. However, our results do not explain why the modes observed appear always to be aligned with the axis of the magnetic field. Finally, the estimated perturbed eigenfunctions contain strong components of spherical harmonics that differ from those of the original unperturbed modes.     

Deformation and patterns of sedimentation,South San Clemente Basin,California Borderland

Seismic reflection profiling in the South San Clemente Basin and the southern portion of the San Diego Trough has revealed at least six sedimentary units exhibiting varying stages of deformation. Four of the units are interpreted to be marine turbidites supplied by adjacent submarine canyons. Sediments comprising the Descanso Plain and correlative material within the South San Clemente Basin are attributed to a southerly source (Banda Canyon), while the more recent Quaternary turbidites from Coronado Canyon filled the southern San Diego Trough and then spilled over into South San Clemente Basin. The relatively high but intermittent rates of sedimentation, together with shifting sources and areas of deposition, have resulted in sedimentary units that were emplaced in comparatively short episodes but which were subjected to relatively continuous tectonic activity. Consequently, the sedimentary layers of each unit appear uniformly affected by deformation which increases in successively older units.     

Trace element distribution in mineral separates of the Allende inclusions and their genetic implications

Concentrations of the REE, Sc, Co, Fe, Zn, Ir, Na and Cr were determined by instrumental neutron activation and mass spectrometric isotope dilution analysis for mineral separates of the coarseand fine-grained types (group I and II of Martin and Mason's classification) of the Allende inclusions.These data, combined with data on mineral/liquid partition coefficients, oxygen isotope distributions and diffusion calculations, suggest the following: (1) Minerals in the coarse-grained inclusions (group I) crystallized in a closed system with respect to refractory elements. On the other hand, differences in oxygen isotope distributions among minerals preclude a totally molten stage in the history of the inclusion. Group I inclusions were formed by rapid condensation (either to liquid or solid) in a supercooled solar nebula; extrasolar pyroxene and spinel dust were included but not melted in the condensing inclusions, thus preserving their extrasolar oxygen isotope composition. REE were distributed by diffusion during the subsequent heating at subsolidus temperatures; because oxygen diffuses much more slowly at these temperatures, the oxygen isotope anomalies were preserved. (2) The fine-grained (group II) inclusions were also formed by condensation from a super-cooled nebular gas; however, REE-rich clinopyroxene and spinel were formed early and REE-poor sodalite and nepheline were formed later and mechanically mixed with clinopyroxene and spinel to form the inclusions. The REE patterns of the bulk inclusions and the mineral separates are fractionated, indicating that REE abundances in the gaseous phase were already fractionated at the time of condensation of the minerals. (3) Pre-existing Mg isotope anomalies in the coarse-grained inclusions must have been erased during the heating stage thus resetting the ²⁶Al-²⁶Mg chronometer.     

The age,origin, and volcanological significance of the Y-5 ash layer in the Mediterranean

The Y-5 ash is the most widespread layer in deep-sea sediments from the eastern Mediterranean. This ash layer was previously correlated with the Citara-Serrara tuff on Ischia Island and dated at approximately 25,000 yr B.P. New data on the glass chemistry of the Y-5 ash and pyroclastic deposits from the Neopolitan volcanic province suggest that the layer is correlative with the large-volume Campanian ignimbrite and not with the deposit from Ischia Island. The volume of the Y-5 ash is approximately 65 km³ which is comparable in magnitude to the volume of the Campanian ignimbrite. An interpolated age of approximately 38,000 yr B.P. is estimated based on sedimentation rates derived from δ¹⁸O stratigraphy. There is a discrepancy between this estimate and previously reported radiocarbon ages which range from 24,000 to 35,000 yr B.P. We propose that the “Campanian tuff ash layer” should be adopted as the full stratigraphic name for the Y-5 ash. The deep-sea ash layer is divisible into two units in proximal localities, probably correlating with two major phases of the eruption: plinian and ignimbrite.     

Mg2+ removal in the system Mg2+—amorphous SiO2—H2O by adsorption and Mg-hydroxysilicate precipitation

Two chemical processes can remove Mg²⁺ from suspensions containing amorphous silica (am-SiO₂) at low temperatures: adsorption and precipitation of a Mg-hydroxysilicate resembling sepiolite. Mg²⁺ removal from am-SiO₂ suspensions was investigated, and the relative role of the two removal processes evaluated, as a function of: pH, ionic strength, Mg²⁺ concentration, and temperature.The extent of Mg²⁺ adsorption onto am-SiO₂ decreases with increasing NaCl concentration due to displacement of Mg²⁺ by Na⁺. At NaCl concentrations of 0.05 M and above, adsorption occurs only at pH values above 8.5, where rapid dissolution of am-SiO₂ gives rise to high concentrations of dissolved silica, resulting in supersaturation with respect to sepiolite. Removal of Mg²⁺, at concentrations of 40 to 650 μM, from am-SiO₂ suspensions in 0.70 M NaCl at 25 °C occurs at pH 9.0 and above. Experiments show that under these conditions adsorption and Mg-hydroxysilicate precipitation remove Mg²⁺ at similar rates. For 0.05 M Mg²⁺, at 0.70 M ionic strength and 25 °C, measurable Mg²⁺ removal occurs down to ca. pH 7.5 but is primarily due to Mg-hydroxysilicate precipitation. For the same solution conditions at 5°C, Mg²⁺ removal occurs above pH 8.0 and is primarily due to adsorption.Assuming that increasing pressure does not greatly enhance adsorption, Mg²⁺ adsorption onto am-SiO₂ is an insignificant process in sea water. The surface charge of pristine am-SiO₂ in sea water is primarily controlled by interactions with Na⁺. The principal reaction between Mg²⁺ and am-SiO₂ in marine sediments is sepiolite precipitation.The age distribution of sepiolite in siliceous pelagic sediments is influenced by temperatures of bottom waters and by geothermal gradients.