During many years the analysis of some geophysical results of Charles Darwin was being carried out in Department. Darwin has connected almost 200 years ago results of catastrophic earthquakes with vertical movement of a surface of the Earth. Usually this movement less horizontal movement and its influence on destruction of cities is not considered. Earthquake hazard assessment studies were focused usually on the horizontal ground motion. Effects of the strong vertical motion were not, practically, discussed. The margins of safety against gravity-induced static vertical forces in constructed buildings usually provide adequate resistance to dynamic forces induced by the vertical acceleration during an earthquake. However, the earthquake in Christchurch is an example of the vertical seismic shock . The earthquake magnitude was rather small - nearby 6.3. However, the result was catastrophic. The same took place in 1835. It allowed to Darwin to formulate a few great ideas. Charles Darwin has explained qualitatively results of an interaction of huge seismic waves with volcanoes and the nature of volcanism and seismicity of our planet. These important data of Charles Darwin became very actual recently. It is possible to tell also the same about tsunami and extreme ocean waves described by Charles Darwin. Therefore this data were analyzed using modern mechanics, mathematics and physics in Department. In particular, the theory of catastrophic waves was developed based on Darwin's data. The theory tried to explain occurrence, evolution and distribution the catastrophic waves in various natural systems, since atoms, oceans, surfaces of the Earth and up to the very early Universe. Some results of the research were published in prestigious magazines. Later they were presented in two books devoted to Charles Darwin's anniversary (2009). Last from them was published in Russian (2011). We give here key ideas of this research which is a part of interdisciplinary researches of Department. Some ideas are discussed. Not less important purpose is very short historical review of some researches of Darwin. In particular, we underline Darwin' priority in the formulation of the bases of Dynamics Earth.
Probabilistic Structural Fire Engineering (PSFE) has been introduced to overcome the limitations of current conventional approaches used for the design of fire-exposed structures. Current structural fire design investigates worst-case fire scenarios and include multiple thermal and structural analyses. PSFE permits buildings to be designed to a level of life safety or economic loss that may occur in future fire events with the help of a probabilistic approach. This thesis presents modifications to the adoption of a Performance-Based Earthquake Engineering (PBEE) framework in Probabilistic Structural Fire Engineering (PSFE). The probabilistic approach runs through a series of interrelationships between different variables, and successive convolution integrals of these interrelationships result in probabilities of different measures. The process starts with the definition of a fire severity measure (FSM), which best relates fire hazard intensity with structural response. It is identified by satisfying efficiency and sufficiency criteria as described by the PBEE framework. The relationship between a fire hazard and corresponding structural response is established by analysis methods. One method that has been used to quantify this relationship in PSFE is Incremental Fire Analysis (IFA). The existing IFA approach produces unrealistic fire scenarios, as fire profiles may be scaled to wide ranges of fire severity levels, which may not physically represent any real fires. Two new techniques are introduced in this thesis to limit extensive scaling. In order to obtain an annual rate of exceedance of fire hazard and structural response for an office building, an occurrence model and an attenuation model for office fires are generated for both Christchurch city and New Zealand. The results show that Christchurch city is 15% less likely to experience fires that have the potential to cause structural failures in comparison to all of New Zealand. In establishing better predictive relationships between fires and structural response, cumulative incident radiation (a fire hazard property) is found to be the most appropriate fire severity measure. This research brings together existing research on various sources of uncertainty in probabilistic structural fire engineering, such as elements affecting post-flashover fire development factors (fuel load, ventilation, surface lining and compartment geometry), fire models, analysis methods and structural reliability. Epistemic uncertainty and aleatory uncertainty are investigated in the thesis by examining the uncertainty associated with modelling and the factors that influence post-flashover development of fires. A survey of 12 buildings in Christchurch in combination with recent surveys in New Zealand produced new statistical data on post-flashover development factors in office buildings in New Zealand. The effects of these parameters on temperature-time profiles are evaluated. The effects of epistemic uncertainty due to fire models in the estimation of structural response is also calculated. Parametric fires are found to have large uncertainty in the prediction of post-flashover fires, while the BFD curves have large uncertainties in prediction of structural response. These uncertainties need to be incorporated into failure probability calculations. Uncertainty in structural modelling shows that the choices that are made during modelling have a large influence on realistic predictions of structural response.
On 14 November 2016, the Mw 7.8 Kaikōura earthquake caused widespread damage along the east coast of the South Island, New Zealand. Kaikōura town itself was isolated from the rest of the country by landslides blocking off major roads. While impacts from the Kaikōura earthquake on large, urban population centres have been generally well documented, this thesis aims to fill gaps in academic knowledge regarding small rural towns. This thesis investigates what, where and when critical infrastructure and lifeline service disruption occurred following the 2016 Kaikōura earthquake in a selection of small towns, and how the communities in these areas adapted to disruption. Following a robust review of literature and news media, four small rural towns were selected from North Canterbury (Culverden & Waiau) and Marlborough (Seddon & Ward) in the South Island, New Zealand. Semi-structured interview sessions with a special focus on these towns were held with infrastructure managers, emergency response and recovery officials, and organisation leaders with experience or expertise in the 2016 Kaikōura earthquake. Findings were supplemented with emergency management situation reports to produce hazard maps and infrastructure exposure maps. A more detailed analysis was conducted for Waiau involving interdependence analyses and a level of service timeline for select lifeline services. The earthquake impacted roads by blocking them with landslides, debris and surface rupture. Bridges where shaken off their abutments, breaking infrastructure links such as fibre landlines as they went. Water supplies and other forms of infrastructure relied heavily on the level of service of roads, as rough rural terrain left few alternatives. Adapting to an artificial loss of road service, some Waiau locals created their own detour around a road cordon in order to get home to family and farms. Performance of dwellings was tied to socioeconomic factors as much as proximity to the epicentre. Farmers who lost water access pulled out fences to allow stock to drink from rivers. Socioeconomic differences between farmland and township residents also contributed to resilience variations between the towns assessed in this study. Understanding how small rural towns respond and adapt to disaster allows emergency management officials and policy to be well informed and flexible with planning for multiple size classes of towns.
Depending on their nature and severity, disasters can create large volumes of debris and waste. Waste volumes from a single event can be the equivalent of many times the annual waste generation rate of the affected community. These volumes can overwhelm existing solid waste management facilities and personnel. Mismanagement of disaster waste can affect both the response and long term recovery of a disaster affected area. Previous research into disaster waste management has been either context specific or event specific, making it difficult to transfer lessons from one disaster event to another. The aim of this research is to develop a systems understanding of disaster waste management and in turn develop context- and disaster-transferrable decision-making guidance for emergency and waste managers. To research this complex and multi-disciplinary problem, a multi-hazard, multi-context, multi-case study approach was adopted. The research focussed on five major disaster events: 2011 Christchurch earthquake, 2009 Victorian Bushfires, 2009 Samoan tsunami, 2009 L’Aquila earthquake and 2005 Hurricane Katrina. The first stage of the analysis involved the development of a set of ‘disaster & disaster waste’ impact indicators. The indicators demonstrate a method by which disaster managers, planners and researchers can simplify the very large spectra of possible disaster impacts, into some key decision-drivers which will likely influence post-disaster management requirements. The second stage of the research was to develop a set of criteria to represent the desirable environmental, economic, social and recovery effects of a successful disaster waste management system. These criteria were used to assess the effectiveness of the disaster waste management approaches for the case studies. The third stage of the research was the cross-case analysis. Six main elements of disaster waste management systems were identified and analysed. These were: strategic management, funding mechanisms, operational management, environmental and human health risk management, and legislation and regulation. Within each of these system elements, key decision-making guidance (linked to the ‘disaster & disaster waste’ indicators) and management principles were developed. The ‘disaster & disaster waste’ impact indicators, the effects assessment criteria and management principles have all been developed so that they can be practically applied to disaster waste management planning and response in the future.
High-quality ground motion records are required for engineering applications including response history analysis, seismic hazard development, and validation of physics-based ground motion simulations. However, the determination of whether a ground motion record is high-quality is poorly handled by automation with mathematical functions and can become prohibitive if done manually. Machine learning applications are well-suited to this problem, and a previous feed-forward neural network was developed (Bellagamba et al. 2019) to determine high-quality records from small crustal events in the Canterbury and Wellington regions for simulation validation. This prior work was however limited by the omission of moderate-to-large magnitude events and those from other tectonic environments, as well as a lack of explicit determination of the minimum usable frequency of the ground motion. To address these shortcomings, an updated neural network was developed to predict the quality of ground motion records for all magnitudes and all tectonic sources—active shallow crustal, subduction intraslab, and subduction interface—in New Zealand. The predictive performance of the previous feed-forward neural network was matched by the neural network in the domain of small crustal records, and this level of predictive performance is now extended to all source magnitudes and types in New Zealand making the neural network applicable to global ground motion databases. Furthermore, the neural network provides quality and minimum usable frequency predictions for each of the three orthogonal components of a record which may then be mapped into a binary quality decision or otherwise applied as desired. This framework provides flexibility for the end user to predict high-quality records with various acceptability thresholds allowing for this neural network to be used in a range of applications.
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.
Geologic phenomena produced by earthquake shaking, including rockfalls and liquefaction features, provide important information on the intensity and spatiotemporal distribution of earthquake ground motions. The study of rockfall and liquefaction features produced in contemporary well- instrumented earthquakes increases our knowledge of how natural and anthropogenic environments respond to earthquakes and improves our ability to deduce seismologic information from analogous pre-contemporary (paleo-) geologic features. The study of contemporary and paleo- rockfall and liquefaction features enables improved forecasting of environmental responses to future earthquakes. In this thesis I utilize a combination of field and imagery-based mapping, trenching, stratigraphy, and numerical dating techniques to understand the nature and timing of rockfalls (and hillslope sedimentation) and liquefaction in the eastern South Island of New Zealand, and to examine the influence that anthropogenic activity has had on the geologic expressions of earthquake phenomena. At Rapaki (Banks Peninsula, NZ), field and imagery-based mapping, statistical analysis and numerical modeling was conducted on rockfall boulders triggered by the fatal 2011 Christchurch earthquakes (n=285) and compared with newly identified prehistoric (Holocene and Pleistocene) boulders (n=1049) deposited on the same hillslope. A significant population of modern boulders (n=26) travelled farther downslope (>150 m) than their most-travelled prehistoric counterparts, causing extensive damage to residential dwellings at the foot of the hillslope. Replication of prehistoric boulder distributions using 3-dimensional rigid body numerical models requires the application of a drag-coefficient, attributed to moderate to dense slope vegetation, to account for their spatial distribution. Radiocarbon dating provides evidence for 17th to early 20th century deforestation at the study site during Polynesian and European colonization and after emplacement of prehistoric rockfalls. Anthropocene deforestation enabled modern rockfalls to exceed the limits of their prehistoric predecessors, highlighting a shift in the geologic expression of rockfalls due to anthropogenic activity. Optical and radiocarbon dating of loessic hillslope sediments in New Zealand’s South Island is used to constrain the timing of prehistoric rockfalls and associated seismic events, and quantify spatial and temporal patterns of hillslope sedimentation including responses to seismic and anthropogenic forcing. Luminescence ages from loessic sediments constrain timing of boulder emplacement to between ~3.0 and ~12.5 ka, well before the arrival of Polynesians (ca AD 1280) and Europeans (ca AD 1800) in New Zealand, and suggest loess accumulation was continuing at the study site until 12-13 ka. Large (>5 m3) prehistoric rockfall boulders preserve an important record of Holocene hillslope sedimentation by creating local traps for sediment aggradation and upbuilding soil formation. Sediment accumulation rates increased considerably (>~10 factor increase) following human arrival and associated anthropogenic burning of hillslope vegetation. New numerical ages are presented to place the evolution of loess-mantled hillslopes in New Zealand’s South Island into a longer temporal framework and highlight the roles of earthquakes and humans on hillslope surface process. Extensive field mapping and characterization for 1733 individual prehistoric rockfall boulders was conducted at Rapaki and another Banks Peninsula site, Purau, to understand their origin, frequency, and spatial and volumetric distributions. Boulder characteristics and distributions were compared to 421 boulders deposited at the same sites during the 2010-2011 Canterbury earthquake sequence. Prehistoric boulders at Rapaki and Purau are comprised of two dominant lithofacies types: volcanic breccia and massive (coherent) lava basalt. Volcanic breccia boulders are found in greatest abundance (64-73% of total mapped rockfall) and volume (~90-96% of total rockfall) at both locations and exclusively comprise the largest boulders with the longest runout distances that pose the greatest hazard to life and property. This study highlights the primary influence that volcanic lithofacies architecture has on rockfall hazard. The influence of anthropogenic modifications on the surface and subsurface geologic expression of contemporary liquefaction created during the 2010-2011 Canterbury earthquake sequence (CES) in eastern Christchurch is examined. Trench observations indicate that anthropogenic fill layer boundaries and the composition/texture of discretely placed fill layers play an important role in absorbing fluidized sand/silt and controlling the subsurface architecture of preserved liquefaction features. Surface liquefaction morphologies (i.e. sand blows and linear sand blow arrays) display alignment with existing utility lines and utility excavations (and perforated pipes) provided conduits for liquefaction ejecta during the CES. No evidence of pre-CES liquefaction was identified within the anthropogenic fill layers or underlying native sediment. Radiocarbon dating of charcoal within the youngest native sediment suggests liquefaction has not occurred at the study site for at least the past 750-800 years. The importance of systematically examining the impact of buried infrastructure on channelizing and influencing surface and subsurface liquefaction morphologies is demonstrated. This thesis highlights the importance of using a multi-technique approach for understanding prehistoric and contemporary earthquake phenomena and emphasizes the critical role that humans play in shaping the geologic record and Earth’s surface processes.
Measurement of basement seismic resonance frequencies can elucidate shallow velocity structure, an important factor in earthquake hazard estimation. Ambient noise cross correlation, which is well-suited to studying shallow earth structure, is commonly used to analyze fundamental-mode Rayleigh waves and, increasingly, Love waves. Here we show via multicomponent ambient noise cross correlation that the basement resonance frequency in the Canterbury region of New Zealand can be straightforwardly determined based on the horizontal to vertical amplitude ratio (H/V ratio) of the first higher-mode Rayleigh waves. At periods of 1-3 s, the first higher-mode is evident on the radial-radial cross-correlation functions but almost absent in the vertical-vertical cross-correlation functions, implying longitudinal motion and a high H/V ratio. A one-dimensional regional velocity model incorporating a ~ 1.5 km-thick sedimentary layer fits both the observed H/V ratio and Rayleigh wave group velocity. Similar analysis may enable resonance characteristics of other sedimentary basins to be determined. © 2013. American Geophysical Union. All Rights Reserved.
This dissertation addresses a diverse range of applied aspects in ground motion simulation validation via the response of complex structures. In particular, the following topics are addressed: (i) the investigation of similarity between recorded and simulated ground motions using code-based 3D irregular structural response analysis, (ii) the development of a framework for ground motion simulations validation to identify the cause of differences between paired observed and simulated dataset, and (iii) the illustration of the process of using simulations for seismic performance-based assessment. The application of simulated ground motions is evaluated for utilisation in engineering practice by considering responses of 3D irregular structures. Validation is performed in a code-based context when the NZS1170.5 (NZS1170.5:2004, 2004) provisions are followed for response history analysis. Two real buildings designed by engineers and physically constructed in Christchurch before the 2010-2011 Canterbury earthquake sequence are considered. The responses are compared when the buildings are subjected to 40 scaled recorded and their subsequent simulated ground motions selected from 22 February 2011 Christchurch. The similarity of recorded and simulated responses is examined using statistical methods such as bootstrapping and hypothesis testing to determine whether the differences are statistically significant. The findings demonstrate the applicability of simulated ground motion when the code-based approach is followed in response history analysis. A conceptual framework is developed to link the differences between the structural response subjected to simulated and recorded ground motions to the differences in their corresponding intensity measures. This framework allows the variability to be partitioned into the proportion that can be “explained” by the differences in ground motion intensity measures and the remaining “unexplained” variability that can be attributed to different complexities such as dynamic phasing of multi-mode response, nonlinearity, and torsion. The application of this framework is examined through a hierarchy of structures reflecting a range of complexity from single-degree-of-freedom to 3D multi-degree-of-freedom systems with different materials, dynamic properties, and structural systems. The study results suggest the areas that ground motion simulation should focus on to improve simulations by prioritising the ground motion intensity measures that most clearly account for the discrepancies in simple to complex structural responses. Three approaches are presented to consider recorded or simulated ground motions within the seismic performance-based assessment framework. Considering the applications of ground motions in hazard and response history analyses, different pathways in utilising ground motions in both areas are explored. Recorded ground motions are drawn from a global database (i.e., NGA-West2 Ancheta et al., 2014). The NZ CyberShake dataset is used to obtain simulations. Advanced ground motion selection techniques (i.e., generalized conditional intensity measure, GCIM) are used for ground motion selection at a few intensity levels. The comparison is performed by investigating the response of an example structure (i.e., 12-storey reinforced concrete special moment frame) located in South Island, NZ. Results are compared and contrasted in terms of hazard, groundmotion selection, structural responses, demand hazard, and collapse risk, then, the probable reasons for differences are discussed. The findings from this study highlight the present opportunities and shortcomings in using simulations in risk assessment. i
The Porter's Pass-Amberley Fault Zone (PPAFZ) is a complex zone of anastomosing faults and folds bounding the south-eastern edge of the transition from subducting Pacific Plate to continental collision on the Australia Plate boundary. This study combines mapping of a 2000 km2 zone from the Southern Alps northeast to the coast near Amberley, 40 km north of metropolitan Christchurch, with an analysis of seismicity and a revision of regional seismic hazard. Three structural styles: 1) a western strike-slip, and 2) a more easterly thrust and reverse domain, pass into 3) a northwest verging fold belt on the northern Canterbury Plains, reflecting the structural levels exposed and the evolving west to east propagation. Basal remnants of a Late Cretaceous-Cenozoic, largely marine sedimentary cover sequence are preserved as outliers that unconformably overlie Mesozoic basement (greywacke and argillite of the Torlesse terrain) in the mountains of the PPAFZ and are underlain by a deeply leached zone which is widely preserved. Structure contouring of the unconformity surface indicates maximum, differential uplift of c.2600 m in the southwest, decreasing to c.1200 m in the coastal fold belt to the northeast. Much lower rates (or reversal) of uplift are evident a few kilometres southeast of the PPAFZ range-front escarpment. The youngest elements of the cover sequence are basement-derived conglomerates of Plio-Pleistocene age preserved on the SE margin. The source is more distant than the intervening mountains of the PPAFZ, probably from the Southern Alps, to the west and northwest. The absence of another regional unconformity on Mesozoic basement, older than Pleistocene, indicates that this uplift is post-Pliocene. Late Pleistocene(<100 kyr) differential uplift rates of c.0.5-2.7 m/kyr from uplifted marine terraces at the east coast, and rates of 2.5-3.3 m/kyr for tectonically-induced river-down cutting further west, suggest that uplift commenced locally during the last 1 Ma, and possibly within the last 0.5 Ma, if average rates are assumed to be uniform over time. Analysis of seismicity, recorded during a 10 week regional survey of micro earthquakes in 1990, identified two seismic zones beneath North Canterbury: 1) a sub-horizontal zone of activity restricted to the upper crust (≤12 km); and 2) a seismic zone in the lower crust (below a ceiling of ≤17 km), that broadens vertically to the north and northwest to a depth of c.40 km, with a bottom edge which dips 10°N and 15°NW, respectively. No events were recorded at depths between 12 km and 17 km, which is interpreted as a relatively aseismic, mid-crustal ductile layer. Marked differences (up to 60°) in the trend of strain axes for events above and below the inferred ductile layer are observed only north of the PPAFZ. A fundamental, north-to-south increase in the Wave-length of major geological structures occurs across the PPAFZ, and is interpreted as evidence that the upper crust beneath the Canterbury Plains is coupled to the lower crust, whereas the upper crust further north is not. Most of the recorded micro earthquakes <12 km deep beneath the PPAFZ have strike-slip mechanisms. It is probable that faults splay upward into the thrusts and folds at the surface as an evolving transpression zone in response to deep shear in basement. There have been no historic surface ruptures of the PPAFZ, but the zone has been characterised historically by frequent small earthquakes. Paleoseismic data (dated landslides and surface ruptures) compiled in this study, indicate a return period of 1500-1900 years between the last two M>7-7.5 earthquakes, and 500-700 years have elapsed since the last. The magnitudes of these events are estimated at c.M7.5, which represents a probable maximum magnitude for the PPAFZ. There are insufficient data to determine whether or not the frequency of large earthquakes conforms to a recognised model of behaviour, but comparison of the paleoseismic data with the historic record of smaller earthquakes, suggests that the magnitudes of the largest earthquakes in this zone are not exponentially distributed. A seismicity model for the PPAFZ (Elder et al., 1991) is reviewed, and a b-value of 1.0 is found to be consistent with the newly acquired paleoseismic data. This b-value reduces the predicted frequency of large earthquakes (M≥7.0) in this zone by a factor of 3.5, while retaining a conservative margin that allows for temporal variations in the frequency of large events and the possibility that the geological database is incomplete, suggesting grounds for revising the hazard model for Christchurch.
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.
A zone of active tectonism occurs in mid and north Canterbury, from the Rakaia to the Waipara Rivers, which coincides with seismicity concentrations and several Quaternary surface anomalies and is here defined as the Porters Pass Tectonic Zone. Although parallel to the Marlborough faults to the north, the lack of regional definition suggests this zone is much younger in its inception reflecting a southward movement of the plate rotation vector. The objectives of this study were to map the structures associated with this zone in the segment between the Rakaia and Waimakariri Rivers with detailed analysis concentrated in the upper Kawai Valley. Quaternary offsets on the main lineament of the Porters Pass Fault were traced through the area and evidence for the rate of movement, probable magnitudes and return periods of related seismic events was sought. The basement was found to be complicated by pre-existing deformation structures in Torlesse Group rocks which have been subsequently been re-activated or rotated by recent fault movement probably beginning in the Pleistocene. This phase is dominantly thrusting and uplift has lead to the erosion of most of the overlying sedimentary cover. Remnants of the Cret-Tertiary sediments still remain as fault-bounded packets. Evidence suggests that a change to development of a regional lateral shear associated with the Porters Pass Tectonic Zone transects the thrust system with complex interaction between the older reverse and new strike-slip faults. Offset rates along the segments of the Porters Pass Fault are not well constrained but are believed to be approximately in the range of 11-13 mm/year for at least the last 130,000 years. This rate is similar to other large faults in the Marlborough region. Two earthquake events have been identified and dated at 600 and 2000 years ago, with a magnitude of greater than 6.5. Evidence suggests characteristic earthquakes along the Porters Pass Fault are greater than Magnitude 7. This result has some major ramifications for the expected seismic hazards for nearby Christchurch.
Though generally considered “natural” disasters, cyclones and earthquakes are increasingly being associated with human activities, incubated through urban settlement patterns and the long-term redistribution of natural resources. As society is becoming more urbanized, the risk of human exposure to disasters is also rising. Architecture often reflects the state of society’s health: architectural damage is the first visible sign of emergency, and reconstruction is the final response in the process of recovery. An empirical assessment of architectural projects in post-disaster situations can lead to a deeper understanding of urban societies as they try to rebuild. This thesis offers an alternative perspective on urban disasters by looking at the actions and attitudes of disaster professionals through the lens of architecture, situated in recent events: the 2010 Christchurch earthquake, the 2010 Haiti earthquake, and the 2005 Hurricane Katrina. An empirical, multi-hazard, cross-sectional case study methodology was used, employing grounded theory method to build theory, and a critical constructivist strategy to inform the analysis. By taking an interdisciplinary approach to understanding disasters, this thesis positions architecture as a conduit between two divergent approaches to disaster research: the hazards approach, which studies the disaster cycles from a scientific perspective; and the sociological approach, which studies the socially constructed vulnerabilities that result from disasters, and the elements of social change that accompany such events. Few studies to date have attempted to integrate the multi-disciplinary perspectives that can advance our understanding of societal problems in urban disasters. To bridge this gap, this thesis develops what will be referred to as the “Rittelian framework”—based on the work of UC Berkeley’s architecture professor Horst Rittel (1930-1990). The Rittelian framework uses the language of design to transcend the multiple fields of human endeavor to address the “design problems” in disaster research. The processes by which societal problems are addressed following an urban disaster involve input by professionals from multiple fields—including economics, sociology, medicine, and engineering—but the contribution from architecture has been minimal to date. The main impetus for my doctoral thesis has been the assertion that most of the decisions related to reconstruction are made in the early emergency recovery stages where architects are not involved, but architects’ early contribution is vital to the long-term reconstruction of cities. This precipitated in the critical question: “How does the Rittelian framework contribute to the critical design decisions in modern urban disasters?” Comparative research was undertaken in three case studies of recent disasters in New Orleans (2005), Haiti (2010) and Christchurch (2010), by interviewing 51 individuals who were selected on the basis of employing the Rittelian framework in their humanitarian practice. Contextualizing natural disaster research within the robust methodological framework of architecture and the analytical processes of sociology is the basis for evaluating the research proposition that architectural problem solving is of value in addressing the ‘Wicked Problems’ of disasters. This thesis has found that (1) the nuances of the way disaster agents interpret the notion of “building back better” can influence the extent to which architectural professionals contribute in urban disaster recovery, (2) architectural design can be used to facilitate but also impede critical design decisions, and (3) framing disaster research in terms of design decisions can lead to innovation where least expected. This empirical research demonstrates how the Rittelian framework can inform a wider discussion about post-disaster human settlements, and improve our resilience through disaster research.
Timber has experienced renewed interests as a sustainable building material in recent times. Although traditionally it has been the prime choice for residential construction in New Zealand and some other parts of the world, its use can be increased significantly in the future through a wider range of applications, particularly when adopting engineered wood material, Research has been started on the development of innovative solutions for multi-storey non-residential timber buildings in recent years and this study is part of that initiative. Application of timber in commercial and office spaces posed some challenges with requirements of large column-free spaces. The current construction practice with timber is not properly suited for structures with the aforementioned required characteristics and new type of structures has to be developed for this type of applications. Any new structural system has to have adequate capacity for carry the gravity and lateral loads due to occupancy and the environmental effects. Along with wind loading, one of the major sources of lateral loads is earthquakes. New Zealand, being located in a seismically active region, has significant risk of earthquake hazard specially in the central region of the country and any structure has be designed for the seismic loading appropriate for the locality. There have been some significant developments in precast concrete in terms of solutions for earthquake resistant structures in the last decade. The “Hybrid” concept combining post-tensioning and energy dissipating elements with structural members has been introduced in the late 1990s by the precast concrete industry to achieve moment-resistant connections based on dry jointed ductile connections. Recent research at the University of Canterbury has shown that the concept can be adopted for timber for similar applications. Hybrid timber frames using post-tensioned beams and dissipaters have the potential to allow longer spans and smaller cross sections than other forms of solid timber frames. Buildings with post-tensioned frames and walls can have larger column-free spaces which is a particular advantage for non-residential applications. While other researchers are focusing on whole structural systems, this research concentrated on the analysis and design of individual members and connections between members or between member and foundation. This thesis extends existing knowledge on the seismic behaviour and response of post-tensioned single walls, columns under uni-direction loads and small scale beam-column joint connections into the response and design of post-tensioned coupled walls, columns under bi-directional loading and full-scale beam-column joints, as well as to generate further insight into practical applications of the design concept for subassemblies. Extensive experimental investigation of walls, column and beam-column joints provided valuable confirmation of the satisfactory performance of these systems. In general, they all exhibited almost complete re-centering capacity and significant energy dissipation, without resulting into structural damage. The different configurations tested also demonstrated the flexibility in design and possibilities for applications in practical structures. Based on the experimental results, numerical models were developed and refined from previous literature in precast concrete jointed ductile connections to predict the behaviour of post-tensioned timber subassemblies. The calibrated models also suggest the values of relevant parameters for applications in further analysis and design. Section analyses involving those parameters are performed to develop procedures to calculate moment capacities of the subassemblies. The typical features and geometric configurations the different types of subassemblies are similar with the only major difference in the connection interfaces. With adoption of appropriate values representing the corresponding connection interface and incorporation of the details of geometry and configurations, moment capacities of all the subassemblies can be calculated with the same scheme. That is found to be true for both post-tensioned-only and hybrid specimens and also applied for both uni-directional and bi-directional loading. The common section analysis and moment capacity calculation procedure is applied in the general design approach for subassemblies.
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.
Bulk rock strength is greatly dependent on fracture density, so that reductions in rock strength associated with faulting and fracturing should be reflected by reduced shear coupling and hence S-wave velocity. This study is carried out along the Canterbury rangefront and in Otago. Both lie within the broader plate boundary deformation zone in the South Island of New Zealand. Therefore built structures are often, , located in areas where there are undetected or poorly defined faults with associated rock strength reduction. Where structures are sited near to, or across, such faults or fault-zones, they may sustain both shaking and ground deformation damage during an earthquake. Within this zone, management of seismic hazards needs to be based on accurate identification of the potential fault damage zone including the likely width of off-plane deformation. Lateral S-wave velocity variability provides one method of imaging and locating damage zones and off-plane deformation. This research demonstrates the utility of Multi-Channel Analysis of Surface Waves (MASW) to aid land-use planning in such fault-prone settings. Fundamentally, MASW uses surface wave dispersive characteristics to model a near surface profile of S-wave velocity variability as a proxy for bulk rock strength. The technique can aid fault-zone planning not only by locating and defining the extent of fault-zones, but also by defining within-zone variability that is readily correlated with measurable rock properties applicable to both foundation design and the distribution of surface deformation. The calibration sites presented here have well defined field relationships and known fault-zone exposure close to potential MASW survey sites. They were selected to represent a range of progressively softer lithologies from intact and fractured Torlesse Group basement hard rock (Dalethorpe) through softer Tertiary cover sediments (Boby’s Creek) and Quaternary gravels. This facilitated initial calibration of fracture intensity at a high-velocity-contrast site followed by exploration of the limits of shear zone resolution at lower velocity contrasts. Site models were constructed in AutoCAD in order to demonstrate spatial correlations between S-wave velocity and fault zone features. Site geology was incorporated in the models, along with geomorphology, river profiles, scanline locations and crosshole velocity measurement locations. Spatial data were recorded using a total-station survey. The interpreted MASW survey results are presented as two dimensional snapshot cross-sections of the three dimensional calibration-site models. These show strong correlations between MASW survey velocities and site geology, geomorphology, fluvial profiles and geotechnical parameters and observations. Correlations are particularly pronounced where high velocity contrasts exist, whilst weaker correlations are demonstrated in softer lithologies. Geomorphic correlations suggest that off-plane deformation can be imaged and interpreted in the presence of suitable topographic survey data. A promising new approach to in situ and laboratory soft-rock material and mass characterisation is also presented using a Ramset nail gun. Geotechnical investigations typically involve outcrop and laboratory scale determination of rock mass and material properties such as fracture density and unconfined compressive strength (UCS). This multi-scale approach is espoused by this study, with geotechnical and S-wave velocity data presented at multiple scales, from survey scale sonic velocity measurements, through outcrop scale scanline and crosshole sonic velocity measurements to laboratory scale property determination and sonic velocity measurements. S-wave velocities invariably increased with decreasing scale. These scaling relationships and strategies for dealing with them are investigated and presented. Finally, the MASW technique is applied to a concealed fault on the Taieri Ridge in Macraes Flat, Central Otago. Here, high velocity Otago Schist is faulted against low velocity sheared Tertiary and Quaternary sediments. This site highlights the structural sensitivity of the technique by apparently constraining the location of the principal fault, which had been ambiguous after standard processing of the seismic reflection data. Processing of the Taieri Ridge dataset has further led to the proposal of a novel surface wave imaging technique termed Swept Frequency Imaging (SFI). This inchoate technique apparently images the detailed structure of the fault-zone, and is in agreement with the conventionally-determined fault location and an existing partial trench. Overall, the results are promising and are expected to be supported by further trenching in the near future.
Rock mass defect controlled deep-seated landslides are widespread within the deeply incised landscapes formed in Tertiary soft rock terrain in New Zealand. The basal failure surfaces of deep-seated slope failures are defined by thin, comparatively weak and laterally continuous bedding parallel layers termed critical stratigraphic horizons. These horizons have a sedimentary origin and have typically experienced some prior tectonically induced shear displacement at the time of slope failure. The key controls on the occurrence and form of deep-seated landslides are considered in terms of rock mass defect properties and tectonic and climatic forcing. The selection of two representative catchments (in southern Hawke's Bay and North Canterbury) affected by tectonic and climatic forcing has shown that the spatial and temporal initiation of deep-seated bedrock landslides in New Zealand Tertiary soft rock terrain is a predictable rather than a stochastic process; and that deep-seated landslides as a mass wasting process have a controlling role in landscape evolution in many catchments formed in Tertiary soft rock terrain. The Ella Landslide in North Canterbury is a deep-seated (~85 m) translational block slide that has failed on a 5 - 10 mm thick, kaolinite-rich, pre-sheared critical stratigraphic horizon. The residual strength of this sedimentary horizon, (C'R 2.6 - 2.7 kPa, and Ѳ'R = 16 - 21°), compared to the peak strength of the dominant lithology (C' = 176 kPa, and Ѳ' = 37°) defines a high strength contrast in the succession, and therefore a critical location for the basal failure surface of deep-seated slope failures. The (early to mid Holocene) Ella Landslide debris formed a large landslide dam in the Kate Stream catchment and this has significantly retarded rates of mass wasting in the middle catchment. Numerical stability analysis shows that this slope failure would have most likely required the influence of earthquake induced strong ground motion and the event is tentatively correlated to a Holocene event on the Omihi Fault. The influence of this slope failure is likely to affect the geomorphic development of the catchment on a scale of 10⁴ - 10⁵ years. In deeply incised catchments at the southeastern margin of the Maraetotara Plateau, southern Hawke's Bay, numerous widespread deep-seated landslides have basal failure surfaces defined by critical stratigraphic horizons in the form of thin « 20 mm) tuffaceous beds in the Makara Formation flysch (alternating sandstone and mudstone units). The geometry of deep-seated slope failures is controlled by these regularly spaced (~70 m), very weak critical stratigraphic horizons (C'R 3.8 - 14.2 kPa, and Ѳ'R = 2 - 5°), and regularly spaced (~45 m) and steeply dipping (-50°) critical conjugate joint/fault sets, which act as slide block release surfaces. Numerical stability analysis and historical precedent show that the temporal initiation of deep-seated landslides is directly controlled by short term tectonic forcing in the form of periodic large magnitude earthquakes. Published seismic hazard data shows the recurrence interval of earthquakes producing strong ground motions of 0.35g at the study site is every 150 yrs, however, if subduction thrust events are considered the level of strong ground motion may be much higher. Multiple occurrences of deep-seated slope failure are correlated to failure on the same critical stratigraphic horizon, in some cases in three adjacent catchments. Failure on multiple critical stratigraphic horizons leads to the development of a "stepped" landscape morphology. This slope form will be maintained during successive accelerated stream incision events (controlled by long term tectonic and climatic forcing) for as long as catchments are developing in this specific succession. Rock mass defect controlled deep seated landslides are controlling catchment head progression, landscape evolution and hillslope morphology in the Hawke's Bay study area and this has significant implications for the development of numerical landscape evolution models of landscapes formed in similar strata. Whereas the only known numerical model to consider deep seated landslides as an erosion process (ZSCAPE) considers them as stochastic in time and space, this study shows that this could not be applied to a landscape where the widespread spatial occurrence of deep-seated landslides is controlled by rock mass defects. In both of the study areas for this project, and by implication in many catchments in Tertiary soft rock terrain, deep-seated landslides controlled by rock mass defect strength, spacing and orientation, and tectonic and climatic forcing have an underlying control on landscape evolution. This study quantifies parameters for the development of numerical landscape evolution models that would assess the role of specific parameters, such as uplift rates, incision rates and earthquake recurrence in catchment evolution in Tertiary soft rock terrain.
Prognostic modelling provides an efficient means to analyse the coastal environment and provide effective knowledge for long term urban planning. This paper outlines how the use of SWAN and Xbeach numerical models within the ESRI ArcGIS interface can simulate geomorphological evolution through hydrodynamic forcing for the Greater Christchurch coastal environment. This research followed the data integration techniques of Silva and Taborda (2012) and utilises their beach morphological modelling tool (BeachMM tool). The statutory requirements outlined in the New Zealand Coastal Policy Statement 2010 were examined to determine whether these requirements are currently being complied with when applying the recent sea level rise predictions by the Intergovernmental Panel on Climate Change (2013), and it would appear that it does not meet those requirements. This is because coastal hazard risk has not been thoroughly quantified by the installation of the Canterbury Earthquake Recovery Authority (CERA) residential red zone. However, the Christchurch City Council’s (CCC) flood management area does provide an extent to which managed coastal retreat is a real option. This research assessed the effectiveness of the prognostic models, forecasted a coastline for 100 years from now, and simulated the physical effects of extreme events such as storm surge given these future predictions. The results of this research suggest that progradation will continue to occur along the Christchurch foreshore due to the net sediment flux retaining an onshore direction and the current hydrodynamic activity not being strong enough to move sediment offshore. However, inundation during periods of storm surge poses a risk to human habitation on low lying areas around the Avon-Heathcote Estuary and the Brooklands lagoon similar to the CCC’s flood management area. There are complex interactions at the Waimakariri River mouth with very high rates of accretion and erosion within a small spatial scale due to the river discharge. There is domination of the marine environment over the river system determined by the lack of generation of a distinct river delta, and river channel has not formed within the intertidal zone clearly. The Avon-Heathcote ebb tidal delta aggrades on the innner fan and erodes on the outer fan due to wave domination. The BeachMM tool facilitates the role of spatial and temporal analysis effectively and the efficiency of that performance is determined by the computational operating system.
