An explosive component in a December 2020 Milan earthquake suggests outgassing of deeply recycled carbon

Carbon dragged at sub-arc depths and sequestered in the asthenospheric upper mantle during cold subduction is potentially released after millions of years during the breakup of continental plates. However, it is unclear whether these deep-carbon reservoirs can be locally remobilized on shorter-term timescales. Here we reveal the fate of carbon released during cold subduction by analyzing an anomalously deep earthquake in December 2020 in the lithospheric mantle beneath Milan (Italy), above a deep-carbon reservoir previously imaged in the mantle wedge by geophysical methods. We show that the earthquake source moment tensor includes a major explosive component that we ascribe to carbon-rich melt/fluid migration along upper-mantle shear zones and rapid release of about 17,000 tons of carbon dioxide when ascending melts exit the carbonate stability field. Our results underline the importance of carbon-rich melts at active continental margins for emission budgets and suggest their potential episodic contributions to atmospheric carbon dioxide.

Mapping Main Structures and Related Mineralization of the Arabian Shield (Saudi Arabia) Using Sharp Edge Detector of Transformed Gravity Data

Saudi Arabia covers most of the Arabian Peninsula and is characterized by tectonic regimes ranging from Precambrian to Recent. Using gravity data to produce the lateral boundaries of subsurface density bodies, and edge detection of potential field data, a new subsurface structural map was created to decipher the structural framework controls on the distribution of gold deposits in Saudi Arabia. Moreover, we detected the relationships between major structures and mineral accumulations, thereby simultaneously solving the problem of edge detectors over complex tectonic patterns for both deeper and shallower origins. Analytic signal (ASg), theta map (TM), TDX, and softsign function (SF) filters were applied to gravity data of Saudi Arabia. The results unveil low connectivity along the Najd fault system (NFS) with depth, except perhaps for the central zones along each segment. The central zones are the location of significant gold mineralization, i.e., Fawarah, Gariat Avala, Hamdah, and Ghadarah. Moreover, major fault zones parallel to the Red Sea extend northward from the south, and their connectivity increases with depth and controls numerous gold mines, i.e., Jadmah, Wadi Bidah, Mamilah, and Wadi Leif. These fault zones intersect the NFS in the Midyan Terrane at the northern part of the AS, and their conjugation is suggested to be favorable for gold mineralization. The SF maps revealed the boundary between the Arabian Shield and Arabian Shelf, which comprises major shear zones, implying that most known mineralization sites are linked to post-accretionary structures and are not limited to the Najd fault system (NFS).

Lithofacies and petrochemical characterization of volcano-sedimentary sequence of Chandil Formation around Kharidih-Bareda area, Seraikela-Kharsawan District, Jharkhand: implications for uranium mineralization

Mesoproterozoic Chandil Formation (ca. 1600 Ma) of North Singhbhum Mobile Belt record numerous features of felsic volcaniclastics and felsic to intermediate volcanics preserved in the central sector of the fold belt around Kharidih-Bareda area, Seraikela-Kharsawan district, Jharkhand. The felsic volcanic rocks exhibit flow bands, autoclasts and layering of crystal mushes revealing viscous nature of eruptives. The volcaniclastic sediments comprise of significant proportion of volcanic epiclasts and accidental lithic fragments. These volcaniclastics have been categorized into five prominent lithofacies viz, stratified lapilli tuff, banded tuff, tuff with penecontemporaneous deformation, welded lapilli stones, vitric tuff and volcanic bombs by field and petrographic studies of outcrops and subsurface borehole cores. The welded lapilli tuffs display fiamme and eutaxitic texture. Interlayering of the volcaniclastics, which are most often pyrite-rich, with psamo-pelitic lithology like carbonaceous phyllite, variegated phyllite, quartzite and minor limestone is suggestive of marine euxenic depositional environment. Petrographic study of the volcaniclastics indicated presence of glass shards, garnet phenocrysts, spherules of tremolite, ovoid to lenticular accretionary lapilli along with devitrified glassy material. Compositionally these felsic volcanics and volcaniclastics are rhyodacitic to andesitic in nature with peraluminous to meta aluminous in character. A/CNK values vary from 0.52 to 2.42 in felsic volcanics and from 0.12 to 1.63 in volcaniclastics. Signatures of arc magmatism is indicated by low concentration of HFS elements such as Nb (5-17 ppm), Ga (11-17 ppm) and Y (5-28 ppm). Elevated intrinsic content of uranium (3-8 ppm), Th/U ratio ranging from 1.2 to 13.2, presence of metamict allanite and zircon in volcanics and volcaniclastics reveal their suitability as a prospective source for search of uranium mineralization. The volcanic-volcaniclastic-clastic association of the Chandil Formation provides an ideal situation where provenance and province both are available. Thus, suitable litho-structural locales such as the concealed shear zones sympathetic to the Dalma thrust and South Purulia Shear Zone within the volcano-sedimentary package of Chandil Formation may be targeted as preferable sites for locating concealed uranium mineralization.

The role of geophysics in the construction and stability assessment of tunnels

Geophysics play a vital role in the constructions of any major manmade structures in the world. One of those being the tunnels. In depth understanding of geophysical methods and a lot of information are needed in order to design a tunnel construction project. Comprehensive investigation on the ground condition has to be done before the field preparation study that will determine the stand-up time and the groundwater condition that may disrupt the tunnel construction. For tunnel stability assessment, an integration of geophysical methods is a must in order to obtain the most accurate results. Satellite imaging interpretation emphasizes on the structural tracing of negative lineament while field mapping emphasizes on location of underground seepage and major tectonic structures such as faults, joints and shear zones. Geoelectrical resistivity tomography survey is able to identify the differences in resistivity of Earth’s materials based on the water content inside of them. The best course of remediation could only be chosen once the output from all these studies are made available.

Landslide Susceptibility Assessment using Open-Source Data in the Far Western Nepal Himalaya: Case Studies from Selected Local Level Units

This paper explores openly available geo-spatial and earth observatory data to understand landslide risk in data scarce rural areas of Nepal. It attempts to explore the application of open-source data and analytical models to inform future landslide research. The first step of this procedure starts from the review of global open datasets, literatures and case studies relevant to landslide research. The second step is followed by the case study in one of the mountainous municipalities of Nepal where we tested the identified open-source data and models to produce landslide susceptibility maps. Past studies and experiences show that the major potential sites of landslide in Nepal are highly concentrated in a geologically weak area such as the active fault regions, shear zones, axis of folds and unfavorable setting of lithology. Triggering factors like concentrated precipitation, frequent earthquake phenomenon and haphazard infrastructural development activities in the marginally stable mountain slopes have posed serious issues of landslides mostly through the geologically weak regions. In this context, openly available geo-spatial datasets can provide baseline information for exploring the landslide hazard scenario in the data scarce areas of Nepal. This research has used the available open-source data to produce a landslide susceptibility map of the Bithadchir Rural Municipality in Bajhang District and Budiganga Municipality in Bajura District of the Sudurpaschim Province of Nepal. We used qualitative analysis to evaluate the parameters and assess the susceptibility of landslide; the result was classified into five susceptibility zones: Very High, High, Moderate, Low, and Very Low. Slope and Aspect were identified to be the major determinants for the assessment. This approach is applicable, specifically, for the preliminary investigation in the data scarce region using open data sources. Furthermore, the result can be used to plan and prioritize effective disaster risk reduction strategies.

DEM investigation on strain localization in a dense periodic granular assembly with high coordination number

Strain localization is one of key phenomena which have been studied extensively in geomaterials and for different kinds of materials including metals and polymers. This well-known phenomenon appears when structure/material is closed to failure. Theoretical, experimental, and numerical research have been dedicated to this subject for a long while. In the numerical aspects, strain localization inside the periodic granular assembly has not been well studied in the literature. In this paper, we investigate the occurrence and development of strain localization within a dense cohesive-frictional granular assembly with high coordination number under bi-periodic boundary conditions by Discrete Element Modeling (DEM). The granular assembly is composed of 2D circular disks and subjected to biaxial loading with constant lateral pressure. The results show that the formation of shear bands is of periodic type, consistent with the boundary conditions. This formation has the origins of the irreversible losing of cohesive contacts, viewed as micro-crackings which strongly concentrated in the periodic shear zones. This micromechanical feature is therefore strongly related to the strain localization observed at the sample scale. Finally, we also show that the strain localization is in perfect agreement with the sample’s displacement fluctuation fields.

Complex ⁴⁰Ar/³⁹Ar age spectra from low metamorphic grade rocks: resolving the input of detrital and metamorphic components in a case study from the Delamerian Orogen

Low metamorphic grade rocks contain both detrital minerals and minerals newly grown or partly recrystallised during diagenesis and metamorphism. However, rocks such as these typically yield complex 40Ar/39Ar age spectra that can be difficult to interpret. In this study, we have analysed a suite of variably deformed rocks from a region of low metamorphic grade within the c. 514–490 Ma Delamerian Orogen, South Australia. The samples analysed range from siltstone and shale to phyllite and all contain either muscovite or phengite determined by hyperspectral mineralogical characterisation. Furnace step heating 40Ar/39Ar analysis produced complex apparent age spectra with multiple age components. Using the concept of asymptotes that define minimum and maximum ages for different components, we interpret the age spectra to preserve a range of detrital mineral ages, along with younger components related to either cooling or deformation- induced recrystallisation. Two samples contain Mesoproterozoic detrital age components, up to c. 1170 Ma, while the c. 515 Ma Heatherdale Shale which has both c. 566 Ma and c. 530 Ma detrital components. All samples contain younger lower (younger) asymptotes in the age spectra defined from multiple heating steps that range from c. 476 to c. 460 Ma. One interpretation of these younger ages is that they are caused by post-metamorphic cooling. However, the shape of the age spectra and the degree of deformation in the phyllites suggest the ages may record recrystallisation of detrital minerals and/or new mica growth during deformation. Potentially these c. 476 to c. 460 Ma ages suggest deformation in the upper portion of the orogen was facilitated by movement along regional faults and shear zones up to around 20 million years after the cessation of deformation in the high-metamorphic grade regions of the Delamerian Orogen.

When did the North Anatolian fault reach southern Marmara, Turkey?

We dated syntectonic calcites on fault planes from the southern branch of the western North Anatolian fault (NAF) in northern Turkey using U-Th geochronology. We selected strike-slip faults that are kinematically related to the current regional strain field. The isotopic ages cluster around different periods during the past ~700 k.y. The most prominent cluster peak of 510.5 ± 9.5 ka (1σ) is consistent with the maximum cumulative strike-slip offset data and tectonic plate motions measured by GPS data, highlighting the fact that the present configuration of the NAF in the southern Marmara region started at ca. 500 ka or earlier. These new isotopic ages, combined with previous considerations of regional tectonics, reveal that faulting along the western NAF initiated primarily in the southern Marmara region at least a few hundred thousand years earlier than the timing suggested for the northern branch of the western NAF. This study presents an innovative approach to constrain the timing of initiation of currently active fault segments along the NAF in southern Marmara. U-Th geochronology of fault-hosted calcite thus has a wide application in determining absolute ages of fault episodes in wider shear zones along plate boundaries.

Determination of Abrasiveness in Copper-Gold Sulfide Ores: A Contribution to the Geometallurgical Model of the Sossego Deposit

The geological context of this study is established in the iron oxide-copper-gold (IOCG) deposit of Sossego (Canaã dos Carajás, Brazil), where hydrothermal alterations in shear zones concentrated the metals of interest and added new characteristics to the metavolcanic-sedimentary and granite rocks. The mineral transformation of rocks by hypersaline fluids enriched in metals and silica also modifies some metallurgical properties, such as abrasiveness. Special bench tests on rock drill cores are used in mapping the abrasiveness of rocks, with the Bond abrasion test being more commonly used in the mining industry, but it has a restrictive sampling protocol and mass requirement for geometallurgical studies. As a counterpoint, the test of the Laboratoire Central des Ponts et Chaussées/Central Laboratory of Bridges and Roads (LCPC) requires a smaller amount of fine material and a finer granulometric range. The study on the use of LCPC was implemented in 40 samples, using Bond Ai as a reference. The results showed a strong correlation between both methodologies (R2 = 95%), validating the use of LCPC to quantify abrasiveness in the Sossego mine. It was also possible to classify the most abrasive lithologies.

Whole-lithosphere shear during oblique rifting

Processes controlling the formation of continental whole-lithosphere shear zones are debated, but their existence requires that the lithosphere is mechanically coupled from base to top. We document the formation of a dextral, whole-lithosphere shear zone in the Death Valley region (DVR), southwest United States. Dextral deflections of depth gradients in the lithosphere-asthenosphere boundary and Moho are stacked vertically, defining a 20–50-km-wide, lower lithospheric shear zone with ~60 km of shear. These deflections underlie an upper-crustal fault zone that accrued ~60 km of dextral slip since ca. 8–7 Ma, when we infer that whole-lithosphere shear began. This dextral offset is less than net dextral offset on the upper-crustal fault zone (~90 km, ca. 13–0 Ma) and total upper-crustal extension (~250 km, ca. 16–0 Ma). We show that, before ca. 8–7 Ma, weak middle crust decoupled upper-crustal deformation from deformation in the lower crust and mantle lithosphere. Between 16 and 7 Ma, detachment slip thinned, uplifted, cooled, and thus strengthened the middle crust, which is exposed in metamorphic core complexes collocated with the whole-lithosphere shear zone. Midcrustal strengthening coupled the layered lithosphere vertically and therefore enabled whole-lithosphere dextral shear. Where thick crust exists (as in pre–16 Ma DVR), midcrustal strengthening is probably a necessary condition for whole-lithosphere shear.