A gravity analysis of the Alpine Fault and the DFDP-2 drill site, Whataroa valley, South Westland, South Island, New Zealand

<p>The second phase of drilling into the Alpine Fault (DFDP-2), was completed in the Whataroa River valley, a former glacial valley located in central Westland, South Island, New Zealand. The site is located next to a steep hillside on the hanging-wall, ~1 km southeast of the mapped surface trace of the Alpine Fault. Projection of the hillside suggests a sediment thickness of 100 ± 40 m at the drill site; however, the sediment thickness was approximately double pre-drill estimates. Additionally, the surface expression and shallow geometry of the Alpine Fault in the Whataroa River valley, is not well-defined due to post-glacial burial of the fault zone. This thesis describes a gravity study designed to better constrain sub-surface structure beneath the DFDP-2 drill site and across the Alpine Fault. During this study, 466 new high-precision gravity observations were collected (standard error = 0.015 mGal) and amalgamated with 134 existing gravity stations, yielding comprehensive coverage of gravity data across the study area. A high density of observations was achieved within pre-determined zones, in addition to regional measurements so that residual gravity anomaly maps could be produced. The maps reveal: a negative residual gravity anomaly interpreted as a dextrally-offset glacial channel at least 350-450 m deep; steep localised gravity gradients near the Alpine Fault and DFDP-2 drill site that are interpreted as faulted and/or eroded boundaries; and a negative gravity anomaly adjacent to the DFDP-2 drill site that is interpreted as the deepest point of an over-deepened glacial lake. Gravity models were used to estimate the bedrock-sediment interface geometry near the DFDP-2 drill site and Alpine Fault. Structural inversion of the density boundary next to the drill site suggests either a moderately-dipping reverse fault or sub-vertical erosional wall exists beneath the hillside. Additional constraints on physical properties from direct density measurements or seismic velocity determinations and direct constraints on sediment thickness and layer geometry from seismic surveys will in future allow this new high-precision gravity dataset to be modelled more effectively.</p>

A gravity analysis of the Alpine Fault and the DFDP-2 drill site, Whataroa valley, South Westland, South Island, New Zealand

<p>The second phase of drilling into the Alpine Fault (DFDP-2), was completed in the Whataroa River valley, a former glacial valley located in central Westland, South Island, New Zealand. The site is located next to a steep hillside on the hanging-wall, ~1 km southeast of the mapped surface trace of the Alpine Fault. Projection of the hillside suggests a sediment thickness of 100 ± 40 m at the drill site; however, the sediment thickness was approximately double pre-drill estimates. Additionally, the surface expression and shallow geometry of the Alpine Fault in the Whataroa River valley, is not well-defined due to post-glacial burial of the fault zone. This thesis describes a gravity study designed to better constrain sub-surface structure beneath the DFDP-2 drill site and across the Alpine Fault. During this study, 466 new high-precision gravity observations were collected (standard error = 0.015 mGal) and amalgamated with 134 existing gravity stations, yielding comprehensive coverage of gravity data across the study area. A high density of observations was achieved within pre-determined zones, in addition to regional measurements so that residual gravity anomaly maps could be produced. The maps reveal: a negative residual gravity anomaly interpreted as a dextrally-offset glacial channel at least 350-450 m deep; steep localised gravity gradients near the Alpine Fault and DFDP-2 drill site that are interpreted as faulted and/or eroded boundaries; and a negative gravity anomaly adjacent to the DFDP-2 drill site that is interpreted as the deepest point of an over-deepened glacial lake. Gravity models were used to estimate the bedrock-sediment interface geometry near the DFDP-2 drill site and Alpine Fault. Structural inversion of the density boundary next to the drill site suggests either a moderately-dipping reverse fault or sub-vertical erosional wall exists beneath the hillside. Additional constraints on physical properties from direct density measurements or seismic velocity determinations and direct constraints on sediment thickness and layer geometry from seismic surveys will in future allow this new high-precision gravity dataset to be modelled more effectively.</p>

Improved earthquake detection as a probe for active fault structures in New Zealand's central Southern Alps

<p>This thesis concerns the detection and analysis of micro-seismicity and low-frequency earthquakes in New Zealand's central Southern Alps. We make use of the 6.5 year continuous seismic dataset collected using the Southern Alps Microearthquake Borehole Array (SAMBA), alongside other temporary and permanent seismic deployments nearby. The small station spacing of this deployment allows for high resolution seismic studies near the Alpine Fault, a dextral-transpressive plate boundary fault between the Pacific and Australian plates. Using this dataset we have documented the rst evidence of low-frequency earthquakes on or near the deep extent of the Alpine Fault. By using a network based crosscorrelation detection method we have generated a 3 year catalogue of 14 low-frequency earthquake families. These low-frequency earthquake families locate close to other indicators and models of the deep extent of the Alpine Fault, and we interpret these low-frequency earthquakes to represent shear failure on or near the deep extent of the Alpine Fault. These low-frequency earthquakes highlight a near-continuous background rate of deformation, punctuated by short periods of tremor. We also observe higher rates of low-frequency earthquake generation after large regional earthquakes. The magnitudes of our low-frequency earthquakes range from Mʟ‒0.8‒1.8, and appear to follow an exponential distribution, implying that there might be a characteristic length-scale of failure. We have extended the catalogue of low-frequency earthquake templates using the full 6.5 year dataset and an objective synthetic detection methodology. We developed a new methodology for template detection after other methods failed, or were not feasible. This method employs simple synthetic template events, which, rather than trying to capture all of the complexities of the body waves we try to detect, approximate a simple waveform that does not correlate well with background noise. To undertake this method we have developed a multi-parallel Python package, which is highly portable (we have run this on computers ranging from dual-core, 8GB RAM laptops to a 393 node, 6349 CPU cluster computer) and distributed via an open-source model. This package was run through the 6.5 year dataset on the New Zealand E-Science PAN cluster to e fficiently (<48 hours clock-time) generate a spatially and temporally continuous catalogue of low-frequency earthquake templates. Using this method to detect an initial suite of over 25,000 detections grouped into 600 families we have generated 600 good quality, discrete stacked waveforms for use in further matched-filter detection routines. We have shown that, for templates with both P and S-phase picks, these templates locate near to our previously determined low-frequency earthquake family locations. Using a network matched- filter detection technique we have generated a catalogue of micro-seismicity in a region of low-seismicity near the Whataroa Valley, motivated by the Deep-Fault Drilling Project; Phase-2. We detected 300 earthquakes that include a selection of near-repeating earthquakes. We find that most detected events are not similar enough to be termed repeating. For 106 earthquakes we are able to generate high-precision magnitudes calculated by singular-value decomposition of similar waveforms. We find a high b-value of 1.44 for these earthquakes, with no earthquakes above Mʟ1.6. By generating high precision cross-correlation derived picks for individual detections and employing a double-difference location methodology we show that seismicity does not delineate a single structure; rather we interpret the detected seismicity as temporally-limited earthquake sequences on small asperities adjacent to the Alpine Fault. Focal mechanisms for the best recorded events show dominantly strike-slip mechanisms, with lesser reverse and normal components. During the drilling of the Deep-Fault Drilling Project: Phase-2 borehole we operated a real-time earthquake detection system around the drill-site. This was a multi-national effort involving 16 seismologists in three countries monitoring the automatic detections in shifts. During the 5 month real-time monitoring period we detected and located 493 earthquakes, none of which occurred within 3km of the drill-site, nor required changes to the drilling operations. We undertook this monitoring using open-source software, which employed a standard energy based detection scheme. This thesis has contributed four complementary earthquake catalogues, a further three years of continuous seismic data from the central Southern Alps, and an opensource Python package for detection and analysis of earthquakes using cross-correlation techniques. The characteristics of these catalogues highlight deformation modes on and near one of the world's major strike-slip plate boundaries, both at depth, and at the upper extent of the seismogenic zone.</p>

Improved earthquake detection as a probe for active fault structures in New Zealand's central Southern Alps

<p>This thesis concerns the detection and analysis of micro-seismicity and low-frequency earthquakes in New Zealand's central Southern Alps. We make use of the 6.5 year continuous seismic dataset collected using the Southern Alps Microearthquake Borehole Array (SAMBA), alongside other temporary and permanent seismic deployments nearby. The small station spacing of this deployment allows for high resolution seismic studies near the Alpine Fault, a dextral-transpressive plate boundary fault between the Pacific and Australian plates. Using this dataset we have documented the rst evidence of low-frequency earthquakes on or near the deep extent of the Alpine Fault. By using a network based crosscorrelation detection method we have generated a 3 year catalogue of 14 low-frequency earthquake families. These low-frequency earthquake families locate close to other indicators and models of the deep extent of the Alpine Fault, and we interpret these low-frequency earthquakes to represent shear failure on or near the deep extent of the Alpine Fault. These low-frequency earthquakes highlight a near-continuous background rate of deformation, punctuated by short periods of tremor. We also observe higher rates of low-frequency earthquake generation after large regional earthquakes. The magnitudes of our low-frequency earthquakes range from Mʟ‒0.8‒1.8, and appear to follow an exponential distribution, implying that there might be a characteristic length-scale of failure. We have extended the catalogue of low-frequency earthquake templates using the full 6.5 year dataset and an objective synthetic detection methodology. We developed a new methodology for template detection after other methods failed, or were not feasible. This method employs simple synthetic template events, which, rather than trying to capture all of the complexities of the body waves we try to detect, approximate a simple waveform that does not correlate well with background noise. To undertake this method we have developed a multi-parallel Python package, which is highly portable (we have run this on computers ranging from dual-core, 8GB RAM laptops to a 393 node, 6349 CPU cluster computer) and distributed via an open-source model. This package was run through the 6.5 year dataset on the New Zealand E-Science PAN cluster to e fficiently (<48 hours clock-time) generate a spatially and temporally continuous catalogue of low-frequency earthquake templates. Using this method to detect an initial suite of over 25,000 detections grouped into 600 families we have generated 600 good quality, discrete stacked waveforms for use in further matched-filter detection routines. We have shown that, for templates with both P and S-phase picks, these templates locate near to our previously determined low-frequency earthquake family locations. Using a network matched- filter detection technique we have generated a catalogue of micro-seismicity in a region of low-seismicity near the Whataroa Valley, motivated by the Deep-Fault Drilling Project; Phase-2. We detected 300 earthquakes that include a selection of near-repeating earthquakes. We find that most detected events are not similar enough to be termed repeating. For 106 earthquakes we are able to generate high-precision magnitudes calculated by singular-value decomposition of similar waveforms. We find a high b-value of 1.44 for these earthquakes, with no earthquakes above Mʟ1.6. By generating high precision cross-correlation derived picks for individual detections and employing a double-difference location methodology we show that seismicity does not delineate a single structure; rather we interpret the detected seismicity as temporally-limited earthquake sequences on small asperities adjacent to the Alpine Fault. Focal mechanisms for the best recorded events show dominantly strike-slip mechanisms, with lesser reverse and normal components. During the drilling of the Deep-Fault Drilling Project: Phase-2 borehole we operated a real-time earthquake detection system around the drill-site. This was a multi-national effort involving 16 seismologists in three countries monitoring the automatic detections in shifts. During the 5 month real-time monitoring period we detected and located 493 earthquakes, none of which occurred within 3km of the drill-site, nor required changes to the drilling operations. We undertook this monitoring using open-source software, which employed a standard energy based detection scheme. This thesis has contributed four complementary earthquake catalogues, a further three years of continuous seismic data from the central Southern Alps, and an opensource Python package for detection and analysis of earthquakes using cross-correlation techniques. The characteristics of these catalogues highlight deformation modes on and near one of the world's major strike-slip plate boundaries, both at depth, and at the upper extent of the seismogenic zone.</p>

Marine Palynomorphs from the Plio-Pleistocene interval of the AND-1B Drill-Core McMurdo Sound, Antarctica

<p>The ANDRILL project recovered over 600 m of Plio-Pleistocene sediments within the Ross embayment, Antarctica. These sediments contain a record of local and regional paleoenvironmental conditions and glacial dynamism. They also provide a proxy for ice dynamics of the West Antarctic Ice Sheet (WAIS) during a period when global temperatures were ~3OC higher than modern. This unique record provides an analogue for future global climate change, which is expected to rise by 3OC by the end of the 21st century. Sixty-one samples from the upper 600 m of the AND-1B core were analysed for their palynomorph content yielding 4 to 5380 grains per sample (with an average frequency of 34 grains per gram). Marine palynomorphs including fossil dinoflagellate cysts, acritarchs, and prasinophyte algae were the focus of this study and fluctuations in their abundance and diversity reflect changes in paleoenvironment and glacial dynamics. The upper 600 m can be divided into 6 discrete units based on the palynomorph assemblage: The early-Pliocene (~5.0 – 4.6 Ma. Unit 1) is characterised by relatively high abundances of in situ round brown dinoflagellate cysts, microforaminiferal linings, and Leiosphaeridia, suggesting warmer than modern paleoenvironmental conditions and seasonal ice within the Ross embayment. The WAIS was likely small and highly dynamic during Unit 1. The mid-Pliocene (~4.6 – 3.4 Ma. Unit 2) exhibits relatively high abundances of round brown dinoflagellate cysts, microforaminiferal linings, and scolecodonts. The relatively low abundance of Leiosphaeridia (understood to indicate proximal/seasonal ice) suggests that ice free conditions at the drill site may have existed for up to ~1.2 Ma and that this may be the warmest period recorded in the core. During the warm, mid-Pliocene interval a sudden increase in scolecodonts (fossilized polychaete remains) may give indications into the water depth at the drill site because of their dependence on physical disturbance (decreasing with depth) for population growth. Further study of the scolecodonts is required before confident estimates of water depth can be made. The mid- to late-Pliocene (~3.4 – 2.6 Ma. Units 3, 4 & 5) is characterised by a variable palynomorph assemblage indicating variability in paleoenvironmental conditions, ice cover at the drill site, and ultimately a variable WAIS. A spike in the prasinophyte alga Cymatiosphaera (understood to indicate reduced salinity) at the base of a diatomite unit in the late-Pliocene may be a record of algae thriving in meltwater from the collapse of the WAIS. Further highresolution analysis is needed to help resolve this event. The Quaternary interval (~2.6 Ma and younger. Unit 6) is significantly different from previous units and is dominated by reworked Eocene dinoflagellate cysts and acritarchs (the “Transantarctic Flora”). This interval records a period of significant cooling and glacial expansion and the WAIS likely grew to its modern “polar” state. The WAIS may have undergone several collapses during super-interglacial periods in the Pleistocene but if it did it did not persist in its collapsed state for significant periods of time.</p>

Marine Palynomorphs from the Plio-Pleistocene interval of the AND-1B Drill-Core McMurdo Sound, Antarctica

<p>The ANDRILL project recovered over 600 m of Plio-Pleistocene sediments within the Ross embayment, Antarctica. These sediments contain a record of local and regional paleoenvironmental conditions and glacial dynamism. They also provide a proxy for ice dynamics of the West Antarctic Ice Sheet (WAIS) during a period when global temperatures were ~3OC higher than modern. This unique record provides an analogue for future global climate change, which is expected to rise by 3OC by the end of the 21st century. Sixty-one samples from the upper 600 m of the AND-1B core were analysed for their palynomorph content yielding 4 to 5380 grains per sample (with an average frequency of 34 grains per gram). Marine palynomorphs including fossil dinoflagellate cysts, acritarchs, and prasinophyte algae were the focus of this study and fluctuations in their abundance and diversity reflect changes in paleoenvironment and glacial dynamics. The upper 600 m can be divided into 6 discrete units based on the palynomorph assemblage: The early-Pliocene (~5.0 – 4.6 Ma. Unit 1) is characterised by relatively high abundances of in situ round brown dinoflagellate cysts, microforaminiferal linings, and Leiosphaeridia, suggesting warmer than modern paleoenvironmental conditions and seasonal ice within the Ross embayment. The WAIS was likely small and highly dynamic during Unit 1. The mid-Pliocene (~4.6 – 3.4 Ma. Unit 2) exhibits relatively high abundances of round brown dinoflagellate cysts, microforaminiferal linings, and scolecodonts. The relatively low abundance of Leiosphaeridia (understood to indicate proximal/seasonal ice) suggests that ice free conditions at the drill site may have existed for up to ~1.2 Ma and that this may be the warmest period recorded in the core. During the warm, mid-Pliocene interval a sudden increase in scolecodonts (fossilized polychaete remains) may give indications into the water depth at the drill site because of their dependence on physical disturbance (decreasing with depth) for population growth. Further study of the scolecodonts is required before confident estimates of water depth can be made. The mid- to late-Pliocene (~3.4 – 2.6 Ma. Units 3, 4 & 5) is characterised by a variable palynomorph assemblage indicating variability in paleoenvironmental conditions, ice cover at the drill site, and ultimately a variable WAIS. A spike in the prasinophyte alga Cymatiosphaera (understood to indicate reduced salinity) at the base of a diatomite unit in the late-Pliocene may be a record of algae thriving in meltwater from the collapse of the WAIS. Further highresolution analysis is needed to help resolve this event. The Quaternary interval (~2.6 Ma and younger. Unit 6) is significantly different from previous units and is dominated by reworked Eocene dinoflagellate cysts and acritarchs (the “Transantarctic Flora”). This interval records a period of significant cooling and glacial expansion and the WAIS likely grew to its modern “polar” state. The WAIS may have undergone several collapses during super-interglacial periods in the Pleistocene but if it did it did not persist in its collapsed state for significant periods of time.</p>

Brief communication: Estimating the ice thickness of the Müller Ice Cap using an inversion of the shallow ice approximation

The Müller Ice Cap will soon set the scene for a new drilling project. Therefore, ice thickness estimates are necessary for planning since thickness measurements of the ice cap are sparse. Here, two models are presented and compared, i) a simple inversion of the shallow ice approximation (SIA inversion) by the use of a single radar line in combination with the glacier outline, surface slope, and elevation, and ii) an iterative inverse method using the Parallel Ice Sheet Model (PISM). The two methods mostly agree about a good drill site candidate. However, the new semi-empirical SIA inversion is insensitive to mass balance, computationally fast, and provides better fits.

Ice core drilling on a high-elevation accumulation zone of Trambau Glacier in the Nepal Himalaya

We drilled an 81.2-m-long ice core in the accumulation area (5860 m a.s.l.) of Trambau Glacier in the Rolwaling region during October–November 2019. The drilling operation was conducted with a lightweight electro-mechanical drill system after two reconnaissance fieldworks in 2017 and 2018, during which two shallow firn cores were drilled with a hand auger. The drill system and ice core samples were transported by helicopters at a high elevation of 6000 m a.s.l. A further challenging issue was the ice core transportation between Nepal and Japan, as no regular commercial flight was available for the frozen samples. The addition of dry ice imported from India immediately prior to leaving Nepal allowed the ice core samples to be successfully transported to a cold room in Japan, and remain in a frozen state. Stratigraphic observations during the drilling operation suggest the drill site has been affected by melting and refreezing.

Indication of high basal melting at the EastGRIP drill site on the Northeast Greenland Ice Stream

The accelerated ice flow of ice streams that reach far into the interior of the ice sheets is associated with lubrication of the ice sheet base by basal meltwater. However, the amount of basal melting under the large ice streams – such as the Northeast Greenland Ice Stream (NEGIS) – is largely unknown. In situ measurements of basal melt rates are important from various perspectives as they indicate the heat budget, the hydrological regime and the relative importance of sliding in glacier motion. The few previous estimates of basal melt rates in the NEGIS region were 0.1 m a−1 and more, based on radiostratigraphy methods. These findings raised the question of the heat source, since even an increased geothermal heat flux could not deliver the necessary amount of heat. Here, we present basal melt rates at the recent deep drill site EastGRIP, located in the centre of NEGIS. Within 2 subsequent years, we found basal melt rates of 0.19±0.04 m a−1 that are based on analysis of repeated phase-sensitive radar measurements. In order to quantify the contribution of processes that contribute to melting, we carried out an assessment of the energy balance at the interface and found the subglacial water system to play a key role in facilitating such high melt rates.

Surface velocity and ice thickness of the Müller ice cap, Axel Heiberg Island

Muller ice cap is situated on Axel Heiberg Island in Arctic Canada. It is characterised by a mountanious region separating the ice cap in the east from the outlet glaciers in the west. Research has taken place on the outlet glaciers of the ice cap since 1959, but only limited research has been conducted on the main ice cap, and no full depth ice cores have ever been drilled. The interesting location of the ice cap facing the Arctic Ocean, the chance of finding ice dating back to the Innuitian Ice Sheet, and the fact that no full depth ice cores have been drilled, makes it an obvious place to do so. In order to achieve a long and undisturbed chronology of the ice core, one needs to find a location for the drilling site, where there is a great ice thickness, low surface velocity and little melt. In this project the aim is to make surface velocity maps of the ice cap and estimate the ice thickness to be able to come up with suggestions of possible drill site areas.Surface velocities are calculated using feature tracking of optical satellite images from the Landsat satellites in the period of 2004-2019. A median velocity map of all Landsat 8 velocity maps is used as validation in modelling the ice thickness and in the investigation of possible drill site areas.To estimate the ice thickness various methods are used and are being compared to the ice thickness measured by Operation IceBridge. The first method is an iterative inverse method where the ice sheet model PISM works as a forward model. The model is found to work rather well on the ice cap, with a root mean squared error (RMS) of 138.9 m, but overestimates the ice thickness on the outlet glaciers. The second model uses a simple inversion of the shallow ice approximation. It overestimates the ice thickness in areas with low surface slope, but has a RMS of 131.4 m on the ice cap. The third and and fourth models uses Monte Carlo sampling methods of the shallow ice approximation without and with sliding, respectively. The latter uses an initial ice thickness guess, and the modelled ice thickness was proofed not to differ from that initial guess at all. The RMSs on the icecap of the two models were found to be 132.1 m and 129.9 m, respectively. Finally, the fifth model uses the PISM setup but with an initial geometry defined by the SIA inversion. The RMS on the ice cap is found to be 135.4 m.Based on the median Landsat 8 surface velocity map, the modelled ice thicknesses and the surface elevation from the Arctic Digital Elevation Model, a map of suggested drill site areas is made. The site which fulfilled the criteria the most is located at 526629 m easting and 8866463 m northing in UTM zone 15N. In this site the surface velocity is 1.2 m yr−1, the surface elevation is 1804 m and the modelled ice thicknesses varies from 535-579 m. The melt in this area is estimated to be less than 20 melt days per year based on the backscatter from Sentinel-1.