Hikurangi Subduction Zone Research Articles

Engineering geomorphological reconnaissance of the December 2018 Waimata Valley mud volcano eruption, Gisborne, New Zealand

Gisborne (North Island, NZ) is affected by rainfall-induced landslides, earthquakes and tsunami, as well as mud volcanoes (MVs). The latter form via upward mobilization of Eocene–age sediments, and have not been studied from an engineering geological standpoint, so the 15 December 2018 Waimata Valley MV eruption provided a unique opportunity. The event erupted c. 16 900 m 3 of mud, forming an elevated vent area, and three mudflows propagating north, east and south. Scanning electron microscopy (SEM) identified smectite as the dominant clay in the Sr-rich mud, and Atterberg limits indicate high plasticity. In-situ testing using dynamic cone penetrometer and shear vane (3–168 kPa) revealed wide variability in strength properties with depth, while ring shear values are 11.3–13.5°. A fault extends NW beneath the Waimata Valley MV, coinciding with the pre-existing Arakihi Road MV. The Waimata Valley MV area was subject to uplift and cracking during the September 2016 Te Araroa earthquake ( M w 7.1), which caused increased activity at pre-existing mud volcanoes at that time. Geodetic data for the Gisborne district shows an uplift phase culminating around the December 15 2018 MV eruption, followed by the commencement of a Hikurangi subduction zone ‘slow slip event’. Nevertheless, relationships between tectonics and MV eruptions remain equivocal.

Re-evaluation of New Zealand seismic hazard for geotechnical assessment and design

This paper scrutinises the seismic hazard of New Zealand (NZ) from a geotechnical engineering perspective. The two codified versions of the seismic hazard in NZS1170.5 (structural loading standard) and NZTA Bridge Manual (NZTA-BM) are shown to yield consistently different peak ground acceleration (PGA) hazards throughout NZ. Results from site-specific PSHA for 24 locations in NZ are used to examine key hazard characteristics, including earthquake magnitude and the effects of site conditions (classes) on the PGA hazard. The comparative evaluations show that for most of the locations considered, NZS1170.5 and NZTA-BM overestimate the PGA hazard. However, NZS1170.5, and NZTA-BM in particular, significantly underestimate the PGA hazard for locations that are at a short source-to-site distance from the Hikurangi Subduction Zone (HSZ), and for which HSZ significantly contributes to their hazard. Using the results from this study, an interim PGA hazard is recommended for geotechnical assessment and design in support of the NZ guidelines for geotechnical earthquake engineering practice. The recommended interim PGA hazard is applicable to all site classes without any modification or use of site amplification factors.

Crustal Structure of the Hikurangi Margin From SHIRE Seismic Data and the Relationship Between Forearc Structure and Shallow Megathrust Slip Behavior

AbstractMarine multichannel and wide‐angle seismic data constrain crustal structure along a 530 km margin‐parallel transect of the Hikurangi subduction zone. The subducting Hikurangi Plateau crust (VP 5.0–7.4 km/s) is ∼1 km thicker (11 ± 1 km) and mantle velocities are ∼0.2 km/s higher (VP 8.3–8.5 km/s) beneath south/central Hikurangi relative to north Hikurangi. In the overthrusting plate, an abrupt 0.5 km/s south‐to‐north reduction in forearc wavespeeds occurs in concert with a change in seismic reflection character. We analyze legacy seismic data to show that the forearc transition likely reflects lateral variability in the updip extent of the Torlesse Backstop. Furthermore, we map this unit along‐strike and note a broad correlation between the backstop and down‐dip extent of shallow slow‐slip. We propose that the geological architecture of the overthrusting plate contributes to spatial variability in the location of shallow frictional transitions along the Hikurangi margin, impacting both seismic and tsunami hazard.

Strong Upper‐Plate Heterogeneity at the Hikurangi Subduction Margin (North Island, New Zealand) Imaged by Adjoint Tomography

AbstractWe use earthquake‐based adjoint tomography to invert for three‐dimensional structure of the North Island, New Zealand, and the adjacent Hikurangi subduction zone. The study area, having a shallow depth to the plate interface below the North Island, offers a rare opportunity for imaging material properties at an active subduction zone using land‐based measurements. Starting from an initial model derived using ray tomography, we perform iterative model updates using spectral element and adjoint simulations to fit waveforms with periods ranging from 4–30 s. We perform 28 model updates using an L‐BFGS optimization algorithm, improving data fit and introducing P‐ and S‐wave velocity changes of up to ±30%. Resolution analysis using point spread functions show that our measurements are most sensitive to heterogeneities in the upper 30 km. The most striking velocity changes coincide with areas related to the active Hikurangi subduction zone. Lateral velocity structures in the upper 5 km correlate well with New Zealand geology. The inversion reveals increased along‐strike heterogeneity on the margin. In Cook Strait we observe a low‐velocity zone interpreted as deep sedimentary basins. In the central North Island, low‐velocity anomalies are linked to surface geology, and we relate velocity structures at depth to crustal magmatic activity below the Taupō Volcanic Zone. Our velocity model provides more accurate synthetic seismograms with respect to the initial model, better constrains small (50 km), shallow (15 km) and near‐offshore velocity structures, and improves our understanding of volcanic and tectonic structures related to the active Hikurangi subduction zone.

Segmentation of Shallow Slow Slip Events at the Hikurangi Subduction Zone Explained by Along‐Strike Changes in Fault Geometry and Plate Convergence Rates

Over the last two decades, geodetic and seismic observations have revealed a spectrum of slow earthquakes along the Hikurangi subduction zone in New Zealand. Of those, shallow slow slip events (SSEs) that occur at depths of less than 15 km along the plate interface show a strong along‐strike segmentation in their recurrence intervals, which vary from ∼1 yr from offshore Tolaga Bay in the northeast to ∼5 yr offshore Cape Turnagain ∼300 km to the southwest. To understand the factors that control this segmentation, we conduct numerical simulations of SSEs incorporating laboratory‐derived rate‐and‐state friction laws with both planar and non‐planar fault geometries. We find that a relatively simple model assuming a realistic non‐planar fault geometry reproduces the characteristics of shallow SSEs as constrained by geodetic observations. Our preferred model captures the magnitudes and durations of SSEs, as well as the northward decrease of their recurrence intervals. Our results indicate that the segmentation of SSE recurrence intervals is favored by along‐strike changes in both the plate convergence rate and the downdip width of the SSE source region. Modeled SSEs with longer recurrence intervals concentrate in the southern part of the fault (offshore Cape Turnagain), where the plate convergence rate is lowest and the source region of SSEs is widest due to the shallower slab dip angle. Notably, the observed segmentation of shallow SSEs cannot be reproduced with a simple planar fault model, which indicates that a realistic plate interface is an important factor to account for in modeling SSEs.

Evidence for Deeply Subducted Lower‐Plate Seamounts at the Hikurangi Subduction Margin: Implications for Seismic and Aseismic Behavior

AbstractSeamounts are found at many subduction zones and act as seafloor heterogeneities that affect slip behavior on megathrusts. At the Hikurangi subduction zone offshore the North Island, New Zealand, seamounts have been identified on the incoming Pacific plate and below the accretionary prism, but there is little concrete evidence for seamounts subducted beyond the present‐day coastline. Using a high‐resolution, adjoint tomography‐derived velocity model of the North Island, we identify two high‐velocity anomalies below the East Coast and an intraslab low‐velocity zone up‐dip of one of these anomalies. We interpret the high‐velocity anomalies as previously unidentified, deeply subducted seamounts, and the low‐velocity zone as fluid in the subducting slab. The seamounts are inferred to be 10–30 km wide and on the plate interface at 12–15 km depth. Resolution analysis using point spread functions confirms that these are well‐resolved features. The locations of the two seamounts coincide with bathymetric features whose geometries are consistent with those predicted from analog experiments and numerical simulations of seamount subduction. The spatial characteristics of seismicity and slow slip events near the inferred seamounts agree well with previous numerical modeling predictions of the effects of seamount subduction on megathrust stress and slip. Anomalous geophysical signatures, magnetic anomalies, and swarm seismicity have also been observed previously at one or both seamount locations. We propose that permanent fracturing of the northern Hikurangi upper plate by repeated seamount subduction may be responsible for the dichotomous slow slip behavior observed geodetically, and partly responsible for along‐strike variations in plate coupling on the Hikurangi subduction interface.

Seismic Constraint on Heterogeneous Deformation and Stress State in the Forearc of the Hikurangi Subduction Zone, New Zealand

Abstract The Hikurangi subduction zone (HSZ) is the collisional boundary between the Pacific and Australian tectonic plates along the eastern coast of the North Island of New Zealand. The region is believed to be capable of hosting large megathrust earthquakes and associated tsunamis. Recent studies observe a range of slip behavior along the plate interface, with a sharp contrast between locked and creeping parts of the megathrust along the margin. This work uses teleseismic scattering data (receiver functions [RFs]) recorded at 53 long-running seismograph stations on the North Island of New Zealand to constrain the structure and mechanical properties of the forearc in the HSZ. We observe directional variations in RF phases at P–S converted delay times (i.e., depths) associated with the overlying forearc crust and note a general correlation with spatial variations in plate coupling as well as other geophysical properties. Our results suggest differences in the nature of crustal deformation (and stress state) along the Hikurangi margin, with evidence of clockwise rotation and/or extension in the northern HSZ, where the overriding forearc crust is uncoupled from the subducting Pacific slab.

Origin and diagenetic priming of a potential slow-slip trigger zone in volcaniclastic deposits flanking a seamount on the subducting plate, Hikurangi margin, New Zealand

ABSTRACT Slow-slip megathrust events in the Hikurangi forearc region may be related to seamount subduction. IODP Expedition 375 drilled the flank of Tūranganui Knoll, a Cretaceous seamount east of the subduction zone containing units analogous to intervals associated with subducted seamounts. Detailed investigations of the volcaniclastic Unit V at Site U1520, based on scanned core images, petrographic observations, point-counting, and μXRF mapping, reveal that it is a mafic volcaniclastic succession rich in clay minerals which are products of a multiphased history of diagenetic alteration. This unit is interpreted as resulting from downslope transport of volcaniclastic debris from submarine eruptions. In contrast, most volcanic clasts recovered at the nearby Site U1526, near the top of Tūranganui Knoll, are interpreted as products of subaerial eruptions reworked by wave action. This shallow water sedimentary package experienced minimal clay-mineral alteration and more extensive early carbonate cementation. Differences in the clast origins, depositional settings and diagenetic histories of these volcaniclastic units highlight the heterogeneity of volcanic systems on the Hikurangi Plateau. However, the presence of voluminous hydrous clays within thick, altered flanking volcaniclastic units at Site U1520 provides a ready source for excess pore fluids that may enable slow-slip events along the Hikurangi subduction zone.

Precise Monitoring of Pore Pressure at Boreholes Around Nankai Trough Toward Early Detecting Crustal Deformation

In our recent study, we detected the pore pressure change due to the slow slip event (SSE) in March 2020 at the two borehole stations (C0002 and C0010), where the other borehole (C0006) close to the Nankai Trough seems not because of instrumental drift for the reference pressure on the seafloor to remove non-crustal deformation such as tidal and oceanic fluctuations. To overcome this problem, we use the seafloor pressure gauges of cabled network Dense Oceanfloor Network System for Earthquakes and Tsunamis (DONET) stations nearby boreholes instead of the reference by introducing time lag between them. We confirm that the time lag is explained from superposition of theoretical tide modes. By applying this method to the pore pressure during the SSE, we find pore pressure change at C0006 about 0.6 hPa. We also investigate the impact of seafloor pressure due to ocean fluctuation on the basis of ocean modeling, which suggests that the decrease of effective normal stress from the onset to the termination of the SSE is explained by Kuroshio meander and may promote updip slip migration, and that the increase of effective normal stress for the short-term ocean fluctuation may terminate the SSE as observed in the Hikurangi subduction zone.

Evidence of Seismic Slip on a Large Splay Fault in the Hikurangi Subduction Zone

AbstractThe Hikurangi subduction zone is capable of producing moderate to large earthquakes as well as regularly repeating slow slip events. However, it is unclear what structures host these different slip styles along the margin. Here we address whether splay faults can host seismic slip at shallow (<1 km) depth by investigating the Pāpaku fault, sampled during International Ocean Discovery Program Expedition 375. We use biomarker thermal maturity to search for evidence of frictional heating within turbiditic sediments of the Pāpaku fault. Four zones of localized high temperature are found near the top of the fault zone, which are interpreted to be zones of localized seismic slip. Thermal modeling shows that the most likely maximum displacement on the shallow Pāpaku fault during each event was 14–17 m. Our results demonstrate that the Pāpaku fault, and potentially other splay faults along the margin, host coseismic slip and have the potential to produce large tsunami (e.g., runup heights of >1 m as observed in the 1947 Poverty and Tolaga Bay earthquakes.

Nested regional-global seismic tomography and precise earthquake relocation along the Hikurangi subduction zone, New Zealand

SUMMARY Seismic and geodetic examinations of the Hikurangi subduction zone (HSZ) indicate a remarkably diverse and complex system. Here, we investigate the 3-D P-wave velocity structure of the HSZ by applying an iterative, nested regional-global tomographic algorithm. The new model reveals enhanced details of seismic variations along the HSZ. We also relocate over 57 000 earthquakes using this newly developed 3-D model and then further improve the relative locations for 75 per cent of the seismicity using waveform cross-correlation. Double seismic zone characteristics, including occurrence, depth distribution and thickness change along the strike of the HSZ. An aseismic but fast Vp zone separates the upper and lower planes of seismicity in the southern and northern North Island. The upper plane of seismicity correlates with low Vp zones below the slab interface, indicating fluid-rich channels formed on top and/or within a dehydrated crust. A broad low Vp zone is resolved in the lower part of the subducting slab that could indicate hydrous mineral breakdown in the slab mantle. In the northern North Island and southern North Island, the lower plane of seismicity mostly correlates with the top of these low Vp zones. The comparison between the thermal model and the lower plane of seismicity in the northern North Island supports dehydration in the lower part of the slab. The mantle wedge of the Taupo volcanic zone (TVZ) is characterized by a low velocity zone underlying the volcanic front (fluid-driven partial melting), a fast velocity anomaly in the forearc mantle (a stagnant cold nose) and an underlying low velocity zone within the slab (fluids from dehydration). These arc-related anomalies are the strongest beneath the central TVZ with known extensive volcanism. The shallow seismicity (<40 km depth) correlates with geological terranes in the overlying plate. The aseismic impermeable terranes, such as the Rakaia terrane, may affect the fluid transport at the plate interface and seismicity in the overlying plate, which is consistent with previous studies. The deep slow slip events (25–60 km depths) mapped in the Kaimanawa, Manawatu and Kapiti regions coincide with low Vp anomalies. These new insights on the structure along the HSZ highlight the change in the locus of seismicity and dehydration at depth that is governed by significant variations in spatial and probably temporal attributes of subduction zone processes.

Paleotsunamis on the Southern Hikurangi Subduction Zone, New Zealand, Show Regular Recurrence of Large Subduction Earthquakes

Abstract Geological records of subduction earthquakes, essential for seismic and tsunami hazard assessment, are difficult to obtain at transitional plate boundaries, because upper-plate fault earthquake deformation can mask the subduction zone signal. Here, we examine unusual shell layers within a paleolagoon at Lake Grassmere, at the transition zone between the Hikurangi subduction zone and the Marlborough fault system. Based on biostratigraphic and sedimentological analyses, we interpret the shell layers as tsunami deposits. These are dated at 2145–1837 and 1505–1283 yr B.P., and the most likely source of these tsunamis was ruptures of the southern Hikurangi subduction interface. Identification of these two large earthquakes brings the total record of southern Hikurangi subduction earthquakes to four in the past 2000 yr. For the first time, it is possible to obtain a geologically constrained recurrence interval for the southern Hikurangi subduction zone. We calculate a recurrence interval of 500 yr (335–655 yr, 95% confidence interval) and a coefficient of variation of 0.27 (0.0–0.47, 95% confidence interval). The probability of a large subduction earthquake on the southern Hikurangi subduction zone is 26% within the next 50 yr. We find no consistent temporal relationship between subduction earthquakes and large earthquakes on upper-plate faults.

Fluid pressurisation and earthquake propagation in the Hikurangi subduction zone

In subduction zones, seismic slip at shallow crustal depths can lead to the generation of tsunamis. Large slip displacements during tsunamogenic earthquakes are attributed to the low coseismic shear strength of the fluid-saturated and non-lithified clay-rich fault rocks. However, because of experimental challenges in confining these materials, the physical processes responsible for the coseismic reduction in fault shear strength are poorly understood. Using a novel experimental setup, we measured pore fluid pressure during simulated seismic slip in clay-rich materials sampled from the deep oceanic drilling of the Pāpaku thrust (Hikurangi subduction zone, New Zealand). Here, we show that at seismic velocity, shear-induced dilatancy is followed by pressurisation of fluids. The thermal and mechanical pressurisation of fluids, enhanced by the low permeability of the fault, reduces the energy required to propagate earthquake rupture. We suggest that fluid-saturated clay-rich sediments, occurring at shallow depth in subduction zones, can promote earthquake rupture propagation and slip because of their low permeability and tendency to pressurise when sheared at seismic slip velocities.

Water Depth Dependence of Long‐Range Correlation in Nontidal Variations in Seafloor Pressure

AbstractIsolating the source of nontidal oceanographic noise in seafloor pressure data is critical for improving the use of these data for seafloor geodetic applications. Residuals between nearby bottom pressure records have typically been used to remove the nontidal components, as these are largely common‐mode. To evaluate the similarities between pairs of observed bottom pressure records at a range of water depths, we calculate the standard deviations of the time series of residuals between data from all site pairs, installed in the Hikurangi subduction zone offshore New Zealand. We find that the magnitude of the standard deviation depends more on relative water depth than the distance between sites. This confirms the result of previous studies from Cascadia that nontidal components are more similar along isobaths even if the distance between sites is large. Furthermore, in order to reduce noises, the required depth difference between site pairs also varies with site depths.

Controls on the origin and evolution of deep-ocean trench-axial channels

The type and volume of sediment entering subduction zones affects the style of plate-boundary deformation and thus sedimentary and tectonic cycles. Because submarine channels significantly increase the transport efficiency of turbidity currents, their presence or absence in subduction trenches is a primary control on trench fill. To date, comprehensive architectural characterization of trench-axial channels has not been possible, undermining efforts to identify the factors controlling their initiation and evolution. Here, we describe the evolution of the Hikurangi Channel, which traverses the Hikurangi Trench, offshore New Zealand. Analysis of two- and three-dimensional seismic data reveals that the channel was present only during the last ~3.5 m.y. of the ~27 m.y. of the trench’s existence; its inception and propagation resulted from increased sediment supply to the trench following amplified hinterland exhumation. To test if the controls on the evolution of the Hikurangi Channel are universal, multivariate statistical analysis of the geomorphology of subduction trenches globally is used to investigate the formative conditions of axial channels in modern trenches. Terrigenous sediment supply and thickness of sediment cover in a trench are the dominant controls; subsidiary factors such as trench length and rugosity also contribute to the conditions necessary for trench-axial channel development. Axial channels regulate sediment distribution in trenches, and this varies temporally and spatially as a channel propagates along a trench. The presence of a trench-axial channel affects plate-boundary mechanics and has implications for the style of subduction-margin deformation.

Seismicity and velocity structure in the vicinity of repeating slow slip earthquakes, northern Hikurangi subduction zone, New Zealand

Using an 11-month (April 2012 to March 2013) deployment of Ocean Bottom Seismographs along the northern Hikurangi subduction zone, New Zealand, we located 2005 earthquakes (M 0.7-4.5) and obtained detailed 3-D Vp and Vs models for a region where a spectrum of subduction fault slip behavior is observed. Our analysis over a 120 × 110 km2 onshore-offshore area reveals high Vp/Vs (>1.8) in the subducting 10 km thick Hikurangi Plateau oceanic crust below an area of repeating slow slip events (SSEs), that are assumed to occur on the plate megathrust. Low Vp/Vs is mapped in a 15 km wide domain along the coastline, offshore Gisborne, marked by low Qs/Qp. Seismicity is absent from the core of high Vp/Vs anomalies (∼2.0) and within the band with lower Vp/Vs, but earthquakes are distributed vertically and are widespread in adjacent high Vp/Vs (∼1.85-1.95) regions along the trend of the coastline. Down-dip and along-strike variation of the Vp/Vs structure of the subducting Hikurangi Plateau could be the result of dehydration process caused by reactivated bending-related faults. Zones of vertical extensional fault–fracture in the subducting crust may act as conduits for the discharge of overpressured fluid and episodic hydrofracturing. During our deployment a large SSE (Mw6.8) occurred in February 2013. In response, seismicity became active in central Hawke Bay and at the down-dip edge of a previously mapped high-amplitude reflectivity zone. Microseismicity appeared synchronous with the slip propagation of the SSEs and may have been initiated by the slip. Our results hold implications for understanding of the physical environments conducive to slow slip generation and represent an important baseline for future seismicity work in the region.

Near Trench 3D Seismic Attenuation Offshore Northern Hikurangi Subduction Margin, North Island, New Zealand

AbstractWe image seismic attenuation near the Hikurangi trench offshore New Zealand, using ocean bottom and land‐based seismometers, revealing high attenuation above a recurring shallow slow‐slip event and within the subducting Hikurangi Plateau. The Hikurangi subduction margin east of the North Island, New Zealand is the site of frequent shallow slow slip events. Overpressured fluids are hypothesized to lead to slow slip at shallow depths close to the oceanic trench. Seismic attenuation, energy loss of seismic waves, can be used to detect high temperatures, melt, the presence of fluids, and fractures. We use local earthquake P‐ and S‐waves from 180 earthquakes to invert for t*, and subsequently invert for Qp and Qs, offshore the North Island directly above the area of slow slip. We image Qp and Qs to ∼25 km depth, increasing resolution of previously identified coastal low Q (100–300), and finding a new region of even higher attenuation (Qp and Qs < 50–100) directly above the shallow slow slip event of 2014–2015, beneath the offshore seismic array. This highest attenuation is downdip of a subducting seamount, and is spatially correlated with a high seismic reflectivity zone and Vp/Vs > 1.85, all of which provide evidence for the presence of fluids. The Qp and Qs is low at the trench (<50–100) and in the subducting plate (100–200), suggesting that seismic wave scattering due to faults, fractures, and the inherent heterogeneous composition of the Hikurangi Plateau, a large igneous province, plays a role in seismic attenuation.

Distribution, microphysical properties, and tectonic controls of deformation bands in the Miocene subduction wedge (Whakataki Formation) of the Hikurangi subduction zone

Abstract. We analyse deformation bands related to horizontal contraction with an intermittent period of horizontal extension in Miocene turbidites of the Whakataki Formation south of Castlepoint, Wairarapa, North Island, New Zealand. In the Whakataki Formation, three sets of cataclastic deformation bands are identified: (1) normal-sense compactional shear bands (CSBs), (2) reverse-sense CSBs, and (3) reverse-sense shear-enhanced compaction bands (SECBs). During extension, CSBs are associated with normal faults. When propagating through clay-rich interbeds, extensional bands are characterised by clay smear and grain size reduction. During contraction, sandstone-dominated sequences host SECBs, and rare CSBs, that are generally distributed in pervasive patterns. A quantitative spacing analysis shows that most outcrops are characterised by mixed spatial distributions of deformation bands, interpreted as a consequence of overprint due to progressive deformation or distinct multiple generations of deformation bands from different deformation phases. As many deformation bands are parallel to adjacent juvenile normal faults and reverse faults, bands are likely precursors to faults. With progressive deformation, the linkage of distributed deformation bands across sedimentary beds occurs to form through-going faults. During this process, bands associated with the wall-, tip-, and interaction-damage zones overprint earlier distributions resulting in complex spatial patterns. Regularly spaced bands are pervasively distributed when far away from faults. Microstructural analysis shows that all deformation bands form by inelastic pore collapse and grain crushing with an absolute reduction in porosity relative to the host rock between 5 % and 14 %. Hence, deformation bands likely act as fluid flow barriers. Faults and their associated damage zones exhibit a spacing of 9 m on the scale of 10 km and are more commonly observed in areas characterised by higher mudstone-to-sandstone ratios. As a result, extensive clay smear is common in these faults, enhancing the sealing capacity of faults. Therefore, the formation of deformation bands and faults leads to progressive flow compartmentalisation from the scale of 9 m down to about 10 cm – the typical spacing of distributed, regularly spaced deformation bands.

Improving Detectability of Seafloor Deformation From Bottom Pressure Observations Using Numerical Ocean Models

We investigated ocean bottom pressure (OBP) observation data at six plate subduction zones around the Pacific Ocean. The six regions included the Hikurangi Trough, the Nankai Trough, the Japan Trench, the Aleutian Trench, the Cascadia Subduction Zone, and the Chile Trench. For the sake of improving the detectability of seafloor deformation using OBP observations, we used numerical ocean models to represent realistic oceanic variations, and subtracted them from the observed OBP data. The numerical ocean models included four ocean general circulation models (OGCMs) of HYCOM, GLORYS, ECCO2, and JCOPE2M, and a single-layer ocean model (SOM). The OGCMs are mainly driven by the wind forcing. The SOM is driven by wind and/or atmospheric pressure loading. The modeled OBP was subtracted from the observed OBP data, and root-mean-square (RMS) amplitudes of the residual OBP variations at a period of 3–90 days were evaluated by the respective regions and by the respective numerical ocean models. The OGCMs and SOM driven by wind alone (SOMw) contributed to 5–27% RMS reduction in the residual OBP. When SOM driven by atmospheric pressure alone (SOMp) was added to the modeled OBP, residual RMS amplitudes were additionally reduced by 2–15%. This indicates that the atmospheric pressure is necessary to explain substantial amounts of observed OBP variations at the period. The residual RMS amplitudes were 1.0–1.7 hPa when SOMp was added. The RMS reduction was relatively effective as 16–42% at the Hikurangi Trough, the Nankai Trough, and the Japan Trench. The residual RMS amplitudes were relatively small as 1.0–1.1 hPa at the Nankai Trough and the Chile Trench. These results were discussed with previous studies that had identified slow slips using OBP observations. We discussed on further accurate OBP modeling, and on improving detectability of seafloor deformation using OBP observation arrays.

GEBCO and ETOPO1 gridded datasets for GMT based cartographic Mapping of Hikurangi, Puysegur and Hjort Trenches, New Zealand

The study focused on the comparative analysis of the submarine geomorphology of three oceanic trenches: Hikurangi Trench (HkT), Puysegur Trench (PT) and Hjort Trench (HjT), New Zealand region, Pacific Ocean. HjT is characterized by an oblique subduction zone. Unique regional tectonic setting consist in two subduction zones: northern (Hikurangi margin) and southern (Puysegur margin), connected by oblique continental collision along the Alpine Fault, South Island. This cause variations in the geomorphic structure of the trenches. PT/HjT subduction is highly oblique (dextral) and directed southwards. Hikurangi subduction is directed northwestwards. South Island is caught in between by the “subduction scissor”. Methodology is based on GMT (The Generic Mapping Tools) for mapping, plotting and modelling. Mapping includes visualized geophysical, tectonic and geological settings of the trenches, based on sequential use of GMT modules. Data include GEBCO, ETOPO1, EGM96. Comparative histogram equalization of topographic grids (equalized, normalized, quadratic) was done by module ’grdhisteq’, automated cross-sectioning – by ’grdtrack’. Results shown that HjT has a symmetric shape form with comparative gradients on both western and eastern slopes. HkT has a trough-like flat wide bottom, steeper gradient slope on the North Island flank. PT has an asymmetric V-form with steep gradient on the eastern slopes and gentler western slope corresponding to the relatively gentle slope of a subducting plate and steeper slope of an upper one. HkT has shallower depths < 2,500 m, PT is <-6,000 m. The deepest values > 6,000 m for HjT. The surrounding relief of the HjT presents the most uneven terrain with gentle slope oceanward, and a steep slope on the eastern flank for PT, surrounded by complex submarine relief along the Macquarie Arc. Data distribution for the HkT demonstrates almost equal pattern for the depths from -600 m to ₋2,600 m. PT has a bimodal data distribution with 2 peaks: 1) -4,250 to -4,500 m (18%); 2) -2,250 to -3,000 m, < 7,5%. The second peak corresponds to the Macquarie Arc. Data distribution for HjT is classic bell-shaped with a clear peak at -3,250 to -3,500 m. The asymmetry of the trenches resulted in geomorphic shape of HkT, PT and HjT affected by geologic processes.