The Master of Engineering Management Project was sponsored by the Canterbury Earthquake Recovery Authority (CERA) and consisted of two phases: The first was an analysis of existing information detailing the effects of hazardous natural events on Canterbury Lifeline Utilities in the past 15 years. The aim of this “Lessons Learned” project was to produce an analysis report that identified key themes from the research, gaps in the existing data and to provide recommendations from these “Lessons Learned.” The Second phase was the development of a practical “Disaster Mitigation Guideline” that outlined lessons in the field of Emergency Sanitation. This research would build upon the first stage and would draw from international reference to develop a guideline that has practical implementation possibilities throughout the world.
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.
Landslides are significant hazards, especially in seismically-active mountainous regions, where shaking amplified by steep topography can result in widespread landsliding. These landslides present not only an acute hazard, but a chronic hazard that can last years-to-decades after the initial earthquake, causing recurring impacts. The Mw 7.8 Kaikōura earthquake caused more than 20,000 landslides throughout North Canterbury and resulted in significant damage to nationally significant infrastructure in the coastal transport corridor (CTC), isolating Kaikōura from the rest of New Zealand. In the years following, ongoing landsliding triggered by intense rainfall exacerbated the impacts and slowed the recovery process. However, while there is significant research on co-seismic landslides and their initial impacts in New Zealand, little research has explored the evolution of co-seismic landslides and how this hazard changes over time. This research maps landslides annually between 2013 and 2021 to evaluate the changes in pre-earthquake, co-seismic and post-earthquake rates of landsliding to determine how landslide hazard has changed over this time. In particular, the research explores how the number, area, and spatial distribution of landslides has changed since the earthquake, and whether post-earthquake mitigation works have in any way affected the long-term landslide hazard. Mapping of landslides was undertaken using open-source, medium resolution Landsat-8 and Sentinel-2 satellite imagery, with landslides identified visually and mapped as single polygons that capture both the source zone and deposit. Three study areas with differing levels of post-earthquake mitigation are compared: (i) the northern CTC, where the majority of mitigation was in the form of active debris removal; (ii) the southern CTC, where mitigation was primarily via passive protection measures; and (iii) Mount Fyffe, which has had no mitigation works since the earthquake. The results show that despite similar initial impacts during the earthquake, the rate of recovery in terms of landslide rates varies substantially across the three study areas. In Mount Fyffe, the number and area of landslides could take 45 and 22 years from 2021 respectively to return to pre-earthquake levels at the current rate. Comparatively, in the CTC, it could take just 5 years and 3-4 years from 2021 respectively. Notably, the fastest recovery in terms of landslide rates in the CTC was primarily located directly along the transport network, whereas what little recovery did occur in Mount Fyffe appeared to follow no particular pattern. Importantly, recovery rates in the northern CTC were notably higher than in the southern CTC, despite greater co-seismic impacts in the former. Combined, these results suggest the active, debris removal mitigation undertaken in the northern CTC may have had the effect of dramatically reducing the time for landslide rates to return to pre-earthquake levels. The role of slope angle and slope aspect were explored to evaluate if these observations could be driven by local differences in topography. The Mount Fyffe study area has higher slope angles than the CTC as a whole and landslides predominantly occurred on slightly steeper slopes than in the CTC. This may have contributed to the longer recovery times for landsliding in Mount Fyffe due to greater gravitational instability, however the observed variations are minor compared to the differences in recovery rates. In terms of slope aspect, landslides in Mount Fyffe preferentially occurred on north- and south-facing slopes whereas landslides in the CTC preferred the east- and south-facing slopes. The potential role of these differences in landslide recovery remains unclear but may be related to the propagation direction of the earthquake and the tracking direction of post-earthquake ex-tropical cyclones. Finally, landslides in the CTC are observed to be moving further away from the transport network and the number of landslides impacting the CTC decreased significantly since the earthquake. Nevertheless, the potential for further landslide reactivation remains. Therefore, despite the recovery in the CTC, it is clear that there is still risk of the transport network being impacted by further landsliding, at least for the next 3-5 yrs.
Farming and urban regions are impacted by earthquake disasters in different ways, and feature a range of often different recovery requirements. In New Zealand, and elsewhere, most earthquake impact and recovery research is urban focused. This creates a research deficit that can lead to the application of well-researched urban recovery strategies in rural areas to suboptimal effect. To begin to reduce this deficit, in-depth case studies of the earthquake impacts and recovery of three New Zealand farms severely impacted by the 14th November 2016, M7.8 Hurunui-Kaikōura earthquake were conducted. The initial earthquake, its aftershocks and coseismic hazards (e.g., landslides, liquefaction, surface rupture) affected much of North Canterbury, Marlborough and the Wellington area. The three case study farms were chosen to broadly represent the main types of farming and topography in the Hurunui District in North Canterbury. The farms were directly and indirectly impacted by earthquakes and related hazards. On-farm infrastructure (e.g., woolsheds, homesteads) and essential services (e.g., water, power), frequently sourced from distributed networks, were severely impacted. The earthquake occurred after two years of regional drought had already stressed farm systems and farmers to restructuring or breaking point. Cascading interlinked hazards stemming from the earthquakes and coseismic hazards continued to disrupt earthquake recovery over a year after the initial earthquake. Semi-structured interviews with the farmers were conducted nine and fourteen months after the initial earthquake to capture the timeline of on-going impacts and recovery. Analysis of both geological hazard data and interview data resulted in the identification of key factors influencing farm level earthquake impact and recovery. These include pre-existing conditions (e.g., drought); farm-specific variations in recovery timelines; and resilience strategies for farm recovery resources. The earthquake recovery process presented all three farms with opportunities to change their business plans and adapt to mitigate on-going and future risk.
This report provides an initial overview and gap analysis of the multi-hazards interactions that might affect fluvial and pluvial flooding (FPF) hazard in the Ōpāwaho Heathcote catchment. As per the terms of reference, this report focuses on a one-way analysis of the potential effects of multi-hazards on FPF hazard, as opposed to a more complex multi-way analysis of interactions between all hazards. We examined the relationship between FPF hazard and hazards associated with the phenomena of tsunamis; coastal erosion; coastal inundation; groundwater; earthquakes; and mass movements. Tsunamis: Modelling research indicates the worst-case tsunami scenarios potentially affecting the Ōpāwaho Heathcote catchment are far field. Under low probability, high impact tsunami scenarios waves could travel into Pegasus Bay and the Avon-Heathcote Estuary Ihutai, reaching the mouth and lower reaches of the Heathcote catchment and river, potentially inundating and eroding shorelines in sub-catchments 1 to 5, and temporarily blocking fluvial drainage more extensively. Any flooding infrastructure or management actions implemented in the area of tsunami inundation would ideally be resilient to tsunami-induced inundation and erosion. Model results currently available are a first estimate of potential tsunami inundation under contemporary sea and land level conditions. In terms of future large tsunami events, these models likely underestimate effects in riverside sub-catchments, as well as effects under future sea level, shoreline and other conditions. Also of significance when considering different FPF management structures, it is important to be mindful that certain types of flood structures can ‘trap’ inundating water coming from ocean directions, leading to longer flood durations and salinization issues. Coastal erosion: Model predictions indicate that sub-catchments 1 to 3 could potentially be affected by coastal erosion by the timescale of 2065, with sub-catchments 1-6 predicted to be potentially affected by coastal erosion by the time scale of 2115. In addition, the predicted open coast effects of this hazard should not be ignored since any significant changes in the New Brighton Spit open coast would affect erosion rates and exposure of the landward estuary margins, including the shorelines of the Ōpāwaho Heathcote catchment. Any FPF flooding infrastructure or management activities planned for the potentially affected sub-catchments needs to recognise the possibility of coastal erosion, and to have a planned response to the predicted potential shoreline translation. Coastal inundation: Model predictions indicate coastal inundation hazards could potentially affect sub-catchments 1 to 8 by 2065, with a greater area and depth of inundation possible for these same sub-catchments by 2115. Low-lying areas of the Ōpāwaho Heathcote catchment and river channel that discharge into the estuary are highly vulnerable to coastal inundation since elevated ocean and estuary water levels can block the drainage of inland systems, compounding FPF hazards. Coastal inundation can overwhelm stormwater and other drainage network components, and render river dredging options ineffective at best, flood enhancing at worst. A distinction can be made between coastal inundation and coastal erosion in terms of the potential impacts on affected land and assets, including flood infrastructure, and the implications for acceptance, adaptation, mitigation, and/or modification options. That is, responding to inundation could include structural and/or building elevation solutions, since unlike erosion, inundation does not necessarily mean the loss of land. Groundwater: Groundwater levels are of significant but variable concern when examining flooding hazards and management options in the Ōpāwaho Heathcote catchment due to variability in soils, topographies, elevations and proximities to riverine and estuarine surface waterbodies. Much of the Canterbury Plains part of the Ōpāwaho Heathcote catchment has a water table that is at a median depth of <1m from the surface (with actual depth below surface varying seasonally, inter-annually and during extreme meteorological events), though the water table depth rapidly shifts to >6m below the surface in the upper Plains part of the catchment (sub-catchments 13 to 15). Parts of Waltham/Linwood (sub-catchments 5 & 6) and Spreydon (sub-catchment 10) have extensive areas with a particularly high water table, as do sub-catchments 18, 19 and 20 south of the river. In all of the sub-catchments where groundwater depth below surface is shallow, it is necessary to be mindful of cascading effects on liquefaction hazard during earthquake events, including earthquake-induced drainage network and stormwater infrastructure damage. In turn, subsidence induced by liquefaction and other earthquake processes during the CES directly affected groundwater depth below surface across large parts of the central Ōpāwaho Heathcote catchment. The estuary margin of the catchment also faces increasing future challenges with sea level rise, which has the potential to elevate groundwater levels in these areas, compounding existing liquefaction and other earthquake associated multi-hazards. Any increases in subsurface runoff due to drainage system, development or climate changes are also of concern for the loess covered hill slopes due to the potential to enhance mass movement hazards. Earthquakes: Earthquake associated vertical ground displacement and liquefaction have historically affected, or are in future predicted to affect, all Ōpāwaho Heathcote sub-catchments. During the CES, these phenomena induced a significant cascades of changes in the city’s drainage systems, including: extensive vertical displacement and liquefaction induced damage to stormwater ‘greyware’, reducing functionality of the stormwater system; damage to the wastewater system which temporarily lowered groundwater levels and increased stormwater drainage via the wastewater network on the one hand, creating a pollution multi-hazard for FPF on the other hand; liquefaction and vertical displacement induced river channel changes affected drainage capacities; subsidence induced losses in soakage and infiltration capacities; changes occurred in topographic drainage conductivity; estuary subsidence (mainly around the Ōtākaro Avon rivermouth) increased both FPF and coastal inundation hazards; estuary bed uplift (severe around the Ōpāwaho Heathcote margins), reduced tidal prisms and increased bed friction, producing an overall reduction the waterbody’s capacity to efficiently flush catchment floodwaters to sea; and changes in estuarine and riverine ecosystems. All such possible effects need to be considered when evaluating present and future capacities of the Ōpāwaho Heathcote catchment FPF management systems. These phenomena are particularly of concern in the Ōpāwaho Heathcote catchment since stormwater networks must deal with constraints imposed by stream and river channels (past and present), estuarine shorelines and complex hill topography. Mass movements: Mass movements are primarily a risk in the Port Hills areas of the Ōpāwaho Heathcote catchment (sub-catchments 1, 2, 7, 9, 11, 16, 21), though there are one or two small but susceptible areas on the banks of the Ōpāwaho Heathcote River. Mass movements in the form of rockfalls and debris flows occurred on the Port Hills during the CES, resulting in building damage, fatalities and evacuations. Evidence has also been found of earthquake-triggered tunnel gully collapsesin all Port Hill Valleys. Follow-on effects of these mass movements are likely to occur in major future FPF and other hazard events. Of note, elevated groundwater levels, coastal inundation, earthquakes (including liquefaction and other effects), and mass movement exhibit the most extensive levels of multi-hazard interaction with FPF hazard. Further, all of the analysed multi-hazard interactions except earthquakes were found to consistently produce increases in the FPF hazard. The implications of these analyses are that multihazard interactions generally enhance the FPF hazard in the Ōpāwaho Heathcote catchment. Hence, management plans which exclude adjustments for multi-hazard interactions are likely to underestimate the FPF hazard in numerous different ways. In conclusion, although only a one-way analysis of the potential effects of selected multi-hazards on FPF hazard, this review highlights that the Ōpāwaho Heathcote catchment is an inherently multi- hazard prone environment. The implications of the interactions and process linkages revealed in this report are that several significant multi-hazard influences and process interactions must be taken into account in order to design a resilient FPF hazard management strategy.
With sea level rise (SLR) fast becoming one of the most pressing matters for governments worldwide, there has been mass amounts of research done on the impacts of SLR. However, these studies have largely focussed on the ways that SLR will impact both the natural and built environment, along with how the risk to low-lying coastal communities can be mitigated, while the inevitable impacts that this will have on mental well-being has been understudied. This research has attempted to determine the ways in which SLR can impact the mental well-being of those living in a low-lying coastal community, along with how these impacts could be mitigated while remaining adaptable to future environmental change. This was done through conducting an in-depth literature review to understand current SLR projections, the key components of mental well-being and how SLR can influence changes to mental well-being. This literature review then shaped a questionnaire which was distributed to residents of the New Brighton coastline. This questionnaire asked respondents how they interact with the local environment, how much they know about SLR and its associated hazards, whether SLR causes any level of stress or worry along with how respondents feel that these impacts could be mitigated. This research found that SLR impacts the mental well-being of those living in low-lying coastal communities through various methods: firstly, the respondents perceived risk to SLR and its associated hazards, which was found to be influenced by the suburbs that respondents live in, their knowledge of SLR, their main sources of information and the prior experience of the Canterbury Earthquake Sequence (CES). Secondly, the financial aspects of SLR were also found to be drivers of stress or worry, with depreciating property values and rising insurance premiums being frequently noted by respondents. It was found that the majority of respondents agreed that being involved in and informed of the protection process, having more readable and accurate information, and an increased engagement with community events and greenspaces would help to reduce the stress or worry caused by SLR, while remaining adaptable to future environmental change.
Lake Taupō in New Zealand is associated with frequent unrest and small to moderate eruptions. It presents a high consequence risk scenario with immense potential for destruction to the community and the surrounding environment. Unrest associated with eruptions may also trigger earthquakes. While it is challenging to educate people about the hazards and risks associated with multiple eruptive scenarios, effective education of students can lead to better mitigation strategies and risk reduction. Digital resources with user-directed outcomes have been successfully used to teach action oriented skills relevant for communication during volcanic crisis [4]. However, the use of choose your own adventure strategies to enhance low probability risk literacy for Secondary school outreach has not been fully explored. To investigate how digital narrative storytelling can mediate caldera risk literacy, a module “The Kid who cried Supervolcano” will be introduced in two secondary school classrooms in Christchurch and Rotorua. The module highlights four learning objectives: (a) Super-volcanoes are beautiful but can be dangerous (b) earthquake (unrest) activity is normal for super-volcanoes (c) Small eruptions are possible from super-volcanoes and can be dangerous in our lifetimes (d) Super-eruptions are unlikely in our lifetimes. Students will create their digital narrative using the platform Elementari (www.elementari.io). The findings from this study will provide clear understanding of students’ understanding of risk perceptions of volcanic eruption scenarios and associated hazards and inform the design of educational resources geared towards caldera risk literacy.
This article argues that active coordination of research engagement after disasters has the potential to maximize research opportunities, improve research quality, increase end-user engagement, and manage escalating research activity to mitigate ethical risks posed to impacted populations. The focus is on the coordination of research activity after the 22nd February 2011 Mw6.2 Christchurch earthquake by the then newly-formed national research consortium, the Natural Hazards Research Platform, which included a social science research moratorium during the declared state of national emergency. Decisions defining this organisation’s functional and structural parameters are analyzed to identify lessons concerning the need for systematic approaches to the management of post disaster research, in collaboration with the response effort. Other lessons include the importance of involving an existing, broadly-based research consortium, ensuring that this consortium's coordination role is fully integrated into emergency management structures, and ensuring that all aspects of decision-making processes are transparent and easily accessed.
Tsunami events including the 2004 Indian Ocean Tsunami and the 2011 Tohoku Earthquake and Tsunami confirmed the need for Pacific-wide comprehensive risk mitigation and effective tsunami evacuation planning. New Zealand is highly exposed to tsunamis and continues to invest in tsunami risk awareness, readiness and response across the emergency management and science sectors. Evacuation is a vital risk reduction strategy for preventing tsunami casualties. Understanding how people respond to warnings and natural cues is an important element to improving evacuation modelling techniques. The relative rarity of tsunami events locally in Canterbury and also globally, means there is limited knowledge on tsunami evacuation behaviour, and tsunami evacuation planning has been largely informed by hurricane evacuations. This research aims to address this gap by analysing evacuation behaviour and movements of Kaikōura and Southshore/New Brighton (coastal suburb of Christchurch) residents following the 2016 Kaikōura earthquake. Stage 1 of the research is engaging with both these communities and relevant hazard management agencies, using a survey and community workshops to understand real-event evacuation behaviour during the 2016 Kaikōura earthquake and subsequent tsunami evacuations. The second stage is using the findings from stage 1 to inform an agent-based tsunami evacuation model, which is an approach that simulates of the movement of people during an evacuation response. This method improves on other evacuation modelling approaches to estimate evacuation times due to better representation of local population characteristics. The information provided by the communities will inform rules and interactions such as traffic congestion, evacuation delay times and routes taken to develop realistic tsunami evacuation models. This will allow emergency managers to more effectively prepare communities for future tsunami events, and will highlight recommended actions to increase the safety and efficiency of future tsunami evacuations.
Recent tsunami events have highlighted the importance of effective tsunami risk management strategies (including land-use planning, structural and natural mitigation, warning systems, education and evacuation planning). However, the rarity of tsunami means that empirical data concerning reactions to tsunami warnings and evacuation behaviour is rare when compared to findings for evacuations from other hazards. More knowledge is required to document the full evacuation process, including responses to warnings, pre-evacuation actions, evacuation dynamics, and the return home. 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 the timings of evacuations, rather than behaviours of those evacuating. In this research, survey data was gathered from coastal communities in Banks Peninsula and Christchurch, New Zealand, relating to behaviours and actions during the November 14th 2016 Kaikōura earthquake tsunami. Survey questions asked about immediate actions following the earthquake shaking, reactions to tsunami warnings, pre-evacuation actions, evacuation dynamics and details on congestion. This data was analysed to characterise trends and identify factors that influenced evacuation actions and behaviour, and was further used to develop a realistic evacuation model prototype to evaluate the capacity of the roading network in Banks Peninsula during a tsunami evacuation. The evacuation model incorporated tsunami risk management strategies that have been implemented by local authorities, and exposure and vulnerability data, alongside the empirical data collected from the survey. This research enhances knowledge of tsunami evacuation behaviour and reactions to tsunami warnings, and can be used to refine evacuation planning to ensure that people can evacuate efficiently, thereby reducing their tsunami exposure and personal risk.
The Canterbury region of New Zealand experienced four earthquakes greater than MW 6.0 between September 2010 and December 2011. This study employs system dynamics as well as hazard, recovery and organisational literature and brings together data collected via surveys, case studies and interviews with organisations affected by the earthquakes. This is to show how systemic interactions and interdependencies within and between industry and geographic sectors affect their recovery post-disaster. The industry sectors in the study are: construction for its role in the rebuild, information and communication technology which is a regional high-growth industry, trucking for logistics, critical infrastructure, fast moving consumer goods (e.g. supermarkets) and hospitality to track recovery through non-discretionary and discretionary spend respectively. Also in the study are three urban centres including the region’s largest Central Business District, which has been inaccessible since the earthquake of 22 February 2011 to the time of writing in February 2013. This work also highlights how earthquake effects propagated between sectors and how sectors collaborated to mitigate difficulties such as product demand instability. Other interacting factors are identified that influence the recovery trajectories of the different industry sectors. These are resource availability, insurance payments, aid from central government, and timely and quality recovery information. This work demonstrates that in recovering from disaster it is crucial for organisations to identify what interacting factors could affect their operations. Also of importance are efforts to reduce the organisation’s vulnerability and increase their resilience to future crises and in day-to-day operations. Lastly, the multi-disciplinary approach to understanding the recovery and resilience of organisations and industry sectors after disaster, leads to a better understanding of effects as well as more effective recovery policy.
The research is funded by Callaghan Innovation (grant number MAIN1901/PROP-69059-FELLOW-MAIN) and the Ministry of Transport New Zealand in partnership with Mainfreight Limited. Need – The freight industry is facing challenges related to climate change, including natural hazards and carbon emissions. These challenges impact the efficiency of freight networks, increase costs, and negatively affect delivery times. To address these challenges, freight logistics modelling should consider multiple variables, such as natural hazards, sustainability, and emission reduction strategies. Freight operations are complex, involving various factors that contribute to randomness, such as the volume of freight being transported, the location of customers, and truck routes. Conventional methods have limitations in simulating a large number of variables. Hence, there is a need to develop a method that can incorporate multiple variables and support freight sustainable development. Method - A minimal viable model (MVM) method was proposed to elicit tacit information from industrial clients for building a minimally sufficient simulation model at the early modelling stages. The discrete-event simulation (DES) method was applied using Arena® software to create simulation models for the Auckland and Christchurch corridor, including regional pick-up and delivery (PUD) models, Christchurch city delivery models, and linehaul models. Stochastic variables in freight operations such as consignment attributes, customer locations, and truck routes were incorporated in the simulation. The geographic information system (GIS) software ArcGIS Pro® was used to identify and analyse industrial data. The results obtained from the GIS software were applied to create DES models. Life cycle assessment (LCA) models were developed for both diesel and battery electric (BE) trucks to compare their life cycle greenhouse gas (GHG) emissions and total cost of ownership (TCO) and support GHG emissions reduction. The line-haul model also included natural hazards in several scenarios, and the simulation was used to forecast the stock level of Auckland and Christchurch depots in response to each corresponding scenario. Results – DES is a powerful technique that can be employed to simulate and evaluate freight operations that exhibit high levels of variability, such as regional pickup and delivery (PUD) and linehaul. Through DES, it becomes possible to analyse multiple factors within freight operations, including transportation modes, routes, scheduling, and processing times, thereby offering valuable insights into the performance, efficiency, and reliability of the system. In addition, GIS is a useful tool for analysing and visualizing spatial data in freight operations. This is exemplified by their ability to simulate the travelling salesman problem (TSP) and conduct cluster analysis. Consequently, the integration of GIS into DES modelling is essential for improving the accuracy and reliability of freight operations analysis. The outcomes of the simulation were utilised to evaluate the ecological impact of freight transport by performing emission calculations and generating low-carbon scenarios to identify approaches for reducing the carbon footprint. LCA models were developed based on simulation results. Results showed that battery-electric trucks (BE) produced more greenhouse gas (GHG) emissions in the cradle phase due to battery manufacturing but substantially less GHG emissions in the use phase because of New Zealand's mostly renewable energy sources. While the transition to BE could significantly reduce emissions, the financial aspect is not compelling, as the total cost of ownership (TCO) for the BE truck was about the same for ten years, despite a higher capital investment for the BE. Moreover, external incentives are necessary to justify a shift to BE trucks. By using simulation methods, the effectiveness of response plans for natural hazards can be evaluated, and the system's vulnerabilities can be identified and mitigated to minimize the risk of disruption. Simulation models can also be utilized to simulate adaptation plans to enhance the system's resilience to natural disasters. Novel contributions – The study employed a combination of DES and GIS methods to incorporate a large number of stochastic variables and driver’s decisions into freight logistics modelling. Various realistic operational scenarios were simulated, including customer clustering and PUD truck allocation. This showed that complex pickup and delivery routes with high daily variability can be represented using a model of roads and intersections. Geographic regions of high customer density, along with high daily variability could be represented by a two-tier architecture. The method could also identify delivery runs for a whole city, which has potential usefulness in market expansion to new territories. In addition, a model was developed to address carbon emissions and total cost of ownership of battery electric trucks. This showed that the transition was not straightforward because the economics were not compelling, and that policy interventions – a variety were suggested - could be necessary to encourage the transition to decarbonised freight transport. A model was developed to represent the effect of natural disasters – such as earthquake and climate change – on road travel and detour times in the line haul freight context for New Zealand. From this it was possible to predict the effects on stock levels for a variety of disruption scenarios (ferry interruption, road detours). Results indicated that some centres rather than others may face higher pressure and longer-term disturbance after the disaster subsided. Remedies including coastal shipping were modelled and shown to have the potential to limit the adverse effects. A philosophical contribution was the development of a methodology to adapt the agile method into the modelling process. This has the potential to improve the clarification of client objectives and the validity of the resulting model.
Gravelly soils’ liquefaction has been documented since early 19th century with however the focus being sand-silts mixture – coarse documentation of this topic, that gravels do in fact liquefy was only acknowledged as an observation. With time, we have been impacted by earthquakes, EQ causing more damage to our property: life and environment-natural and built. In this realm of EQ related-damage the ground or soils in general act as buffer between the epicentre and the structures at a concerned site. Further, in this area, upon the eventual acknowledgement of liquefaction of soils as a problem, massive efforts were undertaken to understand its mechanics, what causes and thereby how to mitigate its ill-effects. Down that lane, gravelly soils’ liquefaction was another milestone covered in early 20th century, and thus regarded as a problematic subject. This being a fairly recent acknowledgement, efforts have initiated in this direction (or area of particle size under consideration-gravels>2mm), with this research outputs intended to complement that research for the betterment of our understanding of what’s happening and how shall we best address it, given the circumstances: socio (life) - environment (structures) - economic (cost or cost-“effectiveness’) and of course political (our “willingness” to want to address the problem). Case histories from at least 29 earthquakes worldwide have indicated that liquefaction can occur in gravelly soils (both in natural deposits and manmade reclamations) inducing large ground deformation and causing severe damage to civil infrastructures. However, the evaluation of the liquefaction resistance of gravelly soils remains to be a major challenge in geotechnical earthquake engineering. To date, laboratory tests aimed at evaluating the liquefaction resistance of gravelly soils are still very limited, as compared to the large body of investigations carried out on assessing the liquefaction resistance of sandy soils. While there is a general agreement that the liquefaction resistance of gravelly soils can be as low as that of clean sands, previous studies suggested that the liquefaction behaviour of gravelly soils is significantly affected by two key factors, namely relative density (Dr) and gravel content (Gc). While it is clear that the liquefaction resistance of gravels increases with the increasing Dr, there are inconclusive and/or contradictory results regarding the effect of Gc on the liquefaction resistance of gravelly soils. Aimed at addressing this important topic, an investigation is being currently carried out by researchers at the University of Canterbury, UC. As a first step, a series of undrained cyclic triaxial tests were conducted on selected sand-gravel mixtures (SGMs), and inter-grain state framework concepts such as the equivalent and skeleton void ratios were used to describe the joint effects of Gc and Dr on the liquefaction resistance of SGMs. Following such experimental effort, this study is aimed at providing new and useful insights, by developing a critical state-based method combined with the inter-grain state framework to uniquely describe the liquefaction resistance of gravelly soils. To do so, a series of monotonic drained triaxial tests will be carried out on selected SGMs. The outcomes of this study, combined with those obtained to date by UC researchers, will greatly contribute to the expansion of a worldwide assessment database, and also towards the development of a reliable liquefaction triggering procedure for characterising the liquefaction potential of gravelly soils, which is of paramount importance not only for the New Zealand context, but worldwide. This will make it possible for practising engineers to identify liquefiable gravelly soils in advance and make sound recommendations to minimise the impact of such hazards on land, and civil infrastructures.
Natural hazard disasters often have large area-wide impacts, which can cause adverse stress-related mental health outcomes in exposed populations. As a result, increased treatment-seeking may be observed, which puts a strain on the limited public health care resources particularly in the aftermath of a disaster. It is therefore important for public health care planners to know whom to target, but also where and when to initiate intervention programs that promote emotional wellbeing and prevent the development of mental disorders after catastrophic events. A large body of literature assesses factors that predict and mitigate disaster-related mental disorders at various time periods, but the spatial component has rarely been investigated in disaster mental health research. This thesis uses spatial and spatio-temporal analysis techniques to examine when and where higher and lower than expected mood and anxiety symptom treatments occurred in the severely affected Christchurch urban area (New Zealand) after the 2010/11 Canterbury earthquakes. High-risk groups are identified and a possible relationship between exposure to the earthquakes and their physical impacts and mood and anxiety symptom treatments is assessed. The main research aim is to test the hypothesis that more severely affected Christchurch residents were more likely to show mood and anxiety symptoms when seeking treatment than less affected ones, in essence, testing for a dose-response relationship. The data consisted of mood and anxiety symptom treatment information from the New Zealand Ministry of Health’s administrative databases and demographic information from the National Health Index (NHI) register, when combined built a unique and rich source for identifying publically funded stress-related treatments for mood and anxiety symptoms in almost the whole population of the study area. The Christchurch urban area within the Christchurch City Council (CCC) boundary was the area of interest in which spatial variations in these treatments were assessed. Spatial and spatio-temporal analyses were done by applying retrospective space-time and spatial variation in temporal trends analysis using SaTScan™ software, and Bayesian hierarchical modelling techniques for disease mapping using WinBUGS software. The thesis identified an overall earthquake-exposure effect on mood and anxiety symptom treatments among Christchurch residents in the context of the earthquakes as they experienced stronger increases in the risk of being treated especially shortly after the catastrophic 2011 Christchurch earthquake compared to the rest of New Zealand. High-risk groups included females, elderly, children and those with a pre-existing mental illness with elderly and children especially at-risk in the context of the earthquakes. Looking at the spatio-temporal distribution of mood and anxiety symptom treatments in the Christchurch urban area, a high rates cluster ranging from the severely affected central city to the southeast was found post-disaster. Analysing residential exposure to various earthquake impacts found that living in closer proximity to more affected areas was identified as a risk factor for mood and anxiety symptom treatments, which largely confirms a dose-response relationship between level of affectedness and mood and anxiety symptom treatments. However, little changes in the spatial distribution of mood and anxiety symptom treatments occurred in the Christchurch urban area over time indicating that these results may have been biased by pre-existing spatial disparities. Additionally, the post-disaster mobility activity from severely affected eastern to the generally less affected western and northern parts of the city seemed to have played an important role as the strongest increases in treatment rates occurred in less affected northern areas of the city, whereas the severely affected eastern areas tended to show the lowest increases. An investigation into the different effects of mobility confirmed that within-city movers and temporary relocatees were generally more likely to receive care or treatment for mood or anxiety symptoms, but moving within the city was identified as a protective factor over time. In contrast, moving out of the city from minor, moderately or severely damaged plain areas of the city, which are generally less affluent than Port Hills areas, was identified as a risk factor in the second year post-disaster. Moreover, residents from less damaged plain areas of the city showed a decrease in the likelihood of receiving care or treatment for mood or anxiety symptoms compared to those from undamaged plain areas over time, which also contradicts a possible dose-response relationship. Finally, the effects of the social and physical environment, as well as community resilience on mood and anxiety symptom treatments among long-term stayers from Christchurch communities indicate an exacerbation of pre-existing mood and anxiety symptom treatment disparities in the city, whereas exposure to ‘felt’ earthquake intensities did not show a statistically significant effect. The findings of this thesis highlight the complex relationship between different levels of exposure to a severe natural disaster and adverse mental health outcomes in a severely affected region. It is one of the few studies that have access to area-wide health and impact information, are able to do a pre-disaster / post-disaster comparison and track their sample population to apply spatial and spatio-temporal analysis techniques for exposure assessment. Thus, this thesis enhances knowledge about the spatio-temporal distribution of adverse mental health outcomes in the context of a severe natural disaster and informs public health care planners, not only about high-risk groups, but also where and when to target health interventions. The results indicate that such programs should broadly target residents living in more affected areas as they are likely to face daily hardship by living in a disrupted environment and may have already been the most vulnerable ones before the disaster. Special attention should be focussed on women, elderly, children and people with pre-existing mental illnesses as they are most likely to receive care or treatment for stress-related mental health symptoms. Moreover, permanent relocatees from affected areas and temporarily relocatees shortly after the disaster may need special attention as they face additional stressors due to the relocation that may lead to the development of adverse mental health outcomes needing treatment.