Modeling of source parameters of the 15 December 2015 Deogarh earthquake of M<Subscript>w</Subscript> 4.0

We herein present source parameters and focal mechanism of a rare cratonic upper crustal earthquake of M_w4.0, which occurred at 8 km depth (centroid depth) below a region near Deogarh, Jharkhand. For our study, we used broadband waveform data from a seismic network of 15 three-component seismographs in the eastern Indian craton. The average seismic moment, moment magnitude and source radius are estimated to be 1.1 × 10¹⁵ N-m, 4.0 and 180.6 m, respectively. The high average stress drop of 14.27 MPa could be attributed to its lower-crustal origin. The mean corner frequency is calculated to be 4.1 Hz. To study the source mechanism, we perform a deviatoric constrained full waveform moment tensor inversion of multiple point sources on the band-passed (0.06 – 0.14 Hz) broadband displacement data of the Deogarh event, using ISOLA software. The best fit is obtained for the source at 8 km centroid depth, with a moment magnitude 3.7, and a right-lateral strike-slip mechanism with strike 162°, dip 72° and rake 169°. The P-axis orients N24°E, which is parallel to the direction of the absolute plate motion direction of Indian plate, while T-axis orients E-W, which is parallel to the strike of the pre-existing Damodar Graben (DG) of Gondwana age. The occurrence of this earthquake is attributed to the neotectonic reactivation of a fault associated with the E-W trending DG shear zone.     

Stability assessment of slide zones in Lesser Himalayan part of Yamunotri pilgrimage route,Uttarakhand, India

The present study deals with the slope stability analysis and geotechnical assessment of a part of pilgrimage route to one of the holy shrines of India, i.e. Yamunotri. The route also embraces a proposed site for 204 m high concrete gravity dam across River Yamuna near Lakhwar village with the aim of generation of 300 MW power. Several slide zones were identified and based on the discontinuity orientation, structural features and debris materials, they are recognised as planar, wedge or circular failure types. The morphological dimensions, structural data, orientations and geotechnical parameters of circular failure slides within weathered quartzites, phyllites and shales were evaluated by extensive field work and by laboratory tests for their stability analysis. Internal mechanisms, cohesion and angle of internal friction that resist shear stress in slope materials, obtained from direct shear test are showing minor variation due to relatively consistent grain size distribution and mineralogical composition. General slope is about 40°–43° with sparse vegetation. Materials in and around slide zones are sands with appreciable amount of fines falling in SP–SM category as per Unified Soil Classification System (USCS), except Niste B slide which has clean sands lying in SP group. Factor of safety, computed by requisite parameters of strength, soil and slope properties in circular failure charts, varies from 1.02 to 1.23 in dry conditions while it reduces below unity with increasing saturation, representing stable conditions in dry conditions but with seepage and saturation along the cracks and discontinuities during rainfall make them unstable. Presence of steep slopes, proximity to stream channels and significant weathered and jointed area are causative factors in the route with rainfall and road widening as major triggers initiating the failure.     

Experimental seismic shear force amplification in scaled RC cantilever shear walls with base irregularities

RC structural slender walls under large seismic excitation are expected to reach base moment capacity mainly affected by the first vibration mode. However, the base shear could be affected by higher modes once yielding in flexure has occurred, which might result in base shear underestimation in linear design. In this work, an experimental program is carried out on five RC rectangular walls 1:10 scaled. All five specimens considered irregularities at base, common in construction and one specimen did not consider shear reinforcement or boundary detailing. Tests are carried on a unidirectional shaking table and excitation is based on two Chile earthquake records with different intensities. Damage is concentrated at the wall base for all specimens; primary due to flexure with some participation of shear. For one of the records an average amplification of 1.3 is obtained, and a decrease in height of the resultant equivalent lateral force closes to 0.4 h_w. By increasing the intensity of the input record, amplification grows to an average of 1.7, while it decreases drastically when subjected to input records with low frequency content. No significant difference is observed in shear amplification in specimens with a base central opening, nor with setback, even though the cracking and failure mode was different for such specimens. Ductility demand shows no correlation when two different earthquakes are considered, whereas the frequency content and Arias intensity (I_a) of the input record directly affected the shear amplification.     

On the nature of intraplate earthquakes

Continental intraplate regions are characterized by uniform stresses over thousands of kilometers. Local stresses, with wavelengths of tens to hundreds of kilometers can accumulate at inhomogeneities lying within these regional fields. A variety of geological structures, herein called local stress concentrators (LSCs), act as elastic inhomogeneities. The temporal buildup of stress depends on the particular structure and its geometrical relationship with the regional stress field. The interaction of the local and the regional stress fields can result in the rotation of the latter over wavelengths of tens to hundreds of kilometers. This rotation can be detected by direct measurement or from seismicity data. Intraplate earthquakes (IPEs) result when the local stresses become comparable with their regional counterparts, i.e., hundreds of megapascals. Globally, most of the seismic energy release associated with IPEs occurs within old rifts which contain LSCs most favorable for stress buildup by stress inversion. Of the various LSCs, stepover en echelon faults are associated the largest IPEs. In low tectonic strain rate regions, IPEs are associated with larger stress drops. With the availability of a variety of LSCs, there is generally an absence of repeat earthquakes. Instead, successive earthquakes occur on different structures, leading to the observation of “roaming” earthquakes. These observations suggest a need for a reevaluation of seismic hazard estimation techniques. This study addresses some of these facets of the nature of IPEs with global examples, including a unique, detailed seismicity and geodetic data set collected in a dozen years following the 2001 M 7.7 Bhuj earthquake in western India.     

Non-Linear Theory and Power-Law Models for Information Integration and Mineral Resources Quantitative Assessments

Singular physical or chemical processes may result in anomalous amounts of energy release or mass accumulation that, generally, are confined to narrow intervals in space or time. Singularity is a property of different types of non-linear natural processes including cloud formation, rainfall, hurricanes, flooding, landslides, earthquakes, wildfires, and mineralization. The end products of these non-linear processes can be modeled as fractals or multifractals. Hydrothermal processes in the Earth’s crust can result in ore deposits characterized by high concentrations of metals with fractal or multifractal properties. Here we show that the non-linear properties of the end products of singular mineralization processes can be applied for prediction of undiscovered mineral deposits and for quantitative mineral resource assessment, whether for mineral exploration or for regional, national and global planning for mineral resource utilization. In addition to the general theory and framework for the non-linear mineral resources assessment, this paper focuses on several power-law models proposed for characterizing non-linear properties of mineralization and for geoinformation extraction and integration. The theories, methods, and computer system discussed in this paper were validated using a case study dealing with hydrothermal Au mineral potential in southern Nova Scotia, Canada.     

Heat Transport Processes in the Earth’s Crust

The heat of the Earth derives from internal and external sources. A heat balance shows that most of the heat provided by external sources is re-emitted by long-wavelength heat radiation and that the dominant internal sources are original heat and heat generated by decay of unstable radioactive isotopes. Understanding of the thermal regime of the Earth requires appreciation of properties and mechanisms for heat generation, storage, and transport. Both experimental and indirect methods are available for inferring the corresponding rock properties. Heat conduction is the dominant transport process in the Earth’s crust, except for settings where appreciable fluid flow provides a mechanism for heat advection. For most crustal and mantle rocks, heat radiation becomes significant only at temperatures above 1200°C.

中国"东部高原"东北部黏土矿物特征研究

通过对"东部高原"东北部地区中侏罗统—白垩系中黏土矿物分析,结果显示黏土矿物的组合及含量变化能够较好反映当时的气候、物源以及古地理变化;同时,伊利石结晶度不具有随地层深度发生变化的规律。黏土矿物组合以及伊利石结晶度在东部高原隆升与消亡前后,发生了比较明显的变化,这从另一方面反映了东部高原的存在。     

∼3,850 Ma tonalites in the Nuuk region,Greenland: geochemistry and their reworking within an Eoarchaean gneiss complex

The Eoarchaean (>3,600 Ma) Itsaq Gneiss Complex of southern West Greenland is dominated by polyphase orthogneisses with a complex Archaean tectonothermal history. Some of the orthogneisses have c. 3,850 Ma zircons, and they vary from rare single phase metatonalites to more common complexly banded migmatites. This is due to heterogeneous strain, in situ anatexis and granitic veining superimposed during younger tectonothermal events. In the single-phase tonalites with c. 3,850 Ma zircon, oscillatory-zoned prismatic zircon is all 3,850 Ma old, but shows patchy ancient loss of radiogenic Pb. SHRIMP spot analyses and laser ablation ICP-MS depth profiling show that thin (usually < 10 μm) younger (3,660–3,590 Ma and Neoarchaean) shells of lower Th/U metamorphic zircon are present on these 3,850 Ma zircons. Several samples with this simple zircon population occur on islands near Akilia. In contrast, migmatites usually contain more complex zircon populations, with often more than one generation of igneous zircon present. Additional zircon dating of banded gneisses across the Complex shows that samples with c. 3,850 Ma igneous zircon are not just a phenomenon restricted to Akilia and adjacent islands. For example, migmatites from Itilleq (c. 65 km from Akilia) contain variable amounts of oscillatory-zoned 3,850 Ma and 3,650 Ma zircon, interpreted, respectively, as the rock age and the time of crustal melting under Eoarchaean metamorphism. With only 110–140 ppm Zr in the tonalites and likely magmatic temperatures of >850°C, zircon solubility–melt composition relationships show that they were only one-third saturated in zircon. Any zircon entrained in the precursor magmas would thus have been highly soluble. Combined with the cathodoluminesence imaging, this demonstrates that the c. 3,850 Ma oscillatory zoned zircon crystallised out of the melt and hence gives a magmatic age. Thus the rare well-preserved tonalites and palaeosome in migmatites testify that c. 3,850 Ma quartzo–feldspathic rocks are a widespread (but probably minor) component in the Itsaq Gneiss Complex. C. 3,850 Ma zircon with negative Eu anomalies (showing growth in felsic systems) also occurs as detrital grains in rare c. 3,800 Ma metaquartzites and as inherited grains in some 3,660 Ma granites (sensu stricto). These demonstrate that still more c. 3,850 Ma rocks were present, but were recycled into Eoarchaean sediments and crustally derived granites. The major and trace element characteristics (e.g. LREE enrichment, HREE depletion, low MgO) of the best-preserved c. 3,850 Ma rocks are typical of Archaean TTG suites, and thus argue for crust formation processes involving important contributions from melting of hydrated mafic crust to the earliest Archaean. Five c. 3,850 Ma tonalites were selected as the best preserved on the basis of field criteria and zircon petrology. Four of these samples have overlapping initial ɛ_Nd (3,850 Ma) values from +2.9 to +3.6± 0.5, with the fourth having a slightly lower value of +0.6. These data provide additional evidence for a markedly LREE-depleted early terrestrial mantle reservoir. The role of c. 3,850 Ma crust should be considered in interpreting isotope signatures of the younger (3,800–3,600 Ma) rocks of the Itsaq Gneiss Complex. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Citizens as sensors: the world of volunteered geography

In recent months there has been an explosion of interest in using the Web to create, assemble, and disseminate geographic information provided voluntarily by individuals. Sites such as Wikimapia and OpenStreetMap are empowering citizens to create a global patchwork of geographic information, while Google Earth and other virtual globes are encouraging volunteers to develop interesting applications using their own data. I review this phenomenon, and examine associated issues: what drives people to do this, how accurate are the results, will they threaten individual privacy, and how can they augment more conventional sources? I compare this new phenomenon to more traditional citizen science and the role of the amateur in geographic observation.     

A case study to verify methods for estimating transmissivity distributions along boreholes

Two methods, both using data from fixed-interval-length hydraulic tests and geological mapping on different scales, have been utilised to estimate transmissivity distributions of individual fractures. Individual fracture properties are of importance in estimating transport of water and solutes; they influence the spreading of grout when sealing tunnels. One of the methods is non-parametric, while the other is based on combinatorics and the multiplication principle. The study verified the usefulness of both methods for estimating the probability of conductive and non-conductive fractures, and the transmissivity distributions. Data for estimating transmissivity distributions originate from flow logging of a borehole in Sweden. Fracture frequency data are from core logging and images using a borehole image processing system (BIPS). Estimated transmissivity distributions, based on flow log data from 3-m sections (test scale length, L: 3 m, and step length, dL: 3 m), were compared with transmissivities obtained directly from detailed, overlapping flow logging which identify individual conductive fractures (test scale, L: 0.5 m, step length, dL: 0.1 m). The latter data are considered to be an acceptable approximation of the individual fracture transmissivities; the agreement, compared with the calculated transmissivity distributions based on 3-m section data, is good for both methods.

---

Resumen Se han utilizado dos métodos, ambos procesando datos de pruebas hidráulicas de longitud e intervalo fijo y mapas geológicos en distintas escalas, para estimar las distribuciones de transmisividad de fracturas individuales. Las propiedades de fracturas individuales son de importancia para estimar el transporte de agua y solutos; las propiedades influyen en la propagación de lechada al sellar túneles. Uno de los métodos es no paramétrico mientras que el otro se basa en el principio de multiplicación y combinatorio. El estudio verifica la utilidad de ambos métodos para estimar la probabilidad de fracturas conductivas y no conductivas y las distribuciones de transmisividad. Los datos para estimar las distribuciones de transmisividad se originan del registro de flujo de un pozo en Suecia. Los datos de frecuencia de fracturas provienen de registros de núcleo e imágenes usando un Sistema de Procesamiento de Imágenes de Pozo (BIPS). Se compararon las distribuciones estimadas de transmisividad basadas en datos del registro de flujo en secciones de 3 m (longitud de escala de la prueba L: 3 m, y escalón de longitud dL: 3 m), con transmisividades obtenidas directamente de registros detallados de flujo superpuesto el cual identifica fracturas conductivas individuales (escala de la prueba, L: 0.5 m, escalón de longitud, dL: 0.1 m). Los últimos datos se consideran una aproximación aceptable de las transmisividades de fracturas individuales; la consistencia, comparada con las distribuciones de transmisividad calculada basadas en datos de sección de 3 m, es buena para ambos métodos.

---

Résumé Deux méthodes ont été testées afin d’estimer la distribution des transmissivités de fractures unitaires; toutes deux utilisent des données issues de tests hydrauliques effectués à intervalles réguliers et de cartographies géologiques à plusieurs échelles. La connaissance des propriétés des fractures unitaires est importante pour estimer le transport d’eau et de solutés; elles influencent en outre l’étalement du coulis lors du scellement des tunnels. La première méthode est non-paramétrique, et la seconde repose sur l’analyse combinatoire et le principe de multiplication. Cette étude a vérifié l’utilité des deux méthodes pour estimer la probabilité d’occurrence de fractures conductrices et non-conductrices, et la distribution des transmissivités. Les distributions des transmissivités ont été estimées à partir de diagraphies effectuées sur un forage situé en Suède. Les densités de fracturation proviennent de carottages et d’imagerie BIPS (Borehole Image Processing System). Les distributions estimées des transmissivités, basées sur des diagraphies de flux effectuées sur des sections de 3 m (échelle du test L:3 m; pas dL:3 m), ont été comparées aux transmissivités obtenues sur des diagraphies de flux détaillées et se chevauchant, identifiant les fractures conductrices unitaires (échelle du test L: 0.5 m ; pas dL: 0.1 m). Il est admis que ces dernières données représentent une approximation satisfaisante des transmissivités des fractures unitaires; les deux méthodes sont en accord avec les distributions des transmissivités calculées à partir des données acquises sur des sections de 3 m.