Search

Chlorine 36 (36Cl, T1/2 = 301,000 years) is a radionuclide with natural and anthropogenic origin that can be rejected during decommissioning of nuclear power plants or chronically during recycling of used nuclear fuels. Once emitted into the atmosphere, chlorine 36 (gas and particles) can be transferred to the soil and vegetation cover by dry and wet deposition. However, quantitative constraints of these deposits are very scarce. Because of its relatively high mobility in the geosphere and its high bioavailability, chlorine 36 fate in the environment should be studied for environmental and human impact assessments.The aim of this study is therefore to develop dry depositional models for the gaseous and particulate fractions of chlorine 36. The model used for the parameterization of gaseous chlorine 36 dry deposition on grass is the << Big-leaf >> model based on the electrical analogy (Seinfield 1985). It was adapted from the Bah (2020) model developed for iodine. For particulate chlorine 36, the Damay-Pellerin (2017) model was used. This model requiring the diameter of particulate chlorine 36 to be known, a sample of aerosol was taken in the chlorine 36 plume to determine on which particle size it is fixed. The sampling was performed on a low pressure impinger (LPI, DEKATI) of 13 levels for particles with diameters between 30nm and 10μm. In order to obtain model input data, meteorological and micrometeorological parameters were measured continuously at the IRSN La-Hague technical platform (PTILH, France). The PTILH is located 2 km north of Orano La-Hague plant which emits small quantities of chlorine 36. An ultrasonic anemometer (Young 81000V) installed at a height of 4,5 m from the ground measured the direction (°), the velocity u (m/s) and the air friction velocity u* (m/s), the sensible heat flux H, the Monin-Obukhov length (L) and the atmospheric stability (1/L). A weather station (Spectrum Watchdog, series 2000) measured temperature Ts (°C), relative humidity RH (%), dew point (°C), global radiation SR (Wat/m2) and photosynthetic active radiation PAR (μM/m2.s). Sampling campaigns of 2 weeks were also conducted at the PTILH site to determine experimental depositional rates of chlorine 36. Chlorine 36 measurements were carried out by acceleration mass spectrometry at CEREGE-LN2C (France).The particle size distribution of the aerosol sample from chlorine 36 plume shows two peaks, a main one at 2μm and a second one at 10μm, typical distribution of marine aerosol. The models yield chlorine 36 depositional velocities between 6,7.10-3 and 10-2 m/s for the particulate fraction, and between 5.10-3 and 1,3.10-2 m/s for the gaseous one. The total dry depositional velocities (particles and gas) calculated from the model are less than one order of magnitude than the ones obtained experimentally.

The Mw 6.4 Petrinja earthquake (2020, Croatia) is among the strongest continental earthquakes that occurred in Eastern Europe for decades. In such low-strain contexts, the sparse terrestrial-monitoring (few seismic and geodetic stations) of rare but strong earthquakes often prevents a detailed analysis of their seismic source. Here, we take advantage of > 160 geodetic benchmarks and optical image correlation to obtain a dense near-field coverage of the coseismic surface displacements. The geodetic dataset is obtained by repeated measurements of benchmark networks designed for civilian purposes and constitutes a unique dataset of coseismic displacements in the near-field of the fault. The optical image correlation is based on pre-earthquake (December 2017) WorldView and post-earthquake (February 2021) Pleiades satellite images with a 50 cm resolution. We also complete these displacement fields with unwrapped coseismic interferograms based on Sentinel-1 products, except in the near field affected by decorrelation. These displacement fields are consistent and thus suitable for modeling the slip distribution of the Petrinja earthquake. The elastic inversion of the geodetic benchmarks revealed interesting characteristics of this event: the rupture occurred on a near-vertical strike-slip fault, at a shallow depth (< 10 km), with significant slip reaching the surface. It also suggests that the deformation was partly accommodated by a subparallel strand 2.5 km from the main source northward. The aim of this research is to improve the source model of Petrinja 2020earthquake sequence, with a joint inversion of the geodetic benchmarks, optical image correlation, and InSAR data. Nevertheless, the comparison of the geodetic and coseismic offsets measurement on the field, shows that > 70% of the slip is likely distributed at the surface. Moreover, the coseismic strain maps derived from the unique benchmark data set helped us to identify zones where deformation appears distributed. Finally, the new data raises questions about whether such moderate earthquakes are accompanied by subsurface off-fault deformation or residual elastic strain.

Measuring 36Cl cosmogenic nuclides on exposed bedrock fault scarps has now been used in several places in the Mediterranean to retrieve ages of the fault seismic history (e.g. Mechernich et al. 2022, Iezzi et al. 2021 and Cowie et al. 2017).In Central Apennines, around the Fucino basin, at least 15 36Cl sampling sites were analyzed in previous studies to interpret the 36Cl data as seismic history or slip-rates. Several codes (e.g., Beck et al. 2018, Shlagenhauf et al. 2010) were used as a basis for solving 36Cl production equations to calculate the 36Cl concentration resulting from bedrock scarp exhumation history. Some codes included an MCMC routine to retrieve the seismic histories the closest to the dataset. The main differences between the various codes lie in: 1-the fault history prior to exhumation, 2-the parameters previous authors decided to inverse (as an example, mean density of the colluvium is inversed in Beck et al. 2018 but not in Tesson et al. 2019) and 3-the a priori distribution of those parameters (for instance, the time between two earthquakes follows an inverse gaussian distribution for Beck et al. 2018 but a uniform distribution for Tesson et al. 2019). I have compared the various codes and run them on the same dataset (one site at Campo Felice, one site at Roccapreturo and one site at Magnola) and found that retrieved seismic histories are similar, although the estimation of uncertainties differs.Moreover, all previous cited codes run under Matlab or Fortran. Fortran codes have the advantage of fast computing time but could be cumbersomeI here propose a new code, adapted from Tesson et al. 2019, in the more accessible and widely used Python language. The inferred pre-exposure is also inversed and is a function of the height of the fault cumulative escarpment. The parameters considered are the number of events, ages of event, the associated slips, the long term slip rate, the quiescence and the pre-exposure and their optimal evaluation is done with a MCMC algorithm provided by Numpyro (Du Phan et al. 2019).Using this new code, we have reanalyzed the 15 36Cl sites around the Fucino and, through a gaussian mixture algorithm, checked the hypothesis of common periods of activity throughout all the Fucino basin.

Slope failures and deep seated landslides are usually considered as the most efficient processes for hillslope erosion in active orogens. Erosion in the central Himalaya (Nepal) confirms such assertion (Morin et al., 2018), with in addition the probable predominance of the very large landslides in the erosive budget of the range (Marc et al., 2019).In this contribution, we report geological evidence for a giant rockslide that occurred around 1190 AD (14C burial age, 36Cl exposure age, and IRSL dating) in the Annapurna massif (central Nepal), and filled a wide glacial cirque (Sabche cirque) and the Seti valley farther downstream with a rock-avalanche deposit up to 1km thick, made of finely-crushed breccia. This rockslide, which involves a total rock volume of ~23 km3, decapitated the paleo-Annapurna IV, a paleo-summit culminating likely above 8000 m of altitude.Such giant rockslides are rare but not uncommon in central Himalaya with former description of pluri-kilometric rockslides involving up to ten cubic kilometres mass wasting (e.g. Weidinger et al., 2002). They have major implications for landscape evolution since they could represent the main mode of erosion of the high glaciated peaks, leading to the sudden reduction of ridge crest elevation by several hundred meters.They also have major implications on the fluvial network and how it responds to such massive and sudden supply of sediments. Downstream of the Seti river, the Pokhara Basin is filled by ~5 km3 of mostly conglomeratic sediment emplaced by fluvial, hyperconcentrated, or turbulent, sediment-laden flows. Extensive 14C dating of organic fragments found in the fine-grained units of the Pokhara sediments (Schwanghart et al., 2016; this study) provide robust constraints on the timing of aggradation. Onset of aggradation around 1200 AD (as well as the calcareous clasts whose only possible origin is the Sabche cirque) indicates that the conglomerates are the result of the active erosion of the Sabche rock-avalanche breccia. The aggradation lasted for about a century at an average rate of 1m/yr. During that period, the erosion of the fine-grained breccia has been extremely efficient and rapid: of the initial 27 km3 of rockslide debris, less than 10% can still be observed today in the Sabche cirque. During that period, the rate of sediment yield delivered by the Seti river was comparable to the highest measured values of post-volcanic eruption sediment transport. It led to an overwhelming content of carbonate-rich material in Narayani river sediments up to 300 km downstream, in the Narayani megafan in the middle of the Gangetic plain. Compared to the Narayani basin (30,000 km2), which exports annually ~0.05km3/yr of sediment, post-collapse erosion of the Sabche breccias (only 0.2% of the Narayani total basin area) would have increased this annual flux by a factor of 3 over 100-150 years.This particular example illustrates how the erosional/sediment routing system in active mountain range, as well as the landscapes, can be durably affected by one single extreme event.

Europe has experienced over the last years earthquakes of moderate magnitude (Mw 5-6), yet destructive, reminding us of the seismogenic potential of slowly deforming regions. Among them, the 2020 Mw 6.4 Petrinja earthquake ruptured the Petrinja-Pokupsko Fault (PPKF) in Central Croatia, about 50-km southeast of Zagreb, a region in which the caracterisation of seismogenic faults had been insufficiently studied before that event. Understanding the strain accommodation through time and space is critical for accurate assessment of the regional seismic hazard.Using field observations and high-resolution topographical data derived from airborne LiDAR (~10 cm resolution) and tri-stereo satellite images (Pléiades, resolution 50 cm), we accurately mapped the fault trace, underlined at several sites by geomorphic markers such as valleys, terrace risers, and alluvial fans that have recorded cumulative displacements ranging from 5 to > 50 m and potentially up to ~180 m. Along the studied section, our fault mapping is composed of a clear NW-SE-trending 10-km-long strand between Donja and Cepelis, and of 1-4-km-long right-stepping segments marked by a non-negligible vertical component. The southern strand is composed of 2-3 sub-parallel segments that accommodate the deformation within a < 500 m wide fault zone.We have identified several sites on the main southern strand where offsets have been accurately measured and where displaced markers have been sampled for cosmogenic nuclide exposure dating and radiocarbon datings. This will allow to estimate the slip-rate for this fault at different sites and over several time spans.The mapped fault appears very discontinuous with the deformation absorbed by a series of small fault sections rather than on a single fault strand. This likely reflects a recent transpressive deformation, with immature faults, in agreement with the source parameter of the 2020 Petrinja earthquake derived from seismology.Finaly, the 2020 coseismic surface ruptures affected the northern section of the PPKF, while the mapped cumulative displacements appears more prominent along the southern section. A better knowledge of the seismic history of this entire fault system is thus crucial for seismic hazard assessment of this area.

Long-term fault escarpments are built by the accumulation of individual earthquakes producing incremental surface displacements on the fault releasing crustal tectonic loading. Cumulative escarpment studies have revealed a spatial slip variability along active faults as well as a temporal variability with the alternation of phases of intense seismic activity over a short period of time followed by long periods of quiescence. Understanding this spatial and temporal slip variability on individual faults and over a complex fault system provide a better knowledge of co-seismic rupture extents, essential for estimating past earthquakes magnitude and for seismic hazard assessment.Up to now, most studies have focused on a timeframe over few seismic cycles, making it difficult to apprehend the rupture barriers persistence and cumulative slip distribution. Here, we aim at quantifying the slip variability over several timescales ranging from a few months to a few million years on the same fault.Our study focusses on the ~50 km-long Liri fault, SW of the Fucino basin. The fault is located at the contact between Cretaceous limestone and patches of Quaternary deposits locally convering Mio-Pliocene flysch sediments. Detailed mapping of the fault trace on high-resolution Digital Elevation Model (DEM) from UAV-acquired images, Pleiades images and Lidar together with field observations revealed changes in the morphological expression of the fault north and south of an important wind gap located at Capistrello. To the north, the faut trace is ~16 km-long located on the eastern side of ~2km-wide limestone ridge, reaching ~1300m asl elevation. Two bends in the fault trace, made of ~5km long segments, can be observed with the fault strike varying between N115° and N140°. In this northern section, the fault scarp appears subtle and we did not observe Quaternary deposits on the hanging wall. In the 30 km-long section, south of Capistrello, the cumulative scarp composed of numerous splays is evidenced by a sharp trace, offsetting several morphological surfaces and associated Quaternary sediment packages. Three major bends are observed in this section of the fault, separating 10 to 30 km-long segments striking between N110° and N160°. An alluvial surface offset by ~14 m of cumulative displacement was dated at ~35kyr using 36Cl cosmogenic exposure dating suggesting a minimum slip rate of 0.4 mm/yr. Other morphological markers that have accumulated displacement between ~10 and 70 m-high have also been sampled for 36Cl cosmogenic exposure dating. Moreover, we excavated two small trenches at the base of the fault scarp within the Quaternary deposits affected by the fault revealing 3 rupture-surfacing earthquakes over the last 2500 yr, the last one occurring after 1226 CE. We will present those results and will discuss how the displacement varies along the fault both in time and space.

Fault displacement can be a source of hazard for critical infrastructures located in the nearby of a capable fault. This issue is usually addressed with zonation and avoidance strategies, but sometime the facilities have not this option. An alternative approach to assess likelihood of exceeding a certain level of displacement for pre-existing infrastructures is the Probabilistic Fault Displacement Hazard Analysis. Different empirical approaches have been proposed since the early 2000s to assess the probability of occurrence and the probability of exceedance of certain values of displacement, for both Primary and Distributed faulting, starting from the fault parameters. We propose the methodological approach used to gain the needed parameters and the results of the PFDHA applied to the Anghiari Fault, a poorly constrained NE-dipping segmented normal fault located in the Upper Tiber Valley (Italy) and belonging to the well-known Altotiberina low-angle normal fault system.In order to constrain the fault geometry and to select sites suitable for paleoseismologic trenching we performed geological survey, morphotectonic analysis and geophysical investigations. To assess the capability of the fault and its rate of activity we carried out a paleoseimic campaign, investigating several segments of the Anghiari fault. To obtain a multiscale evaluation of the fault slip rate, we collected samples to date paleosurfaces displaced by the fault with the cosmogenic nuclides methodology. At the end we performed the PFDHA obtaining curves and maps of hazard for both primary and distributed faulting, managing the uncertainties through various rupture scenario involving different fault segment.

The Italian Apennines is a slowly deforming area, despite not properly being an intraplate region. This is particularly true for the Northern Apennines, where~1ka) recurrence intervals. This makes earthquake geology a fundamental tool for characterizing the seismic hazard.The Anghiari fault is a 11 km-long segmented NE-dipping normal fault bounding the western side of the Upper Tiber Valley (Northern Apennines, Italy), and belonging to the well-known Altotiberina low-angle normal fault system. Here, we provide unprecedented evidence of the Holocene activity of the Anghiari fault through geological, geophysical and palaeoseismological investigations. The fault is composed of at least two nearly parallel splays. One splay runs at the base of the Pleistocene Anghiari ridge, downfaulting the late Quaternary alluvial deposits of the Tiber Valley against Middle Pleistocene continental deposits. The other splay is located within the Middle Pleistocene units of the Anghiari ridge. We focus on the latter. Detailed geomorphological analysis, geological mapping and near-surface geophysics, enabled us to select two sites for palaeoseismological trenching. Radiocarbon dating of faulted sediments provides constraints for late Holocene and historical surface faulting events significantly contributing to the estimation of the seismic hazard in the region.

Chlorine-36 (36Cl) is currently the only in situ cosmogenic nuclide applicable in carbonates, Ca- and K-rich feldspars and aphyric silicate rocks. Because the production reactions of 36Cl are more numerous and complex than those of other cosmogenic nuclides (e.g. 10Be, 3He), comprehensive and user-friendly calculators are essential for routine application of 36Cl to Earth surface research questions. However, the existing, commonly used calculators are affected by several drawbacks, such as oversimplification of the muon production model, or the absence of choice regarding production rate parameters. Here, we present the novel 36Cl exposure age calculation tool in development, “CREp 36Cl dev”, which we have implemented as a new functionality in the existing cosmic-ray-exposure program CREp for 10Be and 3He exposure ages. Taking advantage of the most recently developed scaling, muon and low-energy neutron models, “CREp 36Cl dev” allows calculating surface exposure ages, providing several choices between scaling and muon models and production parameters. “CREp 36Cl dev” can be accessed at https://crep-dev.otelo.univ-lorraine.fr/#/.

Active faults play a major role in relief building, partly through the accumulation of vertical co seismic displacement during major earthquakes. Triangular facets are geomorphic features recording normal fault cumulative displacements on relatively long time scales (10-100ka). To unravel the relationships between the rate at which slip accumulates on a fault scarp and the long-term evolution of triangular facets, we have to acquire quantitative datasets on normal fault slip rates at various timescales and rates of erosion of the facets.Here we present a study on facet build-up over 10-100 ka time range in the central Apennines in Italy. The normal fault systems that control the present tectonic activity of the range are very well studied with numerous detailed paleoseismological records. We focus on the Magnola-Velino fault system which displays well preserved triangular facets and accurate chronological constraints on the 10-20m high fault scarp located at the base of the facets. We combine high resolution morphometric analysis (gullies steepness, facets slope and others), using Digital Elevation Models derived from Pléiades imagery and a new dataset for cosmogenic nuclides concentration (36Cl) including 54 bedrock samples on 9 gullies and facets profiles above the scarps. Magnola-Velino fault system is 20-25km long, and morphometric parameters such as steepness index display a systematic evolution along strike. First measured 36Cl concentrations, on Magnola, range from 6 to 50 x104 at/g on gullies and 50 to 150 x104 at/g on facets, with systematic variations along strike. We document the joint evolution of geochronological data and morphological parameters for this fault system and compare it with existing constraints on long-term slip rate.

We present a database of field data for active faults in the central Apennines, Italy, including trace, fault and main fault locations with activity and location certainties, and slip-rate, slip-vector and surface geometry data. As advances occur in our capability to create more detailed fault-based hazard models, depending on the availability of primary data and observations, it is desirable that such data can be organized in a way that is easily understood and incorporated into present and future models. The database structure presented herein aims to assist this process. We recommend stating what observations have led to different location and activity certainty and presenting slip-rate data with point location coordinates of where the data were collected with the time periods over which they were calculated. Such data reporting allows more complete uncertainty analyses in hazard and risk modelling. The data and maps are available as kmz, kml, and geopackage files with the data presented in spreadsheet files and the map coordinates as txt files. The files are available at: 10.1594/PANGAEA.922582., Measurement(s) active geological fault Technology Type(s) digital curation Factor Type(s) location Sample Characteristic - Environment geological fault Sample Characteristic - Location Central Apennines Machine-accessible metadata file describing the reported data: 10.6084/m9.figshare.14102021

Long-term landscape evolution is controlled by tectonic and climatic forcing acting through surface processes. Rivers are the main drivers of continental denudation because they set the base level of most hillslopes. The mechanisms of fluvial incision are thus a key focus in geomorphological research and require accurate representation and models. River incision is often modeled with a stream power model (SPM) based on the along-stream evolution of drainage area and channel elevation gradient but can also incorporate more complex processes such as threshold effects and statistical discharge distributions, which are fundamental features of river dynamics. Despite their importance in quantitative geomorphology, such model formulations have been confronted with field data only in a limited number of cases. Here we investigate the behavior of stochastic-threshold incision models across the southeastern margin of the French Massif Central, which is characterized by significant relief and the regular occurrence of high-discharge events. Our study is based on a new dataset combining measurements of discharge variability from gauging stations, denudation rates from 34 basins from 10Be cosmogenic radionuclide (CRN) concentration measurements in river sediments, morphometric analysis of river long profiles, and field observations. This new dataset is used for a systematic investigation of various formulations of the SPM and to discuss the importance of incision thresholds. Denudation rates across the SE margin of the Massif Central are in the 20–120 mm kyr−1 (equivalent to mm/ka in the figures) range, and they positively correlate with slope and precipitation. However, the relationship with the steepness index is complex and supports the importance of taking into account spatial variations in parameters (D50, discharge variability k, runoff) controlling the SPM. Overall, the range of denudation rate across the margin can mainly be explained using a simple version of the SPM accounting for spatially heterogeneous runoff. More complex formulations including stochastic discharge and incision thresholds yield poorer performances unless the spatial variations in bedload characteristics controlling incision thresholds are taken into account. Our results highlight the importance of the hypotheses used for such a threshold in SPM application to field studies and notably the impact of actual constraints on bedload size.

(IF 4.45: Q1); International audience; The emerged coral reef terrace sequence at Cape Laundi, on the north coast of Sumba Island (Indonesia), with at least 18 successive strandlines, remains poorly dated in spite of numerous previous data. The age discrepancies within these coral reef terraces (CRTs) were previously explained by their polycyclic nature, triggered by marine erosion and reoccupation of old coral colonies by new ones. This study aims at highlighting these processes, as well as the continental denudation that participates in the partial stripping of the thin superficial coral reef layer overlying the pre-existing surface, exhuming older coral colonies. For this purpose, we use a combined analysis of 36Cl cosmogenic concentrations, new 230Th/U ages, and previous dating in order to quantify denudation rates affecting the sequence and to highlight the role of marine erosion in reworking the lowest CRT surface. Our results demonstrate that 1) the lowermost CRT is composite, i.e., formed by different reefal limestone units constructed and eroded during successive highstands of the last interglacial, 2) following the last deglaciation, this CRT has been subjected again to coastal erosion and reoccupation during the Mid Holocene highstand, 3) its distal edge is affected by the current marine erosion and shows denudation rates higher by one to two orders of magnitude (from 279 ± 0.4 to 581 ± 0.4 mm ka−1) than the continental denudation values of higher CRTs (14.7 ± 8.3 mm ka−1 on average), 4) at the scale of a single CRT surface, variations in continental denudation rates are caused by epikarstification roughness, and 5) the distal edges have the highest continental denudation rate due to diffusion and regressive erosion produced by the runoff occurring along the steep downward cliff.

International audience; The footwall of the surface rupturing Paganica normal fault, the source of the 2009 L'Aquila earthquake (Mw 6.1) in the Central Apennines (Italy), was investigated using integrated geological and geomorphological approaches. The aim was to constrain the active tectonics by studying the Raiale River that orthogonally crosscuts the fault trace, where it provides a useful geomorphological marker of long-term fluvial incision and footwall uplift. Using morphostratigraphy and paleomagnetic analysis, the Plio-Pleistocene morphotectonic evolution of the area was reconstructed, comprising an ancient continental basin and paleolandforms that predate the footwall incision. Starting from the Late Early Pleistocene-Middle Pleistocene, fluvial dissection was mainly due to marked river downcutting triggered by significant activity of the Paganica Fault, which caused progressive base-level lowering. The Raiale River downcutting formed five Middle-Late Pleistocene fluvial terraces, that, along with absolute Optically Stimulated Luminescence (OSL) dating, allowed the identification of paleolongitudinal profiles with a diverging downstream configuration. Terrace dating yielded a minimum incision rate of 0.25 ± 0.02 mm/a, which only partially compensates the footwall uplift and can thus be considered as a minimum value for the Paganica Fault throw rate, which could reach up to ~0.45 mm/a. In parallel, using terrestrial cosmogenic nuclides, a denudation rate of 0.02-0.04 mm/a was measured on the summit of the footwall block. This denudation is in keeping with the drainage incision, suggesting a non-steady state for the fault footwall topography and a dominance of relief growth. Last, the analysis of the modern Raiale River longitudinal profile denoted an ungraded status, with two main knickzones that we interpret as transient forms due to tectonic perturbations, likely triggered by activity of the Paganica Fault during the end Early Pleistocene and the Late Pleistocene. Considering the 2009 L'Aquila earthquake coseismic rupture, we observe that the younger transience on the Raiale River longitudinal profile, if it is of tectonic origin, could have only been produced by much larger seismic events (i.e., Mw > 6.5) than those documented in the area by paleoseismological investigations. The collective results confirmed that in the Central Apennines, conditions of dynamic equilibrium are often not met, and that the persistence of transient perturbations induced by tectonics should be accounted for.

International audience

International audience

The Petrinja-Pokupsko fault-system is a NW-SE right-lateral fault system that ruptured during the 29 December 2020 Mw 6.4 earthquake (~40km south-east of Zagreb, Croatia). Field analysis revealed opening of cracks and offsets of several centimeters (3 to 40 cm) along a ~20 km long fault zone extending from the Kupa river (in the northwest) to the Petrinjčica river (in the southeast). Optical image correlation based on WorldView satellite images has been used to document the first-order near-field rupture signal. The pre-event image was acquired on 7th December 2017, and the post-event image on 15th January 2021. The first results indicate a right-lateral displacement of ≈75 cm with a small (6 earthquakes. Since the fault extends farther NW and SE, from the Vukomeričke Gorice hills to Mount Kozara (Bosnia), for a total length of about 100 km, it could generate potentially larger events. It is also noteworthy that the 2020 Petrinja event occurred only 9 months after the Zagreb March 2020 (Mw 5.3) earthquake. This event occurred on an ENE-WSW-trending thrust fault, broadly orthogonal to the right-lateral Petrinja-Pokupsko fault system, ~45 km north of Petrinja, and raises the prospect of potential interplay between strike-slip and thrust faults in moderate strain-rate intra-plate settings. To address this problem, future works will aim at constraining the geometry of this fault network and its seismogenic potential.

The evolution of continental relief results from the combined action of tectonic and climatic forcings. These processes do not act continuously but often through punctual events (earthquakes, major floods, landslides) whose integrated action over time (100 Kyr to Myr) leads to the formation of landscapes. The distribution of these extreme events is often described by statistical functions involving power-law relationships between frequency and magnitude, which, coupled with the non-linearity of the geomorphological response and threshold effects for the activation of erosion agents, leads to a complex and often poorly understood relief dynamics.Studying the influence of discharge variability helps to better constrain river incision and long-term relief evolution. The south-eastern margin of the Massif Central (France) is a very interesting target for such investigations because it presents episodes of very intense precipitation focused on the relief resulting in marked differences in the statistical discharges distributions across the landscape. Some theoretical river incision models incorporate such variability (Lague et al., 2005) but they have been confronted with real data only in a limited number of cases (DiBiase et al., 2011; Scherler et al., 2017; Campfort et al., 2020). Here we test these models in the Massif Central area and in particular on Cévennes, Ardèche and Margeride mountains by quantifying denudation rates using cosmogenic nuclides (10Be), characterizing discharges variability and performing morphological analysis on longitudinal rivers profiles.The analysis of 326 river gauging stations allow us to observe a strong gradient in discharge variability from the external SE border to the interior of the Massif Central. The 10Be concentrations measured from river sediments in 36 catchments imply a large variation of denudation rates between 29 mm/kyr and 126 mm/kyr. We compare these denudation rates with the spatial distribution of mean annual precipitations, local relief, slope and concavity index, and also integrate all the observations in the frame of a stochastic threshold incision model. Our results confirm the complex model predictions of non-linear relationships between mean denudation rates and the channel steepness index and their dependence on hydrological variability and run-off.key-words : extreme events, stochastic threshold incision model, denudation rates, discharge variability, morphometric parameters, Massif Central

The 29 December 2020, Mw 6.4 Petrinja earthquake nucleated at a depth of ~10 km in the Sisak-Moslavina County in northern Croatia, ~6 km WSW of the Petrinja town. Focal mechanisms, aftershocks distribution, and preliminary Sentinel-1 InSAR interferogram suggest that the NW-SE right-lateral strike-slip Pokupsko-Petrinja fault was the source of this event.The Croatian Geological Survey, joined by a European team of earthquake geologists from France, Slovenia and Italy, performed a prompt systematic survey of the area to map the surface effects of the earthquake. The field survey was guided by geological maps, preliminary morphotectonic mapping based on 1:5,000 topographical maps and InSAR interferogram. Locally, field mapping was aided by drone survey.We mapped unambiguous evidence of surface faulting at several sites between Župić to the NW and Hrastovica to the SE, in the central part of the Pokupsko-Petrinja fault, for a total length of ~6.5 km. This is probably a minimum length since several portions of the fault have not been explored yet, and in part crossing forbidden uncleared minefields. Surface faulting was observed on anthropic features (roads, walls) and on Quaternary sediments (soft colluvium and alluvium) and Miocene bedrock (calcarenites). The observed ruptures strike mostly NW-SE, with evidences of strike-slip right-lateral displacement and zones of extension (opening) or contraction (small pressure ridges, moletracks) atlocal bends of the rupture trace. Those ruptures are interpreted as evidences of coseismic surface faulting (primary effects) as they affect the morphology independently from the slope direction. Ground failures due to gravitational sliding and liquefaction occurrences were also observed, mapped and interpreted as secondary effects (see Amoroso et al., and Vukovski et al., this session). SE of Križ, the rupture broke a water pipeline with a right-lateral offset of several centimetres. Measured right-lateral net displacement varies from a few centimetres up to ~35 cm. A portion of the maximum measured displacement could be due to afterlisp, as it was mapped several days after the main shock. Hybrid surface ruptures (shear plus opening and liquefaction), striking SW-NE, with cm-size left-lateral strike-slip offsets were mapped on the northern side of the Petrinja town, ~3 km NE of the main fault.Overall, the rupture zone appears discontinuous. Several factors might be inferred to explain this pattern such as incomplete mapping of the rupture, inherited structural discontinuities within the Pokupsko-Petrinja fault system, or specific mechanical properties of the Neogene-Quaternary strata

International audience; Drawing upon a socio-constructivist perspective, this paper aims to gain insight into how two communities in the Italian Apennines region - differentially exposed to severe earthquakes in the past - remember, understand and plan for seismic phenomena. Avezzano was completely devastated in 1915, while the last significant event in Sulmona dates to 1706. However, both communities have indirectly experienced recent serious seismic events (e.g., L'Aquila earthquake in 2009 and Amatrice earthquake in 2016) as well as directly experienced lesser local tremblors. Interviews with citizens (N = 37) and stakeholders (N = 18) were conducted in the two cities and content-analysed with the support of NVivo software. The results show substantial differences in how respondents understand seismic events, access this information and prepare preventive strategies depending on their place of residence and role. The Avezzano community shares a more detached view, while the Sulmona community gives a more concerned picture. Moreover, the stakeholders tend to present an optimistic scenario, while citizens express a more critical viewpoint. Overall, the interview contents may be organised around two main oppositions: fatalism vs. empowerment in the face of an earthquake threat; and the will to forget about this risk vs. the need to remember it. The positioning on these two oppositions is defined by psychological, social and structural characteristics of individuals, and give rise to different representations of seismic risk and management. Such representations have practical implications for how the issue is addressed and faced in the communities' everyday life.

The earthquake with magnitude ML=6.2 that occurred on 29th December 2020 has caused significant material damage to objects and infrastructure in the towns of Petrinja, Sisak,Glina and the surrounding area. According to the satellite interferometry data, the coseismic and postseismic deformation area covers around 500 square kilometers. The existing geodetic benchmarks have been set in the affected towns, and their coordinates have been determined based on previous GPS campaigns. The GPS network was set up and adjusted at the State Geodetic Administration's request for geodetic monitoring of infrastructure and cadastral projects. These points are not primarily intended for high accuracy measurements at the level of a few millimeters, so their accuracy and the absolute shift concerning geodynamic processes in the region should be taken into account. Nevertheless, the data obtained by their observation after the earthquake can provide valuable information about the horizontal and vertical displacements with a certain level of confidence. The field survey has detected disappearance of a large number of benchmarks and some valuable information has been lost. Still, 58 points were found and observed and it has been concluded that 52 points are reliable and can be used for future research. Because the network of benchmarks is not developed in rural areas, there is a gap in the distribution of benchmarks in affected area. Therefore, the additional data was collected using the benchmarks established for the engineering and cadastral projects and studies. From a total of 67 points that have been found and observed, 42 points will be used. Along with the data collected in urban areas, there will be a total of 94 benchmarks. The accuracy of the geodetic benchmark measurements is at the centimeter level, while the values of deformation are at the level of a few decimeters. Therefore, the obtained data can be used to better assess the displacement recorded during the 29 December 2020 event. In the future, field research will focus on finding additional benchmarks to reach a better spatial distribution.

In the Central Apennines (Italy), up to now, no absolute dating directly based on the moraines has been carried out to constrain glacial oscillation. However, climatic constrains are often used in the Central Apennine to estimate long term (> 10 ka) fault slip rate. In addition slip rate assessments based on offset morphotectonic markers on the main branches of fault systems and encompassing several seismic cycles (> 10 ka) are sparse. This is particularly true for the Monte Vettore-Monte Bove fault system which triggered the 2016-2017 seismic sequence. We thus provide new assessment for the vertical slip rates along the Mt Vettore-Mt Bove fault system. Offset measurements were made using a 5-cm resolution DEM obtained through a drone survey and constrain a fault scarp height of 15.5 ± 1.4 m and a cumulative offset of 32-40.5 m. Samples were collected from the Valle Lunga terminal moraine at 1710 m asl and yield 36Cl exposure ages of 12.7 + 2.2/-1.9 ka while the flat, abraded surface located on top of the tectonic scarp yield 36Cl exposure ages of 23.4 + 5.3/-4.3 ka. Assuming the offset started to accumulate when climate conditions allow its preservation, thus once the surface was abandoned, we constrain a vertical slip rate of 1.2 ± 0.2 mm/yr along the master branch of the Mt Vettore normal fault. This rate is higher than the ones previously obtained from trenches along secondary splays of the Mt Vettore-Mt Bove and on the Norcia fault systems. Besides, the yielded chronology for the last glacial maximum in that area at ~23 ka is in good agreement with the timing previously proposed for the LGM in the Apennines.

River terraces are geomorphological markers recording deformation. Here, we use four strath and fill river terraces along the Kizilirmak River in Turkey to unravel the deformation along the convex arc formed by the central North Anatolian Fault (NAF), a continental transform fault. 10Be, 26Al, 36Cl cosmogenic exposure ages of T3 and T2 strath terraces constrain their formation at 83 ± 15 ka and 67 ± 7 ka, respectively. We evidence frost-cracking during humid cold periods that bought younger carbonate cobbles on both terraces. T1, a younger fill terrace was probably emplaced during the MIS3 (30–40 ka) aggradational period or during the last glacial–interglacial transition. T0, the most recent fill terrace, was incised shortly before the 1668 earthquake based on 14C dating. It records a cumulated 14 ± 2 m offset linked to the 1943 and 1668 earthquakes. T3 shows a maximum offset of ~ 845 m and constrains a 10 mm/year geological slip rate that is lower than Holocene slip rates. It suggests temporal change in slip rates along the NAF. T2 and T3 also evidenced an uplift of 1 mm/year induced by transpressive deformation accommodated close to the NAF. Compared to the 0.28 mm/year obtained to the north, a larger portion of the shortening in the Central Pontides is accommodated close to the driving plate boundary. We also evidenced high cosmogenic-based erosion rates in Pontides during the Holocene, and even higher rates during T3 planation. Erosion rates were combined with present-day relief to infer a first-order minimal shortening of 12–16 km using simple mass balance principles in the Central Pontides. This long-term shortening is larger than the one induced by the Anatolian rotation suggesting a far field effect of the collision.

International audience; The aim of the Fault2SHA European Seismological Commission Working Group Central Apennines laboratory is to enhance the use of geological data in fault-based seismic hazard and risk assessment and to promote synergies between data providers (earthquake geologists), end-users and decision-makers. Here we use the Fault2SHA Central Apennines Database where geologic data are provided in the form of characterized fault traces, grouped into faults and master faults, with individual slip rate estimates. The proposed methodology first derives slip rate profiles for each master fault. Master faults are then divided into distinct sections of length comparable to the seismogenic depth to allow consideration of variable slip rates along master faults and the exploration of multi-fault ruptures in the computations. The methodology further allows exploration of epistemic uncertainties documented in the database (e.g. master fault definition, slip rates) as well as additional parameters required to characterize the seismogenic potential of fault sources (e.g. 3D fault geometries). To illustrate the power of the methodology, in this paper we consider only one branch of the uncertainties affecting each step of the computation procedure. The resulting hazard and typological risk maps allow both data providers and end-users (1) to visualize the faults that threaten specific localities the most, (2) to appreciate the density of observations used for the computation of slip rate profiles, and (3) interrogate the degree of confidence on the fault parameters documented in the database (activity and location certainty). Finally, closing the loop, the methodology highlights priorities for future geological investigations in terms of where improvements in the density of data within the database would lead to the greatest decreases in epistemic uncertainties in the hazard and risk calculations. Key to this new generation of fault-based seismic hazard and risk methodology are the user-friendly open source codes provided with this publication, documenting, step-by-step, the link between the geological database and the relative contribution of each section to seismic hazard and risk at specific localities.

The interaction between sedimentation/erosion and faulting represents one of the most intriguing topics in landscape and tectonics evolution. Only few studies have been able to document the feedback between faulting and sedimentary loading from field observations. Here, we focus on how sediment loading/unloading influences the dynamics of fault systems in the Fucino basin, in the Central Apennines (Italy). The Fucino basin represents a remarkable case study with respect to the other main extensional basins in the Apennines because of its large dimension, square shape, significant sediment thickness, and its endorheic nature throughout its evolution. We present a detailed structural and geomorphologic analysis of the Fucino basin and its surroundings, investigating the kinematic and geometry of each main fault strand. The slickenlines analysis reveals multiple families of slip-vectors and timing of activity, suggesting a change in extension slip-direction from N240° to N200° during middle Pleistocene. Using a local isostatic model, we estimate that up to the 30% of the vertical geological displacement of the faults, which overall ranges from 0.5 to 2.5km, is related to the sediment loading/unloading. We demonstrate a positive feedback between sedimentation and faulting which may also lead to a reorganization in fault kinematics related to a significant increase in vertical stress. We propose a conceptual model for the permanent endorheic configuration of the Fucino basin, which includes the effect of sediment loading.

International audience; The Mw 6.5 Norcia earthquake occurred on 30 October 2016, along the Mt Vettore fault (Central Apennines, Italy), it was the largest earthquake of the 2016-2017 seismic sequence that started 2 months earlier with the Mw 6.0 Amatrice earthquake (24 August). To detect potential slow slip during the sequence, we produced Interferometric Synthetic Aperture Radar (InSAR) time series using 12- to 6-day repeat cycles of Sentinel-1A/1B images. Time series indicates that centimeter-scale surface displacements took place during the 10 weeks following the Norcia earthquake. Two areas of subsidence are detected: one in the Castelluccio basin (hanging wall of the Mt Vettore fault) and one in the southern extent of the Norcia earthquake surface rupture, near an inherited thrust. Poroelastic and viscoelastic models are unable to explain these displacements. In the Castelluccio basin, the displacement reaches 13.2 +/- 1.4 mm in the ascending line of sight on 6 January 2017. South of the Norcia earthquake surface rupture (a zone between the Norcia and Amatrice earthquakes), the postseismic surface displacements affect a smaller area but reach 35.5 +/- 1.7 mm in ascending line of sight by January 2017 and follow a logarithmic temporal decay consistent with postseismic afterslip. Our analysis suggests that the structurally complex area located south of the Norcia rupture (30 October) is characterized by a conditionally stable frictional regime. This geometrical and frictional barrier likely halted rupture propagation during the Amatrice (24 August) and Norcia (30 October) earthquakes at shallow depth (

Item S1 (additional information on the modeling) and Item S2, (additional information on the sampling, analysis, and obtained dataset).

International audience; In 2016, the Mt. Vettore-Mt. Bove normal fault system (VBFS) broke during three earthquakes (Mw 6.0, Mw 5.9, and Mw 6.5), associated with clear coseismic ruptures. Based on high-resolution topography and geological field data, we determined the displacements of the VBFS. The distributions of the coseismic and postglacial displacements exhibit similar asymmetric shapes, suggesting self-similar slip profiles over the last 18 kyr. The highest displacement during the 2016 events is localized on the southern Mt. Vettore segments, the southern tip of the VBFS, which also showed the maximum throw of 32 m over the last 18 kyr yielding to a very fast throw rate of 1.6 ± 0.5 mm/yr. Assuming a constant throw rate, the geological displacement we determined suggests an inception age of 200–250 kyr for the surface rupturing faults that ruptured in 2016 on the Mt. Vettore sector. This value can be interpreted as the minimum age for the onset of those faults in the Mt. Vettore and might result from a southward shift of VBFS activity with a change of the fault system pattern. We infer that the observed slip distribution, maximum at Mt. Vettore and tapering-off southward, might be due to the VBFS lengthening processes towards the Laga Fault and its interaction with pre-existing structure such as Olevano-Antrodoco-Sibillini thrust that might play a key role in controlling the VBFS's evolution.

Active deformation in the Central Apennines is mostly accommodated by NW-SE normal faults systems that produce moderate to large earthquakes at shallow depth. Recent examples include the 1915 Mw≈7 Avezzano earthquake (Fucino basin) and the 2009 Mw=6.1 L’Aquila earthquake (Aterno basin) which were both associated with major loss of life and massive damage to buildings and infrastructure. Here, we study the 40-km-long Ovindoli – Piano di Pezza – Campo Felice – Monti d’Ocre (OPCM) fault system, a major NNW-SSE system that potentially links the Fucino and the Aterno fault-systems. The OPCM exhibits linear and arcuate sections with four main segments and borders the eastern margin of the Aterno basin. Paleo-earthquake rupture data on the Piano di Pezza (PPF) and Campo Felice (CFF) faults exhibit remarkable synchronicity with the Fucino fault system, with the most recent surface-rupturing earthquake likely occurring in the XIVth century. In order to better understand the relationships between and earthquake rupture scenarios, we focus on the basin geometry and fault surface expression of the Piano di Pezza fault combining geomorphology and subsurface geophysics. We map the fault trace with unprecedented detail using terrestrial laser scanner surveys and quantify surface deformation affecting alluvial fans as well as glacial moraines. We obtain a mean vertical offset of 2.5 m +/- 0.3 m for the most recent features, well in agreement with paleoseismological data. Furthermore, we document slip distributions at different time scales along strike with a maximum value at the connection between the PPF and the OF. Beneath the scarp, geophysical data reveal a complex faulting geometry with several parallel strands and two minor blind splays. Electrical resistivity tomography images show a cumulative vertical offset of ~ 15 m affecting an interface attributed to the Last Glacial Maximum and confirm the high vertical slip across the fault zone. Gravimetric anomalies across the basin also indicate the sedimentary fill has recorded a maximum finite cumulative throw of the PdP fault system of 110-140 m. This suggests a maximum vertical slip rate of 0.2-0.3 mm/year since the Pleistocene, which contrasts with the high post-LGM slip rate estimated from trenches. Overall, our observations suggest that the arcuate PPF originally formed as a reverse fault during the Mio-Pliocene compressive stage and is now reactivated as an extensional horsetail-like feature by ruptures along a major strike-slip fault (OF). This finding points to the PPF as mostly built through ruptures along the OF leaking onto an inherited structure. The time-varying slip rates may also denote an episodic behavior marked by short periods of high seismic activity (a few centuries) and long intervals of seismic quiescence (a few millennia). Furthermore, possible earthquake rupture scenarios along the OPCM may encompass the whole OPCM fault system (cumulative length ca. 40 km) or rupture termination along the PPF (cumulative length ca. 15-20 km) with significantly different impacts over the populated Fucino and Aterno basins.

A driving motivator in many active tectonics studies is to learn more about the recurrence large and potentially destructive earthquakes, providing the means to assess the respective fault’s future seismic behavior. Doing so requires long records of earthquake recurrence. The lack of sufficiently long instrumental seismic records (that would be best suited for this task) has led to the development of other approaches that may constrain the recurrence of surface rupturing earthquakes along individual faults. These approaches take different forms, depending on the specific tectonic and geographic conditions of an investigated region.For example, around the Mediterranean Sea, we frequently find bedrock scarps along normal faults. Assuming that bedrock (i.e., fault free-face) exposure is caused by the occurrence of sub-sequent large earthquakes, we may measure certain rock properties to constrain the time and size of past earthquakes as well as the fault’s geologic slip-rate. A now-classic example in this regard is the measurement of 36Cl concentrations along exposed fault scarps in limestones.For the presented study, we looked at another property of the exposed fault free-face, namely its morphologic roughness. We aim to identify whether fault free-face roughness contains information to constrain earthquake occurrence and fault slip-rates following the assumption that sub-sequent exposure to the elements and sub-areal erosional conditions may leave a signal in how rough (or smooth) the fault free-face is (assuming a somewhat uniform pre-exposure roughness). Here, we present observations of fault free-face surface roughness for the Mt. Vettore fault (last ruptured in 2016) and the Rocca Preturo fault (The underlying models of fault free-face morphology were generated using the Structure-from-Motion approach and a large suite of unregistered optical images.). Employing different metrics to quantify morphologic roughness, we were indeed able to observe a) an increase in surface roughness with fault scarp height (i.e., longer exposure to sub-areal erosion causes higher roughness), and b) distinct (rather than gradual) changes in surface roughness, suggesting a correlation to individual exposure events such as earthquakes. Hence, fault free-face morphology of bedrock faults may serve as an additional metric to reconstruct earthquake recurrence patterns.

The interaction between sedimentation/erosion and faulting represents one of the most intriguing topics in landscape and tectonics evolution. Recently, several studies attempted to unravel this issue but only few of them have been able to document the feedback between faulting and sedimentary loading. Here, we focus on how the sediment loading/unloading influences the dynamic of the faults system taking as study case the Fucino Basin in Central Apennines (Italy). The Fucino Basin represents a remarkable case study with respect to the other main extensional basins in the Apennines, because of its large dimension, rectangular shape, significant sediment thickness and more important, its endorheic nature throughout its evolution.We present a detailed structural analysis all around the basin, investigating the kinematic and geometry of each main fault strand. The slickensides analysis reveals multiples families of slip-vectors and timing of activity which suggest a changing from N240° to N200° occurred during middle-Pleistocene. Moreover, using a simple isostatic model, we estimate that up to the 30% of the total geological displacement of the faults, which overall ranges from 0.5 to 3.5 km, is related to the sediments loading/unloading. Then, we demonstrate a positive feedback between sedimentation and faulting which may also lead to a re-organization in fault-kinematic related to a significant increasing in the vertical stress. Finally, we propose a conceptual model to support the permanent endorheic configuration of the Fucino basin, mainly related to a fault-slip increasing and kinematic changing due to the sediment loading.

In Central Italy more than 393 thousands people live in villages and towns located at less than 5 km distance from a known, mapped, active fault, capable of generating Mw>6 earthquake. Improving seismic risk estimates in such places requires the use of (i) informative databases of active faults and (ii) the implementation of appropriate building-codes. The current level of knowledge regarding activity of active faults in Central Italy has been stored in a recently compiled database (160 slip rates estimates for 88 faults). Given the complex nature of fault ruptures, we adopted a multi-fault rupture approach (SHERIFS) that accounts for both individual ruptures and multi-fault complex ruptures, involving more than one seismogenic fault section. Our earthquake rupture forecast model includes 1249 possible combinations of fault ruptures with lengths ranging from 7 to 42 km. Slip rates and associated errors are used to estimate recurrences of the ruptures assuming a Gutenberg-Richter frequency-magnitudedistribution. The computed distribution is validated against the CPTI15 catalogue.The multi-fault model approach and a seismogenic area approach are used to estimate damages based on published typological fragility curves for typical building classes derived from 30 years of data in Italy (Rota et al., 2006) assuming earthquake occurrence for the faults follows a Poisson time-independent process. Two fragility curves are considered here: one for reinforced concrete designed according to seismic regulations and one for masonry with irregular layout and without tie rods and tie beams, a typical typology for the region. Expected levels of damage for 150 villages and towns in Central Italy are computed for all damage states considering a 50 years risk target period. Results obtained with the fault approach show a much higher variability of the estimated risk depending on the location of the village/town w.r.t. the fault system and the hanging-wall/footwall location. The probability of collapse in 50 years for a typical masonry building ranges between 0.01 and 0.07 in the fault approach and 0.01 and 0.04 for the area approach. For both approaches, the probability of collapse for reinforced concrete buildings is ~90 % less than that for typical masonry structures. Even if this can be considered obvious, it must be underlined that most buildings in Italy were built before 1975 (before the first applicative decree of the seismic Italian law No. 64 of 1974). Thanks to the availability of the detailed database of active faults a strategy to prioritize resources for seismic risk reduction could be adopted. Rota, M., Penna, A. & Strobbia, C. (2006). Typological fragility curves from Italian earthquake damage data. First European Conference on Earthquake Engineering and Seismology Geneva, Switzerland, 3-8 September 2006 Paper Number: 386

Aix-Marseille University launched the VirtuaField project whose objective is to integrate DOMs into a VR application to provide students with a pedagogical tool enabling learning field practice.Indeed, students have few occasions to train in the field during their academic curricula. Field trips are expensive and require a complex logistics. Nowadays, the photogrammetry or LIDAR techniques allow geoscientists to obtain High-Resolution 3D representations of outcrop geometries and textures, often termed as Digital Outcrop Models (DOM). DOMs are already used as pedagogical supports for practical exercises on computers such as fault throw or seismic occurrence calculation, or modelling 3D geological structures from outcrop interpretations. However, these exercises do not cover all required skills to gain autonomy and consistence in the field, such as the pertinent observation sampling. The computer engines are not convenient support for that task because the visualization, although in 3D, still depends on a 2D screen and does not preserve the 1:1 scale, which is of paramount importance for Geoscience interpretations.The Virtual Reality (VR) technique is the ultimate way to provide a full 3D view, which can preserve the 1:1 scale, while benefiting from the numerical nature of the support (DOMs, DEM).First prototypes were provided by the VR2Planets company from the case study of La Fare les Oliviers (SE France), which shows diffuse fractures and fracture corridors, in addition to sedimentological and geomorphological structures. The prototypes have been tested in training experiences with volunteer students. Surveys have been performed in order to obtain feedbacks from students on the ability of the VirtuaField application to gain field skills, but also on the more pertinent way to design the pedagogical tools. The synthesis of these feedbacks will be presented as well as a first outline of the pedagogical guidelines on using VR tools for educational purposes.

Correction for ‘Mass-shift mode to quantify low level 129I in environmental samples by ICP-MS/MS’ by Coralie Carrier et al., J. Anal. At. Spectrom., 2022, 37, 1309–1317, https://doi.org/10.1039/D2JA00128D.

International audience; The dynamics of carbonate-dominated landscapes is still debated, and first-order questions pertaining to the respective importance of chemical and mechanical weathering processes or the influence of climatic variability on denudation rates are largely open. To address these issues, we investigate the denudation pattern across a climatic gradient in Provence, SE France. We focused on a series of major carbonate mountain ranges with elevation spanning 200 to 2000 m and associated with >2-fold variations in mean annual temperature and precipitation. We collected 25 samples from 10 different mountain ranges and quantified denudation rates through the measurement of their Cl-36 concentrations in order to explore tectonic and climatic controls on these carbonate landscapes evolution. We observe that most of our samples experience denudation rates in the 30-60 mm/ka range (averaged over 5 to 40 ka), with no evident dependence on variations of mean annual temperature or precipitation. One range (Lure mountain) stands out as a significant outlier, with much higher denudation rates (>100 mm/ka) than adjacent sites (Ventoux mountain) at the same elevation and experiencing similar present-day environmental conditions. We explore in details the climatic, tectonic and geomorphological contexts of these two ranges in order to understand the reasons for this apparent 3-fold difference in denudation rates. A regionally-documented major Holocene hillslope destabilization event appears to be a likely candidate to explain the observed lower Cl-36 concentrations at the top of the Lure mountain. At longer wavelengths, our regional dataset shows no clear control of modem precipitation gradients on denudation rates, consistent with previous observations in landscapes underlain by other lithologies.

International audience; Key-points Folded alluvial surfaces and balanced cross-section constrain post-LGM and Neogene shortening-rates Decrease in shortening-rate and change in compression direction are inferred at 1-2 Ma Kinematic reorganization is best explained by a shift from Nubia-driven to Adria-driven convergence Abstract The eastern Southern Alps are located at the northern tip of the Adria microplate, which imposes 2.0-2.5-mm/yr of N-S-convergence relative to stable Eurasia. We map surface evidence of recent folding/faulting in this area from a 5-m-DEM. In the eastern part of the belt, observations reveal a 30-km-wide zone of active folding composed of at least five growing anticlines. The most recent ones warped the post-glacial alluvial surface by

International audience; In this study, we present a new glacial chronology based on 20 in situ-produced 36Cl-based cosmic ray exposure datings from moraine boulders and bedrock from the Throne of Zeus (TZ) and Megala Kazania (MK) cirques on Mount Olympus. The 36Cl derived ages of glacial landforms range from 15.6 ± 2.0 to 0.64 ± 0.08 ka, spanning the Late-glacial and the Holocene. The Late-glacial, recorded in both cirques, is partitioned in three distinct phases (LG1-3): an initial phase of moraine stabilization at 15.5 ± 2.0 ka with subsequent deglaciation starting at ∼14 ka (LG1), followed by a shift to marginal conditions for glaciation at 13.5 ± 2.0 ka (LG2), sustained by large amounts of wind-blown snow, despite regional warming. Glacial conditions returned at 12.5 ± 1.5 ka (LG3) and were characterized by low air temperatures and glacier shrinking. The Holocene glacial phases (HOL1-3) are recorded only in the MK cirque, likely due to its topographic attributes. An early Holocene glacier stillstand (HOL1) at 9.6 ± 1.1 ka follows the regional temperatures recovery. No glacier activity is observed during the mid-Holocene. The Late Holocene glacier expansions, include a moraine stabilization phase (HOL2) at 2.5 ± 0.3 ka, during wet conditions and solar insolation minima, while (HOL3) corresponds to the early part of the Little Ice Age (0.64 ± 0.08ka). Our glacial chronology is coherent with glacial chronologies from several cirques along the northeast Mediterranean mountains and in pace with numerous proxies from terrestrial and marine systems from the north Aegean Sea.

International audience; Landscape evolution from the Early 1st millennium BCE to the mid-1st millennium CE is poorly documented around major archaeological sites in Crete. In a previous publication, the general landscape configuration in the vicinity of ancient Phaistos was reconstructed using a palaeoenvironmental approach, from the Proto-Palatial period (ca. 2000 BCE) to the Late-Proto Geometric period (ca. 8th cent. BCE). However, the physiography of the landscape, its hydrology and vegetation history remained uncertain for the later archaeological periods. In the present study, additional radiocarbon dates (8) together with pollen, mollusc and sedimentological analyses (CM diagram) were conducted on previously documented sediment cores. These new results enable us to reconstruct in greater detail the landscape history from the Early Archaic period to Late Byzantine times. The results indicate the continuous presence of swampland from the Proto-Geometric period (10th cent. BCE) probably until the initial stages of the Classical period (5th cent. BCE). Subsequently, during the Classical and Hellenistic periods, there was a short interval of alluvial input of terriginous sediments (not exceeding two centuries in duration) which is directly linked with the complete drying up and drainage of the swampland. We address the issue of the possible climatic origin of this abrupt hydrological change, especially in relation to regional climate change and the sedimentary history of adjacent rivers and streams. Tectonic activity in the area is also an important factor and can be invoked as a potential environmental influence. Anthropogenic factors are also considered, even though there is no direct archaeological evidence of drainage in the western Messara Plain during the Archaic and Classical periods. Finally, from Roman times to the Early Byzantine period, floodplain development prevailed in the area and ponds formed locally, in particular from the Late Hellenistic to Early Byzantine periods; this was related to the climatic conditions of the Roman Warm period. Pollen analysis reveals an open forested landscape during the time interval under investigation, within which domesticated plants such as Olea (olive) were present. However, the representation of Olea decreases continuously from the Late Geometric period to Byzantine times, probably indicating much lower intensity of land use than during Minoan times and possibly also related to the generally colder climatic conditions in Crete from the 8th cent. BCE until the 1st Cent. CE. During this latter interval, there is also the first pollen analytical evidence of Vitis sp.

An Mw 6.1, devastating earthquake, on April 6, 2009, struck the Middle Aterno Valley (Abruzzi Apennines, Italy) due to the activation of a poorly known normal fault system. Structural analysis of the fault population and investigation of the relationships with the Quaternary continental deposits through integrated field and laboratory techniques were conducted in order to reconstruct the long-term, tectono-sedimentary evolution of the basin and hypothesize the size of the fault segment. A polyphasic evolution of the Middle Aterno Valley is characterized by a conjugate, ∼E-W and ∼N S-striking fault system, during the early stage of basin development, and by a dip-slip, NW-striking fault system in a later phase. The old conjugate fault system controlled the generation of the largest sedimentary traps in the area and is responsible for the horst and graben structures within the basin. During the Early Pleistocene the E-W and N S system reactivated with dip-slip kinematics. This gave rise to intra-basin bedrock highs and a significant syn-tectonic deposition, causing variable thickness and hiatuses of the continental infill. Subsequently, since the end of the Early Pleistocene, with the inception of the NW-striking fault system, several NW-strands linked into longer splays and their activity migrated toward a leading segment affecting the Paganica-San Demetrio basin: the Paganica-San Demetrio fault alignment. The findings from this work constrain and are consistent with the subsurface basin geometry inferred from previous geophysical investigations. Notably, two major elements of the ∼E-W and ∼N S-striking faults likely act as transfer to the nearby stepping active fault systems or form the boundaries, as geometric complexities, that limit the Paganica-San Demetrio fault segment overall length to 19 ± 3 km. The resulting size of the leading fault segment is coherent with the extent of the 6 April 2009 L'Aquila earthquake causative fault. The positive match between the geologic long-term and coseismic images of the 2009 seismogenic fault highlights that the comprehensive reconstruction of the deformation history offers a unique contribution to the understanding faults seismic potential.