Land cover change information in urban areas supports decision makers in dealing with public policy planning and resource management. Remote sensing has been demonstrated as an efficient and accurate way to monitor land cover change over large extents. The Canterbury Earthquake Sequence (CES) caused massive damage in Christchurch, New Zealand and resulted in significant land cover change over a short time period. This study combined two types of remote sensing data, aerial imagery (RGB) and LiDAR, as the basis for quantifying land cover change in Christchurch between 2011 – 2015, a period corresponding to the five years immediately following the 22 February 2011 earthquake, which was part of the CES. An object based image analysis (OBIA) approach was adopted to classify the aerial imagery and LiDAR data into seven land cover types (bare land, building, grass, shadow, tree and water). The OBIA approach consisted of two steps, image segmentation and object classification. For the first step, this study used multi-level segmentation to better segment objects. For the second step, the random forest (RF) classifier was used to assign a land cover type to each object defined by the segmentation. Overall classification accuracies for 2011 and 2015 were 94.0% and 94.32%, respectively. Based on the classification result, land cover changes between 2011 and 2015 were then analysed. Significant increases were found in road and tree cover, while the land cover types that decreased were bare land, grass, roof, water. To better understand the reasons for those changes, land cover transitions were calculated. Canopy growth, seasonal differences and forest plantation establishment were the main reasons for tree cover increase. Redevelopment after the earthquake was the main reason for road area growth. By comparing the spatial distribution of these transitions, this study also identified Halswell and Wigram as the fastest developing suburbs in Christchurch. These results provided quantitative information for the effects of CES, with respect to land cover change. They allow for a better understanding for the current land cover status of Christchurch. Among those land cover changes, the significant increase in tree cover aroused particularly interest as urban forests benefit citizens via ecosystem services, including health, social, economic, and environmental benefits. Therefore, this study firstly calculated the percentages of tree cover in Christchurch’s fifteen wards in order to provide a general idea of tree cover change in the city extent. Following this, an automatic individual tree detection and crown delineation (ITCD) was undertaken to determine the feasibility of automated tree counting. The accuracies of the proposed approach ranged between 56.47% and 92.11% in thirty different sample plots, with an overall accuracy of 75.60%. Such varied accuracies were later found to be caused by the fixed tree detection window size and misclassifications from the land cover classification that affected the boundary of the CHM. Due to the large variability in accuracy, tree counting was not undertaken city-wide for both time periods. However, directions for further study for ITCD in Christchurch could be exploring ITCD approaches with variable window size or optimizing the classification approach to focus more on producing highly accurate CHMs.
As a global phenomenon, many cities are undergoing urban renewal to accommodate rapid growth in urban population. However, urban renewal can struggle to balance social, economic, and environmental outcomes, whereby economic outcomes are often primarily considered by developers. This has important implications for urban forests, which have previously been shown to be negatively affected by development activities. Urban forests serve the purpose of providing ecosystem services and thus are beneficial to human wellbeing. Better understanding the effect of urban renewal on city trees may help improve urban forest outcomes via effective management and policy strategies, thereby maximising ecosystem service provision and human wellbeing. Though the relationship between certain aspects of development and urban forests has received consideration in previous literature, little research has focused on how the complete property redevelopment cycle affects urban forest dynamics over time. This research provides an opportunity to gain a comprehensive understanding of the effect of residential property redevelopment on urban forest dynamics, at a range of spatial scales, in Christchurch, New Zealand following a series of major earthquakes which occurred in 2010 – 2011. One consequence of the earthquakes is the redevelopment of thousands of properties over a relatively short time-frame. The research quantifies changes in canopy cover city-wide, as well as, tree removal, retention, and planting on individual residential properties. Moreover, the research identifies the underlying reasons for these dynamics, by exploring the roles of socio-economic and demographic factors, the spatial relationships between trees and other infrastructure, and finally, the attitudes of residential property owners. To quantify the effect of property redevelopment on canopy cover change in Christchurch, this research delineated tree canopy cover city-wide in 2011 and again in 2015. An object-based image analysis (OBIA) technique was applied to aerial imagery and LiDAR data acquired at both time steps, in order to estimate city-wide canopy cover for 2011 and 2015. Changes in tree canopy cover between 2011 and 2015 were then spatially quantified. Tree canopy cover change was also calculated for all meshblocks (a relatively fine-scale geographic boundary) in Christchurch. The results show a relatively small magnitude of tree canopy cover loss, city-wide, from 10.8% to 10.3% between 2011 and 2015, but a statistically significant change in mean tree canopy cover across all the meshblocks. Tree canopy cover losses were more likely to occur in meshblocks containing properties that underwent a complete redevelopment cycle, but the loss was insensitive to the density of redevelopment within meshblocks. To explore property-scale individual tree dynamics, a mixed-methods approach was used, combining questionnaire data and remote sensing analysis. A mail-based questionnaire was delivered to residential properties to collect resident and household data; 450 residential properties (321 redeveloped, 129 non- redeveloped) returned valid questionnaires and were identified as analysis subjects. Subsequently, 2,422 tree removals and 4,544 tree retentions were identified within the 450 properties; this was done by manually delineating individual tree crowns, based on aerial imagery and LiDAR data, and visually comparing the presence or absence of these trees between 2011 and 2015. The tree removal rate on redeveloped properties (44.0%) was over three times greater than on non-redeveloped properties (13.5%) and the average canopy cover loss on redeveloped properties (52.2%) was significantly greater than on non-redeveloped properties (18.8%). A classification tree (CT) analysis was used to model individual tree dynamics (i.e. tree removal, tree retention) and candidate explanatory variables (i.e. resident and household, economic, land cover, and spatial variables). The results indicate that the model including land cover, spatial, and economic variables had the best predicting ability for individual tree dynamics (accuracy = 73.4%). Relatively small trees were more likely to be removed, while trees with large crowns were more likely to be retained. Trees were most likely to be removed from redeveloped properties with capital values lower than NZ$1,060,000 if they were within 1.4 m of the boundary of a redeveloped building. Conversely, trees were most likely to be retained if they were on a property that was not redeveloped. The analysis suggested that the resident and household factors included as potential explanatory variables did not influence tree removal or retention. To conduct a further exploration of the relationship between resident attitudes and actions towards trees on redeveloped versus non-redeveloped properties, this research also asked the landowners from the 450 properties that returned mail questionnaires to indicate their attitudes towards tree management (i.e. tree removal, tree retention, and tree planting) on their properties. The results show that residents from redeveloped properties were more likely to remove and/or plant trees, while residents from non- redeveloped properties were more likely to retain existing trees. A principal component analysis (PCA) was used to explore resident attitudes towards tree management. The results of the PCA show that residents identified ecosystem disservices (e.g. leaf litter, root damage to infrastructure) as common reasons for tree removal; however, they also noted ecosystem services as important reasons for both tree planting and tree retention on their properties. Moreover, the reasons for tree removal and tree planting varied based on whether residents’ property had been redeveloped. Most tree removal occurred on redeveloped properties because trees were in conflict with redevelopment, but occurred on non- redeveloped properties because of perceived poor tree health. Residents from redeveloped properties were more likely to plant trees due to being aesthetically pleasing or to replace trees removed during redevelopment. Overall, this research adds to, and complements, the existing literature on the effects of residential property redevelopment on urban forest dynamics. The findings of this research provide empirical support for developing specific legislation or policies about urban forest management during residential property redevelopment. The results also imply that urban foresters should enhance public education on the ecosystem services provided by urban forests and thus minimise the potential for tree removal when undertaking property redevelopment.
This is a joint Resilience Framework undertaken by the Electrical, Computer and Software Engineering Department of the University of Auckland in association with West Power and Orion networks and partially funded by the New Zealand National Science Challenge and QuakeCoRE. The Energy- Communication research group nearly accomplished two different researches focusing on both asset resilience and system resilience. Asset resilience research which covers underground cables system in Christchurch region is entitled “2010-2011 Canterbury Earthquake Sequence Impact on 11KV Underground Cables” and system resilience research which covers electricity distribution and communication system in West Coast region is entitled “NZ Electricity Distribution Network Resilience Assessment and Restoration Models following Major Natural Disturbance“. As the fourth milestone of the aforementioned research project, the latest outcome of both projects has been socialised with the stakeholders during the Cigre NZ 2019 Forum.
Welcome to the Recover newsletter Issue 4 from the Marine Ecology Research Group (MERG) of the University of Canterbury. Recover is designed to keep you updated on our MBIE-funded earthquake recovery project called RECOVER (Reef Ecology, Coastal Values & Earthquake Recovery). This 4th instalment covers recent work on seaweed recovery in the subtidal zone, ecological engineering in Waikoau / Lyell Creek, and a sneak preview of drone survey results!
Welcome to the Recover newsletter Issue 5 from the Marine Ecology Research Group (MERG) of the University of Canterbury. Recover is designed to keep you updated on our MBIE-funded earthquake recovery project called RECOVER (Reef Ecology, Coastal Values & Earthquake Recovery). This 5th instalment covers the question of how much of the coast uplifted how much, recent lab work on seaweed responses to stressors, and more on our drone survey work to quantify earthquake impacts and recovery along 130 km of coastline in the intertidal zone!
Results from a series of 1D seismic effective stress analyses of natural soil deposits from Christchurch are summarized. The analysed soil columns include sites whose performance during the 2010-2011 Canterbury earthquakes varied significantly, from no liquefaction manifestation at the ground surface to very severe liquefaction, in which case a large area of the site was covered by thick soil ejecta. Key soil profile characteristics and response mechanisms affecting the severity of surface liquefaction manifestation and subsequent damage are explored. The influence of shaking intensity on the triggering and contribution of these mechanisms is also discussed. Careful examination of the results highlights the importance of considering the deposit as a whole, i.e. a system of layers, including interactions between layers in the dynamic response and through pore water pressure redistribution and water flow.
The latest two great earthquake sequences; 2010- 2011 Canterbury Earthquake and 2016 Kaikoura Earthquake, necessitate a better understanding of the New Zealand seismic hazard condition for new building design and detailed assessment of existing buildings. It is important to note, however, that the New Zealand seismic hazard map in NZS 1170.5.2004 is generalised in effort to cover all of New Zealand and limited to a earthquake database prior to 2001. This is “common” that site-specific studies typically provide spectral accelerations different to those shown on the national map (Z values in NZS 1170.5:2004); and sometimes even lower. Moreover, Section 5.2 of Module 1 of the Earthquake Geotechnical Engineering Practice series provide the guidelines to perform site- specific studies.
We measure the longer-term effect of a major earthquake on the local economy, using night-time light intensity measured from space, and investigate whether insurance claim payments for damaged residential property affected the local recovery process. We focus on the destructive Canterbury Earthquake Sequence (CES) 2010 -2011 as our case study. Uniquely for this event, more than 95% of residential housing units were covered by insurance, but insurance payments were staggered over 5 years, enabling us to identify their local impact. We find that night-time luminosity can capture the process of recovery and describe the recovery’s determinants. We also find that insurance payments contributed significantly to the process of economic recovery after the earthquake, but delayed payments were less affective and cash settlement of claims were more effective than insurance-managed repairs in contributing to local recovery.
A major lesson from the 2011 Christchurch earthquake was the apparent lack of ductility of some lightly reinforced concrete (RC) wall structures. In particular, the structural behaviour of the critical wall in the Gallery Apartments building demonstrated that the inelastic deformation capacity of a structure, as well as potentially brittle failure of the reinforcement, is dependent on the level of bond deterioration between reinforcement and surrounding concrete that occurs under seismic loading. This paper presents the findings of an experimental study on bond behaviour between deformed reinforcing bars and the surrounding concrete. Bond strength and relative bond slip was evaluated using 75 pull-out tests under monotonic and cyclic loading. Variations of the experiments include the loading rate, loading history, concrete strength (25 to 70 MPa), concrete age, cover thickness, bar diameter (16 and 20 mm), embedded length, and the position of the embedded bond region within the specimen (deep within or close to free surface). Select test results are presented with inferred implications for RC structures.
Research following the 2010-2011 Canterbury earthquakes investigated the minimum vertical reinforcement required in RC walls to generate well distributed cracking in the plastic hinge region. However, the influence of the loading sequence and rate has not been fully addressed. The new minimum vertical reinforcement limits in NZS 3101:2006 (Amendment 3) include consideration of the material strengths under dynamic load rates, but these provisions have not been validated at a member or system level. A series of tests were conducted on RC prisms to investigate the effect of loading rate and sequence on the local behaviour of RC members. Fifteen axially loaded RC prisms with the designs representing the end region of RC walls were tested under various loading rates to cover the range of pseudo-static and earthquake loading scenarios. These tests will provide substantial data for understanding the local behaviour of RC members, including hysteretic load-deformation behaviour, crack patterns, failure mode, steel strain, strain rate and ductility. Recommendations will be made regarding the effect of loading rate and reinforcement content on the cracking behaviour and ductility of RC members.
The Avon-Heathcote Estuary is of significant value to Christchurch due to its high productivity, biotic diversity, proximity to the city, and its cultural, recreational and aesthetic qualities. Nonetheless, it has been subjected to decades of degradation from sewage wastewater discharges and encroaching urban development. The result was a eutrophied estuary, high in nitrogen, affected by large blooms of nuisance macroalgae and covered by degraded sediments. In March 2010, treated wastewater was diverted from the estuary to a site 3 km offshore. This quickly reduced water nitrogen by 90% within the estuary and, within months, there was reduced production of macroalgae. However, a series of earthquakes beginning in September 2010 brought massive changes: tilting of the estuary, changes in channels and water flow, and a huge influx of liquefied sediments that covered up to 65% of the estuary floor. Water nitrogen increased due to damage to sewage infrastructure and the diversion pipeline being turned off. Together, these drastically altered the estuarine ecosystem. My study involves three laboratory and five in situ experiments that investigate the base of the food chain and responses of benthic microalgae to earthquake-driven sediment and nutrient changes. It was predicted that the new sediments would be coarser and less contaminated with organic matter and nutrients than the old sediments, would have decreased microalgal biomass, and would prevent invertebrate grazing and bioturbation activities. It was believed that microalgal biomass would become similar across new and old sediments types as the unstable new sediments were resuspended and distributed over the old sediments. Contact cores of the sediment were taken at three sites, across a eutrophication gradient, monthly from September 2011 to March 2012. Extracted chlorophyll a pigments showed that microalgal biomass was generally lower on new liquefied sediments compared to old sediments, although there was considerable site to site variation, with the highly eutrophic sites being the most affected by the emergence of the new sediments. Grazer experiments showed that invertebrates had both positive and negative site-specific effects on microalgal biomass depending on their identity. At one site, new sediments facilitated grazing by Amphibola crenata, whereas at another site, new sediments did not alter the direct and indirect effects of invertebrates (Nicon aestuariensis, Macropthalmus hirtipes, and A. crenata) on microalgae. From nutrient addition experiments it was clear that benthic microalgae were able to use nutrients from within both old and new sediments equally. This implied that microalgae were reducing legacy nutrients in both sediments, and that they are an important buffer against eutrophication. Therefore, in tandem with the wastewater diversion, they could underpin much of the recovery of the estuary. Overall, the new sediments were less favourable for benthic microalgal growth and recolonisation, but were less contaminated than old sediments at highly eutrophic sites. Because the new sediments were less contaminated than the old sediments, they could help return the estuary to a noneutrophic state. However, if the new sediments, which are less favourable for microalgal growth, disperse over the old sediments at highly eutrophic sites, they could become contaminated and interfere with estuarine recovery. Therefore, recovery of microalgal communities and the estuary was expected to be generally long, but variable and site-specific, with the least eutrophic sites recovering quickly, and the most eutrophic sites taking years to return to a pre-earthquake and non-eutrophied state. changes in channels and water flow, and a huge influx of liquefied sediments that covered up to 65% of the estuary floor. Water nitrogen increased due to damage to sewage infrastructure and the diversion pipeline being turned off. Together, these drastically altered the estuarine ecosystem. My study involves three laboratory and five in situ experiments that investigate the base of the food chain and responses of benthic microalgae to earthquake-driven sedimen tand nutrient changes. It was predicted that the new sediments would be coarser and less contaminated with organic matter and nutrients than the old sediments, would have decreased microalgal biomass, and would prevent invertebrate grazing and bioturbation activities. It was believed that microalgal biomass would become similar across new and old sediments types as the unstable new sediments were resuspended and distributed over the old sediments. Contact cores of the sediment were taken at three sites, across a eutrophication gradient, monthly from September 2011 to March 2012. Extracted chlorophyll a pigments showed that microalgal biomass was generally lower on new liquefied sediments compared to old sediments, although there was considerable site to site variation, with the highly eutrophic sites being the most affected by the emergence of the new sediments. Grazer experiments showed that invertebrates had both positive and negative site-specific effects on microalgal biomass depending on their identity. At one site, new sediments facilitated grazing by Amphibola crenata, whereas at another site, new sediments did not alter the direct and indirect effects of invertebrates (Nicon aestuariensis, Macropthalmus hirtipes, and A. crenata) on microalgae. From nutrient addition experiments it was clear that benthic microalgae were able to use nutrients from within both old and new sediments equally. This implied that microalgae were reducing legacy nutrients in both sediments, and that they are
Nowadays the telecommunication systems’ performance has a substantial impact on our lifestyle. Their operationality becomes even more substantial in a post-disaster scenario when these services are used in civil protection and emergency plans, as well as for the restoration of all the other critical infrastructure. Despite the relevance of loss of functionality of telecommunication networks on seismic resilience, studies on their performance assessment are few in the literature. The telecommunication system is a distributed network made up of several components (i.e. ducts, utility holes, cabinets, major and local exchanges). Given that these networks cover a large geographical area, they can be easily subjected to the effects of a seismic event, either the ground shaking itself, or co-seismic events such as liquefaction and landslides. In this paper, an analysis of the data collected after the 2010-2011 Canterbury Earthquake Sequence (CES) and the 2016 Kaikoura Earthquake in New Zealand is conducted. Analysing these data, information gaps are critically identified regarding physical and functional failures of the telecommunication components, the timeline of repair/reconstruction activities and service recovery, geotechnical tests and land planning maps. Indeed, if these missing data were presented, they could aid the assessment of the seismic resilience. Thus, practical improvements in the post-disaster collection from both a network and organisational viewpoints are proposed through consultation of national and international researchers and highly experienced asset managers from Chorus. Finally, an outline of future studies which could guide towards a more resilient seismic performance of the telecommunication network is presented.
This poster presents work to date on ground motion simulation validation and inversion for the Canterbury, New Zealand region. Recent developments have focused on the collection of different earthquake sources and the verification of the SPECFEM3D software package in forward and inverse simulations. SPECFEM3D is an open source software package which simulates seismic wave propagation and performs adjoint tomography based upon the spectral-element method. Figure 2: Fence diagrams of shear wave velocities highlighting the salient features of the (a) 1D Canterbury velocity model, and (b) 3D Canterbury velocity model. Figure 5: Seismic sources and strong motion stations in the South Island of New Zealand, and corresponding ray paths of observed ground motions. Figure 3: Domain used for the 19th October 2010 Mw 4.8 case study event including the location of the seismic source and strong motion stations. By understanding the predictive and inversion capabilities of SPECFEM3D, the current 3D Canterbury Velocity Model can be iteratively improved to better predict the observed ground motions. This is achieved by minimizing the misfit between observed and simulated ground motions using the built-in optimization algorithm. Figure 1 shows the Canterbury Velocity Model domain considered including the locations of small-to-moderate Mw events [3-4.5], strong motion stations, and ray paths of observed ground motions. The area covered by the ray paths essentially indicates the area of the model which will be most affected by the waveform inversion. The seismic sources used in the ground motion simulations are centroid moment tensor solutions obtained from GeoNet. All earthquake ruptures are modelled as point sources with a Gaussian source time function. The minimum Mw limit is enforced to ensure good signal-to-noise ratio and well constrained source parameters. The maximum Mw limit is enforced to ensure the point source approximation is valid and to minimize off-fault nonlinear effects.
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.
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.
In September 2010 and February 2011, the Canterbury region experienced devastating earthquakes with an estimated economic cost of over NZ$40 billion (Parker and Steenkamp, 2012; Timar et al., 2014; Potter et al., 2015). The insurance market played an important role in rebuilding the Canterbury region after the earthquakes. Homeowners, insurance and reinsurance markets and New Zealand government agencies faced a difficult task to manage the rebuild process. From an empirical and theoretic research viewpoint, the Christchurch disaster calls for an assessment of how the insurance market deals with such disasters in the future. Previous studies have investigated market responses to losses in global catastrophes by focusing on the insurance supply-side. This study investigates both demand-side and supply-side insurance market responses to the Christchurch earthquakes. Despite the fact that New Zealand is prone to seismic activities, there are scant previous studies in the area of earthquake insurance. This study does offer a unique opportunity to examine and document the New Zealand insurance market response to catastrophe risk, providing results critical for understanding market responses after major loss events in general. A review of previous studies shows higher premiums suppress demand, but how higher premiums and a higher probability of risk affect demand is still largely unknown. According to previous studies, the supply of disaster coverage is curtailed unless the market is subsidised, however, there is still unsettled discussion on why demand decreases with time from the previous disaster even when the supply of coverage is subsidised by the government. Natural disaster risks pose a set of challenges for insurance market players because of substantial ambiguity associated with the probability of such events occurring and high spatial correlation of catastrophe losses. Private insurance market inefficiencies due to high premiums and spatially concentrated risks calls for government intervention in the provision of natural disaster insurance to avert situations of noninsurance and underinsurance. Political economy considerations make it more likely for government support to be called for if many people are uninsured than if few people are uninsured. However, emergency assistance for property owners after catastrophe events can encourage most property owners to not buy insurance against natural disaster and develop adverse selection behaviour, generating larger future risks for homeowners and governments. On the demand-side, this study has developed an intertemporal model to examine how demand for insurance changes post-catastrophe, and how to model it theoretically. In this intertemporal model, insurance can be sought in two sequential periods of time, and at the second period, it is known whether or not a loss event happened in period one. The results show that period one demand for insurance increases relative to the standard single period model when the second period is taken into consideration, period two insurance demand is higher post-loss, higher than both the period one demand and the period two demand without a period one loss. To investigate policyholders experience from the demand-side perspective, a total of 1600 survey questionnaires were administered, and responses from 254 participants received representing a 16 percent response rate. Survey data was gathered from four institutions in Canterbury and is probably not representative of the entire population. The results of the survey show that the change from full replacement value policy to nominated replacement value policy is a key determinant of the direction of change in the level of insurance coverage after the earthquakes. The earthquakes also highlighted the plight of those who were underinsured, prompting policyholders to update their insurance coverage to reflect the estimated cost of re-building their property. The survey has added further evidence to the existing literature, such as those who have had a recent experience with disaster loss report increased risk perception if a similar event happens in future with females reporting a higher risk perception than males. Of the demographic variables, only gender has a relationship with changes in household cover. On the supply-side, this study has built a risk-based pricing model suitable to generate a competitive premium rate for natural disaster insurance cover. Using illustrative data from the Christchurch Red-zone suburbs, the model generates competitive premium rates for catastrophe risk. When the proposed model incorporates the new RMS high-definition New Zealand Earthquake Model, for example, insurers can find the model useful to identify losses at a granular level so as to calculate the competitive premium. This study observes that the key to the success of the New Zealand dual insurance system despite the high prevalence of catastrophe losses are; firstly the EQC’s flat-rate pricing structure keeps private insurance premiums affordable and very high nationwide homeowner take-up rates of natural disaster insurance. Secondly, private insurers and the EQC have an elaborate reinsurance arrangement in place. By efficiently transferring risk to the reinsurer, the cost of writing primary insurance is considerably reduced ultimately expanding primary insurance capacity and supply of insurance coverage.
The Mw 6.2 February 22nd 2011 Christchurch earthquake (and others in the 2010-2011 Canterbury sequence) provided a unique opportunity to study the devastating effects of earthquakes first-hand and learn from them for future engineering applications. All major events in the Canterbury earthquake sequence caused widespread liquefaction throughout Christchurch’s eastern suburbs, particularly extensive and severe during the February 22nd event. Along large stretches of the Avon River banks (and to a lesser extent along the Heathcote) significant lateral spreading occurred, affecting bridges and the infrastructure they support. The first stage of this research involved conducting detailed field reconnaissance to document liquefaction and lateral spreading-induced damage to several case study bridges along the Avon River. The case study bridges cover a range of ages and construction types but all are reinforced concrete structures which have relatively short, stiff decks. These factors combined led to a characteristic deformation mechanism involving deck-pinning and abutment back-rotation with consequent damage to the abutment piles and slumping of the approaches. The second stage of the research involved using pseudo-static analysis, a simplified seismic modelling tool, to analyse two of the bridges. An advantage of pseudo-static analysis over more complicated modelling methods is that it uses conventional geotechnical data in its inputs, such as SPT blowcount and CPT cone resistance and local friction. Pseudo-static analysis can also be applied without excessive computational power or specialised knowledge, yet it has been shown to capture the basic mechanisms of pile behaviour. Single pile and whole bridge models were constructed for each bridge, and both cyclic and lateral spreading phases of loading were investigated. Parametric studies were carried out which varied the values of key parameters to identify their influence on pile response, and computed displacements and damages were compared with observations made in the field. It was shown that pseudo-static analysis was able to capture the characteristic damage mechanisms observed in the field, however the treatment of key parameters affecting pile response is of primary importance. Recommendations were made concerning the treatment of these governing parameters controlling pile response. In this way the future application of pseudo-static analysis as a tool for analysing and designing bridge pile foundations in liquefying and laterally spreading soils is enhanced.
Several concrete cladding panels were damaged during the 2011 Christchurch Earthquakes in New Zealand. Damage included partial collapse of panels, rupture of joint sealants, cracking and corner crushing. Installation errors, faulty connections and inadequate detailing were also contributing factors to the damage. In New Zealand, two main issues are considered in order to accommodate story drifts in the design of precast cladding panels: 1) drift compatibility of tieback or push-pull connections and 2) drift compatibility of corner joints. Tieback connections restrain the panels in the out-of-plane direction while allowing in-plane translation with respect to the building frame. Tieback connections are either in the form of slots or oversized holes or ductile rods usually located at the top of the panels. Bearing connections are also provided at the bottom of panels to transfer gravity loads. At the corners of a building, a vertical joint gap, usually filled with sealants, is provided between the two panels on the two orthogonal sides to accommodate the relative movement. In cases where the joint gap is not sufficient to accommodate the relative movements, panels can collide, generating large forces and the likely failure of the connections. On the other hand, large gaps are aesthetically unpleasing. The current design standards appear to recognize these issues but then leave most of the design and detailing to the discretion of the designers. In the installation phase, the alignment of panels is one of the main challenges faced by installers (and/or contractors). Many prefer temporary props to guide, adjust and hold the panels in place whilst the bearing connections are welded. Moreover, heat generated from extensive welding can twist the steel components inducing undesirable local stresses in the panels. Therefore, the installation phase itself is time-consuming, costly and prone to errors. This paper investigates the performance of a novel panel system that is designed to accommodate lateral inter-story drift through a ‘rocking’ motion. In order to gauge the feasibility of the system, six 2m high precast concrete panels within a single-story steel frame structure have been tested under increasing levels of lateral cyclic drift at the University of Canterbury, New Zealand. Three different panel configurations are tested: 1) a panel with return cover and a flat panel at a corner under unidirectional loading, 2) Two adjacent flat panels under unidirectional loading, and 3) Two flat panels at another oblique corner under bidirectional loading. A vertical seismic joint of 25 mm, filled with one-stage joint sealant, is provided between two of the panels. The test results show the ability of the panels with ‘rocking’ connection details to accommodate larger lateral drifts whilst allowing for smaller vertical joints between panels at corners, quick alignment and easy placement of panels without involving extensive welding on site.
Worldwide turbidity is a huge concern for the health of aquatic ecosystems. Human activities on the land such as construction, deforestation, agriculture, and mining all have impacts on the amount of particulate solids that enter the world’s waterways. These particulate solids can pose a number of risks to aquatic life, but primary among them is the turbidity that they create in the water column. The way suspended solids interact with light creates cloudiness in the water which interferes with the vision, and visually mediated behaviours of aquatic organisms, particularly fish. The Avon-Heathcote estuary of Christchurch, New Zealand, is one such body of water that is subject to tremendous variation in turbidity, no doubt exacerbated by the destruction of Christchurch in the 2010 and 2011 earthquakes, as well as the subsequent ongoing rebuild. The yellow eyed mullet, Aldrichetta Forsteri, is one species that is common with the estuary, and uses it as a habitat for breeding. Though very common throughout New Zealand, and even a part of the catch of commercial fisheries, the yellow eyed mullet is a largely unstudied organism, with virtually no published scientific enquiry based on the species. The present work assesses how several behaviours of the yellow eyed mullet are effected by acute turbidity at 10, 50, 90, 130 and 170 NTU, finding that: 1) The optomotor response of mullet to 2.5 mm stripes drops to insignificant levels between 10 and 50 NTU, 2) The swimming activity of the yellow eyed mullet is highest at 10 NTU and drops to a significantly lower level at higher turbidities, 3) The grouping behaviour of small groups of yellow eyed mullet are unchanged by increasing turbidity levels, 4) that yellow eyed mullet do not exhibit significantly different behavioural response to a simulated predator at any of the tested turbidities, and 5) that yellow eyed mullet to do significantly alter their oxygen consumption during exposure to the turbidities in an increasing series. The results presented in these studies indicate that turbidites above 50 NTU pose a significant risk to the lifestyle of the yellow eyed mullet, potentially impacting their ability to perceive their surroundings, feed, school, and avoid predation. Future work has a lot of ground to cover to more precisely determine the relationship between yellow eyed mullet behaviour and physiology, and the turbidity of their environment. In particular, future work should focus more closely on the turbidities between 10 and 50 NTU, as well as looking to field work to see what the predominant predators of the mullet are, and specifically whether turbidity increases or decreases the risk of mullet being subject to avian predation. There is also considerable scope for studies on the effects of chronic turbidity upon mullet, which will add understand to the predicament of escalating turbidity and its effects upon this common and yet mysterious native fish.
“much of what we know about leadership is today redundant because it is literally designed for a different operating model, a different context, a different time” (Pascale, Sternin, & Sternin, p. 4). This thesis describes a project that was designed with a focus on exploring ways to enhance leadership capacity in non-government organisations operating in Christchurch, New Zealand. It included 20 CEOs, directors and managers from organisations that cover a range of settings, including education, recreation, and residential and community therapeutic support; all working with adolescents. The project involved the creation of a peer-supported professional learning community that operated for 14 months; the design and facilitation of which was informed by the Appreciative Inquiry principles of positive focus and collaboration. At the completion of the research project in February 2010, the leaders decided to continue their collective processes as a self-managing and sustaining professional network that has grown and in 2014 is still flourishing under the title LYNGO (Leaders of Youth focussed NGOs). Two compelling findings emerged from this research project. The first of these relates to efficacy of a complexity thinking framework to inform the actions of these leaders. The leaders in this project described the complexity thinking framework as the most relevant, resonant and dynamic approach that they encountered throughout the research project. As such this thesis explores this complexity thinking informed leadership in detail as the leaders participating in this project believed it offers an opportune alternative to more traditional forms of positional leadership and organisational approaches. This exploration is more than simply a rationale for complexity thinking but an iterative in-depth exploration of ‘complexity leadership in action’ which in Chapter 6 elaborates on detailed leadership tools and frameworks for creating the conditions for self-organisation and emergence. The second compelling finding relates to efficacy of Appreciative Inquiry as an emergent research and development process for leadership learning. In particular the adoption of two key principles; positive focus and inclusivity were beneficial in guiding the responsive leadership learning process that resulted in a professional learning community that exhibited high engagement and sustainability. Additionally, the findings suggest that complexity thinking not only acts as a contemporary framework for adaptive leadership of organisations as stated above; but that complexity thinking has much to offer as a framework for understanding leadership development processes through the application of Appreciative Inquiry (AI)-based principles. A consideration of the components associated with complexity thinking has promise for innovation and creativity in the development of leaders and also in the creation of networks of learning. This thesis concludes by suggesting that leaders focus on creating hybrid organisations, ones which leverage the strengths (and minimise the limitations) of self-organising complexity-informed organisational processes, while at the same time retaining many of the strengths of more traditional organisational management structures. This approach is applied anecdotally to the place where this study was situated: the post-earthquake recovery of Christchurch, New Zealand.
This thesis focuses on the role of legal preparedness for managing large-scale urban disasters in Aotearoa New Zealand. It uses the Auckland Volcanic Field as a case study to answer the question: ‘is New Zealand’s current legal framework prepared to respond to and recover from a large-scale urban disaster?’. The Auckland Volcanic Field was chosen as the main case study because a future eruption is a low likelihood, high-impact event that New Zealand is going to have to manage in the future. Case studies are a key feature of this thesis as both New Zealand based and overseas examples are used to explore the role of legal preparedness by identifying and investigating a range of legal issues that need to be addressed in advance of a future Auckland Volcanic Field eruption. Of particular interest is the impact of legal preparedness for the recovery phase. The New Zealand case studies include; Canterbury earthquake sequence 2010-2011, the Kaikōura earthquake 2016, the Auckland flooding 2018, and the North Island Severe Weather event 2023, which encompasses both the Auckland Anniversary weekend flooding and Cyclone Gabrielle. As New Zealand has not experienced a large-scale urban volcanic eruption, overseas examples are explored to provide insights into the legal issues that are volcano specific. The overseas volcanic case studies cover eruptions in Heimaey (Iceland), the Soufrière Hills (Montserrat and the Grenadines), La Soufrière (St Vincent) and Tungurahua (Ecuador). New Zealand’s past experiences highlight a trend for introducing post-event legal frameworks to manage recovery. Consequently, the current disaster management system is not prioritising legal preparedness and instead is choosing to rely on exceptional powers. Unsurprisingly, the introduction of new post-event recovery frameworks has repercussions. In New Zealand, new post-event legal frameworks are introduced swiftly under urgency, they contain broad unstructured decision-making powers, and are often flawed. As these exceptional new frameworks sit outside the ‘normal’ legal frameworks, they in effect create a parallel “shadow system”. Based on the evidence explored in this thesis it does not appear that Auckland’s current disaster management framework is prepared to deal with a large-scale urban event caused by an Auckland Volcanic Field eruption. Following this conclusion, it is the submission of this thesis that New Zealand’s current legal framework is not prepared to respond to and recover from a large-scale urban disaster. To become legally prepared, New Zealand needs to consider the legal tools required to manage large-scale urban disasters in advance. This will prevent the creation of a legal vacuum in the aftermath of disasters and the need for new recovery frameworks. Adopting a new attitude will require a change in approach towards legal preparedness which prioritises it, rather than sidelining it. This may also require changes within New Zealand’s disaster management system including the introduction of a formal monitoring mechanism, which will support and prioritise legal preparedness. This thesis has shown that not legally preparing for future disasters is a choice which carries significant consequences. None of these consequences are inevitable when managing large-scale disasters, however they are inevitable when frameworks are not legally prepared in advance. To not legally prepare, is to prepare to fail and thus create a disaster by choice.
Structures of the Lowry Peaks Range - Waikari Valley district are complex. The majority comprise three members of a predominantly WSW -ENE striking major northwards-directed, leading edge imbricate thrust system, with associated angular, asymmetric fault-propagation folds. This system forms anomalously within a large NESW trending belt of structures characterising the entire east coast of north Canterbury, both onshore and offshore and terminates westwards against N-S striking, east facing fold-fault zone. The objectives of this study address the origin, geometry and kinematics of the interaction between these diversely trending systems. Stratigraphy and small-scale structures denote three periods of deformation, namely: i) Middle Cretaceous deformation of the basement rocks, ii) weak Middle Oligocene deformation associated with the inception of the plate boundary through the South Island, and iii) major Pliocene - Recent deformation that formed the majority of the above-mentioned structures. Stress tensor analyses within competent basement and limestone cover rocks suggest two sets of sub-horizontal compression, NE-SW and NW-SE, the former likely to relate to a localised earlier period of deformation, now overprinted by the latter. NW-SE oriented sub-horizontal compression correlates well with results from other parts of north Canterbury. The result of NW-SE compression on the W-E to WSW-ENE striking structures is a large component of oblique motion, which is manifest in four ways: i) movement on two, differently oriented splays rather than a single fault strand, ii) the development of a sinuous trace for a number of the major folds, whereby the ends are oriented normal to the compression direction, the centres parallel to the strike of the faults, iii) the development of a number of cross-folds, striking NNE-SSW and iv) the apparently recent development of a strike-slip component on at least one of the major thrust faults. The origin of the W-E, or WSW-ENE striking structures may be reactivation of Late Cretaceous faults, stratigraphic evidence for the existence of a "structural high" (the Hurunui High) over the majority of the area in the Late Cretaceous to Early Eocene times suggests the formation of a W-E trending horst structure, with a corresponding asymmetric graben to the south. The junction of WSW-ENE trending structures with N-S trending structures to the west centres on an alluvial-filled depression, Waikari Flat, into which the structures of the WSW-ENE trending imbricate thrust system plunge, locally curling to the SW at their ends to link with N-S trending structures to the south. Roof thrusting on two orientations, W-E and N-S, towards to SE is currently occurring above these structures. Currently the area is not highly seismically active, although a magnitude ~6.4 Ms earthquake in historic times has been recorded. The effects of tectonics on the drainage of the area does suggest that the majority of the systems, are still potentially active, albeit moving at a comparatively slow rate. The majority of the recent motion appears to be concentrated on the roof-thrusting occurring in Waikari Flat, and uplift along the Lowry Peaks Fault System. Increasing amounts of secondary movement on back-thrusts and cross fractures is also implied for western ends of the major imbricate thrust system. In contrast, the southern-most fault system appears to be largely sustaining dextral strike-slip motion, with some local folding in central portions.
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.
The Canterbury Earthquake Sequence (CES) of 2010-2011 produced large seismic moments up to Mw 7.1. These large, near-to-surface (<15 km) ruptures triggered >6,000 rockfall boulders on the Port Hills of Christchurch, many of which impacted houses and affected the livelihoods of people within the impacted area. From these disastrous and unpredicted natural events a need arose to be able to assess the areas affected by rockfall events in the future, where it is known that a rockfall is possible from a specific source outcrop but the potential boulder runout and dynamics are not understood. The distribution of rockfall deposits is largely constrained by the physical properties and processes of the boulder and its motion such as block density, shape and size, block velocity, bounce height, impact and rebound angle, as well as the properties of the substrate. Numerical rockfall models go some way to accounting for all the complex factors in an algorithm, commonly parameterised in a user interface where site-specific effects can be calibrated. Calibration of these algorithms requires thorough field checks and often experimental practises. The purpose of this project, which began immediately following the most destructive rupture of the CES (February 22, 2011), is to collate data to characterise boulder falls, and to use this information, supplemented by a set of anthropogenic boulder fall data, to perform an in-depth calibration of the three-dimensional numerical rockfall model RAMMS::Rockfall. The thesis covers the following topics: • Use of field data to calibrate RAMMS. Boulder impact trails in the loess-colluvium soils at Rapaki Bay have been used to estimate ranges of boulder velocities and bounce heights. RAMMS results replicate field data closely; it is concluded that the model is appropriate for analysing the earthquake-triggered boulder trails at Rapaki Bay, and that it can be usefully applied to rockfall trajectory and hazard assessment at this and similar sites elsewhere. • Detailed analysis of dynamic rockfall processes, interpreted from recorded boulder rolling experiments, and compared to RAMMS simulated results at the same site. Recorded rotational and translational velocities of a particular boulder show that the boulder behaves logically and dynamically on impact with different substrate types. Simulations show that seasonal changes in soil moisture alter rockfall dynamics and runout predictions within RAMMS, and adjustments are made to the calibration to reflect this; suggesting that in hazard analysis a rockfall model should be calibrated to dry rather than wet soil conditions to anticipate the most serious outcome. • Verifying the model calibration for a separate site on the Port Hills. The results of the RAMMS simulations show the effectiveness of calibration against a real data set, as well as the effectiveness of vegetation as a rockfall barrier/retardant. The results of simulations are compared using hazard maps, where the maximum runouts match well the mapped CES fallen boulder maximum runouts. The results of the simulations in terms of frequency distribution of deposit locations on the slope are also compared with those of the CES data, using the shadow angle tool to apportion slope zones. These results also replicate real field data well. Results show that a maximum runout envelope can be mapped, as well as frequency distribution of deposited boulders for hazard (and thus risk) analysis purposes. The accuracy of the rockfall runout envelope and frequency distribution can be improved by comprehensive vegetation and substrate mapping. The topics above define the scope of the project, limiting the focus to rockfall processes on the Port Hills, and implications for model calibration for the wider scientific community. The results provide a useful rockfall analysis methodology with a defensible and replicable calibration process, that has the potential to be applied to other lithologies and substrates. Its applications include a method of analysis for the selection and positioning of rockfall countermeasure design; site safety assessment for scaling and demolition works; and risk analysis and land planning for future construction in Christchurch.
One of the most controversial issues highlighted by the 2010-2011 Christchurch earthquake series and more recently the 2016 Kaikoura earthquake, has been the evident difficulty and lack of knowledge and guidelines for: a) evaluation of the residual capacity damaged buildings to sustain future aftershocks; b) selection and implementation of a series of reliable repairing techniques to bring back the structure to a condition substantially the same as prior to the earthquake; and c) predicting the cost (or cost-effectiveness) of such repair intervention, when compared to fully replacement costs while accounting for potential aftershocks in the near future. As a result of such complexity and uncertainty (i.e., risk), in combination with the possibility (unique in New Zealand when compared to most of the seismic-prone countries) to rely on financial support from the insurance companies, many modern buildings, in a number exceeding typical expectations from past experiences at an international level, have ended up being demolished. This has resulted in additional time and indirect losses prior to the full reconstruction, as well as in an increase in uncertainty on the actual relocation of the investment. This research project provides the main end-users and stakeholders (practitioner engineers, owners, local and government authorities, insurers, and regulatory agencies) with comprehensive evidence-based information to assess the residual capacity of damage reinforced concrete buildings, and to evaluate the feasibility of repairing techniques, in order to support their delicate decision-making process of repair vs. demolition or replacement. Literature review on effectiveness of epoxy injection repairs, as well as experimental tests on full-scale beam-column joints shows that repaired specimens have a reduced initial stiffness compared with the undamaged specimen, with no apparent strength reduction, sometimes exhibiting higher displacement ductility capacities. Although the bond between the steel and concrete is only partially restored, it still allows the repaired specimen to dissipate at least the same amount of hysteretic energy. Experimental tests on buildings subjected to earthquake loading demonstrate that even for severe damage levels, the ability of the epoxy injection to restore the initial stiffness of the structure is significant. Literature review on damage assessment and repair guidelines suggests that there is consensus within the international community that concrete elements with cracks less than 0.2 mm wide only require cosmetic repairs; epoxy injection repairs of cracks less and 2.0 mm wide and concrete patching of spalled cover concrete (i.e., minor to moderate damage) is an appropiate repair strategy; and for severe damaged components (e.g., cracks greater than 2.0 mm wide, crushing of the concrete core, buckling of the longitudinal reinforcement) local replacement of steel and/or concrete in addition to epoxy crack injection is more appropriate. In terms of expected cracking patterns, non-linear finite element investigations on well-designed reinforced concrete beam-to-column joints, have shown that lower number of cracks but with wider openings are expected to occur for larger compressive concrete strength, f’c, and lower reinforcement content, ρs. It was also observed that the tensile concrete strength, ft, strongly affects the expected cracking pattern in the beam-column joints, the latter being more uniformly distributed for lower ft values. Strain rate effects do not seem to play an important role on the cracking pattern. However, small variations in the cracking pattern were observed for low reinforcement content as it approaches to the minimum required as per NZS 3101:2006. Simple equations are proposed in this research project to relate the maximum and residual crack widths with the steel strain at peak displacement, with or without axial load. A literature review on fracture of reinforcing steel due to low-cycle fatigue, including recent research using steel manufactured per New Zealand standards is also presented. Experimental results describing the influence of the cyclic effect on the ultimate strain capacity of the steel are also discussed, and preliminary equations to account for that effect are proposed. A literature review on the current practice to assess the seismic residual capacity of structures is also presented. The various factors affecting the residual fatigue life at a component level (i.e., plastic hinge) of well-designed reinforced concrete frames are discussed, and equations to quantify each of them are proposed, as well as a methodology to incorporate them into a full displacement-based procedure for pre-earthquake and post-earthquake seismic assessment.
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.
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.
The Eastern Humps and Leader faults, situated in the Mount Stewart Range in North Canterbury, are two of the ≥17 faults which ruptured during the 2016 MW7.8 Kaikōura Earthquake. The earthquake produced complex, intersecting ground ruptures of these faults and the co-seismic uplift of the Mount Stewart Range. This thesis aims to determine how these two faults accommodated deformation during the 2016 earthquake and how they interact with each other and with pre-existing geological structures. In addition, it aims to establish the most likely subsurface geometry of the fault complex across the Mount Stewart Range, and to investigate the paleoseismic history of the Leader Fault. The Eastern Humps Fault strikes ~240° and dips 80° to 60° to the northwest and accommodated right- lateral – reverse-slip, with up to 4 m horizontal and 2 m vertical displacement in the 2016 earthquake. The strike of the Leader Fault varies from ~155 to ~300°, and dips ~30 to ~80° to the west/northwest, and mainly accommodated left-lateral – reverse-slip of up to 3.5 m horizontal and 3.5 m vertical slip in the 2016 earthquake. On both the Eastern Humps and Leader faults the slip is variable along strike, with areas of low total displacement and areas where horizontal and vertical displacement are negatively correlated. Fault traces with low total displacement reflect the presence of off-fault (distributed) displacement which is not being captured with field measurements. The negative correlation of horizontal and vertical displacement likely indicates a degree of slip partitioning during the 2016 earthquake on both the Eastern Humps and Leader faults. The Eastern Humps and Leader faults have a complex, interdependent relationship with the local bedrock geology. The Humps Fault appears to be a primary driver of ongoing folding and deformation of the local Mendip Syncline and folding of the Mount Stewart Range, which probably began prior to, or synchronous with, initial rupture of The Humps Fault. The Leader Fault appears to use existing lithological weaknesses in the Cretaceous-Cenozoic bedrock stratigraphy to rupture to the surface. This largely accounts for the strong variability on the strike and dip of the Leader Fault, as the geometry of the surface ruptures tend to reflect the strike and dip of the geological strata which it is rupturing through. The Leader Fault may also accommodate some degree of flexural slip in the Cenozoic cover sequence of the Mendip Syncline, contributing to the ongoing growth of the fold. The similarity between topography and uplift profiles from the 2016 earthquake suggest that growth of the Mount Stewart Range has been primarily driven by multiple (>500) discrete earthquakes that rupture The Humps and Leader faults. The spatial distribution of surface displacements across the Mount Stewart Range is more symmetrical than would be expected if uplift is driven primarily by The Humps and Leader faults alone. Elastic dislocation forward models were used to model potential sub-surface geometries and the resulting patterns of deformation compared to photogrammetry-derived surface displacements. Results show a slight preference for models with a steeply southeast-dipping blind fault, coincident with a zone of seismicity at depth, as a ‘backthrust’ to The Humps and Leader faults. This inferred Mount Stewart Fault accommodated contractional strain during the 2016 earthquake and contributes to the ongoing uplift of the Mount Stewart Range with a component of folding. Right-lateral and reverse shear stress change on the Hope Fault was also modelled using Coulomb 3.3 software to examine whether slip on The Humps and Leader faults could transfer enough stress onto the Hope Fault to trigger through-going rupture. Results indicate that during the 2016 earthquake right-lateral shear and reverse stress only increased on the Hope Fault in small areas to the west of the Leader Fault, and similar ruptures would be unlikely to trigger eastward propagating rupture unless the Hope Fault was close to failure prior to the earthquake. Paleoseismic trenches were excavated on the Leader Fault at four locations from 2018 to 2020, revealing near surface (< 4m depth) contractional deformation of Holocene stratigraphy. Three of the trench locations uncovered clear evidence for rupture of the Leader Fault prior to 2016, with fault displacement of near surface stratigraphy being greater than displacement recorded during the 2016 earthquake. Radiocarbon dating of in-situ organic material from two trenches indicate a date of the penultimate earthquake on the Leader Fault within the past 1000 years. This date is consistent with The Humps and Leader faults having ruptured simultaneously in the past, and with multi-fault ruptures involving The Humps, Leader, Hundalee and Stone Jug faults having occurred prior to the 2016 Kaikōura earthquake. Overall, the results contribute to an improved understanding of the Kaikōura earthquake and highlight the importance of detailed structural and paleoseismic investigations in determining controls on earthquake ‘complexity’.
Coastal margins are exposed to rising sea levels that present challenging circumstances for natural resource management. This study investigates a rare example of tectonic displacement caused by earthquakes that generated rapid sea-level change in a tidal lagoon system typical of many worldwide. This thesis begins by evaluating the coastal squeeze effects caused by interactions between relative sea-level (RSL) rise and the built environment of Christchurch, New Zealand, and also examples of release from similar effects in areas of uplift where land reclamations were already present. Quantification of area gains and losses demonstrated the importance of natural lagoon expansion into areas of suitable elevation under conditions of RSL rise and showed that they may be necessary to offset coastal squeeze losses experienced elsewhere. Implications of these spatial effects include the need to provide accommodation space for natural ecosystems under RSL rise, yet other land-uses are likely to be present in the areas required. Consequently, the resilience of these environments depends on facilitating transitions between human land-uses either proactively or in response to disaster events. Principles illustrated by co-seismic sea-level change are generally applicable to climate change adaptation due to the similarity of inundation effects. Furthermore, they highlight the potential role of non-climatic factors in determining the overall trajectory of change. Chapter 2 quantifies impacts on riparian wetland ecosystems over an eight year period post- quake. Coastal wetlands were overwhelmed by RSL rise and recovery trajectories were surprisingly slow. Four risk factors were identified from the observed changes: 1) the encroachment of anthropogenic land-uses, 2) connectivity losses between areas of suitable elevation, 3) the disproportionate effect of larger wetland vulnerabilities, and 4) the need to protect new areas to address the future movement of ecosystems. Chapter 3 evaluates the unique context of shoreline management on a barrier sandspit under sea-level rise. A linked scenario approach was used to evaluate changes on the open coast and estuarine shorelines simultaneously and consider combined effects. The results show dune loss from a third of the study area using a sea-level rise scenario of 1 m over 100 years and with continuation of current land-uses. Increased exposure to natural hazards and accompanying demand for seawalls is a likely consequence unless natural alternatives can be progressed. In contrast, an example of managed retreat following earthquake-induced subsidence of the backshore presents a new opportunity to restart saltmarsh accretion processes seaward of coastal defences with the potential to reverse decades of degradation and build sea-level rise resilience. Considering both shorelines simultaneously highlights the existence of pinch-points from opposing forces that result in small land volumes above the tidal range. Societal adaptation is delicately poised between the paradigms of resisting or accommodating nature and challenged by the long perimeter and confined nature of the sandspit feature. The remaining chapters address the potential for salinity effects caused by tidal prism changes with a focus on the conservation of īnanga (Galaxias maculatus), a culturally important fish that supports New Zealand‘s whitebait fishery. Methodologies were developed to test the hypothesis that RSL changes would drive a shift in the distribution of spawning sites with implications for their management. Chapter 4 describes a new practical methodology for quantifying the total productivity and spatiotemporal variability of spawning sites at catchment scale. Chapter 5 describes the novel use of artificial habitats as a detection tools to help overcome field survey limitations in degraded environments where egg mortality can be high. The results showed that RSL changes resulted in major shifts in spawning locations and these were associated with new patterns of vulnerability due to the continuation of pre-disturbance land-uses. Unexpected findings includes an improved understanding of the spatial relationship between salinity and spawning habitat, and identification of an invasive plant species as important spawning habitat, both with practical management implications. To conclude, the design of legal protection mechanisms was evaluated in relation to the observed habitat shifts and with a focus on two new planning initiatives that identified relatively large protected areas (PAs) in the lower river corridors. Although the larger PAs were better able to accommodate the observed habitat shifts inefficiencies were also apparent due to spatial disparities between PA boundaries and the values requiring protection. To reduce unnecessary trade-offs with other land-uses, PAs of sufficient size to cover the observable spatiotemporal variability and coupled with adaptive capacity to address future change may offer a high effectiveness from a network of smaller PAs. The latter may be informed by both monitoring and modelling of future shifts and these are expected to include upstream habitat migration driven by the identified salinity relationships and eustatic sea-level rise. The thesis concludes with a summary of the knowledge gained from this research that can assist the development of a new paradigm of environmental sustainability incorporating conservation and climate change adaptation. Several promising directions for future research identified within this project are also discussed.
This thesis addresses the topic of local bond behaviour in RC structures. The mechanism of bond refers to the composite action between deformed steel reinforcing bars and the surrounding concrete. Bond behaviour is an open research topic with a wide scope, particularly because bond it is such a fundamental concept to structural engineers. However, despite many bond-related research findings having wide applications, the primary contribution of this research is an experimental evaluation of the prominent features of local bond behaviour and the associated implications for the seismic performance of RC structures. The findings presented in this thesis attempt to address some structural engineering recommendations made by the Canterbury Earthquakes Royal Commission following the 2010-2011 Canterbury (New Zealand) earthquake sequence. A chapter of this thesis discusses the structural behaviour of flexure-dominated RC wall structures with an insufficient quantity of longitudinal reinforcement, among other in situ conditions, that causes material damage to predominantly occur at a single crack plane. In this particular case, the extent of concrete damage and bond deterioration adjacent to the crack plane will influence the ductility capacity that is effectively provided by the reinforcing steel. As a consequence of these in situ conditions, some lightly reinforced wall buildings in Christchurch lost their structural integrity due to brittle fracture of the longitudinal reinforcement. With these concerning post-earthquake observations in mind, there is the underlying intention that this thesis presents experimental evidence of bond behaviour that allows structural engineers to re-assess their confidence levels for the ability of lightly reinforced concrete structures to achieve the life-safety seismic performance objective the ultimate limit state. Three chapters of this thesis are devoted to the experimental work that was conducted as the main contribution of this research. Critical details of the experimental design, bond testing method and test programme are reported. The bond stress-slip relationship was studied through 75 bond pull-out tests. In order to measure the maximum local bond strength, all bond tests were carried out on deformed reinforcing bars that did not yield as the embedded bond length was relatively short. Bond test results have been presented in two separate chapters in which 48 monotonic bond tests and 27 cyclic bond tests are presented. Permutations of the experiments include the loading rate, cyclic loading history, concrete strength (25 to 70 MPa), concrete age, cover thickness, bar diameter (16 and 20 mm), embedded length, and position of the embedded bond region within the specimen (close or far away to the free surface). The parametric study showed that the concrete strength significantly influences the maximum bond strength and that it is reasonable to normalise the bond stress by the square-root of the concrete compressive strength, √(f'c). The generalised monotonic bond behaviour is described within. An important outcome of the research is that the measured bond strength and stiffness was higher than stated by the bond stress-slip relationship in the fib Model Code 2010. To account for these observed differences, an alternative model is proposed for the local monotonic bond stress-slip relationship. Cyclic bond tests showed a significant proportion of the total bond degradation occurs after the loading cycle in the peak bond strength range, which is when bond slip has exceeded 0.5 mm. Subsequent loading to constant slip values showed a linear relationship between the amount of bond strength degradation and the log of the number of cycles that were applied. To a greater extent, the cyclic bond deterioration depends on the bond slip range, regardless of whether the applied load cycling is half- or fully-reversed. The observed bond deterioration and hysteretic energy dissipated during cyclic loading was found to agree reasonably well between these cyclic tests with different loading protocols. The cyclic bond deterioration was also found to be reasonably consistent exponential damage models found in the literature. This research concluded that the deformed reinforcing bars used in NZ construction, embedded in moderate to high strength concrete, are able to develop high local bond stresses that are mobilised by a small amount of local bond slip. Although the relative rib geometry was not varied within this experimental programme, a general conclusion of this thesis is that deformed bars currently available in NZ have a relative rib bearing area that is comparatively higher than the test bars used in previous international research. From the parametric study it was found that the maximum monotonic bond strength is significant enhanced by dynamic loading rates. Experimental evidence of high bond strength and initial bond stiffness generally suggests that only a small amount of local bond slip that can occur when the deformed test bar was subjected to large tension forces. Minimal bond slip and bond damage limits the effective yielding length that is available for the reinforcing steel to distribute inelastic material strains. Consequently, the potential for brittle fracture of the reinforcement may be a more problematic and widespread issue than is apparent to structural engineers. This research has provided information that improve the reliability of engineering predictions (with respect to ductility capacity) of maximum crack widths and the extent of bond deterioration that might occur in RC structures during seismic actions.
