Evidence for Proterozoic Collision from Airborne Magnetic and Gravity Studies in Southern Granulite Terrain, India and Signatures of Recent Tectonic Activity in the Palghat Gap

The composite airborne total intensity map of the Southern Granulite Terrain (SGT) at an average elevation of 7000' (≈ 2100 m) shows bands of bipolar regional magnetic anomalies parallel to the structural trends suggesting the distribution of mafic/ultramafic rocks that are controlled by regional structures/shear zones and thrusts in this region. The spectrum and the apparent susceptibility map computed from the observed airborne magnetic anomalies provide bands of high susceptibility zones in the upper crust associated with known shear zones/thrusts such as Transition Zone, Moyar-Bhavani and Palghat-Cauvery Shear Zones (MBSZ and PCSZ). The quantitative modelling of magnetic anomalies across Transition Zone, MBSZ and PCSZ suggest the presence of mafic rocks of susceptibility (1.5-4.0 × 10⁻³ CGS units) in upper crust from 8-10 km extending up to about 21-22 km, which may represent the level of Curie point geotherm as indicated by high upper mantle heat flow in this section.Two sets of paired gravity anomalies in SGT and their modelling with seismic constraints suggest gravity highs and lows to be caused by high density mafic rocks along Transition Zone and Cauvery Shear Zone (CSZ) in the upper crust at depth of 6-8 km and crustal thickening of 45-46 km south of them, respectively. High susceptibility and high density rocks (2.8 g/cm³) along these shear zones supported by high velocity, high conductivity and tectonic settings suggest lower crustal mafic/ultramafic granulite rocks thrusted along them. These signatures with lower crustal rocks of metamorphic ages of 2.6-2.5 Ga north of PCSZ and Neoproterozoic period (0.6-0.5 Ga) south of it suggest that the SGT represents mosaic of accreted crust due to compression and thrusting. These observations along with N-verging thrusts and dipping reflectors from Dharwar Craton to SGT suggest two stages of N-S directed compression: (i) between Dharwar Craton and northern block of SGT during 2.6-2.5 Ga with Transition Zone and Moyar Shear towards the west as thrust, and (ii) between northern and southern blocks of SGT with CSZ as collision zone and PCSZ as thrust during Neoproterozoic period (0.6-0.5 Ga). The latter event may even represent just a compressive phase without any collision related to Pan-African event. The proposed sutures in both these cases separate gravity highs and lows of paired gravity anomalies towards north and south, respectively. The magnetic anomalies and causative sources related to Moyar Shear, MBSZ and PCSZ join with those due to Transition Zone, Mettur and Gangavalli Shears in their eastern parts, respectively to form an arcuate-shaped diffused collision zone during 2.6-2.5 Ga.Most of the Proterozoic collision zones are highlands/plateaus but the CSZ also known as the Palghat Gap represents a low lying strip of 80-100 km width, which however, appears to be related to recent tectonic activities as indicated by high upper mantle heat flow and thin crust in this section. It is supported by low density, low velocity and high conductive layer under CSZ and seismic activity in this region as observed in case of passive rift valleys. They may be caused by asthenospheric upwarping along pre-existing faults/thrusts (MBSZ and PCSZ) due to plate tectonic forces after the collision of Indian and Eurasian plates since Miocene time.     

Experimental technique of calibration of symmetrical air pollution models

Based on the inherent property of symmetry of air pollution models, a Symmetrical Air Pollution Model Index (SAPMI) has been developed to calibrate the accuracy of predictions made by such models, where the initial quantity of release at the source is not known. For exact prediction the value of SAPMI should be equal to 1. If the predicted values are overestimating then SAPMI is <1 and if it is underestimating then SAPMI is >1. Specific design for the layout of receptors has been suggested as a requirement for the calibration experiments. SAPMI is applicable for all variations of symmetrical air pollution dispersion models.     

Hydrochemistry of reservoirs of Damodar River basin, India: weathering processes and water quality assessment

Water samples collected from the six reservoirs of Damodar River basin in pre- and post-monsoon, have been analysed, to study the major ion chemistry and the weathering and geochemical processes controlling the water composition. Ca, Na and HCO₃ dominate the chemical composition of the reservoir water. The seasonal data shows a minimum concentration of most of the ions in post-monsoon and a maximum concentration in pre-monsoon seasons, reflecting the concentrating effects due to elevated temperature and increased evaporation during the low water level period of the pre-monsoon season. Water chemistry of the reservoirs strongly reflects the dominance of continental weathering aided by atmospheric and anthropogenic activities in the catchment area. Higher concentration of SO₄ and TDS in Panchet, Durgapur and Tenughat reservoirs indicate mining and anthropogenic impact on water quality. The high contribution of (Ca+Mg) to the total cations, high concentration of dissolved silica, relatively high (Na+K)/TZ⁺ ratio (0.3) and low equivalent ratio of (Ca+Mg)/(Na+K) suggests combined influence of carbonate and silicate weathering. Kaolinite is the possible mineral that is in equilibrium with the water, implying that the chemistry of reservoir water favours kaolinite formation. The calculated values of SAR, RSC and sodium percentage indicate the ‘excellent to good quality’ of water for irrigation uses.     

Pull-apart origin of the Satpura Gondwana basin, central India

The Gondwana basins of peninsular India are traditionally considered as extensional-rift basins due to the overwhelming evidence of fault-controlled synsedimentary subsidence. These basins indeed originated under a bulk extensional tectonic regime, due to failure of the attenuated crust along pre-existing zones of weakness inherited from Precambrian structural fabrics. However, disposition of the basins and their structural architecture indicate that the kinematics of all the basins cannot be extensional. To maintain kinematic compatibility with other basins as well as the bulk lateral extension, some basins ought to be of strike-slip origin. The disposition, shape and structural architecture of the Satpura basin, central India suggest that the basin could be a pull-apart basin that developed above a releasing jog of a left-stepping strike-slip fault system defined by the Son-Narmada south fault and Tapti north fault in consequence to sinistral displacement along WSW-ENE. Development of a sedimentary basin under the above-mentioned kinematic condition was simulated in model experiments with sandpack. The shape, relative size, stratigraphic and structural architecture of the experimental basin tally with that of the Satpura basin. The experimental results also provide insights into the tectono-sedimentary evolution of the Satpura basin in particular and pull-apart basins in general.     

Influence of fluids on deep crustal Jabalpur earthquake of 21, May 1997: Geophysical evidences

We present the geophysical evidences on the role of fluids for generation of the lower crustal Jabalpur earthquake (21 May 1997, mb 6.0, Mw 5.8), in the mid-continental fracture zone of the Indian Peninsular Shield. With a focal depth of 35 km, it indicates a high angled (< 62^∘ enclosed with maximum principal stress direction) reverse fault with small component of left-lateral strike slip in the lower crust. The Son-Narmada-Tapti (SONATA) magalineament, during the past two centuries, has experienced about 25 moderate to strong earthquakes; two of which namely the Son Valley (1927, M 6.5) and Jabalpur (21 May 1997) were disastrous. Historical earthquakes and recent earthquake swarms indicate a moderate to high seismicity in SONATA belt that is due to high strain accumulation, flexuring of the crust and neotectonic movements of the faults in the rift zones. By analyzing geophysical parameters such as Zero-Free air-based (ZFb) gravity anomalies (∼ −10 to –30 mGals), heat flow values (45–47 mWm⁻²), magneto-telluric values (1- Ohm m), strain rate (1.5 × 10⁻⁸) and failure stress conditions, we identify plausible causative factors for the occurrence of lower crustal earthquake in this region Fluids, due to dehydration of serpentinite in the lower crust, are suggested to be present in the earthquake source zone. The estimated pore-fluid factor for the Jabalpur earthquake (λ _v ) is 0.95. The diffusion of pore-pressure relaxation, represented as pressure perturbation generated by coseismic stress change was seen in the form of swarm activity two years prior to the Jabalpur earthquake. We suggest the existence of a deep pre-fractured zone with low shear stress (τ = 15–18 MPa) that indicates the presence of fluid filled fractured mafic material in the felsic crust, in critical state of unstable failure condition, and also fluid driven migration of swarm activity before the Jabalpur earthquake.     

A note on the lithosphere thickness and heat flow density of the Indian Craton from MAGSAT data

The Curie depth map of India compiled from MAGSAT data has been used for preparing the lithosphere thickness and the surface heat flow density maps of the Indian Craton, utilizing the concept of magnetothermometry. The lithosphere thicknesses of the major Indian geological units/provinces, as obtained from the prepared map, are found to be in reasonably good agreement with the previously published values for these regions. Also, the surface heat flow density values obtained from the prepared maps closely follow the previously published results. The maps are useful in providing first order estimates of lithosphere thickness and surface heat flow density of the important geological units/provinces of India.     

Operational mesoscale atmospheric dispersion prediction using a parallel computing cluster

An operational atmospheric dispersion prediction system is implemented on a cluster supercomputer for Online Emergency Response at the Kalpakkam nuclear site. This numerical system constitutes a parallel version of a nested grid meso-scale meteorological model MM5 coupled to a random walk particle dispersion model FLEXPART. The system provides 48-hour forecast of the local weather and radioactive plume dispersion due to hypothetical airborne releases in a range of 100 km around the site. The parallel code was implemented on different cluster configurations like distributed and shared memory systems. A 16-node dual Xeon distributed memory gigabit ethernet cluster has been found sufficient for operational applications. The runtime of a triple nested domain MM5 is about 4 h for a 24 h forecast. The system had been operated continuously for a few months and results were ported on the IMSc home page. Initial and periodic boundary condition data for MM5 are provided by NCMRWF, New Delhi. An alternative source is found to be NCEP, USA. These two sources provide the input data to the operational models at different spatial and temporal resolutions using different assimilation methods. A comparative study on the results of forecast is presented using these two data sources for present operational use. Improvement is noticed in rainfall forecasts that used NCEP data, probably because of its high spatial and temporal resolution     

Nonlinear Response of Single Piles in Sand Subjected to Lateral Loads using k hmax Approach

This paper presents results of analysis of full-scale pile load test data of 14 piles embedded in either loose or medium dense sands. The analysis was performed using two methods, p–y curve approach and a more recently developed k_hmax approach. Comparison of the results obtained using both the methods is also presented. A step-by-step analysis procedure is presented for predicting lateral load deflection response of single piles in sand using the k_hmax approach. The results presented show that the k_hmax approach has promise over the p–y curve approach because of its simplicity and the fact that it provides upper- and lower-bound curves, which are valuable guides to making engineering decisions. For loose sands, a new range of k_hmax values is recommended to better predict the lateral load–deflection response of single piles.     

Local thickening of the Cascadia forearc crust and the origin of seismic reflectors in the uppermost mantle

Seismic reflection profiles from three different surveys of the Cascadia forearc are interpreted using P wave velocities and relocated hypocentres, which were both derived from the first arrival travel time inversion of wide-angle seismic data and local earthquakes. The subduction decollement, which is characterized beneath the continental shelf by a reflection of 0.5 s duration, can be traced landward into a large duplex structure in the lower forearc crust near southern Vancouver Island. Beneath Vancouver Island, the roof thrust of the duplex is revealed by a 5–12 km thick zone, identified previously as the E reflectors, and the floor thrust is defined by a short duration reflection from a < 2-km-thick interface at the top of the subducting plate. We show that another zone of reflectors exists east of Vancouver Island that is approximately 8 km thick, and identified as the D reflectors. These overlie the E reflectors; together the two zones define the landward part of the duplex. The combined zones reach depths as great as 50 km. The duplex structure extends for more than 120 km perpendicular to the margin, has an along-strike extent of 80 km, and at depths between 30 km and 50 km the duplex structure correlates with a region of anomalously deep seismicity, where velocities are less than 7000 m s^− 1. We suggest that these relatively low velocities indicate the presence of either crustal rocks from the oceanic plate that have been underplated to the continent or crustal rocks from the forearc that have been transported downward by subduction erosion. The absence of seismicity from within the E reflectors implies that they are significantly weaker than the overlying crust, and the reflectors may be a zone of active ductile shear. In contrast, seismicity in parts of the D reflectors can be interpreted to mean that ductile shearing no longer occurs in the landward part of the duplex. Merging of the D and E reflectors at 42–46 km depth creates reflectivity in the uppermost mantle with a vertical thickness of at least 15 km. We suggest that pervasive reflectivity in the upper mantle elsewhere beneath Puget Sound and the Strait of Georgia arises from similar shear zones.     

Modelling Transport of Linearly Sorbing Solutes in a Single Fracture: Asymptotic Behavior of Solute Velocity and Dispersivity

The asymptotic behavior of the solute velocity and dispersivity for a system of parallel fractures with matrix diffusion is made using numerical modeling and theoretical analyses. The study is limited to linearly sorbing solutes with a constant continuous source boundary condition. Expressions are provided for solute velocity and effective dispersivity in terms of fracture porosity during asymptotic stage using spatial moment analyses. The importance of matrix porosity and fracture porosity on solute velocity as well as the relationship governing effective dispersivity and fracture porosity is discussed for both non-reactive and linearly sorbing solutes. By using a dimensionless effective dispersivity parameter it is shown that the relationship between the fracture porosity and dimensionless effective dispersivity is linear for non-reactive solutes. It is also shown that this holds true for the linearly sorbing solutes with the same proportionality constant.