New Zealand has experienced several strong earthquakes in its history. While an earthquake cannot be prevented from occurring, planning can reduce its consequences when it does occur. This dissertation research examines various aspects of disaster risk management policy in Aotearoa New Zealand. Chapter 2 develops a method to rank and prioritise high-rise buildings for seismic retrofitting in Wellington, the earthquake-prone capital city of New Zealand. These buildings pose risks to Wellington’s long-term seismic resilience that are of clear concern to current and future policymakers. The prioritization strategy we propose, based on multi-criteria decision analysis (MCDA) methods, considers a variety of data on each building, including not only its structural characteristics, but also its location, its economic value to the city, and its social importance to the community around it. The study demonstrates how different measures, within four general criteria – life safety, geo-spatial location of the building, its economic role, and its socio-cultural role – can be operationalized into a viable framework for determining retrofitting/demolition policy priorities. Chapter 3 and chapter 4 analyse the Residential Red Zone (RRR) program that was implemented in Christchurch after the 2011 earthquake. In the program, approximately 8,000 homeowners were told that their homes were no longer permittable, and they were bought by the government (through the Canterbury Earthquake Recovery Authority). Chapter 3 examines the subjective wellbeing of the RRR residents (around 16000 people) after they were forced to move. We consider three indicators of subjective wellbeing: quality of life, stress, and emotional wellbeing. We found that demographic factors, health conditions, and the type of government compensation the residents accepted, were all significant determinants of the wellbeing of the Red Zone residents. More social relations, better financial circumstances, and the perception of better government communication were also all associated positively with a higher quality of life, less stress, and higher emotional wellbeing. Chapter 4 concentrates on the impact of this managed retreat program on RRR residents’ income. We use individual-level comprehensive, administrative, panel data from Canterbury, and difference in difference evaluation method to explore the effects of displacement on Red Zone residential residents. We found that compared to non-relocated neighbours, the displaced people experience a significant initial decrease in their wages and salaries, and their total income. The impacts vary with time spent in the Red Zone and when they moved away. Wages and salaries of those who were red-zoned and moved in 2011 were reduced by 8%, and 5.4% for those who moved in 2012. Females faced greater decreases in wages and salaries, and total income, than males. There were no discernible impacts of the relocation on people’s self-employment income.
Earthquakes are insured only with public sector involvement in high-income countries where the risk of earthquakes is perceived to be high. The proto-typical examples of this public sector involvement are the public earthquake insurance schemes in California, Japan, and New Zealand (NZ). Each of these insurance programs is structured differently, and the purpose of this paper is to examine these differences using a concrete case-study, the sequence of earthquakes that occurred in the Christchurch, New Zealand, in 2011. This event turned out to have been the most heavily insured earthquake event in history. We examine what would have been the outcome of the earthquakes had the system of insurance in NZ been different. In particular, we focus on the public earthquake insurance programs in California (the California Earthquake Authority - CEA), and in Japan (Japanese Earthquake Reinsurance - JER). Overall, the aggregate cost to the public insurer in NZ was $NZ 11.1 billion in its response to the earthquakes. If a similar-sized disaster event had occurred in Japan and California, homeowners would have received $NZ 2.5 billion and $NZ 1.4 billion from the JER and CEA, respectively. We further describe the spatial and distributive patterns of these different scenarios.
The last seven years have seen southern New Zealand a ected by several large and damaging earthquakes: the moment magnitude (MW) 7.8 Dusky Sound earthquake on 15 July 2009, the MW 7.1 Dar eld (Canterbury) earthquake on 4 September 2010, and most notably the MW 6.2 Christchurch earthquake on 22 February 2011 and the protracted aftershock sequence. In this thesis, we address the postseismic displacement produced by these earthquakes using methods of satellite-based geodetic measurement, known as Interferometric Synthetic Aperture Radar (InSAR) and Global Positioning System (GPS), and computational modelling. We observe several ground displacement features in the Canterbury and Fiordland regions during three periods: 1) Following the Dusky Sound earthquake; 2) Following the Dar eld earthquake and prior to the Christchurch earthquake; and 3) Following the Christchurch earthquake until February 2015. The ground displacement associated with postseismic motion following the Dusky Sound earthquake has been measured by continuous and campaign GPS data acquired in August 2009, in conjunction with Di erential Interferometric Synthetic Aperture Radar (DInSAR) observations. We use an afterslip model, estimated by temporal inversion of geodetic data, with combined viscoelastic rebound model to account for the observed spatio-temporal patterns of displacement. The two postseismic processes together induce a signi cant displacement corresponding to principal extensional and contractual strain rates of the order of 10⁻⁷ and 10⁻⁸ yr⁻¹ respectively, across most of the southern South Island. We also analyse observed postseismic displacement following the Dusky Sound earthquake using a new inversion approach in order to describe afterslip in an elasticviscoelastic medium. We develop a mathematical framework, namely the "Iterative Decoupling of Afterslip and Viscoelastic rebound (IDAV)" method, with which to invert temporally dense and spatially sparse geodetic observations. We examine the IDAV method using both numerical and analytical simulations of Green's functions. For the post-Dar eld time interval, postseismic signals are measured within approximately one month of the mainshock. The dataset used for the post-Dar eld displacement spans the region surrounding previously unrecognised faults that ruptured during the mainshock. Poroelastic rebound in a multi-layered half-space and dilatancy recovery at shallow depths provide a satisfactory t with the observations. For the post-Christchurch interval, campaign GPS data acquired in February 2012 to February 2015 in four successive epochs and 66 TerraSAR-X (TSX) SAR acquisitions in descending orbits between March 2011 and May 2014 reveal approximately three years of postseismic displacement. We detect movement away from the satellite of ~ 3 mm/yr in Christchurch and a gradient of displacement of ~ 4 mm/yr across a lineament extending from the westernmost end of the Western Christchurch Fault towards the eastern end of the Greendale East Fault. The postseismic signals following the Christchurch earthquake are mainly accounted for by afterslip models on the subsurface lineament and nearby faults.
The Mѡ=7.1 Darfield (Canterbury) earthquake struck on 4 September 2010, approximately 45 km west of Christchurch, New Zealand. It revealed a previously unknown fault (the Greendale fault) and caused billions of dollars of damage due to high peak ground velocities and extensive liquefaction. It also triggered the Mw=6.3 Christchurch earthquake on 22 February 2011, which caused further damage and the loss of 185 lives. The objective of this research was to determine the relationship between stress and seismic properties in a seismically active region using manually-picked P and S wave arrival times from the aftershock sequence between 8 September 2010-13 January 2011 to estimate shear-wave splitting (SWS) parameters, VP =VS-ratios, anisotropy (delay-time tomography), focal mechanisms, and tectonic stress on the Canterbury plains. The maximum horizontal stress direction was highly consistent in the plains, with an average value of SHmax=116 18 . However, the estimates showed variation in SHmax near the fault, with one estimate rotating by as much as 30° counter-clockwise. This suggests heterogeneity of stress at the fault, though the cause remains unclear. Orientations of the principal stresses predominantly indicate a strike-slip regime, but there are possible thrust regimes to the west and north/east of the fault. The SWS fast directions (ø) on the plains show alignment with SHmax at the majority of stations, indicating stress controlled anisotropy. However, structural effects appear more dominant in the neighbouring regions of the Southern Alps and Banks Peninsula.
Sea level rise is one consequence of Earth’s changing climate. Century-long tide gauge records show that global-mean sea-level rise reached 11-16 cm during the twentieth century at a mean rate of 1.2 mm/y. Today, the average rate of global-mean sea-level rise is higher at 3-4 mm/y and is expected to increase in the future. This represents a hazard to low elevation coastal zones worldwide. Yet, before global sea level projections can be used to characterise future coastal flood hazard at a local scale, the effects of tectonics (and other processes) that drive vertical land motion (VLM) must be considered. VLM is defined as the vertical velocity (uplift or subsidence) of the solid surface with respect to the centre of Earth. In this study, new VLM maps are generated over coastal strips in New Zealand, using Sentinel-1 InSAR and GNSS data.In New Zealand, measuring VLM using InSAR on naturally vegetated or agricultural land is difficult due to signal decorrelation. Along the rural Bay of Plenty coastal strip, I use a persistent-scatterer approach to generate a VLM map from both east-looking ascending and west-looking descending Sentinel-1 data between 2015-2021. Using time-series data over the same time period from a dense network of 20 GNSS sensors, I tie InSAR-derived line-of-sight velocity to the 2014 ITRF reference frame. I test two different methods for measuring VLM and compare the results against GNSS vertical velocity along the Bay of Plenty coast. Best results are achieved by first removing the interpolated horizontal GNSS velocity field from each of the InSAR datasets, before averaging the two VLM estimates. Measured VLM is between -3 and 3 mm/y, with negative values (subsidence) occurring within the low-lying Rangitāiki Plain and Ōpōtiki valley, and uplift across the elevated region west of Matatā.This thesis integrates geomorphological, geological, and historical levelling VLM records with modern satellite datasets to assess VLM across timescales ranging from 10 to 100,000 years at Matatā. Uplift rate has been variable through time, with average uplift over the last 300,000 years of 1 mm/y, 4.5 mm/y since 1720 years, 2 mm/y between 1950-1978, and 10 mm/y between 2004-2011. Previous modelling has shown that the best fit to the 2004-2011 rapid uplift rates is an inflating magmatic source at ~10 km depth beneath Matatā. To reconcile all data, I present a VLM model that consists of short-lived periods (7 years) of rapid uplift (10 mm/y), separated by longer periods (30 years) of lower background uplift (3 mm/y). The episodic nature of VLM at Matatā likely reflects short-lived periods of magmatic intrusion. Episodic VLM characterised by large rates of uplift (10 mm/y) has been seen at Taupō volcano, and other volcanic centers globally. It has been 12 years since the end of the last intrusion episode; this modelling suggest one may expect to observe increased uplift rates at Matatā in the coming decades. Densely populated urban coastal strips are most at risk from the effects of relative sea-level rise. At the same time, anthropogenic activities associated with urbanization, such as groundwater withdrawal, and land reclamation can lead to local land subsidence (LLS), further exacerbating the risk to urban infrastructure. LLS refers to subsidence relative to nearby land area assumed to be stable. In this thesis, I create the first high-resolution (10 m) maps of LLS at six urban coastal strips in New Zealand, with a combined length of 285 km, using Sentinel-1 InSAR data between 2018-2021. This analysis reveals 89% of urban coastal strips are subsiding at rates of -0.5 mm/y or greater, and 11% is subsiding at higher rates of -3.0 mm/y or greater. On average, subsidence is -0.6 to -2.9 mm/y higher at the coastal strip, compared to inland areas occupied by GNSS stations. This analysis also documents highly-localised hotspots of LLS, with subsidence rates of up to -15 mm/y. In Christchurch, rapid and localised subsidence (-8 mm/y) is observed within coastal suburbs New Brighton and Southshore. In most cities, the highest subsidence rates occur on land reclaimed in the early-late twentieth century, and in areas built on Holocene sediment. Time-series analysis of LLS at sites of reclaimed land shows both linear and non-linear rates of deformation over time periods of up to 6-8 years. This thesis highlights the variable exposure to relative sea-level rise of New Zealand coastal strips, and demonstrates that in many cases current rates of VLM should be expected to continue for the next few decades.
