Beach ridge stratigraphy can provide an important record of both sustained coastal progradation and responses to events such as extreme storms, as well as evidence of earthquake induced sediment pulses. This study is a stratigraphic investigation of the late Holocene mixed sand gravel (MSG) beach ridge plain on the Canterbury coast, New Zealand. The subsurface was imaged along a 370 m shore-normal transect using 100 and 200 MHz ground penetrating radar (GPR) antennae, and cored to sample sediment textures. Results show that, seaward of a back-barrier lagoon, the Pegasus Bay beach ridge plain prograded almost uniformly, under conditions of relatively stable sea level. Nearshore sediment supply appears to have created a sustained sediment surplus, perhaps as a result of post-seismic sediment pulses, resulting in a flat, morphologically featureless beach ridge plain. Evidence of a high magnitude storm provides an exception, with an estimated event return period in excess of 100 years. Evidence from the GPR sequence combined with modern process observations from MSG beaches indicates that a paleo storm initially created a washover fan into the back-barrier lagoon, with a large amount of sediment simultaneously moved off the beach face into the nearshore. This erosion event resulted in a topographic depression still evident today. In the subsequent recovery period, sediment was reworked by swash onto the beach as a sequence of berm deposit laminations, creating an elevated beach ridge that also has a modern-day topographic signature. As sediment supply returned to normal, and under conditions of falling sea level, a beach ridge progradation sequence accumulated seaward of the storm feature out to the modern-day beach as a large flat, uniform progradation plain. This study highlights the importance of extreme storm events and earthquake pulses on MSG coastlines in triggering high volume beach ridge formation during the subsequent recovery period.
The majority of Christchurch’s stormwater has historically been discharged untreated directly into urban surface waterways. These receiving waterways have become adversely affected by the contaminants carried in the stormwater, particularly sediment and heavy metals. An event-based contaminant load model was developed to identify the distribution and magnitude of contaminant loads entering the waterway, as well as to assess the reduction in TSS and heavy metal loads that can be achieved by various stormwater management options. The GIS-Excel based model estimates contaminant loads from an individual storm event based on different contributing impervious surfaces and key rainfall characteristics (rainfall intensity, duration, pH and antecedent dry days). It then calculates contaminant reduction loads that could be achieved through source reduction (e.g. green roofs, repainting) as well as from treatment (e.g. raingardens, wet ponds) applied to different surfaces within the catchment. This model differs from other annual load models as it is event-based and accounts for storm characteristics in its calculation of contaminant loads. Christchurch is a valuable case setting due the unique opportunity for retrofitting improved stormwater management in the post-earthquake rebuild. It is anticipated that this modelling approach could later be adapted for use in other urban settings outside of Christchurch.
<strong>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.</strong>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.
Predictive modelling provides an efficient means to analyse the coastal environment and generate knowledge for long term urban planning. In this study, the numerical models SWAN and XBeach were incorporated into the ESRI ArcGIS interface by means of the BeachMMtool. This was applied to the Greater Christchurch coastal environment to simulate geomorphological evolution through hydrodynamic forcing. Simulations were performed using the recent sea level rise predictions by the Intergovernmental Panel on Climate Change (2013) to determine whether the statutory requirements outlined in the New Zealand Coastal Policy Statement 2010 are consistent with central, regional and district designations. Our results indicate that current land use zoning in Greater Christchurch is not consistent with these predictions. This is because coastal hazard risk has not been thoroughly quantified during the process of installing the Canterbury Earthquake Recovery Authority residential red zone. However, the Christchurch City Council’s flood management area does provide an extent to which managed coastal retreat is a real option. The results of this research suggest that progradation will continue to occur along the Christchurch foreshore due to the net sediment flux retaining an onshore direction and the current hydrodynamic activity not being strong enough to move sediment offshore. However, inundation during periods of storm surge poses a risk to human habitation on low lying areas around the Avon-Heathcote Estuary and the Brooklands lagoon.
The Canterbury earthquake sequence in New Zealand’s South Island induced widespread liquefaction phenomena across the Christchurch urban area on four occasions (4 Sept 2010; 22 Feb; 13 June; 23 Dec 2011), that resulted in widespread ejection of silt and fine sand. This impacted transport networks as well as infiltrated and contaminated the damaged storm water system, making rapid clean-up an immediate post-earthquake priority. In some places the ejecta was contaminated by raw sewage and was readily remobilised in dry windy conditions, creating a long-term health risk to the population. Thousands of residential properties were inundated with liquefaction ejecta, however residents typically lacked the capacity (time or resources) to clean-up without external assistance. The liquefaction silt clean-up response was co-ordinated by the Christchurch City Council and executed by a network of contractors and volunteer groups, including the ‘Farmy-Army’ and the ‘Student-Army’. The duration of clean-up time of residential properties and the road network was approximately 2 months for each of the 3 main liquefaction inducing earthquakes; despite each event producing different volumes of ejecta. Preliminary cost estimates indicate total clean-up costs will be over NZ$25 million. Over 500,000 tonnes of ejecta has been stockpiled at Burwood landfill since the beginning of the Canterbury earthquakes sequence. The liquefaction clean-up experience in Christchurch following the 2010-2011 earthquake sequence has emerged as a valuable case study to support further analysis and research on the coordination, management and costs of large volume deposition of fine grained sediment in urban areas.
