The origin of the Nicaraguan depression

A conspicuous graben extends for 800 kilometers through El Salvador and western Nicaragua to the Caribbean Sea in northeastern Costa Rica. Like the smaller but structurally similar Semangko and Toba Depressions of northern Sumatra, the trough is clearly related to voluminous volcanic eruptions during Late Tertiary time. In the region around Lakes Managua and Nicaragua, where the depression is best defined and reaches its greatest dimension, a thick series of Tertiary sediments and volcanic rocks provides a means of interpreting the Cenozoic history of the region. Following a long period of intermittent volcanic activity and sedimentation, extensive sheets of andesitic and dacitic ignimbrites were erupted during Late Miocene time from fissure sources which appear to have been located near the now-subsided central portion of the graben. Near the coast, ignimbrites flowed across a flat lagunal shore overwhelming and burying the tropical vegetation and finally coming to rest in shallow water. Unusual textures and chaotic mixtures of pumice with sediments and silicified wood characterize the bases of many of the water-laid ignimbrites. Subsidence of the graben does not appear to have occurred concurrently with the ignimbrite eruptions but followed them closely near the end of the Miocene or the beginning of Pliocene time. Subsequent activity has been confined to relatively smaller eruptions from central vents near the boundary faults of the graben.     

Volatile loss from melt inclusions in pyroclasts of differing sizes

We have investigated the loss of H₂O from olivine-hosted melt inclusions (MIs) by designing an experiment using tephra samples that cooled at different rates owing to their different sizes: ash, lapilli, and bomb samples that were deposited on the same day (10/17/74) of the sub-Plinian eruption of Volcán de Fuego in Guatemala. Ion microprobe, laser ablation-ICPMS, and electron probe analyses show that MIs from ash and lapilli record the highest H₂O contents, up to 4.4 wt%. On the other hand, MIs from bombs indicate up to 30 % lower H₂O contents (loss of ~1 wt% H₂O) and 10 % post-entrapment crystallization of olivine. This evidence is consistent with the longer cooling time available for a bomb-sized clast, up to 10 min for a 3–4-cm radius bomb, assuming conductive cooling and the fastest H diffusivities measured in olivine (D~10^?9 to 10^?10 m²/s). On the other hand, several lines of evidence point to some water loss prior to eruption, during magma ascent and degassing in the conduit. Thus, results point to both slower post-eruptive cooling and slower magma ascent affecting MIs from bombs, leading to H₂O loss over the timescale of minutes to hours. The important implication of this study is that a significant portion of the published data on H₂O concentrations in olivine-hosted MIs may reflect unrecognized H₂O loss via diffusion. This work highlights the importance of reporting clast and MI sizes in order to assess diffusive effects and the potential benefit of using water loss as a chronometer of magma ascent.     

Probabilistic domain decomposition for the solution of the two-dimensional magnetotelluric problem

Probabilistic domain decomposition is proposed as a novel method for solving the two-dimensional Maxwell’s equations as used in the magnetotelluric method. The domain is split into non-overlapping sub-domains and the solution on the sub-domain boundaries is obtained by evaluating the stochastic form of the exact solution of Maxwell’s equations by a Monte-Carlo approach. These sub-domains can be naturally chosen by splitting the sub-surface domain into regions of constant (or at least continuous) conductivity. The solution over each sub-domain is obtained by solving Maxwell’s equations in the strong form. The sub-domain solver used for this purpose is a meshless method resting on radial basis function-based finite differences. The method is demonstrated by solving a number of classical magnetotelluric problems, including the quarter-space problem, the block-in-half-space problem and the triangle-in-half-space problem.     

Partitioning of Na between olivine and melt: An experimental study with application to the formation of meteoritic Na2O-rich chondrule glass and refractory forsterite grains

We report experimentally determined 1 atm olivine/melt D_Na partitioning data for low f_O2, a variety of melt compositions and a temperature range of 1325-1522 °C. We demonstrated that high-current electron microprobe analyses (EPMA, I = 500 nA, 600 s on the peak) allow quantitative determination of Na₂O in olivine down to ∼10 μg/g. The mean olivine/melt D_Na from 12 experimental runs is 0.0031 ± 0.0007 (1σ). This is the recommended value for low pressures and a wide range of natural compositions.This result is applied to the problem of the origin of alkalis in chondrules and the formation of chondritic refractory forsterite grains. The data on Semarkona (LL3.0) chondrules show that Na₂O is primordial and was present during olivine crystallization. For refractory forsterite grains from Murchison (CM2), we demonstrate that high CaO contents are not a result of equilibration with Na₂O-rich melts, but require high activities of CaO during their formation.     

Fluid entrapment and reequilibration during subduction and exhumation: A case study from the high-grade Nestos shear zone, Central Rhodope, Greece

A fluid inclusion study on metamorphic minerals of successive growth stages was performed on highly deformed paragneisses from the Nestos Shear Zone at Xanthi (Central Rhodope), in which microdiamonds provide unequivocal evidence for ultrahigh-pressure (UHP) metamorphism. The correlation of fluid inclusion density isochores and fluid inclusion reequilibration textures with geothermobarometric data and the relative chronology of micro- and macro-scale deformation stages allow a better understanding of both the fluid and metamorphic evolution along the P–T–d path. Textural evidence for subduction towards the NE is recorded by the orientation of intragranular NE-oriented fluid inclusion planes and the presence of single, annular fluid inclusion decrepitation textures. These textures occur within quartz “foam” structures enclosed in an earlier generation of garnets with prolate geometries and rarely within recrystallized matrix quartz, and reequilibrated both in composition and density during later stages of exhumation. No fluid inclusions pertaining to the postulated ultrahigh-pressure stage for microdiamond-bearing garnet–kyanite–gneisses have yet been found. The prolate shape of garnets developed during the earliest stages of exhumation that is recorded structurally by (L  S) tectonites, which subsequently accommodated progressive ductile SW shearing and folding up to shallow crustal levels. The majority of matrix kyanite and a later generation of garnet were formed during SW-directed shear under plane-strain conditions. Fluid inclusions entrapped in quartz during this stage of deformation underwent density loss and transformed to almost pure CO₂ inclusions by preferential loss of H₂O. Those inclusions armoured within garnet retained their primary 3-phase H₂O–CO₂ compositions. Reequilibration of fluid inclusions in quartz aggregates is most likely the result of recrystallization along with stress-induced, preferential H₂O leakage along dislocations and planar lattice defects which results in the predominance of CO₂ inclusions with supercritical densities. Carbonic fluid inclusions from adjacent kyanite–corundum-bearing pegmatoids and, the presence of shear-plane-parallel fluid inclusion planes within late quartz boudin structures consisting of pure CO₂-fluid inclusions with negative crystal shapes, bear witness of the latest stage of deformation by NE-directed extensional shear.This study shows that the textures of early fluid inclusions that formed already during the prograde metamorphic path can be preserved and used to derive information about the kinematics of subduction that is difficult to obtain from other sources. The textures of early inclusions, together with later generations of unaltered primary and secondary inclusions in metamorphic index minerals that can be linked to specific deformation stages and even P–T conditions, are a welcome supplement for the reconstruction of a rather detailed P–T–d path.     

Effect of oxygen fugacity on the H2O storage capacity of forsterite in the carbon-saturated systems

High pressure experiments have been performed in the systems Mg₂SiO₄-C-O-H and Mg₂SiO₄-K₂CO₃-C at 6.3 GPa and 1200 to 1600 °C using a split-sphere multi-anvil apparatus. In the Mg₂SiO₄-C-O-H system the composition of fluid was modeled by adding different amounts of water and stearic acid. The fO₂ was controlled by the Mo-MoO₂ or Fe-FeO oxygen buffers. Several experiments in the Mg₂SiO₄-C-O-H system and all experiments in the Mg₂SiO₄-K₂CO₃-C system have been conducted without buffering the fO₂. Forsterite in the system Mg₂SiO₄-K₂CO₃-C does not reveal OH absorption bands in the IR spectra, while forsterite coexisting with carbon-bearing fluid and silicate melt at logfO₂ from FMQ-2 to FMQ-5 (from 2 to 5 log units below fayalite-magnetite-quartz oxygen buffer) contains 800-1850 wt. ppm H₂O. The maximum concentrations were detected at 1400 °C and FMQ-3.5. We observed an increase in the solidus temperature in the system Mg₂SiO₄-C-O-H from 1200 to above 1600 °C with log fO₂ decreasing from FMQ-2 to FMQ-5. The increase of the solidus temperature and the broadening of the stability field of the H₂O-H₂-CH₄ subsolidus fluid phase at 1400-1600 °C explain the high H₂O storage capacity of forsterite relative to that crystallized from carbon-free, oxidized, hydrous, silicic melt. At temperatures above 1400 °C liquidus forsterite precipitated along with diamond from oxidized (FMQ-1) carbonate-silicate melt and from silicate melt dissolving the moderately reduced C-O-H fluid (from FMQ-2 to FMQ-3.5). Formation of diamond was not detected under ultra-reduced conditions (FMQ-5) at 1200-1600 °C. Olivine co-precipitating with diamond from dry carbonate-silicate or hydrous-silicic fluid/melt can provide information on the H₂O contents and speciation of the diamond-forming media in the mantle. The conditions for minimum post-crystallization alteration of olivine and its hydrogen content are discussed.     

Mineral associations and character of isotopically anomalous organic material in the Tagish Lake carbonaceous chondrite

We report a coordinated analytical study of matrix material in the Tagish Lake carbonaceous chondrite in which the same small (?20 μm) fragments were measured by secondary ion mass spectrometry (SIMS), transmission electron microscopy (TEM), energy-dispersive X-ray spectroscopy (EDS), electron energy-loss spectroscopy (EELS), and X-ray absorption near-edge spectroscopy (XANES). SIMS analysis reveals H and N isotopic anomalies (hotspots), ranging from hundreds to thousands of nanometers in size, which are present throughout the fragments. Although the differences in spatial resolution of the SIMS techniques we have used introduce some uncertainty into the exact location of the hotspots, in general, the H and N isotopic anomalies are spatially correlated with C enrichments, suggesting an organic carrier. TEM analysis, enabled by site-specific extraction using a focused-ion-beam scanning-electron microscope, shows that the hotspots contain an amorphous component, Fe-Ni sulfides, serpentine, and mixed-cation carbonates. TEM imaging reveals that the amorphous component occurs in solid and porous forms, EDS indicates that it contains abundant C, and EELS and XANES at the C K edge reveal that it is largely aromatic. This amorphous component is probably macromolecular C, likely the carrier of the isotopic anomalies, and similar to the material extracted from bulk samples as insoluble organic matter. However, given the large sizes of some of the hotspots, the disparity in spatial resolution among the various techniques employed in our study, and the phases with which they are associated, we cannot entirely rule out that some of the isotopic anomalies are carried by inorganic material, e.g., sheet silicates. The isotopic composition of the organic matter points to an initially primitive origin, quite possibly within cold interstellar clouds or the outer reaches of the solar protoplanetary disk. The association of organic material with secondary phases, e.g., serpentine and carbonates, suggests that the organic matter was susceptible to parent-body processing, and thus, isotopic dilution.     

Origin of the northern Indus Fan and Murray Ridge, Northern Arabian Sea: interpretation from seismic and magnetic imaging

The nature and origin of the sediments and crust of the Murray Ridge System and northern Indus Fan are discussed. The uppermost unit consists of Middle Miocene to recent channel–levee complexes typical of submarine fans. This unit is underlain by a second unit composed of hemipelagic to pelagic sediments deposited during the drift phase after the break-up of India–Seychelles–Africa. A predrift sequence of assumed Mesozoic age occurring only as observed above basement ridges is composed of highly consolidated rocks. Different types of the acoustic basement were detected, which reflection seismic pattern, magnetic anomalies and gravity field modeling indicate to be of continental character. The continental crust is extremely thinned in the northern Indus Fan, lacking a typical block-faulted structure. The Indian continent–ocean transition is marked on single MCS profiles by sequences of seaward-dipping reflectors (SDR). In the northwestern Arabian Sea, the Indian plate margin is characterized by several phases of volcanism and deformation revealed from interpretation of multichannel seismic profiles and magnetic anomalies. From this study, thinned continental crust spreads between the northern Murray Ridge System and India underneath the northern Indus Fan.     

Dynamic propagation criteria for catastrophic failure in planar landslides

Quantitative assessment of the risk of submarine landslides is an essential part of the design process for offshore oil and gas developments in deep water, beyond the continental shelf. Landslides may be triggered by a reduction in shear strength of subsea sediments over a given zone, caused for example by seismic activity. Simple criteria are then needed to identify critical conditions whereby the zone of weakness could grow catastrophically to cause a landslide. A number of such criteria have been developed over the last decade, based either on ideas drawn from fracture mechanics, or considering the equilibrium of the initial weakened zone and adjacent process zones of gradually softening material. Accounting for the history of the weak zone initiation is critical for derivation of reliable propagation criteria, in particular considering dynamic effects arising from accumulating kinetic energy of the failing material, which will allow the failure to propagate from a smaller initial zone of weakened sediments. Criteria are developed here for planar conditions, taking full account of such dynamic effects, which are shown to be capable of reducing the critical length of the softened zone by 20% or more compared with criteria based on static conditions. A numerical approach is used to solve the governing dynamic equations for the sliding material, the results from which justify assumptions that allow analytical criteria to be developed for the case where the initial softening occurs instantaneously. The effect of more gradual softening is also explored. Copyright © 2016 John Wiley & Sons, Ltd.