Bedrock geology of DFDP-2B, central Alpine Fault, New Zealand

<p>During the second phase of the Alpine Fault, Deep Fault Drilling Project (DFDP) in the Whataroa River, South Westland, New Zealand, bedrock was encountered in the DFDP-2B borehole from 238.5–893.2 m Measured Depth (MD). Continuous sampling and meso- to microscale characterisation of whole rock cuttings established that, in sequence, the borehole sampled amphibolite facies, Torlesse Composite Terrane-derived schists, protomylonites and mylonites, terminating 200–400 m above an Alpine Fault Principal Slip Zone (PSZ) with a maximum dip of 62°. The most diagnostic structural features of increasing PSZ proximity were the occurrence of shear bands and reduction in mean quartz grain sizes. A change in composition to greater mica:quartz + feldspar, most markedly below c. 700 m MD, is inferred to result from either heterogeneous sampling or a change in lithology related to alteration. Major oxide variations suggest the fault-proximal Alpine Fault alteration zone, as previously defined in DFDP-1 core, was not sampled.</p

Petrophysical, Geochemical, and Hydrological Evidence for Extensive Fracture-Mediated Fluid and Heat Transport in the Alpine Fault's Hanging-Wall Damage Zone

International audienceFault rock assemblages reflect interaction between deformation, stress, temperature, fluid, and chemical regimes on distinct spatial and temporal scales at various positions in the crust. Here we interpret measurements made in the hanging‐wall of the Alpine Fault during the second stage of the Deep Fault Drilling Project (DFDP‐2). We present observational evidence for extensive fracturing and high hanging‐wall hydraulic conductivity (∼10−9 to 10−7 m/s, corresponding to permeability of ∼10−16 to 10−14 m2) extending several hundred meters from the fault's principal slip zone. Mud losses, gas chemistry anomalies, and petrophysical data indicate that a subset of fractures intersected by the borehole are capable of transmitting fluid volumes of several cubic meters on time scales of hours. DFDP‐2 observations and other data suggest that this hydrogeologically active portion of the fault zone in the hanging‐wall is several kilometers wide in the uppermost crust. This finding is consistent with numerical models of earthquake rupture and off‐fault damage. We conclude that the mechanically and hydrogeologically active part of the Alpine Fault is a more dynamic and extensive feature than commonly described in models based on exhumed faults. We propose that the hydrogeologically active damage zone of the Alpine Fault and other large active faults in areas of high topographic relief can be subdivided into an inner zone in which damage is controlled principally by earthquake rupture processes and an outer zone in which damage reflects coseismic shaking, strain accumulation and release on interseismic timescales, and inherited fracturing related to exhumation

Using low-frequency earthquakes to monitor slow tectonic deformation in the central Southern Alps, New Zealand

This thesis involves the study of low-frequency earthquakes (LFEs) in the central Southern Alps. The Alpine Fault is the principal locus of deformation within the Australia–Pacific plate boundary in the South Island of New Zealand and it is late in its typical ∼300-year seismic cycle. Surveying the seismicity associated with slow deformation in the vicinity of the Alpine Fault may provide constraints on the stresses acting on a major transpressive margin prior to an anticipated great (≥M8) earthquake. Here, we use 8 years of data from the Southern Alps Microearthquake Borehole Array (SAMBA) (amongst those, 3 years of data were collected as part of this project) in order to: (1) generate an updated LFE catalogue using an improved matched-filter technique that incorporates phase-weighted stacking; (2) compute LFE focal mechanisms and invert them to infer the crustal stress field on the deep extent of the Alpine Fault; (3) expand the LFE catalogue to cover a wider range of spatial/temporal behaviours; (4) study LFE families’ characteristics to identify periods where slow slip might happen. We first use fourteen primary LFE templates in an iterative matched-filter and stacking routine, which allows the detection of similar signals and produces LFE families sharing common locations. We generate an 8-yr catalogue containing 10,000 LFEs that are combined for each of the 14 LFE families using phase-weighted stacking to produce signals with the highest possible signal-to-noise ratios. We find LFEs to occur almost continuously during the 8-yr study period and we highlight two types of LFE distributions: (1) discrete behaviour with an inter-event time exceeding 2 minutes; (2) burst-like behaviour with an inter-event time below 2 minutes. The discrete events are interpreted as small-scale frequent deformation on the deep extent of the Alpine Fault and the LFE bursts (corresponding in most cases to known episodes of tremor or large regional earthquakes) are interpreted as brief periods of increased slip activity indicative of slow slip. We compute improved non-linear earthquake locations using a 3D velocity model and find LFEs to occur below the seismogenic zone at depths of 17–42 km, on or near the hypothesised deep extent of the Alpine Fault. We then compute the first estimates of LFE focal mechanisms associated with continental faulting. Focal mechanisms, in conjunction with recurrence intervals, are consistent with quasi-continuous shear faulting on the deep extent of the Alpine Fault. We then generate a new catalogue that regroups hundreds of LFE families. This time 638 synthetic LFE waveforms are generated using a 3D grid and used as primary templates in a matched-filter routine. Of those, 529 templates yield enough detections during the first iteration of the matched-filter routine (≥ 500 detections over the 8-yr study period) and are kept for further analysis. We then use the best 25% of correlated events for each LFE family to generate linear stacks which create new LFE templates. From there, we run a second and final iteration of the matched-filter routine with the new LFE templates to obtain our final LFE catalogue. The remaining 529 templates detect between 150 and 1,671 events each totalling 300,996 detections over the 8-yr study period. Of those 529 LFEs, we manage to locate 378 families. Their depths range between 11 and 60 km and LFEs locate mainly in the southern part of the SAMBA network. We finally examine individual LFE family rates and occurrence patterns. They indicate that LFE sources seem to evolve from an episodic or ‘stepped’ to a continuous behaviour with depth. This transition may correspond to an evolution from a stick-slip to a stable-sliding slip regime. Hence, we propose that the distinctive features of LFE occurrence patterns reflect variations in the in-situ stress and frictional conditions at the individual LFE source locations on the Alpine Fault. Finally, we use this new extensive catalogue as a tool for in-depth analyses of the deep central Alpine Fault structure and its slip behaviour. We identify eight episodes of increased LFE activity between 2009 and 2017 and provide time windows for further investigations of tremor and slow slip. We also study the spatial and temporal behaviours of LFEs and find that LFEs with synchronous occurrence patterns tend to be clustered in space. We thus suggest that individual LFE sources form spatially coherent clusters that may represent localised asperities or elastic patches on the deep Alpine Fault interface. We infer that those clusters may have a similar rheological response to tectonic forcing or to potential slow slip events. Eventually, we discover slow (10km/day) and rapid (∼20-25km/h) migrations of LFEs along the Alpine Fault. The slow migration might be controlled by slow slip events themselves while the rapid velocities could be explained by the LFE sources’ intrinsic properties

Using low-frequency earthquakes to monitor slow tectonic deformation in the central Southern Alps, New Zealand

This thesis involves the study of low-frequency earthquakes (LFEs) in the central Southern Alps. The Alpine Fault is the principal locus of deformation within the Australia–Pacific plate boundary in the South Island of New Zealand and it is late in its typical ∼300-year seismic cycle. Surveying the seismicity associated with slow deformation in the vicinity of the Alpine Fault may provide constraints on the stresses acting on a major transpressive margin prior to an anticipated great (≥M8) earthquake. Here, we use 8 years of data from the Southern Alps Microearthquake Borehole Array (SAMBA) (amongst those, 3 years of data were collected as part of this project) in order to: (1) generate an updated LFE catalogue using an improved matched-filter technique that incorporates phase-weighted stacking; (2) compute LFE focal mechanisms and invert them to infer the crustal stress field on the deep extent of the Alpine Fault; (3) expand the LFE catalogue to cover a wider range of spatial/temporal behaviours; (4) study LFE families’ characteristics to identify periods where slow slip might happen.  We first use fourteen primary LFE templates in an iterative matched-filter and stacking routine, which allows the detection of similar signals and produces LFE families sharing common locations. We generate an 8-yr catalogue containing 10,000 LFEs that are combined for each of the 14 LFE families using phase-weighted stacking to produce signals with the highest possible signal-to-noise ratios. We find LFEs to occur almost continuously during the 8-yr study period and we highlight two types of LFE distributions: (1) discrete behaviour with an inter-event time exceeding 2 minutes; (2) burst-like behaviour with an inter-event time below 2 minutes. The discrete events are interpreted as small-scale frequent deformation on the deep extent of the Alpine Fault and the LFE bursts (corresponding in most cases to known episodes of tremor or large regional earthquakes) are interpreted as brief periods of increased slip activity indicative of slow slip. We compute improved non-linear earthquake locations using a 3D velocity model and find LFEs to occur below the seismogenic zone at depths of 17–42 km, on or near the hypothesised deep extent of the Alpine Fault. We then compute the first estimates of LFE focal mechanisms associated with continental faulting. Focal mechanisms, in conjunction with recurrence intervals, are consistent with quasi-continuous shear faulting on the deep extent of the Alpine Fault.  We then generate a new catalogue that regroups hundreds of LFE families. This time 638 synthetic LFE waveforms are generated using a 3D grid and used as primary templates in a matched-filter routine. Of those, 529 templates yield enough detections during the first iteration of the matched-filter routine (≥ 500 detections over the 8-yr study period) and are kept for further analysis. We then use the best 25% of correlated events for each LFE family to generate linear stacks which create new LFE templates. From there, we run a second and final iteration of the matched-filter routine with the new LFE templates to obtain our final LFE catalogue. The remaining 529 templates detect between 150 and 1,671 events each totalling 300,996 detections over the 8-yr study period. Of those 529 LFEs, we manage to locate 378 families. Their depths range between 11 and 60 km and LFEs locate mainly in the southern part of the SAMBA network. We finally examine individual LFE family rates and occurrence patterns. They indicate that LFE sources seem to evolve from an episodic or ‘stepped’ to a continuous behaviour with depth. This transition may correspond to an evolution from a stick-slip to a stable-sliding slip regime. Hence, we propose that the distinctive features of LFE occurrence patterns reflect variations in the in-situ stress and frictional conditions at the individual LFE source locations on the Alpine Fault.  Finally, we use this new extensive catalogue as a tool for in-depth analyses of the deep central Alpine Fault structure and its slip behaviour. We identify eight episodes of increased LFE activity between 2009 and 2017 and provide time windows for further investigations of tremor and slow slip. We also study the spatial and temporal behaviours of LFEs and find that LFEs with synchronous occurrence patterns tend to be clustered in space. We thus suggest that individual LFE sources form spatially coherent clusters that may represent localised asperities or elastic patches on the deep Alpine Fault interface. We infer that those clusters may have a similar rheological response to tectonic forcing or to potential slow slip events. Eventually, we discover slow (10km/day) and rapid (∼20-25km/h) migrations of LFEs along the Alpine Fault. The slow migration might be controlled by slow slip events themselves while the rapid velocities could be explained by the LFE sources’ intrinsic properties.</p

GravProcess: An easy-to-use MATLAB software to process campaign gravity data and evaluate the associated uncertainties

International audienceWe present GravProcess, a set of MATLAB routines to process gravity data from complex campaign surveys and calculate the associated gravity field. Data reduction, analysis, and representation are done using the MATLAB Graphical User Interface Tool, which can be installed on most systems and platforms. Data processing is divided into several steps: (1) Integration of gravity data, station location, and gravity line connection input files; (2) Gravity data reduction applying solid-Earth tide and instrumental drift corrections and, depending on the required processing level, air pressure and oceanic tidal corrections; (3) Automatic network adjustment and alignment to absolute base stations; (4) Free air and terrain corrections to calculate gravity values and anomalies, and to estimate the associated errors. The final step is dedicated to post-processing and includes graphical representations of data and an output text file, which can be used by Geographic Information System software. An example of this processing chain applied to a recent survey in northern Morocco is given and compared with previous available results

Incipient mantle delamination, active tectonics and crustal thickening in Northern Morocco: Insights from gravity data and numerical modeling

International audienceThe Betic-Rif orocline surrounding the Alboran Sea, the westernmost tip of the Mediterranean Sea, accommodates the NW–SE convergence between the Nubia and Eurasia plates. Recent GPS observations indicate a ∼4 mm/yr SW motion of the Rif Mountains, relative to stable Nubia, incompatible with a simple two-plate model. New gravity data acquired in this study define a pronounced negative Bouguer anomaly south of the Rif, interpreted as a ∼40 km-thick crust in a state of non-isostatic equilibrium. We study the correlation between these present-day kinematic and geodynamic processes using a finite-element code to model in 2-D the first-order behavior of a lithosphere affected by a downward normal traction (representing the pull of a high-density body in the upper mantle). We show that intermediate viscosities for the lower crust and uppermost mantle (View the MathML source1021–1022Pas) allow an efficient coupling between the mantle and the base of the brittle crust, thus enabling (1) the conversion of vertical movement, resulting from the downward traction, to horizontal movement and (2) shortening in the brittle upper crust. Our results show that incipient delamination of the Nubian continental lithosphere, linked to slab pull, can explain the present-day abnormal tectonics, contribute to the gravity anomaly observed in northern Morocco, and give insight into recent tectonics in the Western Mediterranean region

Extreme hydrothermal conditions at an active plate-bounding fault

International audienceTemperature and fluid pressure conditions control rock deformation and mineralization on geological faults, and hence the distribution of earthquakes1. Typical intraplate continental crust has hydrostatic fluid pressure and a near-surface thermal gradient of 31 ± 15 degrees Celsius per kilometre2, 3. At temperatures above 300–450 degrees Celsius, usually found at depths greater than 10–15 kilometres, the intra-crystalline plasticity of quartz and feldspar relieves stress by aseismic creep and earthquakes are infrequent. Hydrothermal conditions control the stability of mineral phases and hence frictional–mechanical processes associated with earthquake rupture cycles, but there are few temperature and fluid pressure data from active plate-bounding faults. Here we report results from a borehole drilled into the upper part of the Alpine Fault, which is late in its cycle of stress accumulation and expected to rupture in a magnitude 8 earthquake in the coming decades4, 5. The borehole (depth 893 metres) revealed a pore fluid pressure gradient exceeding 9 ± 1 per cent above hydrostatic levels and an average geothermal gradient of 125 ± 55 degrees Celsius per kilometre within the hanging wall of the fault. These extreme hydrothermal conditions result from rapid fault movement, which transports rock and heat from depth, and topographically driven fluid movement that concentrates heat into valleys. Shear heating may occur within the fault but is not required to explain our observations. Our data and models show that highly anomalous fluid pressure and temperature gradients in the upper part of the seismogenic zone can be created by positive feedbacks between processes of fault slip, rock fracturing and alteration, and landscape development at plate-bounding faults

Bedrock geology of DFDP-2B, central Alpine Fault, New Zealand

During the second phase of the Alpine Fault, Deep Fault Drilling Project (DFDP) in the Whataroa River, South Westland, New Zealand, bedrock was encountered in the DFDP-2B borehole from 238.5–893.2 m Measured Depth (MD). Continuous sampling and meso- to microscale characterisation of whole rock cuttings established that, in sequence, the borehole sampled amphibolite facies, Torlesse Composite Terrane-derived schists, protomylonites and mylonites, terminating 200–400 m above an Alpine Fault Principal Slip Zone (PSZ) with a maximum dip of 62°. The most diagnostic structural features of increasing PSZ proximity were the occurrence of shear bands and reduction in mean quartz grain sizes. A change in composition to greater mica:quartz + feldspar, most markedly below c. 700 m MD, is inferred to result from either heterogeneous sampling or a change in lithology related to alteration. Major oxide variations suggest the fault-proximal Alpine Fault alteration zone, as previously defined in DFDP-1 core, was not sampled

Bedrock geology of DFDP-2B, central Alpine Fault, New Zealand

During the second phase of the Alpine Fault, Deep Fault Drilling Project (DFDP) in the Whataroa River, South Westland, New Zealand, bedrock was encountered in the DFDP-2B borehole from 238.5–893.2 m Measured Depth (MD). Continuous sampling and meso- to microscale characterisation of whole rock cuttings established that, in sequence, the borehole sampled amphibolite facies, Torlesse Composite Terrane-derived schists, protomylonites and mylonites, terminating 200–400 m above an Alpine Fault Principal Slip Zone (PSZ) with a maximum dip of 62°. The most diagnostic structural features of increasing PSZ proximity were the occurrence of shear bands and reduction in mean quartz grain sizes. A change in composition to greater mica:quartz + feldspar, most markedly below c. 700 m MD, is inferred to result from either heterogeneous sampling or a change in lithology related to alteration. Major oxide variations suggest the fault-proximal Alpine Fault alteration zone, as previously defined in DFDP-1 core, was not sampled