Genesis of a new slab tear fault in the Indo-Australian plate,offshore northern Sumatra,Indian ocean

Following the December 2004 and March 2005 major shallow foci inter-plate earthquakes in the north Sumatra region, a slab-tear fault located within the subducting Indian plate ruptured across the West Sunda Trench (WST) within the marginal intra-plate region. Trend, length and movement pattern of this New Tear Fault (NTF) segment is almost identical to another such slab-tear fault mapped previously by Hamilton (1979), located around 160 km south of NTF. Seismic activity along the NTF remained quasi-stable till the end of the year 2011, when an earthquake of magnitude 7.2 occurred on 10.01.2012 just at the tip of NTF, only around ～100 km within the intra-plate domain west of WST. The NTF rupture propagated further towards SSW with the generation of two more large earthquakes on 11.04.2012. The foreshock (10.01.12; M7.2) — mainshock (11.04.12; M 8.6) — aftershock (11.04.12; M 8.2) sequence along with numerous smaller magnitude aftershocks unmistakably define the extension of NTF, a slab-tear fault that results tectonic segmentation of the convergent plate margin. Within the intra-plate domain most earthquakes display consistent left-lateral strike slip mechanism along NNE trending fault plane.     

Tectonic imprints within a granite exposed near Srinagar,Rajasthan, India

Partial melting in the middle to lower crustal level produces melts of granitic composition during orogeny. Thrusts play a vital role in their exhumation after consolidation of these granitic melts. In this paper we focus on one such granite along the eastern margin of the Delhi Fold Belt (DFB) rocks near Srinagar, Rajasthan, India. This is the first report of granite within the area and holds a key stratigraphic position in the entire rock package. The said granite is found to be intrusive to the DFB metasediments as well as their basement popularly known as the Banded Gneissic Complex (BGC). We disentangle the deformation fabrics seen within the granite and associated DFB metasediments, suggesting that subsequent to emplacement and consolidation, the granite has co-folded along with the country rocks. Three deformational events could be identified within the DFB metasediments namely, D_1D, D_2D and D_3D. The peak metamorphism was achieved in the D_1D event. The granite magma is generated and emplaced late syn-kinematic to D_1D and thereafter is deformed by D_2D and D_3D producing D_1G and D_2G structural fabrics. These compressive deformations resulted in the collapse of the basin; the combined package of DFB rocks and the granite was thrusted eastwards over the basement rocks. The tectonic transport direction during thrusting is suggested eastwards from our structural analysis. Transverse faults developed perpendicular to the length of the granite have led to partitioning of the strain thereby showing a heterogeneity in the development of fabric within it.     

Plasma clouds associated with Comet P/Borrelly dust impacts

The NASA DS1 spacecraft encountered Comet P/Borrelly on September 22, 2001 at a distance of ∼2171 km on the sunward side of the comet. The flyby speed was ∼16.5 km s⁻¹. Using high temporal resolution (50 μs) absolute electric field amplitude measurements from a ∼1 m dipole antenna, new features of plasma clouds created by cometary dust impacts have been detected. The pulses have 1/e exponential decays of ∼650 μs duration, exponentially shaped overshoots with rise times of ∼2 ms, and exponential-shaped overshoot decay times of ∼10 ms. Assuming a plasma temperature of 10⁴ K, these pulse features have been explained as plasma cloud space charge effects from the electron, proton and heavy ion portions of the clouds passing the antenna. Complex pulse shapes were also detected. These are believed to be due to either plasma cloud scattering off of the spacecraft, or to secondary impacts. Small electric pulses of duration 10-15 ms of cometary origin were detected but are presently unexplained. The electric component of the plasma wave spectra at closest approach had an f^−2.4 power law shape from 10 Hz to 1 kHz. The electron cyclotron frequency was approximately 1 kHz. One possible explanation of the wave spectrum is that whistler mode waves associated with phase steepened cometary plasma waves are dispersed, leading to the broad spectrum. Finally, based on the present results, a new type of low-cost, large-area dust detector is proposed.     

Exposure of developing countries to sea-level rise and storm surges

An increase in sea surface temperature is strongly evident at all latitudes and in all oceans. The scientific evidence to date suggests that increased sea surface temperature will intensify cyclone activity and heighten storm surges. The paper assesses the exposure of (coastal) developing countries to sea-level rise and the intensification of storm surges. Geographic Information System (GIS) software is used to overlay the best available, spatially-disaggregated global data on critical exposed elements (land, population, GDP, agricultural extent and wetlands) with the inundation zones projected with heightened storm surges and a 1 m sea-level rise. Country-level results indicate a significant increase in exposure of developing countries to these climate-induced changes.     

Modelling the pore fluid diffusion process in aftershock initiation for 2004 Sumatra earthquake: implications for marine geohazard estimation in the Andaman region

The role of fluid injection on the occurrence and migration path for the aftershocks of 2004 Sumatra earthquake (Mw 9.3) and January 2005 Andaman earthquake swarm within the aftershock sequence is investigated here from the viewpoint of pore fluid diffusion process. The Sumatra earthquake created a regionally extensive crustal rupture plane exceeding 1,200 km length below the Andaman Sea. The r–t plots (Shapiro et al. 1997) are constructed for these aftershocks in order to examine the role of poroelastic effects as rupturing progressed with time. Their main results are as follows: the r–t plot corresponding to first 3 h of aftershock activity (when only 44 events of mb ≥ 4.5 originated) reveals that 95% of the data points occurred below the modelled parabola with relatively high D value of 20 m²/s, whereas a significantly low D value of 3.5 m²/s characterises the aftershock activity for the first 24 h (when 420 events of mb ≥ 4.0 occurred). Here, the Coulomb stress was transferred from the main shock with a rapid imposition of normal stress, thus inducing the pore-pressure change that started diminishing almost immediately by fluid diffusion, at a rate, defined by the diminishing D value. The modelling results for fault seismicity at far off distances from the main epicentre are interpreted here as potential indicators for large-scale sub-seabed rupturing—consequent to stress changes induced by bending of the Indian Ocean plate. Bathymetric slopes under the Andaman subduction zone are particularly amenable to sub-marine slides where crustal E–W hinge faults inferred seismically cut across the N–S trending regional thrust and strike-slip faults. Seabed rupturing appears to allow deep-slab hydration in these areas, producing pressure gradients along the normal faults. These features are important since they can herald marine geohazards in the Andaman region.     

Role of transverse tectonics in the Himalayan collision: Further evidences from two contemporary earthquakes

Two contemporary earthquakes originating in the central Himalayan arc and its foredeep (Sikkim earthquake of 18.09.2011, M_w 6.9, h: 10–60 (?) km and Bihar-Nepal earthquake of 20.08.1988, Mw 6.8, h: 57 km) are commonly associated with transverse lineaments/faults traversing the region. Such lineaments/faults form active seismic blocks defining promontories for the advancing Indian Craton. These actually produce conjugate shear faulting pattern suggestive of pervasive crustal interplay deep inside the mountains. Focal mechanism solutions allow inferring that large part of the current convergence across the central Himalayan arc is accommodated by lateral slip. Similar slip also continues unabated in the densely populated foredeep for distances up to several tens of kilometers south of the Main Boundary Thrust (MBT).     

An Improved Approach to Non-Force-Free Coronal Magnetic Field Extrapolation

We develop an approach to deriving the three-dimensional non-force-free coronal magnetic field from vector magnetograms. Based on the principle of minimum dissipation rate, a general non-force-free magnetic field is expressed as the superposition of one potential field and two constant-α (linear) force-free fields. Each is extrapolated from its bottom boundary data, providing the normal component only. The constant-α parameters are distinct and determined by minimizing the deviations between the numerically computed and measured transverse magnetic field at the bottom boundary. The boundary conditions required are at least two layers of vector magnetograms, one at the photospheric level and the other at the chromospheric level, presumably. We apply our approach to a few analytic test cases, especially to two nonlinear force-free cases examined by Schrijver et al. (Solar Phys. 235, 161, 2006). We find that for one case with small α parameters, the quantitative measures of the quality of our result are better than the median values of those from a set of nonlinear force-free methods. The reconstructed magnetic-field configuration is valid up to a vertical height of the transverse scale. For the other cases, the results remain valid to a lower vertical height owing to the limitations of the linear force-free-field solver. Because our method is based on the fast-Fourier-transform algorithm, it is much faster and easy to implement. We discuss the potential usefulness of our method and its limitations.     

CO<Subscript>2</Subscript> content of andesitic melts at graphite-saturated upper mantle conditions with implications for redox state of oceanic basalt source regions and remobilization of reduced carbon from subducted eclogite

We have performed experiments to determine the effects of pressure, temperature and oxygen fugacity on the CO₂ contents in nominally anhydrous andesitic melts at graphite saturation. The andesite composition was specifically chosen to match a low-degree partial melt composition that is generated from MORB-like eclogite in the convective, oceanic upper mantle. Experiments were performed at 1–3 GPa, 1375–1550?°C, and fO₂ of FMQ ?3.2 to FMQ ?2.3 and the resulting experimental glasses were analyzed for CO₂ and H₂O contents using FTIR and SIMS. Experimental results were used to develop a thermodynamic model to predict CO₂ content of nominally anhydrous andesitic melts at graphite saturation. Fitting of experimental data returned thermodynamic parameters for dissolution of CO₂ as molecular CO₂: ln(K ⁰) = ?21.79?±?0.04, ΔV ⁰?=?32.91?±?0.65 cm³mol^?1, ΔH ⁰?=?107?±?21 kJ mol^?1, and dissolution of CO₂ as CO₃ ^2?: ln(K ⁰ ) = ?21.38?±?0.08, ΔV ⁰?=?30.66?±?1.33 cm³ mol^?1, ΔH ⁰?=?42?±?37 kJ mol^?1, where K ⁰ is the equilibrium constant at some reference pressure and temperature, ΔV ⁰ is the volume change of reaction, and ΔH ⁰ is the enthalpy change of reaction. The thermodynamic model was used along with trace element partition coefficients to calculate the CO₂ contents and CO₂/Nb ratios resulting from the mixing of a depleted MORB and the partial melt of a graphite-saturated eclogite. Comparison with natural MORB and OIB data suggests that the CO₂ contents and CO₂/Nb ratios of CO₂-enriched oceanic basalts cannot be produced by mixing with partial melts of graphite-saturated eclogite. Instead, they must be produced by melting of a source containing carbonate. This result places a lower bound on the oxygen fugacity for the source region of these CO₂-enriched basalts, and suggests that fO₂ measurements made on cratonic xenoliths may not be applicable to the convecting upper mantle. CO₂-depleted basalts, on the other hand, are consistent with mixing between depleted MORB and partial melts of a graphite-saturated eclogite. Furthermore, calculations suggest that eclogite can remain saturated in graphite in the convecting upper mantle, acting as a reservoir for C.     

Modeling of shear failure in fractured reservoirs with a porous matrix

A finite volume-based numerical modeling framework using a hierarchical fracture representation (HFR) has been developed to compute flow-induced shear failure. To accurately capture the mechanics near fracture manifolds, discontinuous basis functions are employed which ensure continuity of the displacement gradient across fractures. With these special basis functions, traction and compressive forces on the fracture segment can be calculated without any additional constraints, which is extremely useful for estimating the irreversible displacement along the fracture (slip) based on a constitutive friction law. The method is further extended to include slip-dependent hydraulic aperture change and grid convergent results are obtained. Further, the change in hydraulic aperture is modeled using an asymptotic representation which respects the experimentally observed behavior of pore volume dilation due to shear slip. The model allows the initial rapid increase in hydraulic aperture due to shear slip and asymptotically approaches a finite value after repeated shearing of a fracture segment. This aperture increase is the only feedback for mechanics into the fluid flow for a linear elastic mechanics problem. The same model is also extended to include poroelastic relations between flow and mechanics solver. The grid convergence result in the case of poroelastic flow-mechanics coupling for flow-induced shear failure is also obtained. This proves the robustness of the numerical and analytical modeling of fracture and friction in the extended finite volume method (XFVM) set-up. Finally, a grid convergent result for seismic moment magnitude for single fracture and fracture network with random initial hydraulic and friction properties is also obtained. The b-value, which represents the slope of seismic moment occurrence frequency decay vs seismic moment magnitude, which is approximately constant in a semi-logarithmic plot, is estimated. The numerical method leads to converged b-values for both single fracture and fracture network simulations, as grid and time resolutions are increased. For the resulting linear system, a sequential approach is used, that is, first, the flow and then the mechanics problems are solved. The new modeling framework is very useful to predict seismicity, permeability, and flow evolution in geological reservoirs. This is demonstrated with numerical simulations of enhancing a geothermal system.     

Finite element analysis of axisymmetric flow into partially penetrating wells

The finite element method was used to solve the problem of flow to a partially penetrating strainer well in a leaky confined aquifer under steady and unsteady conditions. Predictions of the model were compared with reported pumping test data. A good agreement was obtained between them. The effects of penetration and screen position on discharge through imperfect well were studied. Results were compared with analytical solutions and presented graphically.