Slip rate of trench‐parallel normal faulting along the Mejillones Fault (Atacama Fault System): Relationships with the northern Chile subduction and implications for seismic hazards

The recent tectonics of the arid northern Chile Andean western forearc is characterized by trench‐parallel normal faults within the Atacama Fault System (AFS). Since the 1995‐Mw 8.1 Antofagasta earthquake, the mechanism driving this recent and localized extension is considered to be associated with the seismic cycle within the subduction zone. Analyzing morphotectonic patterns along these faults allows examining the seismic potential associated with the subduction zone. Using field Digital Elevation Models and in situ‐produced cosmogenic ¹⁰Be, we determined a 0.2 mm/a long‐term slip rate along the Mejillones Fault, one of the most prominent structures within the AFS. This result suggests that the AFS corresponds to slow slip rate faults despite the rapid subduction context. However, the size of coseismic slips observed along the AFS faults suggests that larger subduction earthquakes (Mw > 8.1) may occur episodically in the area.     

Failure in accretionary wedges with the maximum strength theorem: numerical algorithm and 2D validation

The objective is to capture the 3D spatial variation in the failure mode occurring in accretionary wedges and their analog experiments in the laboratory from the sole knowledge of the material strength and the structure geometry. The proposed methodology relies on the maximum strength theorem which is inherited from the kinematic approach of the classical limit analysis. It selects the optimum virtual velocity field which minimizes the tectonic force. These velocity fields are constructed by interpolation thanks to the spatial discretization conducted with ten-noded tetrahedra in 3D and six-noded triangles in 2D. The resulting, discrete optimization problem is first presented emphasizing the dual formalism found most appropriate in the presence of nonlinear strength criteria, such as the Drucker–Prager criterion used in all reported examples. The numerical scheme is first applied to a perfectly triangular 2D wedge. It is known that failure occurs to the back for topographic slope smaller than and to the front for slope larger than a critical slope, defining subcritical and supercritical slope stability conditions, respectively. The failure mode is characterized by the activation of a ramp, its conjugate back thrust, and the partial or complete activation of the décollement. It is shown that the critical slope is captured precisely by the proposed numerical scheme, the ramp, and the back thrust corresponding to regions of localized virtual strain. The influence of the back-wall friction on this critical slope is explored. It is found that the failure mechanism reduces to a thrust rooting at the base of the back wall and the absence of back thrust, for small enough values of the friction angle. This influence is well explained by the Mohr construction and further validated with experimental results with sand, considered as an analog material. 3D applications of the same methodology are presented in a companion paper.     

Geochemical modelling of fluid–rock interactions in the context of the Soultz-sous-Forêts geothermal system

The development of the Enhanced Geothermal System (EGS) at Soultz-sous-Forêts (France) has given to scientists an interesting opportunity for the application of geochemical modelling of water–rock interactions, combining theoretical studies with field and experimental data. The main results of four successive and complementary studies are summarized: geochemical modelling of fluid–rock interactions with prediction of dissolution/precipitation of minerals, feed-back effects on the mineralogy and petrography of the rock (major role of silicates in the geological past and of carbonates in the near future of the exploitation), experimental control of the dynamics of silicates under thermal gradient and relation between the evolution of the petrophysics of the rocks and the heat and mass transfers. The thermal cycle of the fluid, between 200 °C and 65 °C in the geothermal loop, may be responsible for dissolution/precipitation of minerals which modify the porosity and permeability of the granite, as it happened in the geological past, in relation with hydrothermal circulations in the Rhine Graben.     

Hydraulic interactions between fractures and bedding planes in a carbonate aquifer studied by means of experimentally induced water-table fluctuations (Coaraze experimental site, southeastern France)

In aquifers with variable permeabilities, the water exchanges between high- and low-permeability regions are controlled by the hydraulic head gradient. Past studies have addressed this problem mainly considering steady-state hydraulic conditions. To study such exchanges during water-table fluctuations, a spring in southeastern France was equipped with a water-gate that creates 10-m artificial fluctuations of the water table. The water exchanges are discussed with respect to hydrochemical and pressure measurements in the groundwater. With successive water-table fluctuations the mineralization and pH decrease, but the bicarbonate content increases in response to carbon dioxide dissolution. At this scale of single fractures and surroundings, the hydrochemistry allows water flows from low- or high-permeability discontinuities to be discriminated. During hydraulic head fluctuations, the waters from low- and high-permeability discontinuities become mixed. During water-table rise, the low-permeability matrix contributed to the refilling of the permeable faults and to the mixing of the waters. Dynamic flows in the opposite direction to the hydraulic gradient can expel mineralized water, pushing it towards permeable discontinuities. This mechanism could be the cause of the peak in the mineralization recorded in some karst springs at the start of flood events.     

Regional modeling of carbonaceous aerosols over Europe—focus on secondary organic aerosols

In this study, an improved and complete secondary organic aerosols (SOA) chemistry scheme was implemented in the CHIMERE model. The implementation of isoprene chemistry for SOA significantly improves agreement between long series of simulated and observed particulate matter concentrations. While simulated organic carbon concentrations are clearly improved at elevated sites by adding the SOA scheme, time correlation are impaired at low level sites in Portugal, Italy and Slovakia. At several sites a clear underestimation by the CHIMERE model is noticed in wintertime possibly due to missing wood burning emissions as shown in previous modeling studies. In Europe, the CHIMERE model gives yearly average SOA concentrations ranging from 0.5 μg m ^− 3 in the Northern Europe to 4 μg m ^− 3 over forested regions in Spain, France, Germany and Italy. In addition, our work suggests that during the highest fire emission periods, fires can be the dominant source of primary organic carbon over the Mediterranean Basin, but the SOA contribution from fire emissions is low. Isoprene chemistry has a strong impact on SOA formation when using current available kinetic schemes.     

A two‐step process for the reflooding of the Mediterranean after the Messinian Salinity Crisis

The Messinian Salinity Crisis is well known to have resulted from a significant drop of the Mediterranean sea level. Considering both onshore and offshore observations, the subsequent reflooding is generally thought to have been very sudden. We present here offshore seismic evidence from the Gulf of Lions and re‐visited onshore data from Italy and Turkey that lead to a new concept of a two‐step reflooding of the Mediterranean Basin after the Messinian Salinity Crisis. The refilling was first moderate and relatively slow accompanied by transgressive ravinement, and later on very rapid, preserving the subaerial Messinian Erosional Surface. The amplitude of these two successive rises of sea level has been estimated at ≤500 m for the first rise and 600–900 m for the second rise. Evaporites from the central Mediterranean basins appear to have been deposited principally at the beginning of the first step of reflooding. After the second step, which preceeded the Zanclean Global Stratotype Section and Point, successive connections with the Paratethyan Dacic Basin, then the Adriatic foredeep, and finally the Euxinian Basin occurred, as a consequence of the continued global rise in sea level. A complex morphology with sills and sub‐basins led to diachronous events such as the so‐called ‘Lago Mare’.This study helps to distinguish events that were synchronous over the entire Mediterranean realm, such as the two‐step reflooding, from those that were more local and diachronous. In addition, the shoreline that marks the transition between these two steps of reflooding in the Provence Basin provides a remarkable palaeogeographical marker for subsidence studies.     

Late Holocene high precipitation events recorded in lake sediments and catchment geomorphology,Lake Vuoksjávrátje,NW Sweden

In this paper, we show the potential of combining multi‐proxy analysis of lake sediments with catchment geomorphology to better understand palaeoenvironmental changes. Previous studies have suggested that alpine lakes in N Sweden may be highly sensitive to variations in catchment erosion and precipitation, and that this sensitivity may influence ecologically based reconstructions of past temperature changes. We analysed lake sediments covering the last 5100 years from the alpine Lake Vuoksjávrátje in NW Sweden in order to identify different erosional regimes in the lake catchment, sediment sources and lake sedimentary processes, which ultimately affect the palaeoecological record. The measured proxies include elemental geochemistry from XRF core scanning, grain size, sediment accumulation rates, fraction of terrestrial organic carbon and geomorphological mapping, supported by previously published chironomid and total organic carbon data. From the integrated results we identified time intervals when increased intensity of precipitation altered sedimentation and lake catchment erosional processes. The most prominent event occurred c. 2900 cal. a BP and is interpreted to be the result of excessive precipitation, possibly related to the climatic shift towards cooler and wetter conditions referred to as the 2.8 ka event.     

Geomorphologic evidence for Plio‐Quaternary shortening in the southern Neuquén basin (40°S,Argentina)

A morphostructural analysis of a Pliocene flood basalt formation in the southern Neuquén basin (40°S) shows evidence of contractional deformation less than 3.5 Ma ago. This formation exhibits a general dip towards the south‐east, with relict outcrops located 100 m higher than the main source volcano, which suggests a local tilting of the lava flow. This tilt has been brought about by Plio‐Quaternary reactivation of the eastern border of the Sañico Massif along two thrusts that offset the lava flow. Another long‐wavelength bulge in the southern part of the lava flow unit indicates a possible Pliocene uplift of the North Patagonian Massif. These results provide new evidence of continuing shortening in the Neuquén basin during the Plio‐Quaternary, challenging the hypothesis that an extensional regime has existed since the end of the Miocene in this basin.     

Thermal stresses in borehole heat exchangers

The backfilling materials of borehole heat exchangers (BHE), particularly the grout material, must provide a suitable thermal contact and ensure durability to the induced thermal stresses because of the heat loading. In this paper, the thermal stresses that occurred in BHEs because of heat injection or extraction is investigated with an analytical solution of a hollow cylinder model that is adapted for time‐dependent heat loading, the geometry of a BHE, and the thermo‐mechanical properties of surrounding ground conditions. Firstly, the hollow cylinder model is solved with the considered boundary conditions in 2D plane stress. Secondly, the temperature differences at the inner and outer circles of the cylinder are evaluated with the heat line source models for continuous and discontinuous loading to observe the impact of the heat loading schedule. The developed analytical solution for thermal stress investigation is validated with numerical models. It is demonstrated that the analytical solutions agree well with numerical results for two types of BHE configurations (co‐axial and single U‐shaped pipes). Furthermore, the calculated maximum stresses are compared with the tensile strength of grout materials obtained from Brazilian tests. It is predicted that the thermal contraction of the grout, partially constrained by the surrounding rock, generates tensile stresses that may lead to cracking in the BHE. According to the results, the stiffness of rock has a primary role on the developed tensile stresses, and the relationship between the thermal conductivity of the ground and of the grout induces a proportional impact on the magnitude of thermal stresses. Copyright © 2015 John Wiley & Sons, Ltd.