Recent global tsunami events have highlighted the importance of effective tsunami risk management strategies (including land-use planning, structural and natural defences, warning systems, education and evacuation measures). However, the rarity of tsunami means that empirical data concerning reactions to tsunami warnings and tsunami evacuation behaviour is rare when compared to findings about evacuations to avoid other sources of hazard. To date empirical research into tsunami evacuations has focused on evacuation rates, rather than other aspects of the evacuation process. More knowledge is required about responses to warnings, pre-evacuation actions, evacuation dynamics and the return home after evacuations. Tsunami evacuation modelling has the potential to inform evidence-based tsunami risk planning and response. However to date tsunami evacuation models have largely focused on timings of evacuations, rather than evacuation behaviours. This Masters research uses a New Zealand case study to reduce both of these knowledge gaps. Qualitative survey data was gathered from populations across coastal communities in Banks Peninsula and Christchurch, New Zealand, required to evacuate due to the tsunami generated by the November 14th 2016 Kaikōura Earthquake. Survey questions asked about reactions to tsunami warnings, actions taken prior to evacuating and movements during the 2016 tsunami evacuation. This data was analysed to characterise trends and identify factors that influenced evacuation actions and behaviour. Finally, it was used to develop an evacuation model for Banks Peninsula. Where appropriate, the modelling inputs were informed by the survey data. Three key findings were identified from the results of the evacuation behaviour survey. Although 38% of the total survey respondents identified the earthquake shaking as a natural cue for the tsunami, most relied on receiving official warnings, including sirens, to prompt evacuations. Respondents sought further official information to inform their evacuation decisions, with 39% of respondents delaying their evacuation in order to do so. Finally, 96% of total respondents evacuated by car. This led to congestion, particularly in more densely populated Christchurch city suburbs. Prior to this research, evacuation modelling had not been completed for Banks Peninsula. The results of the modelling showed that if evacuees know how to respond to tsunami warnings and where and how to evacuate, there are no issues. However, if there are poor conditions, including if people do not evacuate immediately, if there are issues with the roading network, or if people do not know where or how to evacuate, evacuation times increase with there being more bottlenecks leading out of the evacuation zones. The results of this thesis highlight the importance of effective tsunami education and evacuation planning. Reducing exposure to tsunami risk through prompt evacuation relies on knowledge of how to interpret tsunami warnings, and when, where and how to evacuate. Recommendations from this research outline the need for public education and engagement, and the incorporation of evacuation signage, information boards and evacuation drills. Overall these findings provide more comprehensive picture of tsunami evacuation behaviour and decision making based on empirical data from a recent evacuation, which can be used to improve tsunami risk management strategies. This empirical data can also be used to inform evacuation modelling to improve the accuracy and realism of the evacuation models.
Rapid, reliable information on earthquake-affected structures' current damage/health conditions and predicting what would happen to these structures under future seismic events play a vital role in accelerating post-event evaluations, leading to optimized on-time decisions. Such rapid and informative post-event evaluations are crucial for earthquake-prone areas, where each earthquake can potentially trigger a series of significant aftershocks, endangering the community's health and wealth by further damaging the already-affected structures. Such reliable post-earthquake evaluations can provide information to decide whether an affected structure is safe to stay in operation, thus saving many lives. Furthermore, they can lead to more optimal recovery plans, thus saving costs and time. The inherent deficiency of visual-based post-earthquake evaluations and the importance of structural health monitoring (SHM) methods and SHM instrumentation have been highlighted within this thesis, using two earthquake-affected structures in New Zealand: 1) the Canterbury Television (CTV) building, Christchurch; 2) the Bank of New Zealand (BNZ) building, Wellington. For the first time, this thesis verifies the theoretically- and experimentally validated hysteresis loop analysis (HLA) SHM method for the real-world instrumented structure of the BNZ building, which was damaged severely due to three earthquakes. Results indicate the HLA-SHM method can accurately estimate elastic stiffness degradation for this reinforced concrete (RC) pinched structure across the three earthquakes, which remained unseen until after the third seismic event. Furthermore, the HLA results help investigate the pinching effects on the BNZ building's seismic response. This thesis introduces a novel digital clone modelling method based on the robust and accurate SHM results delivered by the HLA method for physical parameters of the monitored structure and basis functions predicting the changes of these physical parameters due to future earthquake excitations. Contrary to artificial intelligence (AI) based predictive methods with black-box designs, the proposed predictive method is entirely mechanics-based with an explicitly-understandable design, making them more trusted and explicable to stakeholders engaging in post-earthquake evaluations, such as building owners and insurance firms. The proposed digital clone modelling framework is validated using the BNZ building and an experimental RC test structure damaged severely due to three successive shake-table excitations. In both structures, structural damage intensifies the pinching effects in hysteresis responses. Results show the basis functions identified from the HLA-SHM results for both structures under Event 1 can online estimate structural damage due to subsequent Events 2-3 from the measured structural responses, making them valuable tool for rapid warning systems. Moreover, the digital twins derived for these two structures under Event 1 can successfully predict structural responses and damage under Events 2-3, which can be integrated with the incremental dynamic analysis (IDA) method to assess structural collapse and its financial risks. Furthermore, it enables multi-step IDA to evaluate earthquake series' impacts on structures. Overall, this thesis develops an efficient method for providing reliable information on earthquake-affected structures' current and future status during or immediately after an earthquake, considerably guaranteeing safety. Significant validation is implemented against both experimental and real data of RC structures, which thus clearly indicate the accurate predictive performance of this HLA-based method.
This thesis documents the development and demonstration of an assessment method for analysing earthquake-related damage to concrete waste water gravity pipes in Christchurch, New Zealand, following the 2010-2011 Canterbury Earthquake Sequence (CES). The method is intended to be internationally adaptable to assist territorial local authorities with improving lifelines infrastructure disaster impact assessment and improvements in resilience. This is achieved through the provision of high-resolution, localised damage data, which demonstrate earthquake impacts along the pipe length. The insights gained will assist decision making and the prioritisation of resources following earthquake events to quickly and efficiently restore network function and reduce community impacts. The method involved obtaining a selection of 55 reinforced concrete gravity waste water pipes with available Closed-Circuit Television (CCTV) inspection footage filmed before and after the CES. The pipes were assessed by reviewing the recordings, and damage was mapped to the nearest metre along the pipe length using Geographic Information Systems. An established, systematic coding process was used for reporting the nature and severity of the observed damage, and to differentiate between pre-existing and new damage resulting from the CES. The damage items were overlaid with geospatial data such as Light Detection and Ranging (LiDAR)-derived ground deformation data, Liquefaction Resistance Index data and seismic ground motion data (Peak Ground acceleration and Peak Ground Velocity) to identify potential relationships between these parameters and pipe performance. Initial assessment outcomes for the pipe selection revealed that main pipe joints and lateral connections were more vulnerable than the pipe body during a seismic event. Smaller diameter pipes may also be more vulnerable than larger pipes during a seismic event. Obvious differential ground movement resulted in increased local damage observations in many cases, however this was not obvious for all pipes. Pipes with older installation ages exhibited more overall damage prior to a seismic event, which is likely attributable to increased chemical and biological deterioration. However, no evidence was found relating pipe age to performance during a seismic event. No evidence was found linking levels of pre-CES damage in a pipe with subsequent seismic performance, and seismic performance with liquefaction resistance or magnitude of seismic ground motion. The results reported are of limited application due to the small demonstration sample size, but reveal the additional level of detail and insight possible using the method presented in this thesis over existing assessment methods, especially in relation to high resolution variations along the length of the pipe such as localised ground deformations evidenced by LiDAR. The results may be improved by studying a larger and more diverse sample pool, automating data collection and input processes in order to improve efficiency and consider additional input such as pipe dip and cumulative damage over a large distance. The method is dependent on comprehensive and accurate pre-event CCTV assessments and LIDAR data so that post-event data could be compared. It is proposed that local territorial authorities should prioritise acquiring this information as a first important step towards improving the seismic resilience of a gravity waste water pipe network.
Liquefaction affects late Holocene, loose packed and water saturated sediment subjected to cyclical shear stress. Liquefaction features in the geological record are important off-fault markers that inform about the occurrence of moderate to large earthquakes (> 5 Mw). The study of contemporary liquefaction features provides a better understanding of where to find past (paleo) liquefaction features, which, if identified and dated, can provide information on the occurrence, magnitude and timing of past earthquakes. This is particularly important in areas with blind active faults. The extensive liquefaction caused by the 2010-2011 Canterbury Earthquake Sequence (CES) gave the geoscience community the opportunity to study the liquefaction process in different settings (alluvial, coastal and estuarine), investigating different aspects (e.g. geospatial correlation with landforms, thresholds for peak ground acceleration, resilience of infrastructures), and to collect a wealth geospatial dataset in the broad region of the Canterbury Plains. The research presented in this dissertation examines the sedimentary architecture of two environments, the alluvial and coastal settings, affected by liquefaction during the CES. The novel aim of this study is to investigate how landform and subsurface sedimentary architecture influence liquefaction and its surface manifestation, to provide knowledge for locating studies of paleoliquefaction in future. Two study cases documented in the alluvial setting showed that liquefaction features affected a crevasse splay and point bar ridges. However, the liquefaction source layer was linked to paleochannel floor deposits below the crevasse splay in the first case, and to the point bar deposits themselves in the second case. This research documents liquefaction features in the coastal dune system of the Canterbury Plains in detail for the first time. In the coastal dune setting the liquefiable layer is near the surface. The pore water pressure is vented easily because the coastal dune soil profile is entirely composed of non-cohesive, very well sorted sandy sediment that weakly resists disturbance from fluidised sediment under pressure. As a consequence, the liquefied flow does not need to find a specific crack through which the sediment is vented at the surface; instead, the liquefied sand finds many closely spaced conduits to vent its excess of pore water pressure. Therefore, in the coastal dune setting it is rare to observe discrete dikes (as they are defined in the alluvial setting), instead A horizon delamination (splitting) and blistering (near surface sills) are more common. The differences in styles of surface venting lead to contrasts in patterns of ejecta in the two environments. Whereas the alluvial environment is characterised by coalesced sand blows forming lineations, the coastal dune environment hosts apparently randomly distributed isolated sand blows often associated with collapse features. Amongst the techniques tested for the first time to investigate liquefaction features are: 3D GPR, which improved the accuracy of the trenching even six years after the liquefaction events; thin section analysis to investigate sediment fabric, which helped to discriminate liquefied sediment from its host sediment, and modern from paleoliquefaction features; a Random Forest classification based on the CES liquefaction map, which was used to test relationships between surface manifestation of liquefaction and topographic parameters. The results from this research will be used to target new study sites for future paleoliquefaction research and thus will improve the earthquake hazard assessment across New Zealand.
In the period between September 2010 and December 2011, Christchurch was shaken by a series of strong earthquakes including the MW7.1 4 September 2010, Mw 6.2 22 February 2011, MW6.2 13 June 2011 and MW6.0 23 December 2011 earthquakes. These earthquakes produced very strong ground motions throughout the city and surrounding areas that resulted in soil liquefaction and lateral spreading causing substantial damage to buildings, infrastructure and the community. The stopbank network along the Kaiapoi and Avon River suffered extensive damage with repairs projected to take several years to complete. This presented an opportunity to undertake a case-study on a regional scale of the effects of liquefaction on a stopbank system. Ultimately, this information can be used to determine simple performance-based concepts that can be applied in practice to improve the resilience of river protection works. The research presented in this thesis draws from data collected following the 4th September 2010 and 22nd February 2011 earthquakes. The stopbank damage is categorised into seven key deformation modes that were interpreted from aerial photographs, consultant reports, damage photographs and site visits. Each deformation mode provides an assessment of the observed mechanism of failure behind liquefaction-induced stopbank damage and the factors that influence a particular style of deformation. The deformation modes have been used to create a severity classification for the whole stopbank system, being ‘no or low damage’ and ‘major or severe damage’, in order to discriminate the indicators and factors that contribute to ‘major to severe damage’ from the factors that contribute to all levels of damage a number of calculated, land damage, stopbank damage and geomorphological parameters were analysed and compared at 178 locations along the Kaiapoi and Avon River stopbank systems. A critical liquefiable layer was present at every location with relatively consistent geotechnical parameters (cone resistance (qc), soil behaviour type (Ic) and Factor of Safety (FoS)) across the study site. In 95% of the cases the critical layer occurred within two times the Height of the Free Face (HFF,). A statistical analysis of the geotechnical factors relating to the critical layer was undertaken in order to find correlations between specific deformation modes and geotechnical factors. It was found that each individual deformation mode involves a complex interplay of factors that are difficult to represent through correlative analysis. There was, however, sufficient data to derive the key factors that have affected the severity of deformation. It was concluded that stopbank damage is directly related to the presence of liquefaction in the ground materials beneath the stopbanks, but is not critical in determining the type or severity of damage, instead it is merely the triggering mechanism. Once liquefaction is triggered it is the gravity-induced deformation that causes the damage rather than the shaking duration. Lateral spreading and specifically the depositional setting was found to be the key aspect in determining the severity and type of deformation along the stopbank system. The presence or absence of abandoned or old river channels and point bar deposits was found to significantly influence the severity and type of deformation. A review of digital elevation models and old maps along the Kaiapoi River found that all of the ‘major to severe’ damage observed occurred within or directly adjacent to an abandoned river channel. Whilst a review of the geomorphology along the Avon River showed that every location within a point bar deposit suffered some form of damage, due to the depositional environment creating a deposit highly susceptible to liquefaction.
The development of Digital City technologies to manage and visualise spatial information has increasingly become a focus of the research community, and application by city authorities. Traditionally, the Geographic Information Systems (GIS) and Building Information Models (BIM) underlying Digital Cities have been used independently. However, integrating GIS and BIM into a single platform provides benefits for project and asset management, and is applicable to a range of issues. One of these benefits is the means to access and analyse large datasets describing the built environment, in order to characterise urban risk from and resilience to natural hazards. The aim of this thesis is to further explore methodologies of integration in two distinct areas. The first, integration through connectivity of heterogeneous datasets where GIS spatial infrastructure data is merged with 3D BIM building data to create a digital twin. Secondly, integration through analysis whereby data from the digital twin are extracted and integrated with computational models. To achieve this, a workflow was developed to identify the required datasets of a digital twin, and develop a process of integrating those datasets through a combination of; semi-autonomous conversion, translation and extension of data; and semantic web and services-based processes. Through use of a designed schema, the data were streamed in a homogenous format in a web-based platform. To demonstrate the value of this workflow with respect to urban risk and resilience, the process was applied to the Taiora: Queen Elizabeth II recreation and sports centre in eastern Christchurch, New Zealand. After integration of as-built GIS and BIM datasets, targeted data extraction was implemented, with outputs tailored for analysis in an infrastructure serviceability loss model, which assessed potable water network performance in the 22nd February 2011 Christchurch Earthquake. Using the same earthquake conditions as the serviceability loss model, performance of infrastructure assets in service at the time of the 22nd February 2011 Christchurch Earthquake was compared to new assets rebuilt at the site, post-earthquake. Due to improved potable water infrastructure resilience resulting from installation of ductile piles, a decrease of 35.5% in the probability of service loss was estimated in the serviceability loss model. To complete the workflow, the results from the external analysis were uploaded to the web-based platform. One of the more significant outcomes from the workflow was the identification of a lack of mandated metadata standards for fittings/valves connecting a building to private laterals. Whilst visually the GIS and BIM data show the building and pipes as connected, the semantic data does not include this connectivity relationship. This has no material impact on the current serviceability loss model as it is not one of the defined parameters. However, a proposed modification to the model would utilise the metadata to further assess the physical connection robustness, and increase the number of variables for estimating probability of service loss. This thesis has made a methodological contribution to urban resilience analysis by demonstrating how readily available up-to-date building and infrastructure data can be integrated, and with tailored extraction from a Digital City platform, be used for disaster impact analysis in an external computational engine, with results in turn imported and visualised in the Digital City platform. The workflow demonstrated that translation and integration of data would be more successful if a regional/national mandate was implemented for the submission of consent documentation in a specified standard BIM format. The results of this thesis have identified that the key to ensuring the success of an integrated tool lies in the initial workflow required to safeguard that all data can be either captured or translated in an interoperable format.
The development of Digital City technologies to manage and visualise spatial information has increasingly become a focus of the research community, and application by city authorities. Traditionally, the Geographic Information Systems (GIS) and Building Information Models (BIM) underlying Digital Cities have been used independently. However, integrating GIS and BIM into a single platform provides benefits for project and asset management, and is applicable to a range of issues. One of these benefits is the means to access and analyse large datasets describing the built environment, in order to characterise urban risk from and resilience to natural hazards. The aim of this thesis is to further explore methodologies of integration in two distinct areas. The first, integration through connectivity of heterogeneous datasets where GIS spatial infrastructure data is merged with 3D BIM building data to create a digital twin. Secondly, integration through analysis whereby data from the digital twin are extracted and integrated with computational models. To achieve this, a workflow was developed to identify the required datasets of a digital twin, and develop a process of integrating those datasets through a combination of; semi-autonomous conversion, translation and extension of data; and semantic web and services-based processes. Through use of a designed schema, the data were streamed in a homogenous format in a web-based platform. To demonstrate the value of this workflow with respect to urban risk and resilience, the process was applied to the Taiora: Queen Elizabeth II recreation and sports centre in eastern Christchurch, New Zealand. After integration of as-built GIS and BIM datasets, targeted data extraction was implemented, with outputs tailored for analysis in an infrastructure serviceability loss model, which assessed potable water network performance in the 22nd February 2011 Christchurch Earthquake. Using the same earthquake conditions as the serviceability loss model, performance of infrastructure assets in service at the time of the 22nd February 2011 Christchurch Earthquake was compared to new assets rebuilt at the site, post-earthquake. Due to improved potable water infrastructure resilience resulting from installation of ductile piles, a decrease of 35.5% in the probability of service loss was estimated in the serviceability loss model. To complete the workflow, the results from the external analysis were uploaded to the web-based platform. One of the more significant outcomes from the workflow was the identification of a lack of mandated metadata standards for fittings/valves connecting a building to private laterals. Whilst visually the GIS and BIM data show the building and pipes as connected, the semantic data does not include this connectivity relationship. This has no material impact on the current serviceability loss model as it is not one of the defined parameters. However, a proposed modification to the model would utilise the metadata to further assess the physical connection robustness, and increase the number of variables for estimating probability of service loss. This thesis has made a methodological contribution to urban resilience analysis by demonstrating how readily available up-to-date building and infrastructure data can be integrated, and with tailored extraction from a Digital City platform, be used for disaster impact analysis in an external computational engine, with results in turn imported and visualised in the Digital City platform. The workflow demonstrated that translation and integration of data would be more successful if a regional/national mandate was implemented for the submission of consent documentation in a specified standard BIM format. The results of this thesis have identified that the key to ensuring the success of an integrated tool lies in the initial workflow required to safeguard that all data can be either captured or translated in an interoperable format.
The Canterbury Earthquake Sequence (CES) of 2010-2011 produced large seismic moments up to Mw 7.1. These large, near-to-surface (<15 km) ruptures triggered >6,000 rockfall boulders on the Port Hills of Christchurch, many of which impacted houses and affected the livelihoods of people within the impacted area. From these disastrous and unpredicted natural events a need arose to be able to assess the areas affected by rockfall events in the future, where it is known that a rockfall is possible from a specific source outcrop but the potential boulder runout and dynamics are not understood. The distribution of rockfall deposits is largely constrained by the physical properties and processes of the boulder and its motion such as block density, shape and size, block velocity, bounce height, impact and rebound angle, as well as the properties of the substrate. Numerical rockfall models go some way to accounting for all the complex factors in an algorithm, commonly parameterised in a user interface where site-specific effects can be calibrated. Calibration of these algorithms requires thorough field checks and often experimental practises. The purpose of this project, which began immediately following the most destructive rupture of the CES (February 22, 2011), is to collate data to characterise boulder falls, and to use this information, supplemented by a set of anthropogenic boulder fall data, to perform an in-depth calibration of the three-dimensional numerical rockfall model RAMMS::Rockfall. The thesis covers the following topics: • Use of field data to calibrate RAMMS. Boulder impact trails in the loess-colluvium soils at Rapaki Bay have been used to estimate ranges of boulder velocities and bounce heights. RAMMS results replicate field data closely; it is concluded that the model is appropriate for analysing the earthquake-triggered boulder trails at Rapaki Bay, and that it can be usefully applied to rockfall trajectory and hazard assessment at this and similar sites elsewhere. • Detailed analysis of dynamic rockfall processes, interpreted from recorded boulder rolling experiments, and compared to RAMMS simulated results at the same site. Recorded rotational and translational velocities of a particular boulder show that the boulder behaves logically and dynamically on impact with different substrate types. Simulations show that seasonal changes in soil moisture alter rockfall dynamics and runout predictions within RAMMS, and adjustments are made to the calibration to reflect this; suggesting that in hazard analysis a rockfall model should be calibrated to dry rather than wet soil conditions to anticipate the most serious outcome. • Verifying the model calibration for a separate site on the Port Hills. The results of the RAMMS simulations show the effectiveness of calibration against a real data set, as well as the effectiveness of vegetation as a rockfall barrier/retardant. The results of simulations are compared using hazard maps, where the maximum runouts match well the mapped CES fallen boulder maximum runouts. The results of the simulations in terms of frequency distribution of deposit locations on the slope are also compared with those of the CES data, using the shadow angle tool to apportion slope zones. These results also replicate real field data well. Results show that a maximum runout envelope can be mapped, as well as frequency distribution of deposited boulders for hazard (and thus risk) analysis purposes. The accuracy of the rockfall runout envelope and frequency distribution can be improved by comprehensive vegetation and substrate mapping. The topics above define the scope of the project, limiting the focus to rockfall processes on the Port Hills, and implications for model calibration for the wider scientific community. The results provide a useful rockfall analysis methodology with a defensible and replicable calibration process, that has the potential to be applied to other lithologies and substrates. Its applications include a method of analysis for the selection and positioning of rockfall countermeasure design; site safety assessment for scaling and demolition works; and risk analysis and land planning for future construction in Christchurch.
Recent surface-rupturing earthquakes in New Zealand have highlighted significant exposure and vulnerability of the road network to fault displacement. Understanding fault displacement hazard and its impact on roads is crucial for mitigating risks and enhancing resilience. There is a need for regional-scale assessments of fault displacement to identify vulnerable areas within the road network for the purposes of planning and prioritising site-specific investigations. This thesis employs updated analysis of data from three historical surface-rupturing earthquakes (Edgecumbe 1987, Darfield 2010, and Kaikoūra 2016) to develop an empirical model that addresses the gap in regional fault displacement hazard analysis. The findings contribute to understanding of • How to use seismic hazard model inputs for regional fault displacement hazard analysis • How faulting type and sediment cover affects the magnitude and spatial distribution of fault displacement • How the distribution of displacement and regional fault displacement hazard is impacted by secondary faulting • The inherent uncertainties and limitations associated with employing an empirical approach at a regional scale • Which sections of New Zealand’s roading network are most susceptible to fault displacement hazard and warrant site-specific investigations • Which regions should prioritise updating emergency management plans to account for post-event disruptions to roading. I used displacement data from the aforementioned historical ruptures to generate displacement versus distance-to-fault curves for various displacement components, fault types, and geological characteristics. Using those relationships and established relationships for along-strike displacement, displacement contours were generated surrounding active faults within the NZ Community Fault Model. Next, I calculated a new measure of 1D strain along roads as well as relative hazard, which integrated 1D strain and normalised slip rate data. Summing these values at the regional level identified areas of heightened relative hazard across New Zealand, and permits an assessment of the susceptibility of road networks using geomorphon land classes as proxies for vulnerability. The results reveal that fault-parallel displacements tend to localise near the fault plane, while vertical and fault-perpendicular displacements sustain over extended distances. Notably, no significant disparities were observed in off-fault displacement between the hanging wall and footwall sides of the fault, or among different surface geology types, potentially attributed to dataset biases. The presence of secondary faulting in the dataset contributes to increased levels of tectonic displacement farther from the fault, highlighting its significance in hazard assessments. Furthermore, fault displacement contours delineate broader zones around dip-slip faults compared to strike-slip faults, with correlations identified between fault length and displacement width. Road ‘strain’ values are higher around dip-slip faults, with notable examples observed in the Westland and Buller Districts. As expected, relative hazard analysis revealed elevated values along faults with high slip rates, notably along the Alpine Fault. A regional-scale analysis of hazard and exposure reveals heightened relative hazard in specific regions, including Wellington, Southern Hawke’s Bay, Central Bay of Plenty, Central West Coast, inland Canterbury, and the Wairau Valley of Marlborough. Notably, the Central West Coast exhibits the highest summed relative hazard value, attributed to the fast-slipping Alpine Fault. The South Island generally experiences greater relative hazard due to larger and faster-slipping faults compared to the North Island, despite having fewer roads. Central regions of New Zealand face heightened risk compared to Southern or Northern regions. Critical road links intersecting high-slipping faults, such as State Highways 6, 73, 1, and 2, necessitate prioritisation for site-specific assessments, emergency management planning and targeted mitigation strategies. Roads intersecting with the Alpine Fault are prone to large parallel displacements, requiring post-quake repair efforts. Mitigation strategies include future road avoidance of nearby faults, modification of road fill and surface material, and acknowledgement of inherent risk, leading to prioritised repair efforts of critical roads post-quake. Implementing these strategies enhances emergency response efforts by improving accessibility to isolated regions following a major surface-rupturing event, facilitating faster supply delivery and evacuation assistance. This thesis contributes to the advancement of understanding fault displacement hazard by introducing a novel regional, empirical approach. The methods and findings highlight the importance of further developing such analyses and extending them to other critical infrastructure types exposed to fault displacement hazard in New Zealand. Enhancing our comprehension of the risks associated with fault displacement hazard offers valuable insights into various mitigation strategies for roading infrastructure and informs emergency response planning, thereby enhancing both national and global infrastructure resilience against geological hazards.
Coastal and river environments are exposed to a number of natural hazards that have the potential to negatively affect both human and natural environments. The purpose of this research is to explain that significant vulnerabilities to seismic hazards exist within coastal and river environments and that coasts and rivers, past and present, have played as significant a role as seismic, engineering or socio-economic factors in determining the impacts and recovery patterns of a city following a seismic hazard event. An interdisciplinary approach was used to investigate the vulnerability of coastal and river areas in the city of Christchurch, New Zealand, following the Canterbury Earthquake Sequence, which began on the 4th of September 2010. This information was used to identify the characteristics of coasts and rivers that make them more susceptible to earthquake induced hazards including liquefaction, lateral spreading, flooding, landslides and rock falls. The findings of this research are applicable to similar coastal and river environments elsewhere in the world where seismic hazards are also of significant concern. An interdisciplinary approach was used to document and analyse the coastal and river related effects of the Canterbury earthquake sequence on Christchurch city in order to derive transferable lessons that can be used to design less vulnerable urban communities and help to predict seismic vulnerabilities in other New Zealand and international urban coastal and river environments for the future. Methods used to document past and present features and earthquake impacts on coasts and rivers in Christchurch included using maps derived from Geographical Information Systems (GIS), photographs, analysis of interviews from coastal, river and engineering experts, and analysis of secondary data on seismicity, liquefaction potential, geology, and planning statutes. The Canterbury earthquake sequence had a significant effect on Christchurch, particularly around rivers and the coast. This was due to the susceptibility of rivers to lateral spreading and the susceptibility of the eastern Christchurch and estuarine environments to liquefaction. The collapse of river banks and the extensive cracking, tilting and subsidence that accompanied liquefaction, lateral spreading and rock falls caused damage to homes, roads, bridges and lifelines. This consequently blocked transportation routes, interrupted electricity and water lines, and damaged structures built in their path. This study found that there are a number of physical features of coastal and river environments from the past and the present that have induced vulnerabilities to earthquake hazards. The types of sediments found beneath eastern Christchurch are unconsolidated fine sands, silts, peats and gravels. Together with the high water tables located beneath the city, these deposits made the area particularly susceptible to liquefaction and liquefaction-induced lateral spreading, when an earthquake of sufficient size shook the ground. It was both past and present coastal and river processes that deposited the types of sediments that are easily liquefied during an earthquake. Eastern Christchurch was once a coastal and marine environment 6000 years ago when the shoreline reached about 6 km inland of its present day location, which deposited fine sand and silts over this area. The region was also exposed to large braided rivers and smaller spring fed rivers, both of which have laid down further fine sediments over the following thousands of years. A significant finding of this study is the recognition that the Canterbury earthquake sequence has exacerbated existing coastal and river hazards and that assessments and monitoring of these changes will be an important component of Christchurch’s future resilience to natural hazards. In addition, patterns of recovery following the Canterbury earthquakes are highlighted to show that coasts and rivers are again vulnerable to earthquakes through their ability to recovery. This city’s capacity to incorporate resilience into the recovery efforts is also highlighted in this study. Coastal and river areas have underlying physical characteristics that make them increasingly vulnerable to the effects of earthquake hazards, which have not typically been perceived as a ‘coastal’ or ‘river’ hazard. These findings enhance scientific and management understanding of the effects that earthquakes can have on coastal and river environments, an area of research that has had modest consideration to date. This understanding is important from a coastal and river hazard management perspective as concerns for increased human development around coastlines and river margins, with a high seismic risk, continue to grow.
One of the most controversial issues highlighted by the 2010-2011 Christchurch earthquake series and more recently the 2016 Kaikoura earthquake, has been the evident difficulty and lack of knowledge and guidelines for: a) evaluation of the residual capacity damaged buildings to sustain future aftershocks; b) selection and implementation of a series of reliable repairing techniques to bring back the structure to a condition substantially the same as prior to the earthquake; and c) predicting the cost (or cost-effectiveness) of such repair intervention, when compared to fully replacement costs while accounting for potential aftershocks in the near future. As a result of such complexity and uncertainty (i.e., risk), in combination with the possibility (unique in New Zealand when compared to most of the seismic-prone countries) to rely on financial support from the insurance companies, many modern buildings, in a number exceeding typical expectations from past experiences at an international level, have ended up being demolished. This has resulted in additional time and indirect losses prior to the full reconstruction, as well as in an increase in uncertainty on the actual relocation of the investment. This research project provides the main end-users and stakeholders (practitioner engineers, owners, local and government authorities, insurers, and regulatory agencies) with comprehensive evidence-based information to assess the residual capacity of damage reinforced concrete buildings, and to evaluate the feasibility of repairing techniques, in order to support their delicate decision-making process of repair vs. demolition or replacement. Literature review on effectiveness of epoxy injection repairs, as well as experimental tests on full-scale beam-column joints shows that repaired specimens have a reduced initial stiffness compared with the undamaged specimen, with no apparent strength reduction, sometimes exhibiting higher displacement ductility capacities. Although the bond between the steel and concrete is only partially restored, it still allows the repaired specimen to dissipate at least the same amount of hysteretic energy. Experimental tests on buildings subjected to earthquake loading demonstrate that even for severe damage levels, the ability of the epoxy injection to restore the initial stiffness of the structure is significant. Literature review on damage assessment and repair guidelines suggests that there is consensus within the international community that concrete elements with cracks less than 0.2 mm wide only require cosmetic repairs; epoxy injection repairs of cracks less and 2.0 mm wide and concrete patching of spalled cover concrete (i.e., minor to moderate damage) is an appropiate repair strategy; and for severe damaged components (e.g., cracks greater than 2.0 mm wide, crushing of the concrete core, buckling of the longitudinal reinforcement) local replacement of steel and/or concrete in addition to epoxy crack injection is more appropriate. In terms of expected cracking patterns, non-linear finite element investigations on well-designed reinforced concrete beam-to-column joints, have shown that lower number of cracks but with wider openings are expected to occur for larger compressive concrete strength, f’c, and lower reinforcement content, ρs. It was also observed that the tensile concrete strength, ft, strongly affects the expected cracking pattern in the beam-column joints, the latter being more uniformly distributed for lower ft values. Strain rate effects do not seem to play an important role on the cracking pattern. However, small variations in the cracking pattern were observed for low reinforcement content as it approaches to the minimum required as per NZS 3101:2006. Simple equations are proposed in this research project to relate the maximum and residual crack widths with the steel strain at peak displacement, with or without axial load. A literature review on fracture of reinforcing steel due to low-cycle fatigue, including recent research using steel manufactured per New Zealand standards is also presented. Experimental results describing the influence of the cyclic effect on the ultimate strain capacity of the steel are also discussed, and preliminary equations to account for that effect are proposed. A literature review on the current practice to assess the seismic residual capacity of structures is also presented. The various factors affecting the residual fatigue life at a component level (i.e., plastic hinge) of well-designed reinforced concrete frames are discussed, and equations to quantify each of them are proposed, as well as a methodology to incorporate them into a full displacement-based procedure for pre-earthquake and post-earthquake seismic assessment.
One of the great challenges facing human systems today is how to prepare for, manage, and adapt successfully to the profound and rapid changes wreaked by disasters. Wellington, New Zealand, is a capital city at significant risk of devastating earthquake and tsunami, potentially requiring mass evacuations with little or short notice. Subsequent hardship and suffering due to widespread property damage and infrastructure failure could cause large areas of the Wellington Region to become uninhabitable for weeks to months. Previous research has shown that positive health and well-being are associated with disaster-resilient outcomes. Preventing adverse outcomes before disaster strikes, through developing strengths-based skill sets in health-protective attitudes and behaviours, is increasingly advocated in disaster research, practise, and management. This study hypothesised that well-being constructs involving an affective heuristic play vital roles in pathways to resilience as proximal determinants of health-protective behaviours. Specifically, this study examined the importance of health-related quality of life and subjective well-being in motivating evacuation preparedness, measured in a community sample (n=695) drawn from the general adult population of Wellington’s isolated eastern suburbs. Using a quantitative epidemiological approach, the study measured the prevalence of key quality of life indicators (physical and mental health, emotional well-being or “Sense of Coherence”, spiritual well-being, social well-being, and life satisfaction) using validated psychometric scales; analysed the strengths of association between these indicators and the level of evacuation preparedness at categorical and continuous levels of measurement; and tested the predictive power of the model to explain the variance in evacuation preparedness activity. This is the first study known to examine multi-dimensional positive health and global well-being as resilient processes for engaging in evacuation preparedness behaviour. A cross-sectional study design and quantitative survey were used to collect self-report data on the study variables; a postal questionnaire was fielded between November 2008 and March 2009 to a sampling frame developed through multi-stage cluster randomisation. The survey response rate was 28.5%, yielding a margin of error of +/- 3.8% with 95% confidence and 80% statistical power to detect a true correlation coefficient of 0.11 or greater. In addition to the primary study variables, data were collected on demographic and ancillary variables relating to contextual factors in the physical environment (risk perception of physical and personal vulnerability to disaster) and the social environment (through the construct of self-determination), and other measures of disaster preparedness. These data are reserved for future analyses. Results of correlational and regression analyses for the primary study variables show that Wellingtonians are highly individualistic in how their well-being influences their preparedness, and a majority are taking inadequate action to build their resilience to future disaster from earthquake- or tsunami-triggered evacuation. At a population level, the conceptual multi-dimensional model of health-related quality of life and global well-being tested in this study shows a positive association with evacuation preparedness at statistically significant levels. However, it must be emphasised that the strength of this relationship is weak, accounting for only 5-7% of the variability in evacuation preparedness. No single dimension of health-related quality of life or well-being stands out as a strong predictor of preparedness. The strongest associations for preparedness are in a positive direction for spiritual well-being, emotional well-being, and life satisfaction; all involve a sense of existential meaningfulness. Spiritual well-being is the only quality of life variable making a statistically significant unique contribution to explaining the variance observed in the regression models. Physical health status is weakly associated with preparedness in a negative direction at a continuous level of measurement. No association was found at statistically significant levels for mental health status and social well-being. These findings indicate that engaging in evacuation preparedness is a very complex, holistic, yet individualised decision-making process, and likely involves highly subjective considerations for what is personally relevant. Gender is not a factor. Those 18-24 years of age are least likely to prepare and evacuation preparedness increases with age. Multidimensional health and global well-being are important constructs to consider in disaster resilience for both pre-event and post-event timeframes. This work indicates a need for promoting self-management of risk and building resilience by incorporating a sense of personal meaning and importance into preparedness actions, and for future research into further understanding preparedness motivations.
The Avon-Heathcote Estuary is of significant value to Christchurch due to its high productivity, biotic diversity, proximity to the city, and its cultural, recreational and aesthetic qualities. Nonetheless, it has been subjected to decades of degradation from sewage wastewater discharges and encroaching urban development. The result was a eutrophied estuary, high in nitrogen, affected by large blooms of nuisance macroalgae and covered by degraded sediments. In March 2010, treated wastewater was diverted from the estuary to a site 3 km offshore. This quickly reduced water nitrogen by 90% within the estuary and, within months, there was reduced production of macroalgae. However, a series of earthquakes beginning in September 2010 brought massive changes: tilting of the estuary, changes in channels and water flow, and a huge influx of liquefied sediments that covered up to 65% of the estuary floor. Water nitrogen increased due to damage to sewage infrastructure and the diversion pipeline being turned off. Together, these drastically altered the estuarine ecosystem. My study involves three laboratory and five in situ experiments that investigate the base of the food chain and responses of benthic microalgae to earthquake-driven sediment and nutrient changes. It was predicted that the new sediments would be coarser and less contaminated with organic matter and nutrients than the old sediments, would have decreased microalgal biomass, and would prevent invertebrate grazing and bioturbation activities. It was believed that microalgal biomass would become similar across new and old sediments types as the unstable new sediments were resuspended and distributed over the old sediments. Contact cores of the sediment were taken at three sites, across a eutrophication gradient, monthly from September 2011 to March 2012. Extracted chlorophyll a pigments showed that microalgal biomass was generally lower on new liquefied sediments compared to old sediments, although there was considerable site to site variation, with the highly eutrophic sites being the most affected by the emergence of the new sediments. Grazer experiments showed that invertebrates had both positive and negative site-specific effects on microalgal biomass depending on their identity. At one site, new sediments facilitated grazing by Amphibola crenata, whereas at another site, new sediments did not alter the direct and indirect effects of invertebrates (Nicon aestuariensis, Macropthalmus hirtipes, and A. crenata) on microalgae. From nutrient addition experiments it was clear that benthic microalgae were able to use nutrients from within both old and new sediments equally. This implied that microalgae were reducing legacy nutrients in both sediments, and that they are an important buffer against eutrophication. Therefore, in tandem with the wastewater diversion, they could underpin much of the recovery of the estuary. Overall, the new sediments were less favourable for benthic microalgal growth and recolonisation, but were less contaminated than old sediments at highly eutrophic sites. Because the new sediments were less contaminated than the old sediments, they could help return the estuary to a noneutrophic state. However, if the new sediments, which are less favourable for microalgal growth, disperse over the old sediments at highly eutrophic sites, they could become contaminated and interfere with estuarine recovery. Therefore, recovery of microalgal communities and the estuary was expected to be generally long, but variable and site-specific, with the least eutrophic sites recovering quickly, and the most eutrophic sites taking years to return to a pre-earthquake and non-eutrophied state. changes in channels and water flow, and a huge influx of liquefied sediments that covered up to 65% of the estuary floor. Water nitrogen increased due to damage to sewage infrastructure and the diversion pipeline being turned off. Together, these drastically altered the estuarine ecosystem. My study involves three laboratory and five in situ experiments that investigate the base of the food chain and responses of benthic microalgae to earthquake-driven sedimen tand nutrient changes. It was predicted that the new sediments would be coarser and less contaminated with organic matter and nutrients than the old sediments, would have decreased microalgal biomass, and would prevent invertebrate grazing and bioturbation activities. It was believed that microalgal biomass would become similar across new and old sediments types as the unstable new sediments were resuspended and distributed over the old sediments. Contact cores of the sediment were taken at three sites, across a eutrophication gradient, monthly from September 2011 to March 2012. Extracted chlorophyll a pigments showed that microalgal biomass was generally lower on new liquefied sediments compared to old sediments, although there was considerable site to site variation, with the highly eutrophic sites being the most affected by the emergence of the new sediments. Grazer experiments showed that invertebrates had both positive and negative site-specific effects on microalgal biomass depending on their identity. At one site, new sediments facilitated grazing by Amphibola crenata, whereas at another site, new sediments did not alter the direct and indirect effects of invertebrates (Nicon aestuariensis, Macropthalmus hirtipes, and A. crenata) on microalgae. From nutrient addition experiments it was clear that benthic microalgae were able to use nutrients from within both old and new sediments equally. This implied that microalgae were reducing legacy nutrients in both sediments, and that they are
In the last century, seismic design has undergone significant advancements. Starting from the initial concept of designing structures to perform elastically during an earthquake, the modern seismic design philosophy allows structures to respond to ground excitations in an inelastic manner, thereby allowing damage in earthquakes that are significantly less intense than the largest possible ground motion at the site of the structure. Current performance-based multi-objective seismic design methods aim to ensure life-safety in large and rare earthquakes, and to limit structural damage in frequent and moderate earthquakes. As a result, not many recently built buildings have collapsed and very few people have been killed in 21st century buildings even in large earthquakes. Nevertheless, the financial losses to the community arising from damage and downtime in these earthquakes have been unacceptably high (for example; reported to be in excess of 40 billion dollars in the recent Canterbury earthquakes). In the aftermath of the huge financial losses incurred in recent earthquakes, public has unabashedly shown their dissatisfaction over the seismic performance of the built infrastructure. As the current capacity design based seismic design approach relies on inelastic response (i.e. ductility) in pre-identified plastic hinges, it encourages structures to damage (and inadvertently to incur loss in the form of repair and downtime). It has now been widely accepted that while designing ductile structural systems according to the modern seismic design concept can largely ensure life-safety during earthquakes, this also causes buildings to undergo substantial damage (and significant financial loss) in moderate earthquakes. In a quest to match the seismic design objectives with public expectations, researchers are exploring how financial loss can be brought into the decision making process of seismic design. This has facilitated conceptual development of loss optimisation seismic design (LOSD), which involves estimating likely financial losses in design level earthquakes and comparing against acceptable levels of loss to make design decisions (Dhakal 2010a). Adoption of loss based approach in seismic design standards will be a big paradigm shift in earthquake engineering, but it is still a long term dream as the quantification of the interrelationships between earthquake intensity, engineering demand parameters, damage measures, and different forms of losses for different types of buildings (and more importantly the simplification of the interrelationship into design friendly forms) will require a long time. Dissecting the cost of modern buildings suggests that the structural components constitute only a minor portion of the total building cost (Taghavi and Miranda 2003). Moreover, recent research on seismic loss assessment has shown that the damage to non-structural elements and building contents contribute dominantly to the total building loss (Bradley et. al. 2009). In an earthquake, buildings can incur losses of three different forms (damage, downtime, and death/injury commonly referred as 3Ds); but all three forms of seismic loss can be expressed in terms of dollars. It is also obvious that the latter two loss forms (i.e. downtime and death/injury) are related to the extent of damage; which, in a building, will not just be constrained to the load bearing (i.e. structural) elements. As observed in recent earthquakes, even the secondary building components (such as ceilings, partitions, facades, windows parapets, chimneys, canopies) and contents can undergo substantial damage, which can lead to all three forms of loss (Dhakal 2010b). Hence, if financial losses are to be minimised during earthquakes, not only the structural systems, but also the non-structural elements (such as partitions, ceilings, glazing, windows etc.) should be designed for earthquake resistance, and valuable contents should be protected against damage during earthquakes. Several innovative building technologies have been (and are being) developed to reduce building damage during earthquakes (Buchanan et. al. 2011). Most of these developments are aimed at reducing damage to the buildings’ structural systems without due attention to their effects on non-structural systems and building contents. For example, the PRESSS system or Damage Avoidance Design concept aims to enable a building’s structural system to meet the required displacement demand by rocking without the structural elements having to deform inelastically; thereby avoiding damage to these elements. However, as this concept does not necessarily reduce the interstory drift or floor acceleration demands, the damage to non-structural elements and contents can still be high. Similarly, the concept of externally bracing/damping building frames reduces the drift demand (and consequently reduces the structural damage and drift sensitive non-structural damage). Nevertheless, the acceleration sensitive non-structural elements and contents will still be very vulnerable to damage as the floor accelerations are not reduced (arguably increased). Therefore, these concepts may not be able to substantially reduce the total financial losses in all types of buildings. Among the emerging building technologies, base isolation looks very promising as it seems to reduce both inter-storey drifts and floor accelerations, thereby reducing the damage to the structural/non-structural components of a building and its contents. Undoubtedly, a base isolated building will incur substantially reduced loss of all three forms (dollars, downtime, death/injury), even during severe earthquakes. However, base isolating a building or applying any other beneficial technology may incur additional initial costs. In order to provide incentives for builders/owners to adopt these loss-minimising technologies, real-estate and insurance industries will have to acknowledge the reduced risk posed by (and enhanced resilience of) such buildings in setting their rental/sale prices and insurance premiums.
Mitigating the cascade of environmental damage caused by the movement of excess reactive nitrogen (N) from land to sea is currently limited by difficulties in precisely and accurately measuring N fluxes due to variable rates of attenuation (denitrification) during transport. This thesis develops the use of the natural abundance isotopic composition of nitrate (δ15N and δ18O of NO₃-) to integrate the spatialtemporal variability inherent to denitrification, creating an empirical framework for evaluating attenuation during land to water NO₃- transfers. This technique is based on the knowledge that denitrifiers kinetically discriminate against 'heavy' forms of both N and oxygen (O), creating a parallel enrichment in isotopes of both species as the reaction progresses. This discrimination can be quantitatively related to NO₃- attenuation by isotopic enrichment factors (εdenit). However, while these principles are understood, use of NO₃- isotopes to quantify denitrification fluxes in non-marine environments has been limited by, 1) poor understanding of εdenit variability, and, 2) difficulty in distinguishing the extent of mixing of isotopically distinct sources from the imprint of denitrification. Through a combination of critical literature analysis, mathematical modelling, mesocosm to field scale experiments, and empirical studies on two river systems over distance and time, these short comings are parametrised and a template for future NO₃- isotope based attenuation measurements outlined. Published εdenit values (n = 169) are collated in the literature analysis presented in Chapter 2. By evaluating these values in the context of known controllers on the denitrification process, it is found that the magnitude of εdenit, for both δ15N and δ18O, is controlled by, 1) biology, 2) mode of transport through the denitrifying zone (diffusion v. advection), and, 3) nitrification (spatial-temporal distance between nitrification and denitrification). Based on the outcomes of this synthesis, the impact of the three factors identified as controlling εdenit are quantified in the context of freshwater systems by combining simple mathematical modelling and lab incubation studies (comparison of natural variation in biological versus physical expression). Biologically-defined εdenit, measured in sediments collected from four sites along a temperate stream and from three tropical submerged paddy fields, varied from -3‰ to -28‰ depending on the site’s antecedent carbon content. Following diffusive transport to aerobic surface water, εdenit was found to become more homogeneous, but also lower, with the strength of the effect controlled primarily by diffusive distance and the rate of denitrification in the sediments. I conclude that, given the variability in fractionation dynamics at all levels, applying a range of εdenit from -2‰ to -10‰ provides more accurate measurements of attenuation than attempting to establish a site-specific value. Applying this understanding of denitrification's fractionation dynamics, four field studies were conducted to measure denitrification/ NO₃- attenuation across diverse terrestrial → freshwater systems. The development of NO₃- isotopic signatures (i.e., the impact of nitrification, biological N fixation, and ammonia volatilisation on the isotopic 'imprint' of denitrification) were evaluated within two key agricultural regions: New Zealand grazed pastures (Chapter 4) and Philippine lowland submerged rice production (Chapter 5). By measuring the isotopic composition of soil ammonium, NO₃- and volatilised ammonia following the bovine urine deposition, it was determined that the isotopic composition of NO₃ - leached from grazed pastures is defined by the balance between nitrification and denitrification, not ammonia volatilisation. Consequently, NO₃- created within pasture systems was predicted to range from +10‰ (δ15N)and -0.9‰ (δ18O) for non-fertilised fields (N limited) to -3‰ (δ15N) and +2‰ (δ18O) for grazed fertilised fields (N saturated). Denitrification was also the dominant determinant of NO₃- signatures in the Philippine rice paddy. Using a site-specific εdenit for the paddy, N inputs versus attenuation were able to be calculated, revealing that >50% of available N in the top 10 cm of soil was denitrified during land preparation, and >80% of available N by two weeks post-transplanting. Intriguingly, this denitrification was driven by rapid NO₃- production via nitrification of newly mineralised N during land preparation activities. Building on the relevant range of εdenit established in Chapters 2 and 3, as well as the soil-zone confirmation that denitrification was the primary determinant of NO₃- isotopic composition, two long-term longitudinal river studies were conducted to assess attenuation during transport. In Chapter 6, impact and recovery dynamics in an urban stream were assessed over six months along a longitudinal impact gradient using measurements of NO₃- dual isotopes, biological populations, and stream chemistry. Within 10 days of the catastrophic Christchurch earthquake, dissolved oxygen in the lowest reaches was <1 mg l⁻¹, in-stream denitrification accelerated (attenuating 40-80% of sewage N), microbial biofilm communities changed, and several benthic invertebrate taxa disappeared. To test the strength of this method for tackling the diffuse, chronic N loading of streams in agricultural regions, two years of longitudinal measurements of NO₃- isotopes were collected. Attenuation was negatively correlated with NO₃- concentration, and was highly dependent on rainfall: 93% of calculated attenuation (20 kg NO₃--N ha⁻¹ y⁻¹) occurred within 48 h of rainfall. The results of these studies demonstrate the power of intense measurements of NO₃- stable isotope for distinguishing temporal and spatial trends in NO₃ - loss pathways, and potentially allow for improved catchment-scale management of agricultural intensification. Overall this work now provides a more cohesive understanding for expanding the use of NO₃- isotopes measurements to generate accurate understandings of the controls on N losses. This information is becoming increasingly important to predict ecosystem response to future changes, such the increasing agricultural intensity needed to meet global food demand, which is occurring synergistically with unpredictable global climate change.