Even though stresses in the crust are triaxial (𝜎1>𝜎2>𝜎3) the overwhelming majority of rock deformation experiments are conducted under axisymmetric (or conventional triaxial) loading (𝜎1>𝜎2=𝜎3). This configuration disregards the effect of 𝜎2 on the physical and deformation properties of rocks, thus complicating and degrading the extrapolation of results to natural crustal conditions. A True Triaxial loading configuration is necessary to overcome this simplification, however, these improvements in addressing real crustal conditions come at a cost, which is the challenging boundary conditions that arise from having six loading rams rather than just two. Two main loading boundary effects can severely impact the stress distribution and failure mechanism of samples deformed in a True Triaxial Apparatus (TTA): 1) the end friction effect caused by the stiffness contrast between the rock sample and the metal loading platens, and 2) the unstressed sample edges resulting from the requirement that loading platens must necessarily be slightly smaller than the rock specimen. Managing and reducing these boundary effects is fundamental for obtaining accurate and representative data from true triaxial experiments, and for the further development of these apparatuses.A novel TTA developed in the UCL Rock & Ice Physics Laboratory was designed to subject cubic or cuboid rock samples to truly triaxial stresses through the independent control of six loading rams. The apparatus is equipped with a confining and pore pressure system that allows for the deformation of saturated samples whilst simultaneously measuring permeability along the three axes. A suite of Finite Element Method (FEM) models was implemented to evaluate the parameters that minimize loading boundary effects in UCL’s TTA for a 50 mm edge-length cubic sample of sandstone. Our results indicate that using aluminum loading platens (𝐸𝑠𝑎𝑚𝑝𝑙𝑒/𝐸𝑝𝑙𝑎𝑡𝑒𝑛=0.47) reduces the end friction effect by a factor of two compared to using steel platens (𝐸𝑠𝑎𝑚𝑝𝑙𝑒/𝐸𝑝𝑙𝑎𝑡𝑒𝑛=0.17). In addition, we find that elevated confining pressure significantly reduces the stress concentration produced by unstressed edges. Specifically, a confining pressure of 10 MPa eliminates tensile stresses at the sample corners. These results are currently being implemented into the experimental protocol and execution in UCL’s TTA in order to ensure that we obtain reliable true triaxial data. However, these observations are generic and could therefore contribute to improved development and operation of true triaxial loading systems generally.

The Taurus-Littrow Valley, location of the Apollo 17 landing site, hosts recent, late-Copernican geomorphological landforms and tectonic structures, namely the Light Mantle avalanche deposit and the Lee-Lincoln lobate scarp. The Light Mantle deposit represents a unique case of a hypermobile avalanche on the Moon (El-Baz 1972; Schmitt et al. 2017). The Lee-Lincoln lobate scarp is the surface expression of a recent thrust fault (Watters et al. 2010), which is considered to be the source of strong seismic shaking throughout Taurus-Littrow Valley (van der Bogert et al. 2012, 2019), and potentially still active (Watters et al. 2019).The Light Mantle represents the only extraterrestrial landslide for which an absolute age is provided (70-110 Ma), thanks to the Apollo 17 returned samples (e.g., Schmitt et al. 2017). Therefore, the Light Mantle deposit can be used as a geomorphological marker and time constraint for surface changes that occurred since its emplacement. By applying the principle of superposition, surface changes superposed on the Light Mantle deposit, and on the slope from which it was generated (the NE-facing slope of the South Massif), must post-date the landslide event. For example, small scale grabens (10-20 m wide) associated with the Lee-Lincoln lobate scarp are found superposed on the Light Mantle unit (Watters et al. 2010). These troughs likely formed less than 50 Ma and are thought to be generated by the flexural bending of the hanging wall (Watters et al. 2010, 2012). Similarly, boulder tracks, whose survival time is estimated to range up to 35 Ma (e.g., Arvidson et al. 1976; Kumar et al. 2019), are found on the NE-facing slope of the South Massif, therefore evidence that boulder falls have occurred after the Light Mantle landslide event.Here, we extend the body of evidence of surface changes that have affected the South Massif since the emplacement of the Light Mantle deposit. We map boulder tracks, areas of disturbed regolith, linear slope structures, and other structures associated with the summit of the South Massif. We identified features (i.e., slope structures oblique to contours, the Nansen Moat and the trough at the NE-base of the Sout Massif) directly related to back-thrust faults associated with the Lee-Lincoln thrust fault, which are re-activating the buried fault that bounds Taurus-Littrow Valley; we identified other features (i.e., crestal graben-like structures, slope structures parallel to contours) that derived from gravitational adjustment following basal slope support removal due to back-thrust faulting. Moreover, the overlapping relationships between the boulder tracks and regolith disturbance suggests that continuous slope deformation has been affecting the NE-facing slope. We attribute the efficiency of the process to repeated ground-shaking perturbation, which maintains the slope in a perpetually unstable state.We conclude that the NE-facing slope of the South Massif has been recently and continuously affected by slope deformation processes. We suggest that the efficiency of these processes is the product of lasting, and perhaps ongoing, effects of activity of the Lee-Lincoln thrust fault, coupled with the influence of the subsurface geometry of the valley inherited from the impact basin formation.

The Hengill region is one of the largest areas of high geothermal gradient and subsurface heat flow in Iceland, and hosts two of its largest geothermal power plants: Hellisheiði and Nesjavellir, which have a combined capacity of 423 MWe and 560 MWth. The Hengill region is located in a unique tectonic setting, characterized by the convergence of three plate boundary segments: the oblique-spreading Reykjanes Peninsula (RP), the orthogonal-spreading rift of the Western Volcanic Zone (WVZ) and the transform, South Iceland Seismic Zone (SISZ). Unlike most tectonic triple junctions, which occur on the ocean floor, the Hengill Triple Junction (HTJ) is exposed above sea level, thus providing a unique opportunity to study the interplay between three plate boundary segments and the local deformation processes occurring at their convergence site. Additionally, the injection of fluids due to on-going geothermal operations enhances the natural tendency of the region for seismic activity, and results in a significant level of induced seismic hazard. Because slip on pre-existing faults is triggered when the applied shear stress surpasses the frictional strength of the fault, regions that are naturally subjected to higher shear stresses are more prone to fault re-activation due to fluid re-injection. Therefore, a spatial variation in tectonic stresses may result in varying induced seismicity potential within a region.The local interplay of the three converging tectonic regimes, and their effect on the stress fields within the triple junction region, has been examined through a combination of regional structural mapping and a numerical model of the plate boundary interactions using the Boundary Element Method (BEM). Large scale structural mapping and analytical models of oblique rifting were used to estimate the degree of rift obliquity for individual fissure swarms. Our results reveal that the transition from the highly oblique rift system of the RP towards the spreading-orthogonal rift of the WVZ is smooth, and manifests as a rotation of the trend of fissures and eruptive ridges, and the strike of normal faults formed in response to the local stress field developed in the HTJ. These results were then correlated with those from the BEM model, which allows us to predict the orientation, relative magnitude, and distribution of stresses within the study area. Finally, the shear stress distribution determined from the BEM model was plotted against the location of both natural and induced seismic events detected in the region over a time span of 26 months, by the COSEISMIQ project (Grigoli et al., 2022). Our results show that seismic events cluster in either at the triple junction or SW of it, within the highly stressed regions of the RP. Furthermore, the seismicity transitions from scattered to non-existent towards the north of the region, where shear stresses also diffuse. The good correlation between the high shear stress regions predicted by the model and distribution of seismicity suggests that this approach may provide a valuable and cost-effective tool for seismic hazard prediction within regions with complex tectonic settings.

<p>Fractures in rock are ubiquitous; from cold dry planetary bodies to the hottest, wettest climates on Earth, and from km-scale tectonic fractures deep in Earth’s crust to microcracks in surficial rocks. Yet, many of these fractures propagate progressively over geologic timescales, making their development complex and enigmatic. Therefore, to measure how fractures have developed in rocks exposed at Earth’s surface over millennia, and how this consequently changes rock physical properties, we collected ten ~25 cm diameter granitic boulders from two sites in the Eastern Sierra, California, USA. The boulders were deposited on the surface of alluvial terraces and fans during geologically instantaneous glacial and alluvial events at different times since about 148ka BP, then the depositional surfaces were subsequently abandoned. The <em>chronosequences </em>of geomorphic surfaces provide a natural laboratory in which rocks of consistent lithology have been exposed to similar environmental conditions for different lengths of time, allowing us to compare rock property evolution on the order of 0 to 10<sup>5</sup> years of environmental exposure; an approach that allows us to better understand and characterize mechanical weathering processes, especially long-term changes in rock fracturing. Note that <em>fresh</em> (time-zero) rocks in this study are represented by boulders found within active channels, and that the measured changes in rocks with longer exposure times are interpreted by comparison with the fresh rocks. Focusing only on similarly sized boulders removes any ambiguities in tectonic and exhumation history that might arise in outcrop samples, thus ensuring that rocks from each site have experienced similar stress conditions; namely those restricted to the environment.</p><p>We performed laboratory measurements on 10 granitic boulders (four from Lundy Canyon, with exposure ages of ~0 to ~148 ka; six from Shepherd Creek, with exposure ages of ~0 to ~117 ka) to quantify how rock physical properties changed as a function of environmental exposure age. We measured key parameters commonly used as proxies for crack damage, including porosity, compressional wave velocity (Vp), and shear wave velocity (Vs). We hypothesize that changes in crack damage are likely to affect rock mechanical properties, so we also measured tensile strength, uniaxial compressive strength (UCS), and Young’s modulus (E). We find that all measured parameters evolve as a function of exposure age, with systematic increases in porosity, and systematic decreases in Vp, Vs, tensile strength, UCS, and E. For example, porosity increases from 0.5 – 1.0 % in the fresh rock to 2.6 – 3.2 % in the oldest rocks. We interpret these changes as reflecting progressive subcritical crack growth that arises due to ubiquitous, but relatively low magnitude, environmental stresses continuously acting on the boulders, as opposed to differences inherited before their erosion from bedrock.</p><p>Apart from demonstrating the importance of environmentally driven cracking in rock weathering, these observations of progressive crack damage accumulation also have significant implications for the interpretation of any measurements made on rocks exposed at Earth’s surface, even if the age of exposure is relatively short compared to the age of the geologic deposit itself.</p>

Long runout landslides are a particular type of mass-wasting phenomena that belongs to the category of surface processes associated with rapid strain rates. The reduction of friction that has to be invoked to explain their high velocities and exceptional travel distance over nearly horizontal surfaces has yet to find satisfactory explanation. Inspired by fault mechanics studies, thermally-activated mechanisms can explain the dynamic frictional strength loss during sliding along the initial failure surface and the early development of velocities higher than expected. However, as slides continue moving along nearly horizontal valley floors, the weakening mechanisms required to sustain their exceptional behaviour are less certain.Long runout landslides are found ubiquitous in our solar system and the slow erosion rates that operate on extraterrestrial planetary bodies allow the preservation of their geomorphological record. The availability of the latest high-resolution imagery of the surface of Mars and the Moon allows to conduct detailed morphometric analysis not so granted on our planet. On the other hand, on Earth, the partial loss of the geomorphological record due to fast erosion rates is compensated by the accessibility of sites that enable us to conduct field work. In order to better understand the mechanisms responsible for the apparent friction weakening we use a comparative planetary geology approach, in the attempt to link the morphology and the internal structures of long runout landslide deposits to the mechanisms involved during the emplacement of such catastrophic events.We focused on the distinctive longitudinal ridges that mark the surface of the landslide deposits. The formation mechanism of longitudinal ridges in long runout landslides has been proposed to require ice, as this low friction material would allow the spreading of the deposit, causing the development of longitudinal ridges by tensile deformation of the slide. However, ice-free laboratory experiments on rapid granular flows have demonstrated that longitudinal ridges can form as a consequence of helicoidal cells that generate from a mechanical instability, which onset requires a rough surface and a velocity threshold to be surpassed. Moreover, such experiments have showed that the wavelength of the longitudinal ridges is always 2 to 3 times the thickness of the flow.We here present the results from three case studies: the 63-km-long Coprates Labe landslide in Valles Marineris on Mars; the 4-km-long El Magnifico landslide in Chile, Earth; and the 50-km-long Tsiolkovskiy crater landslide, at the far side of the Moon. We found that the wavelength of the longitudinal ridges is consistently 2 to 3 times the thickness of the landslide deposit, in agreement with experimental work on rapid granular flows. The recurrence of such scaling relationship suggests a scale- and environment-independent mechanism. We discuss the applicability of high-speed granular flow convection-style mechanisms to long runout landslides and speculate on the existence of an alternative vibration-assisted mechanism.

Optical image correlation (OIC) is a powerful tool for measuring the 3-D near-field surface displacements produced in large earthquakes. The method compares pre- and post-earthquake orthorectified images; shifts between common pixels in the image pair reflects a 2-D (horizontal) offset. The third dimension (vertical displacement) is calculated by differencing the pre- and post-topography, while accounting for the horizontal displacements. Optical image correlation has a sub-pixel detection capability, and can provide information on the displacement field close to fault ruptures (where InSAR typically decorrelates). Small-scale measurements of the distributed damage provide important constraints on the strain distribution within the fault core and the surrounding damage zone, as well as offering insights into the rupture mechanics. OIC is frequently applied to the recent earthquakes where the image footprint is large relative to the rupture extents. However, historical ruptures are documented by aerial photographs which cover a relatively small area. This means that many images are needed to cover the rupture area and all pixels in pre-and post-earthquake images which span the rupture are typically affected by the ground displacement. This creates complications for image co-registration, alignment and correlation of the final mosaics. To address this problem we developed a workflow that automatically generates a DEM (digital elevation model) and an orthorectified image mosaic. The process uses structure-from-motion (SfM) and stereo-matching approaches, and results in precise and accurate registration between the image pairs. We applied this method to the 1959 Hebgen Lake earthquake, SW Montana, U.S. This large (Mw 7.2) intraplate normal event re-activated pre-existing faults north of the Hebgen Lake reservoir and created a complex rupture network. We used 20 pre-earthquake photographs from 1947 and 70 post-earthquake images from 1977 and 1982. The final results show a 3-D displacement localized onto several prominent structures: the Hebgen fault and the Red Canyon fault, consistent with field mapping following the earthquake. A significant vertical offset and a large horizontal NS-component agree well with SW extension on NW-SE-striking normal faults. Additionally, we used fault-perpendicular profiles to explore the along-strike variation in fault displacement and to determine the extent of the off-fault damage. This work demonstrates that the application of OIC techniques to historical earthquakes can provide new information relating to the geometry and displacement of fault ruptures, and isolate the last event from the previously accumulated displacements. Additionally, the method we propose offers potential for the characterisation of historical earthquakes in general, and promises to improve our understanding of rupture behaviour through a statistical analysis of many earthquakes.