Ground deformation monitoring by using the Permanent Scatterers Technique: The example of the Oltrepo Pavese (Lombardia, Italy)

The applicability of the Permanent Scatterers Synthetic Aperture Radar Interferometry (PSInSAR) technique for detecting and monitoring ground displacements was tested in the Oltrepo Pavese territory (Northern Italy, southern Lombardia), which could be representative of similar geological contexts in the Italian Apennines. The study area, which extends for almost 1100 km², is characterized by a complex geological and structural setting and the presence of clay-rich sedimentary formations. These characteristics make the Oltrepo Pavese particularly prone to several geological hazards: shallow and deep landslides, subsidence and swelling/shrinkage of the clayey soils. The PSInSAR technique used in this study overcomes most of the limitations of conventional interferometric approaches by identifying, within the area of interest, a set of “radar benchmarks” (PS), where very precise displacement measurements can be carried out. More than 90,000 PS were identified by processing Synthetic Aperture Radar (SAR) images acquired from 1992 to 2001 by the European Remote Sensing satellites (ERS). The PSInSAR application at a sub-regional scale detected slow ground deformations ranging from + 5 to − 16 mm/year, and resulting from various processes (landslides, swelling/shrinkage of clay soils and water pumping). The PS displacements were analysed by collecting data obtained through geological, geomorphologic field surveys, geotechnical analysis of the soils and the information was integrated within a landslide inventory and the damaged building inventory. Despite the limited number of landslide bodies with PS (7% of the inventoried landslides), the PS data helped to revise the state of activity of several landslides. Furthermore, some previously unknown unstable slopes were detected. Two areas of uplift and two areas of subsidence were identified.     

Permanent Scatterers for landslide investigations: outcomes from the ESA-SLAM project

Within the SLAM project (Service for Landslide Monitoring), launched in 2003 by the European Space Agency (ESA) the Permanent Scatterers (PS) technique, a multi-image interferometric approach, coupled with the interpretation of aerial-photos and optical satellite images, was carried out for landslide investigations. The PS analysis was applied at a regional scale as support for landslide inventory mapping and at local scale for the monitoring of single well-known slope movements. For the integration of the PS measurements within a landslide inventory the Arno river basin (Italy) was chosen as test site for the presence of a high number of mass movements (to date about 300 areas at high landslide risk and more than 27,000 individual landslides mapped by the institutional authorities). About 350 SAR images have been interferometrically processed by means of the PS technique, with the detection of about 600,000 PS. The use of optical images contributed spatial meaning to the point-wise information provided by the PS, making it easier to identify terrain features related to slope instability and the landslide boundaries. Here we describe the employed methodology and its impact in the updating of a preexisting landslide inventory. 6.8% of the total number of landslides were characterized by ground displacement measurements from the PS: 6.1% of already mapped landslides and 0.8% of new unstable areas detected through the PS analysis. Moreover, most of the PS are located in urban areas, showing that the proposed methodology is suitable for landslide mapping in areas with a quite high density of urbanization, but that over vegetated areas it still suffers from the limitations induced by the current space-borne SAR missions (e.g. temporal de-correlation). On the other hand, the use of InSAR for the monitoring of single slow landslides threatening built-up areas has provided satisfactory results, allowing the measurement of superficial deformations with high accuracy on the landslide sectors characterized by a good radar reflectivity and coherence.     

Permanent Scatterers for landslide investigations: outcomes from the ESA-SLAM project

Within the SLAM project (Service for Landslide Monitoring), launched in 2003 by the European Space Agency (ESA) the Permanent Scatterers (PS) technique, a multi-image interferometric approach, coupled with the interpretation of aerial-photos and optical satellite images, was carried out for landslide investigations. The PS analysis was applied at a regional scale as support for landslide inventory mapping and at local scale for the monitoring of single well-known slope movements. For the integration of the PS measurements within a landslide inventory the Arno river basin (Italy) was chosen as test site for the presence of a high number of mass movements (to date about 300 areas at high landslide risk and more than 27,000 individual landslides mapped by the institutional authorities). About 350 SAR images have been interferometrically processed by means of the PS technique, with the detection of about 600,000 PS. The use of optical images contributed spatial meaning to the point-wise information provided by the PS, making it easier to identify terrain features related to slope instability and the landslide boundaries. Here we describe the employed methodology and its impact in the updating of a preexisting landslide inventory. 6.8% of the total number of landslides were characterized by ground displacement measurements from the PS: 6.1% of already mapped landslides and 0.8% of new unstable areas detected through the PS analysis. Moreover, most of the PS are located in urban areas, showing that the proposed methodology is suitable for landslide mapping in areas with a quite high density of urbanization, but that over vegetated areas it still suffers from the limitations induced by the current space-borne SAR missions (e.g. temporal de-correlation). On the other hand, the use of InSAR for the monitoring of single slow landslides threatening built-up areas has provided satisfactory results, allowing the measurement of superficial deformations with high accuracy on the landslide sectors characterized by a good radar reflectivity and coherence.     

Testing the potential of Sentinel-1A TOPS interferometry for the detection and monitoring of landslides at local scale (Veneto Region,Italy)

The recent Sentinel-1 mission, started by the European Space Agency in April 2014, provides the scientific community with new capabilities for the monitoring of the Earth surface. In particular, the Terrain Observation by Progressive Scans imaging technique used in the Interferometric Wide swath acquisition mode permits to acquire data over very wide areas (250 km of swath extension) at 20-m spatial resolution, with 12-day revisit time, making it suitable for ground displacement monitoring applications. With more than 1 year of synthetic aperture radar images available, it is now possible to carry out monitoring activities of slow moving phenomena such as landslides at both regional and local scales. In this work, the potential of Sentinel-1A for the monitoring of shallow (from 2 to 6 m of depth) landslides occurring in the North-Eastern Italian Pre-Alps was tested. Two stacks of Sentinel-1A scenes acquired in both ascending and descending orbits were processed using the Permanent Scatterer Interferometry (PSI) technique. The results, analysed in terms of PS density and quality, were compared with the ERS-1/2 and ENVISAT PSI database available from the Italian National Cartographic Portal to assess the capabilities of Sentinel-1A in detecting and monitoring landslides in respect to the previous satellite missions. The results of this work show the great potential of Sentinel-1A in the continuous monitoring of landslide-prone territories even at local scale. The achievable results can provide information that is useful to delineate the spatial and temporal evolution of landslides and precisely assess their rates of deformation.     

Topographical changes revealed by high-resolution airborne LiDAR data: The 1999 Tsaoling landslide induced by the Chi–Chi earthquake

The 1999 Chi–Chi earthquake triggered the catastrophic Tsaoling landslide in central Taiwan. We mapped the landslide area and estimated the landslide volume, using a high-resolution digital elevation model from airborne LiDAR (Light Detection And Ranging), aerial photographs and topographic maps. The comparison between scar and deposit volumes, about 0.126 km³ and 0.150 km³ respectively, suggests a coseismic volume increase of 19% due to decompaction during landsliding. In July 2003, the scar and deposit volumes were about 0.125 km³ and 0.110 km³ respectively. These estimates suggest that 4 years after the event, the volume of landslide debris removed by river erosion was nearly 0.040 km³. These determinations are confirmed by direct comparison between the most accurate topographic models of the post-landslide period, indicating a very high erosion rate at the local scale (0.01 km³/year) for the deposit area of the landslide. Such a large value highlights the importance of landslide processes for erosion and long-term denudation in the Taiwan mountain belt.     

Monitoring of post-failure landslide deformation by the PS-InSAR technique at Lubietova in Central Slovakia

Interferometric synthetic aperture radar data from ERS and ENVISAT sensors were utilized in the analysis of the post-failure deformations in the area of Lubietova town in Central Slovakia. The catastrophic landslide of 1977 together with surrounding landslides in the Lubietova area were analysed with the help of persistent scatterers (PS) technique in order to evaluate recent and past deformations of the unstable slopes. Although long-term precise geodetic monitoring of the 1977 landslide revealed differential deformations inside the sliding mass, due to the lack of the PS located inside the landside caused by temporal decorrelation, unfortunately, these records could not be directly compared. The adjacent landslides with sufficient number of PS were analysed by transformation of the line of sight displacements recorded by the sensors to the slope vector direction. This procedure allowed identification of the precise boundaries of the actively moving landslide parts and the updating of the landslide inventory for the Lubietova area.     

极具潜力的空间对地观测新技术——合成孔径雷达干涉

“合成孔径雷达干涉(InSAR)”是近十年发展起来的空间对地观测遥感新技术。它具有从覆盖同一地区的星载(或机载)合成孔径雷达复数图像对提取干涉相位图,借助于雷达成像时的姿态数据重建地表三维模型(即数字高程模型)的巨大潜力。尤其是基于多幅雷达复数图像处理的差分干涉技术(D-InSAR)可以用于监测地表形变,精度可达厘米级甚至更高,其监测空间分辨率是前所未有的。介绍了InSAR和D-InSAR的基本原理,对影响干涉结果的一些重要因素做了分析,重点回顾和展望了差分干涉技术在与地表形变有关的地震监测和震后形变测量、地面下沉和山体滑坡、火山运动监测等方面应用的现状和前景。     

Melting of the subcontinental lithospheric mantle by the Emeishan mantle plume; evidence from the basal alkaline basalts in Dongchuan, Yunnan, Southwestern China

The Emeishan continental flood basalt (ECFB) sequence in Dongchuan, SW China comprises a basal tephrite unit overlain by an upper tholeiitic basalt unit. The upper basalts have high TiO₂ contents (3.2–5.2 wt.%), relatively high rare-earth element (REE) concentrations (40 to 60 ppm La, 12.5 to 16.5 ppm Sm, and 3 to 4 ppm Yb), moderate Zr/Nb and Nb/La ratios (9.3–10.2 and 0.6–0.9, respectively) and relatively high _{Nd (t)} values, ranging from − 0.94 to 2.3, and are comparable to the high-Ti ECFB elsewhere. The tephrites have relatively high P₂O₅ (1.3–2.0 wt.%), low REE concentrations (e.g., 17 to 23 ppm La, 4 to 5.3 ppm Sm, and 2 to 3 ppm Yb), high Nb/La (2.0–3.9) ratios, low Zr/Nb ratios (2.3–4.2), and extremely low _{Nd (t)} values (mostly ranging from − 10.6 to − 11.1). The distinct compositional differences between the tephrites and the overlying tholeiitic basalts cannot be explained by either fractional crystallization or crustal contamination of a common parental magma. The tholeiitic basalts formed by partial melting of the Emeishan plume head at a depth where garnet was stable, perhaps > 80 km. We propose that the tephrites were derived from magmas formed when the base of the previously metasomatized, volatile-mineral bearing subcontinental lithospheric mantle was heated by the upwelling mantle plume.     

Perspectives on basaltic magma crystallization and differentiation: Lava-lake blocks erupted at Mauna Loa volcano summit, Hawaii

Explosive eruptions at Mauna Loa summit ejected coarse-grained blocks (free of lava coatings) from Moku'aweoweo caldera. Most are gabbronorites and gabbros that have 0–26 vol.% olivine and 1–29 vol.% oikocrystic orthopyroxene. Some blocks are ferrogabbros and diorites with micrographic matrices, and diorite veins (≤ 2 cm) cross-cut some gabbronorites and gabbros. One block is an open-textured dunite.

The MgO of the gabbronorites and gabbros ranges  7–21 wt.%. Those with MgO > 10 wt.% have some incompatible-element abundances (Zr, Y, REE; positive Eu anomalies) lower than those in Mauna Loa lavas of comparable MgO; gabbros (MgO < 10 wt.%) generally overlap lava compositions. Olivines range Fo_83–58, clinopyroxenes have Mg#s  83–62, and orthopyroxene Mg#s are 84–63 — all evolved beyond the mineral-Mg#s of Mauna Loa lavas. Plagioclase is An_75–50. Ferrogabbro and diorite blocks have  3–5 wt.% MgO (TiO₂ 3.2–5.4%; K₂O 0.8–1.3%; La 16–27 ppm), and a diorite vein is the most evolved (SiO₂ 59%, K₂O 1.5%, La 38 ppm). They have clinopyroxene Mg#s 67–46, and plagioclase An_57–40. The open-textured dunite has olivine  Fo_83.5. Seven isotope ratios are ⁸⁷Sr/⁸⁶Sr 0.70394–0.70374 and ¹⁴³Nd/¹⁴⁴Nd 0.51293–0.51286, and identify the suite as belonging to the Mauna Loa system.

Gabbronorites and gabbros originated in solidification zones of Moku'aweoweo lava lakes where they acquired orthocumulate textures and incompatible-element depletions. These features suggest deeper and slower cooling lakes than the lava lake paradigm, Kilauea Iki, which is basalt and picrite. Clinopyroxene geobarometry suggests crystallization at < 1 kbar P. Highly evolved mineral Mg#s, < 75, are largely explained by cumulus phases exposed to evolving intercumulus liquids causing compositional ‘shifts.’ Ferrogabbro and diorite represent segregation veins from differentiated intercumulus liquids filter pressed into rigid zones of cooling lakes. Clinopyroxene geobarometry suggests < 300 bar P. Open-textured dunite represents olivine-melt mush, precursor to vertical olivine-rich bodies (as in Kilauea Iki). Its Fo_83.5 identifies the most primitive lake magma as  8.3 wt.% MgO. Mass balancing and MELTS show that such a magma could have yielded both ferrogabbro and diorite by ≥ 50% fractional crystallization, but under different fO₂: < FMQ (250 bar) led to diorite, and FMQ (250 bar) yielded ferrogabbro. These segregation veins, documented as similar to those of Kilauea, testify to appreciable volumes of ‘rhyolitic’ liquid forming in oceanic environments. Namely, SiO₂-rich veins are intrinsic to all shields that reached caldera stage to accommodate various-sized cooling, differentiating lava lakes.     

On the applicability of satellite SAR interferometry to landslide hazards detection in hilly areas: a case study of Shuicheng,Guizhou in Southwest China

Landslides - SAR Interferometry (InSAR) has been proven to be effective for measuring landslides deformation. However, the InSAR application of landslide mapping is limited by the blind observation...     

Origin of lamprophyres by the mixing of basic and alkaline melts in magma chamber in Beiya area, western Yunnan, China

Lamprophyres consisting mainly of diopside, phlogopite and K-feldspar formed in the early Tertiary around 60 Ma in the Beiya area and are characterized by low SiO₂ ± 46–50 wt.%), Rb (31–45 ppm) and Sr (225–262 ppm), high Al₂O₃, (11.2–13.1 wt.%), CaO (8.0–8.7 wt.%), MgO (11.5–12.1 wt.%), K₂O(4.9–5.5 wt.%), TiO₂ (2.9–3.3 wt.%) and REE (174–177 ppm), and compatible elements (e.g. Sc, Cr and Ni) and HSF elements (e.g. Th, U, Zr, Nb, Ta, Ti and Y), and low ¹⁴³Nd/¹⁴⁴Nd 0.512372–0.512536, middle ⁸⁷Sr/⁸⁶Sr 0.707322–0.707395, middle ²⁰⁶Pb/²⁰⁴Pb 18.50–18.59, ²⁰⁷Pb/²⁰⁴Pb 15.60–15.65 and ²⁰⁸Pb/²⁰⁴Pb 38.75–38.8. These rocks developed peculiar quartz megacrysts with poly-layer reaction zones, melt inclusions, and partial melted K-feldspar and plagioclase inclusions, and plastic shapes. Important features of these rocks include: (1) hybrid composition of elements, (2) abrupt increase of SiO₂ content of the melt, recorded by zoned diopside, (3) development of sanidine and aegirine-augite reaction zones, (4) alkaline melt and partial melted K-feldspar and plagioclase inclusions, (5) deformed quartz inclusions associated with quartz megacrysts, (6) the presence of quartz megacrysts in plastic shape with their parent melts, (7) the occurrence of olivine, high-MgO ilmenite and spinel inclusions within earlier formed diopside, phlogopite and magnetite. Median ⁸⁷Sr/⁸⁶Sr values between Tertiary alkaline porphyries in the Beiya area and the western Yunnan and Tertiary basalt in the western Yunnan indicate that the Beiya lamprophyre melts were derivative and resulted from the mixing between basic melts that were related to the partial melting of phenocrysts of spinel iherzolite from a mantle source. The alkaline melts originated from partial melting along the Jinshajiang subduction ductile shear zone at the contact between the buried Palaeo-Tethyan oceanic lithosphere and the upper mantle lithosphere. The alkaline melts are composed of 65% sanidine (Or₇₀Ab₂₈An₂) and 35% SiO₂. The melt mixing occurred in magma chambers in the middle-shallow crust at 8–10 km before the derivative lamprophyre melts intruded into the shallow cover in Beiya area. This mixing of basic and alkaline melts might represent a general process for the formation of lamprophyre in the western Yunnan.     

The origin of medium-K ankaramitic arc magmas from Lombok (Sunda arc, Indonesia): Mineral and melt inclusion evidence

High-calcium, nepheline-normative ankaramitic basalts (MgO > 10 wt.%, CaO/Al₂O₃ > 1) from Rinjani volcano, Lombok (Sunda arc, Indonesia) contain phenocrysts of clinopyroxene and olivine (Fo_85–92) with inclusions of spinel (Cr# 58–77) and crystallised melt. Olivine crystals have variable but on average low NiO (0.10–0.23 wt.%) and high CaO (0.22–0.35 wt.%) contents for their forsterite number. The CaO content of Fo_89–91 olivine is negatively correlated with the Al₂O₃ content of enclosed spinel (9–15 wt.%) and positively correlated with the CaO/Al₂O₃ ratios of melt inclusions (0.9–1.5). Major and trace element patterns of melt inclusions are similar to that of the host rock, indicating that the magma could have formed by accumulation of small batches of melt, with compositions similar to the melt inclusions. The liquidus temperature of the magma was  1275 °C, and its oxygen fugacity ≤ FMQ + 2.5. Correlations between K₂O, Zr, Th and LREE in the melt inclusions are interpreted to reflect variable degrees of melting of the source; correlations between Al₂O₃, Na₂O, Y and HREE are influenced by variations in the mineralogy of the source. The melts probably formed from a water-poor, clinopyroxene-rich mantle source.     

Investigating landslides with space-borne Synthetic Aperture Radar (SAR) interferometry

This paper is addressed to readers without advanced knowledge of remote sensing. It illustrates some current and potential uses of satellite Synthetic Aperture Radar interferometry (InSAR) for landslide assessment. Data acquired by SAR systems can provide 3D terrain models and be used to assist in regional scale investigations, e.g. aimed at evaluation of susceptibility of slopes to failure. Under favourable environmental conditions, the innovative Permanent Scatterers (PS) technique, which overcomes several limitations of conventional SAR differential interferometry (DInSAR) applications in landslide studies, is suitable for monitoring slope deformations with millimetric precision. The PS technique combines the wide-area coverage typical of satellite imagery with the capability of providing displacement data relative to individual image pixels. With the currently available radar satellites, however, only very slow ground surface displacements can be reliably detected and measured. The presented case study of a landslide from the Liechtenstein Alps indicates that the most attractive and reliable contribution provided by this remote sensing technique lies in the possibility of (i.) wide-area qualitative distinction between stable and unstable areas and (ii.) qualitative (relative) hazard zonation of large, slow landslides based on the identification of segments characterised by different movement rates. Since only the radar line of sight projection of the displacements can be detected, a quantitative exploitation of the PS data is possible only where sufficient ground truth is available. In site specific or single landslide investigations the PS data can represent a very useful complementary data source with respect to the information acquired through ground based observations and in situ surveying. However, the difficulties associated with the feasibility assessments of the applicability of SAR data to local scale problems, as well as with the interpretation of PS results, require a close collaboration between landslide experts and specialists in advanced processing of radar satellite data. The interpretation of the exact geotechnical significance of small, radar sensed ground surface deformations is challenging, especially where ground truth is lacking. Although any ground deformation is potentially of interest to an engineering geologist, detection of movements in both vertical and horizontal directions is needed in the case of landslides to evaluate slope failure mechanisms. With their high radar viewing angles, however, the current space-borne systems can detect only a fraction of the horizontal component of movement. It is expected that the upcoming SAR dedicated missions with new sensors and different acquisition geometries, combined with the rapid developments in the field of advanced radar data processing, will allow a full 3D reconstruction of deformation data and help to further reduce the current limitations of the PS and similar DInSAR approaches.     

Preconsolidation stress in the Vega Baja and Media areas of the River Segura (SE Spain): Causes and relationship with piezometric level changes

Preconsolidation stress (σ′_p) is the maximum effective stress that a soil has suffered throughout its life. From a geotechnical point of view, preconsolidation stress has a great importance because it separates elastic and reversible deformations from inelastic and only partially irreversible deformations and marks the starting point of high compressibility. This study calculates the preconsolidation stress for 139 undisturbed soil samples from the Vega Baja and Media of the Segura river (SE Spain), using the uniaxial consolidation test and applying the method proposed by Casagrande while using a novel analytical procedure proposed by Gregory et al. [Gregory, A.S., Whalley, W.R., Watts, C.W., Bird, N.R.A., Hallet, P.D., Whitmore, A.P., 2006. Calculation of the compression index and precompression stress from soil compression test data. Soil and Till. Res., 89, 45–57] to avoid subjective interpretations of maximum curvature point. The results show overconsolidation ratio (OCR — the ratio of preconsolidation stress to current natural overburden stress) values for the 10–15 m depth of soil varying from 2 to 14 and maximum preconsolidation stresses above 800 kPa. The main causes of calculated preconsolidation identified are desiccation due to seasonal drying and wetting cycles that have induced additional stresses always lower than 42 kPa for the more superficial samples. Water level decline due to the reduction of recharge suffered by the aquifer system during periods of drought and the uncontrolled withdrawal of water is considered to be the second cause of anomalous OCR values. This second cause induces low stresses to the more superficial layers (lower than 41 kPa) that can reach values higher than 150 kPa for the deeper layers for known water level decreases. In consequence, the soils of the Vega Baja and Media of the Segura river are highly overconsolidated for the first 5 m, decreasing gradually with depth to 10–15 m deep. For samples located deeper than 15 m the soils seem to be underconsolidated, probably due to the existence of confined aquifers that cause deviations from a hydrostatic and linear pore pressure model. This fact has a huge practical significance which implies that deformations affecting superficial layers are lower than those expected for deeper layers for the same load.     

D-INSAR技术在大范围滑坡监测中的应用

对大面积自然山区进行高精度的滑坡灾害监测一直是我国及其它地质灾害频发国家的一个重要课题,尚未得到很好解决。D-INSAR技术是近年来微波遥感发展的一个重要方向,理论上可以以厘米甚至更高的精度对连续大面积的微小地面活动进行监测,被认为是极具应用前景的滑坡灾害监测技术。为此以应用D-INSAR技术对大范围的自然滑坡进行监测和预警为目的,详细地推导了D-INSAR技术监测形变的原理公式并介绍了方法和数据处理流程,在全面回顾该技术在国内外滑坡监测中的应用现状和实例的基础上,分析了D-INSAR技术在滑坡监测中的优势及问题点,并结合最新技术进展提出了解决方案,认为应该综合利用先进地球探测及数字处理技术与D-INSAR技术相结合来克服D-INSAR技术中的问题点。具体介绍了永久散射体（PS）技术、小基线子集（SBAS）技术、GPS技术及基于地面的LISA技术等,并对这几种技术与D-INSAR技术的结合进行了探讨,使其达到对滑坡更有效的监测。     

The origin of basic rocks of the Korosten AMCG complex, Ukrainian shield: Implication of Nd and Sr isotope data

The Korosten complex is a Paleoproterozoic gabbro–anorthosite–rapakivi granite intrusion which was emplaced over a protracted time interval — 1800–1737 Ma. The complex occupies an area of about 12 000 km² in the north-western region of the Ukrainian shield. About 18% of this area is occupied by various mafic rocks (gabbro, leucogabbro, anorthosite) that comprise five rock suites: early anorthositic A₁ (1800–1780 Ma), main anorthositic A₂ (1760 Ma), early gabbroic G₃ (between 1760 and 1758 Ma), late gabbroic G₄ (1758 Ma), and a suite of dykes D₅ (before 1737 Ma). In order to examine the relationships between the various intrusions and to assess possible magmatic sources, Nd and Sr isotopic composition in mafic whole-rock samples were measured. New Sr and Nd isotope measurements combined with literature data for the mafic rocks of the Korosten complex are consistent and enable construction of Rb–Sr and Sm–Nd isochronous regressions that yield the following ages: 1870 ± 310 Ma (Rb–Sr) and 1721 ± 90 Ma (Sm–Nd). These ages are in agreement with those obtained by the U–Pb method on zircons and indicate that both Rb–Sr and Sm–Nd systems have remained closed since the time of crystallisation. In detail, however, measurable differences in isotopic composition of the Korosten mafic rock depending on their suite affiliation were revealed. The oldest, A₁ rocks have lower Sr (⁸⁷Sr/⁸⁶Sr₍₁₇₆₀₎ = 0.70233–0.70288) and higher Nd (εNd₍₁₇₆₀₎ = 1.6–0.9) isotopic composition. The most widespread A₂ anorthosite and leucogabbro display higher Sr and lower Nd isotopic composition: ⁸⁷Sr/⁸⁶Sr₍₁₇₆₀₎ = 0.70362, εNd₍₁₇₆₀₎ varies from 0.2 to − 0.7. The G₃ gabbro–norite has slightly lower εNd₍₁₇₆₀₎ varying from − 0.7 to − 0.9. Finally, G₄ gabbroic rocks show relatively high initial ⁸⁷Sr/⁸⁶Sr (0.70334–0.70336) and the lowest Nd isotopic composition (εNd₍₁₇₆₀₎ varies from − 0.8 to − 1.4) of any of the mafic rocks of the Korosten complex studied to date. On the basis of Sr and Nd isotopic composition we conclude that Korosten initial melts may have inherited their Nd and Sr isotopic characteristics from the lower crust created during the 2.05–1.95 Ga Osnitsk orogeny and 2.0 Ga continental flood basalt event. Indeed, εNd₍₁₇₆₀₎ values in Osnitsk rocks vary from 0.0 to − 1.9 and from 0.2 to 3.4 in flood basalts. We suggest that these rocks being drawn into the upper mantle might melt and give rise to the Korosten initial melts. ⁸⁷Sr/⁸⁶Sr₍₁₇₆₀₎ values also support this interpretation. We suggest that the Sr and Nd isotopic data currently available on mafic rocks of the Korosten complex are consistent with an origin of its primary melts by partial melting of lower crustal material due to downthrusting of the lower crust into upper mantle forced by Paleoproterozoic amalgamation of Sarmatia and Fennoscandia.     

Geochemical and Sr–Nd–Pb isotopic compositions of the Eocene Dölek and Sariçiçek Plutons, Eastern Turkey: Implications for magma interaction in the genesis of high-K calc-alkaline granitoids in a post-collision extensional setting

The major and trace elements and Sr–Nd–Pb isotopes of the host rocks and the mafic microgranular enclaves (MME) gathered from the Dölek and Sariçiçek plutons, Eastern Turkey, were studied to understand the underlying petrogenesis and geodynamic setting. The plutons were emplaced at  43 Ma at shallow depths ( 5 to 9 km) as estimated from Al-in hornblende geobarometry. The host rocks consist of a variety of rock types ranging from diorite to granite (SiO₂ = 56.98–72.67 wt.%; Mg# = 36.8–50.0) populated by MMEs of gabbroic diorite to monzodiorite in composition (SiO₂ = 53.21–60.94 wt.%; Mg# = 44.4–53.5). All the rocks show a high-K calc-alkaline differentiation trend. Chondrite-normalized REE patterns are moderately fractionated and relatively flat [(La/Yb)_N = 5.11 to 8.51]. They display small negative Eu anomalies (Eu/Eu = 0.62 to 0.88), with enrichment of LILE and depletion of HFSE. Initial Nd–Sr isotopic compositions for the host rocks are ε_Nd(43 Ma) = − 0.6 to 0.8 and mostly I_Sr = 0.70482–0.70548. The Nd model ages (T_DM) vary from 0.84 to 0.99 Ga. The Pb isotopic ratios are (²⁰⁶Pb/²⁰⁴Pb) = 18.60–18.65, (²⁰⁷Pb/²⁰⁴Pb) = 15.61–15.66 and (²⁰⁸Pb/²⁰⁴Pb) = 38.69–38.85. Compared with the host rocks, the MMEs are relatively homogeneous in isotopic composition, with I_Sr ranging from 0.70485 to 0.70517, ε_Nd(43 Ma) − 0.1 to 0.8 and with Pb isotopic ratios of (²⁰⁶Pb/²⁰⁴Pb) = 18.58–18.64, (²⁰⁷Pb/²⁰⁴Pb) = 15.60–15.66 and (²⁰⁸Pb/²⁰⁴Pb) = 38.64–38.77. The MMEs have T_DM ranging from 0.86 to 1.36 Ga. The geochemical and isotopic similarities between the MMEs and their host rocks indicate that the enclaves are of mixed origin and are most probably formed by the interaction between the lower crust- and mantle-derived magmas. All the geochemical data, in conjunction with the geodynamic evidence, suggest that a basic magma derived from an enriched subcontinental lithospheric mantle, probably triggered by the upwelling of the asthenophere, and interacted with a crustal melt that originated from the dehydration melting of the mafic lower crust at deep crustal levels. Modeling based on the Sr–Nd isotope data indicates that  77–83% of the subcontinental lithospheric mantle involved in the genesis. Consequently, the interaction process played an important role in the genesis of the hybrid granitoid bodies, which subsequently underwent a fractional crystallization process along with minor amounts of crustal assimilation, en route to the upper crustal levels generating a wide variety of rock types ranging from diorite to granite in an extensional regime.     

How to assess landslide activity and intensity with Persistent Scatterer Interferometry (PSI): the PSI-based matrix approach

We provide a step-by-step analysis and discussion of the ‘PSI-based matrix approach’, a methodology employing ground deformation velocities derived through Persistent Scatterer Interferometry (PSI) for the assessment of the state of activity and intensity of extremely to very slow landslides. Two matrices based on PSI data are designed respectively for landslides already mapped in preexisting inventories and for newly identified phenomena. Conversely, a unique intensity scale is proposed indiscriminately for both. Major influencing factors of the approach are brought to light by the application in the 14 km² area of Verbicaro, in Northern Calabria (Italy). These include lack of PSI data within the landslide boundaries, temporal coverage of the available estimates, and need of field checks as well as the operative procedures to set the activity and intensity thresholds. For the area of Verbicaro, we exploit 1992–2011 PSI data from ERS1/2 and RADARSAT1/2 satellites, projecting them along the maximum slope directions. An activity threshold of ±5 mm/year is determined by applying the average projection factor of local slopes to the PSI data precision. The intensity threshold between extremely and very slow phenomena (16 mm/year) is reduced by ~20 % to account for temporal and spatial averages being applied to attribute representative velocities to each landslide. The methodology allows assessing the state of activity and the intensity for 13 of the 24 landslides premapped in the 2007 inventory and for two newly identified phenomena. Current limitations due to characteristics and spatial coverage of PSI data are critically tackled within the discussion, jointly with respective implications.     

Late Triassic granitoids of the eastern margin of the Tibetan Plateau: Geochronology, petrogenesis and implications for tectonic evolution

Late Triassic granitoids in the Songpan-Garzê Fold Belt (SGFB), on the eastern margin of the Tibetan Plateau, formed at 230 to 220 Ma and can be divided into two groups. Group 1 are high-K calc-alkaline rocks with adakitic affinities (K-adakites), with Sr > 400 ppm, Y < 11 ppm, strongly fractionated REE patterns ((La/Yb)_N = 32–105) and high K₂O/Na₂O (≈ 1). Group 2 are ordinary high-K calc-alkaline I-types with lower Sr (< 400 ppm), higher Y (> 18 ppm) and weakly fractionated REE patterns ((La/Yb)_N < 20). Rocks of both groups have similar negative Eu anomalies (Eu/Eu = 0.50 to 0.94) and initial ⁸⁷Sr/⁸⁶Sr (0.70528 to 0.71086), but group 1 rocks have higher ε_Nd(t) (− 1.01 to − 4.84) than group 2 (− 3.11 to − 6.71). Calculated initial Pb isotope ratios for both groups are: ²⁰⁶Pb/²⁰⁴Pb = 18.343 to 18.627, ²⁰⁷Pb/²⁰⁴Pb = 15.610 to 15.705 and ²⁰⁸Pb/²⁰⁴Pb = 38.269 to 3759. Group 1 magmas were derived through partial melting of thickened and then delaminated TTG-type, eclogitic lower crust, with some contribution from juvenile enriched mantle melts. Group 2 magmas were generated by partial melting of shallower lower crustal rocks. The inferred magma sources of both groups suggest that the basement of the SGFB was similar to the exposed Kangding Complex, and that the SGFB was formed in a similar manner to the South China basement. Here, passive margin crust was greatly thickened and then delaminated, all within a very short time interval ( 20 Myr). Such post-collisional crustal thickening could be the tectonic setting for the generation of many adakitic magmas, especially where there is no spatial and temporal association with subduction.     

Basalt magma genesis and fractionation in collision- and extension-related provinces: A comparison between eastern, central and western Anatolia

Experimental studies of synthetic and natural basalt systems suggest that conditions of magma genesis and fractionation depend fundamentally on mantle temperatures and lithospheric stress fields. In general, compressional settings are more conducive to polybaric fractionation than extensional settings and in this regard, the Anatolian magmatic province offers a natural laboratory for comparing near-coeval basalt eruptions as a function of regional tectonics — compressional (collision-related) régimes dominating in eastern Anatolia and extensional tectonics characterizing a western province related to Aegean Sea opening. Projection of Plio-Quaternary basalt normative compositions from the Western Anatolia Extensional Province (WAEP), the Central Anatolian ‘Ova’ Province (CAOP), and Eastern Anatolia Compressional Province (EACP) are projected onto Ol–Ne–Cpx and Pl–Cpx–Ol planes in the simplified basalt system (Ne–Cpx–Ol–Qz), each showing distinctive liquid lines of descent. WAEP basalts are mostly constrained by low-pressure (< 0.5 GPa) cotectics while CAOP and EACP compositions conform to moderate and/or high-pressure (0.8–3.0 GPa) cotectics. Overall, a quasi-linear shift from moderate and/or high-pressure to low-pressure equilibria matches the westward transgression from compressional east Anatolia to the extensional west Anatolian–Aegean region. Comparison of their respective primary (mantle-equilibrated) magmas–simulated by normalizing their compositions to MgO = 15 wt.% (Mg-15)–with parameterized anhydrous and H₂O-undersaturated experimental melts suggests they segregated from spinel- to garnet-lherzolite mantle facies at pressures between c. 2 and 3 GPa (c. 70–100 km depth) under H₂O-undersaturated conditions. Interpolated potential temperatures (T_p) and lithospheric stretching factors (β) range as follows: (1) eastern Anatolian basalts associated with the Arabian foreland show T_p varying between 1250 and 1400 °C (except for the Karacalidag alkali basalts, south of the Bitlis–Zagros fracture zone, for which T_p ranges up to 1450 °C), for β values of 1.2–1.8. T_p values for central Anatolia (e.g. Sivas) range between 1300 and 1375 °C (except for Karapinar, Egrikuyu and Hasandag, which show < 1150 °C), and β values of 1.3–1.4. For western Anatolian basalts, T_p range mostly between 1250 and 1330 °C, except for a single value for Canakkale of 1400 °C and Kula sample showing T_p < 1200 °C, and β values of 1.3–2.0. Variation of these conditions is as great or greater than that between provinces, although there are clearly significant constraints on the inferred polybaric to low-pressure isobaric fractionation régimes. Covariation of total FeO, TiO₂, La/Yb, Ce/Sm, Zr/Y and Zr/Nb reflects small but significant differences in bulk composition and ambient melt fraction while the covariance of Ce/Sm and Sm/Yb is consistent with the segregation of primitive melts at the spinel- to garnet-lherzolite transition.