At the conclusion of the 2010 and 2011 Canterbury earthquakes more than 5100 homes had been deemed unsafe for habitation. The land and buildings of these were labelled “red zoned” and are too badly damaged for remediation. These homes have been demolished or are destined for demolition. To assist the red zone population to relocate, central government have offered to ‘buy out’ home owners at the Governmental Value (GV) that was last reviewed in 2007. While generous in the economic context at the time, the area affected was the lowest value land and housing in Christchurch and so there is a capital shortfall between the 2007 property value and the cost of relocating to more expensive properties. This shortfall is made worse by increasing present day values since the earthquakes. Red zone residents have had to relocate to the far North and Western extremities of Christchurch, and some chose to move even further to neighbouring towns or cities. The eastern areas and commercial centres close to the red zone are affected as well. They have lost critical mass which has negatively impacted businesses in the catchments of the Red Zone. This thesis aims to repopulate the suburbs most affected by the abandonment of the red zone houses. Because of the relative scarcity of sound building sites in the East and to introduce affordability to these houses, an alternative method of development is required than the existing low density suburban model. Smart medium density design will be tested as an affordable and appropriate means of living. Existing knowledge in this field will be reviewed, an analysis of what East Christchurch’s key characteristics are will occur, and an examination of built works and site investigations will also be conducted. The research finds that at housing densities of 40 units per hectare, the spatial, vehicle, aesthetic needs of East Christchurch can be accommodated. Centralising development is also found to offer better lifestyle choices than the isolated suburbs at the edges of Christchurch, to be more efficient using existing infrastructure, and to place less reliance on cars. Stronger communities are formed from the outset and for a full range of demographics. Eastern affordable housing options are realised and Christchurch’s ever expanding suburban tendencies are addressed. East Christchurch presently displays a gaping scar of devastated houses that ‘The New Eastside’ provides a bandage and a cure for. Displaced and dispossessed Christchurch residents can be re-housed within a new heart for East Christchurch.
A video of a presentation by Professor David Johnston during the fourth plenary of the 2016 People in Disasters Conference. Johnston is a Senior Scientist at GNS Science and Director of the Joint Centre for Disaster Research in the School of Psychology at Massey University. The presentation is titled, "Understanding Immediate Human Behaviour to the 2010-2011 Canterbury Earthquake Sequence, Implications for injury prevention and risk communication".The abstract for the presentation reads as follows: The 2010 and 2011 Canterbury earthquake sequences have given us a unique opportunity to better understand human behaviour during and immediately after an earthquake. On 4 September 2010, a magnitude 7.1 earthquake occurred near Darfield in the Canterbury region of New Zealand. There were no deaths, but several thousand people sustained injuries and sought medical assistance. Less than 6 months later, a magnitude 6.2 earthquake occurred under Christchurch City at 12:51 p.m. on 22 February 2011. A total of 182 people were killed in the first 24 hours and over 7,000 people injured overall. To reduce earthquake casualties in future events, it is important to understand how people behaved during and immediately after the shaking, and how their behaviour exposed them to risk of death or injury. Most previous studies have relied on an analysis of medical records and/or reflective interviews and questionnaire studies. In Canterbury we were able to combine a range of methods to explore earthquake shaking behaviours and the causes of injuries. In New Zealand, the Accident Compensation Corporation (a national health payment scheme run by the government) allowed researchers to access injury data from over 9,500 people from the Darfield (4 September 2010) and Christchurch (22 February 2011 ) earthquakes. The total injury burden was analysed for demography, context of injury, causes of injury, and injury type. From the injury data inferences into human behaviour were derived. We were able to classify the injury context as direct (immediate shaking of the primary earthquake or aftershocks causing unavoidable injuries), and secondary (cause of injury after shaking ceased). A second study examined people's immediate responses to earthquakes in Christchurch New Zealand and compared responses to the 2011 earthquake in Hitachi, Japan. A further study has developed a systematic process and coding scheme to analyse earthquake video footage of human behaviour during strong earthquake shaking. From these studies a number of recommendations for injury prevention and risk communication can be made. In general, improved building codes, strengthening buildings, and securing fittings will reduce future earthquake deaths and injuries. However, the high rate of injuries incurred from undertaking an inappropriate action (e.g. moving around) during or immediately after an earthquake suggests that further education is needed to promote appropriate actions during and after earthquakes. In New Zealand - as in US and worldwide - public education efforts such as the 'Shakeout' exercise are trying to address the behavioural aspects of injury prevention.
<b>New Zealand has experienced several strong earthquakes in its history. While an earthquake cannot be prevented from occurring, planning can reduce its consequences when it does occur. This dissertation research examines various aspects of disaster risk management policy in Aotearoa New Zealand.</b> Chapter 2 develops a method to rank and prioritise high-rise buildings for seismic retrofitting in Wellington, the earthquake-prone capital city of New Zealand. These buildings pose risks to Wellington’s long-term seismic resilience that are of clear concern to current and future policymakers. The prioritization strategy we propose, based on multi-criteria decision analysis (MCDA) methods, considers a variety of data on each building, including not only its structural characteristics, but also its location, its economic value to the city, and its social importance to the community around it. The study demonstrates how different measures, within four general criteria – life safety, geo-spatial location of the building, its economic role, and its socio-cultural role – can be operationalized into a viable framework for determining retrofitting/demolition policy priorities. Chapter 3 and chapter 4 analyse the Residential Red Zone (RRR) program that was implemented in Christchurch after the 2011 earthquake. In the program, approximately 8,000 homeowners were told that their homes were no longer permittable, and they were bought by the government (through the Canterbury Earthquake Recovery Authority). Chapter 3 examines the subjective wellbeing of the RRR residents (around 16000 people) after they were forced to move. We consider three indicators of subjective wellbeing: quality of life, stress, and emotional wellbeing. We found that demographic factors, health conditions, and the type of government compensation the residents accepted, were all significant determinants of the wellbeing of the Red Zone residents. More social relations, better financial circumstances, and the perception of better government communication were also all associated positively with a higher quality of life, less stress, and higher emotional wellbeing. Chapter 4 concentrates on the impact of this managed retreat program on RRR residents’ income. We use individual-level comprehensive, administrative, panel data from Canterbury, and difference in difference evaluation method to explore the effects of displacement on Red Zone residential residents. We found that compared to non-relocated neighbours, the displaced people experience a significant initial decrease in their wages and salaries, and their total income. The impacts vary with time spent in the Red Zone and when they moved away. Wages and salaries of those who were red-zoned and moved in 2011 were reduced by 8%, and 5.4% for those who moved in 2012. Females faced greater decreases in wages and salaries, and total income, than males. There were no discernible impacts of the relocation on people’s self-employment income.
Documenting earthquake-induced ground deformation is significant to assess the characteristics of past and contemporary earthquakes and provide insight into seismic hazard. This study uses airborne light detection and ranging (LiDAR) and conducts multi-disciplinary field techniques to document the surface rupture morphology and evaluate the paleoseismicity and seismic hazard parameters of the Hurunui segment of the Hope Fault in the northern South Island of New Zealand. It also documents and evaluates seismically induced features and ground motion characteristics of the 2010 Darfield and 2011 Christchurch earthquakes in the Port Hills, south of Christchurch. These two studies are linked in that they investigate the near-field coseismic features of large (Mw ~7.1) earthquakes in New Zealand and produce data for evaluating seismic hazards of future earthquakes. In the northern South Island of New Zealand, the Australian-Pacific plate boundary is characterised by strike-slip deformation across the Marlborough Fault System (MFS). The ENE-striking Hope Fault (length: ~230 km) is the youngest and southernmost fault in the MFS, and the second fastest slipping fault in New Zealand. The Hope Fault is a major source of seismic hazard in New Zealand and has ruptured (in-part) historically in the Mw 7.1 1888 Amuri earthquake. In the west, the Hurunui segment of the Hope Fault is covered by beech forest. Hence, its seismic hazard parameters and paleoearthquake chronology were poorly constrained and it was unknown whether the 1888 earthquake ruptured this segment or not and if so, to what extent. Utilising LiDAR and field data, a 29 km-long section of the Hurunui segment of the Hope Fault is mapped. LiDAR-mapping clearly reveals the principal slip zone (PSZ) of the fault and a suite of previously unrecognised structures that form the fault deformation zone (FDZ). FDZ width measurements from 415 locations reveal a spatially-variable, active FDZ up to ~500 m wide with an average width of 200 m. Kinematic analysis of the fault structures shows that the Hurunui segment strikes between 070° and 075° and is optimally oriented for dextral strike-slip within the regional stress field. This implies that the wide FDZ observed is unlikely to result from large-scale fault mis-orientation with respect to regional stresses. The analysis of FDZ width indicates that it increases with increased hanging wall topography and increased topographic relief suggesting that along-strike topographic perturbations to fault geometry and stress states increase fault zone complexity and width. FDZ width also increases where the tips of adjacent PSZ strands locally vary in strike, and where the thickness of alluvial deposits overlying bedrock increases. LiDAR- and photogrammetrically-derived topographic mapping indicates that the boundary between the Hurunui and Hope River segments is characterised by a ~850-m-wide right stepover and a 9º-14° fault bend. Paleoseismic trenching at Hope Shelter site reveals that 6 earthquakes occurred at A.D. 1888, 1740-1840, 1479-1623, 819-1092, 439-551, and 373- 419. These rupture events have a mean recurrence interval of ~298 ± 88 yr and inter-event times ranging from 98 to 595 yrs. The variation in the inter-event times is explained by (1) coalescing rupture overlap from the adjacent Hope River segment on to the Hurunui segment at the study site, (2) temporal clustering of large earthquakes on the Hurunui segment, and/or (3) ‘missing’ rupture events. It appears that the first two options are more plausible to explain the earthquake chronologies and rupture behaviour on the Hurunui segment, given the detailed nature of the geologic and chronologic investigations. This study provides first evidence for coseismic multi-segment ruptures on the Hope Fault by identifying a rupture length of 44-70 km for the 1888 earthquake, which was not confined to the Hope River segment (primary source for the 1888 earthquake). LiDAR data is also used to identify and measure dextral displacements and scarp heights from the PSZ and structures within the FDZ along the Hurunui segment. Reconstruction of large dextrally-offset geomorphic features shows that the vertical component of slip accounts for only ~1% of the horizontal displacements and confirms that the fault is predominantly strike-slip. A strong correlation exists between the dextral displacements and elevations of geomorphic features suggesting the possibility of age correlation between the geomorphic features. A mean single event displacement (SED) of 3.6 ± 0.7 m is determined from interpretation of sets of dextral displacements of ≤ 25 m. Using the available surface age data and the cumulative dextral displacements from Matagouri Flat, McKenzie Fan, Macs Knob and Hope River sites, and the mean SED, a mean slip rate of 12.2 ± 2.4 mm/yr, and a mean recurrence interval of ~320 ± 120 yr, and a potential earthquake magnitude of Mw 7.2 are determined for the Hurunui segment. This study suggests that the fault slip rate has been constant over the last ~15000 yr. Strong ground motions from the 2010 Darfield (Canterbury) earthquake displaced boulders and caused ground damage on some ridge crests in the Port Hills. However, the 2011 Christchurch earthquake neither displaced boulders nor caused ground damage at the same ridge crests. Documentation of locations (~400 m a.s.l.), lateral displacements (8-970 cm), displacement direction (250° ± 20°) of displaced boulders, in addition to their hosting socket geometries (< 1 cm to 50 cm depth), the orientation of the ridges (000°-015°) indicate that boulders have been displaced in the direction of instrumentally recorded transient peak ground horizontal displacements nearby and that the seismic waves have been amplified at the study sites. The co-existence of displaced and non-displaced boulders at proximal sites suggests small-scale ground motion variability and/or varying boulder-ground dynamic interactions relating to shallow phenomena such as variability in soil depth, bedrock fracture density and/or microtopography on the bedrock-soil interface. Shorter shaking duration of the 2011 Christchurch event, differing frequency contents and different source characteristics were all factors that may have contributed to generating circumstances less favourable to boulder displacement in this earthquake. Investigating seismically induced features, fault behaviour, site effects on the rupture behaviour, and site response to the seismic waves provides insights into fault rupture hazards.
Knowledge of past climate variability is essential for understanding present and future climate trends. This study used Halocarpus biformis (pink pine) ring-width chronologies to investigate palaeotemperature history in Westland, New Zealand. The ensuing reconstruction is among the longest palaeoseries produced for New Zealand to date. It is in good agreement with other tree-ring-based records, and with instrumental (both local and hemispheric) data. Thirteen pink pine chronologies were developed. Ring-width measurements were detrended using the Regional Curve Standardisation method to retain as much low-frequency variance as possible. Crossdating revealed the existence of a strong common signal among trees. Inter-site comparison indicated that a common control mechanism affected tree growth not only within sites, but also across sites. To determine whether climate was the main factor that controlled the growth of pink pine in Westland, correlation and response function analyses were employed. Temperature, precipitation and the Southern Oscillation Index were tested for their relationship with tree growth. Mean monthly temperature was identified as the primary growth-limiting factor. Chronologies were positively correlated with temperature over an extended period (5-17 months), and climate response modelling showed that temperature explained 11-60% variance in the tree-ring data. The highest and most stable correlations occurred between tree growth and summer (January-March) temperatures. Tree-ring data from the six sites that contained the strongest temperature signal were combined, and the Westland Regional Chronology (WRC) was developed. The WRC was then used to reconstruct January-March temperatures back to A.D. 1480. The calibration model explained 43% of the variance in temperature, and all calibration and verification tests were passed at high levels of significance. The reconstruction showed that temperatures in Westland have been following a positive trend over the last 520 years. The coolest 25-year period was 1542-1566, while temperatures reached their maximum in 1966-1990. Spectral analysis of the Westland palaeotemperature record revealed cycles at periods of about 3, 5-6, 11, 14, 22, 45 and 125 years. This study also confirmed that climate response is species-dependent. A separate exercise, which compared two species from the same site, demonstrated that while pink pine's growth was mainly influenced by summer temperatures, Libocedrus bidwillii was affected by conditions at the beginning of the growing season. However, the temperature signal in Westland's Libocedrus bidwillii was insufficient to produce a reliable reconstruction. It might be because the climate signal in this species was obscured by disturbances, as was shown in the final section of this project. Frequent growth releases and suppressions implied that Libocedrus bidwillii integrated both major (Alpine Fault earthquakes) and minor (windthrow) disturbances in its ring widths. Pink pine, on the other hand, was not sensitive to disturbance, and was therefore a better indicator of palaeotemperatures in Westland. This research has strengthened the New Zealand network of chronology sites, and confirmed that pink pine has great dendroclimatic value. The last 520 years of temperature fluctuations were reconstructed with a high degree of fidelity - the model developed in this thesis is currently the most accurate estimate of a temperature-growth relationship in the country.
The Avon and Heathcote Rivers, located in the city of Christchurch, New Zealand, are lowland spring-fed rivers linked with the Christchurch Groundwater System. At present, the flow paths and recharge sources to the Christchurch Groundwater System are not fully understood. Study of both the Avon and Heathcote Rivers can provide greater insight into this system. In addition, during the period 2010-2012, Christchurch has experienced large amounts of seismic activity, including a devastating Mw 6.2 aftershock on February 22nd, 2011, which caused widespread damage and loss of life. Associated with these earthquakes was the release of large amounts of water through liquefaction and temporary springs throughout the city. This provided a unique opportunity to study groundwater surface water interactions following a large scale seismic event. Presented herein is the first major geochemical study on the Avon and Heathcote Rivers and the hydrological impact of the February 22, 2011 Christchurch Earthquake. The Avon, Heathcote, and Waimakariri Rivers were sampled in quarterly periods starting in July 2011 and analyzed for stable Isotopes δ¹⁸O, δD, and δ¹³C and major anion composition. In addition, post -earthquake samples were collected over the days immediately following the February 22, 2011 earthquake and analyzed for stable isotopes δ¹⁸O and δD and major anion composition. A variety of analytical methods were used identify the source of the waters in the Avon-Heathcote System and evaluate the effectiveness of stable isotopes as geochemical tracers in the Christchurch Groundwater System. The results of this thesis found that the waters from the Avon and Heathcote Rivers are geochemically the same, originating from groundwater, and exhibit a strong tidal influence within 5km of the Avon-Heathcote Estuary. The surface waters released following the February 22nd, 2011 earthquake were indistinguishable from quarterly samples taken from the Avon and Heathcote Rivers when comparing stable isotopic composition. The anion data suggests the waters released following the February 22nd, 2011 Christchurch Earthquake were sourced primarily from shallow groundwater, and also suggests a presence of urban sewage at some sites. Attempts to estimate recharge sources for the Avon-Heathcote Rivers using published models for the Christchurch Groundwater System yielded results that were not consistent between models. In evaluating the use of geochemical constituents as tracers in the Christchurch Groundwater System, no one isotope could provide a clear resolution, but when used in conjunction, δ¹⁸O, δ¹³C, and DIC, seem to be the most effective tracers. Sample sizes for δ¹³C were too small for a robust evaluation. Variability on the Waimakariri River appears to be greater than previously estimated, which could have significant impacts on geochemical models for the Christchurch Groundwater System. This research demonstrates the value of using multiple geochemical constituents to enrich our understanding of the groundwater surfaces-water interactions and the Christchurch Groundwater System as a whole.
During 2010 and 2011, major earthquakes caused widespread damage and the deaths of 185 people in the city of Christchurch. Damaged school buildings resulted in state intervention which required amendment of the Education Act of 1989, and the development of ‘site sharing agreements’ in undamaged schools to cater for the needs of students whose schools had closed. An effective plan was also developed for student assessment through establishing an earthquake impaired derived grade process. Previous research into traditional explanations of educational inequalities in the United Kingdom, the United States of America, and New Zealand were reviewed through various processes within three educational inputs: the student, the school and the state. Research into the impacts of urban natural disasters on education and education inequalities found literature on post disaster education systems but nothing could be found that included performance data. The impacts of the Canterbury earthquakes on educational inequalities and achievement were analysed over 2009-2012. The baseline year was 2009, the year before the first earthquake, while 2012 is seen as the recovery year as no schools closed due to seismic events and there was no state intervention into the education of the region. National Certificate of Educational Achievement (NCEA) results levels 1-3 from thirty-four secondary schools in the greater Christchurch region were graphed and analysed. Regression analysis indicates; in 2009, educational inequalities existed with a strong positive relationship between a school’s decile rating and NCEA achievement. When schools were grouped into decile rankings (1-10) and their 2010 NCEA levels 1-3 results were compared with the previous year, the percentage of change indicates an overall lower NCEA achievement in 2010 across all deciles, but particularly in lower decile schools. By contrast, when 2011 NCEA results were compared with those of 2009, as a percentage of change, lower decile schools fared better. Non site sharing schools also achieved higher results than site sharing schools. State interventions, had however contributed towards student’s achieving national examinations and entry to university in 2011. When NCEA results for 2012 were compared to 2009 educational inequalities still exist, however in 2012 the positive relationship between decile rating and achievement is marginally weaker than in 2009. Human ethics approval was required to survey one Christchurch secondary school community of students (aged between 12 and 18), teachers and staff, parents and caregivers during October 2011. Participation was voluntary and without incentives, 154 completed questionnaires were received. The Canterbury earthquakes and aftershocks changed the lives of the research participants. This school community was displaced to another school due to the Christchurch earthquake on 22 February 2011. Research results are grouped under four geographical perspectives; spatial impacts, socio-economic impacts, displacement, and health and wellbeing. Further research possibilities include researching the lag effects from the Canterbury earthquakes on school age children.
Capacity design and hierarchy of strength philosophies at the base of modern seismic codes allow inelastic response in case of severe earthquakes and thus, in most traditional systems, damage develops at well-defined locations of reinforced concrete (RC) structures, known as plastic hinges. The 2010 and 2011 Christchurch earthquakes have demonstrated that this philosophy worked as expected. Plastic hinges formed in beams, in coupling beams and at the base of columns and walls. Structures were damaged permanently, but did not collapse. The 2010 and 2011 Christchurch earthquakes also highlighted a critical issue: the reparability of damaged buildings. No methodologies or techniques were available to estimate the level of subsequent earthquakes that RC buildings could still sustain before collapse. No repair techniques capable of restoring the initial condition of buildings were known. Finally, the cost-effectiveness of an eventual repair intervention, when compared with a new building, was unknown. These aspects, added to nuances of New Zealand building owners’ insurance coverage, encouraged the demolition of many buildings. Moreover, there was a perceived strong demand from government and industry to develop techniques for assessing damage to steel reinforcement bars embedded in cracked structural concrete elements. The most common questions were: “Have the steel bars been damaged in correspondence to the concrete cracks?”, “How much plastic deformation have the steel bars undergone?”, and “What is the residual strain capacity of the damaged bars?” Minimally invasive techniques capable of quantifying the level and extent of plastic deformation and residual strain capacity are not yet available. Although some studies had been recently conducted, a validated method is yet to be widely accepted. In this thesis, a least-invasive method for the damage-assessment of steel reinforcement is developed. Based on the information obtained from hardness testing and a single tensile test, it is possible to estimate the mechanical properties of earthquake-damaged rebars. The reduction in the low-cycle fatigue life due to strain ageing is also quantified. The proposed damage assessment methodology is based on empirical relationships between hardness and strain and residual strain capacity. If damage is suspected from in situ measurements, visual inspection or computer analysis, a bar may be removed and more accurate hardness measurements can be obtained using the lab-based Vickers hardness methodology. The Vickers hardness profile of damaged bars is then compared with calibration curves (Vickers hardness versus strain and residual strain capacity) previously developed for similar steel reinforcement bars extracted from undamaged locations. Experimental tests demonstrated that the time- and temperature-dependent strain-ageing phenomenon causes changes in the mechanical properties of plastically deformed steels. In particular, yield strength and hardness increases, whereas ductility decreases. The changes in mechanical properties are quantified and their implications on the hardness method are highlighted. Low-cycle fatigue (LCF) failures of steel reinforcing bars have been observed in laboratory testing and post-earthquake damage inspections. Often, failure might not occur during a first seismic event. However, damage is accumulated and the remaining fatigue life is reduced. Failure might therefore occur in a subsequent seismic event. Although numerous studies exist on the LCF behaviour of steel rebars, no studies had been conducted on the strain-ageing effects on the remaining fatigue life. In this thesis, the reduction in fatigue life due to this phenomenon is determined through a number of experimental tests.
The lateral capacity of a conventional CLT shear wall is often governed by the strength and stiffness of its connections, which do not significantly utilize the in-plane strength of the CLT. Therefore, CLT shear walls are not yet being used efficiently in the construction of mass timber buildings due to a lack of research on high-capacity connections and alternative wall configurations. In this study, cyclic experiments were completed on six full-scale, 5-ply cantilever CLT shear walls with high-capacity hold-downs using mixed angle screws and bolts. All specimens exhibited significantly higher strength and stiffness than previously tested conventional CLT shear walls in the literature. The base connections demonstrated ductile failure modes through yielding of the hold-down connections. Based on the experimental results, numerical models were calibrated to investigate the seismic behaviour of CLT shear walls for prototype buildings of 3 and 6-storeys in Christchurch, NZ. As an alternative to cantilever (single) shear walls, a type of coupled wall with steel link beams between adjacent CLT wall piers was investigated. Effective coupling requires the link beam-to-wall connections to have adequate strength to ensure ductile link beam responses and adequate stiffness to yield the link beams at a relatively low inter-storey drift level. To this end, three beam-to-wall connection types were developed and cyclically tested to investigate their behaviour and feasibility. Based on the test results of the critical connection, a 3-storey, 2/3-scale coupled CLT wall specimen with three steel link beams and mixed angle screwed hold-downs was cyclically tested to evaluate its performance and experimentally validate the system concept. The test results showed a relatively high lateral strength compared to conventional CLT shear walls, as well as a high system ductility ratio of 7.6. Failure of the system was characterised by combined bending and withdrawal of the screws in the mixed angle screw hold-downs, yielding and eventual inelastic buckling of the steel link beams, CLT toe crushing, and local CLT delamination. Following the initial test, the steel link beams, mixed angle screw hold-downs, and damaged CLT regions were repaired, then the wall specimen was re-tested. The repaired wall behaved similarly to the original test and exhibited slightly higher energy dissipation and peak strength, but marginally more rapid strength deterioration under cyclic loading. Several hybrid coupled CLT shear walls were numerically modelled and calibrated based on the results of the coupled wall experiments. Pushover analyses were conducted on a series of configurations to validate a capacity design method for the system and to investigate reasonable parameter values for use in the preliminary design of the system. Additionally, an iterative seismic design method was proposed and used to design sample buildings of 6, 8, and 10-storeys using both nonlinear pushover and nonlinear time history analyses to verify the prototype designs. Results of the sample building analyses demonstrated adequate seismic behaviour and the proposed design parameters were found to be appropriate. In summary, high-capacity CLT shear walls can be used for the resistance of earthquakes by using stronger base connections and coupled wall configurations. The large-scale experimental testing in this study has demonstrated that both cantilever and coupled CLT shear walls are feasible LLRSs which can provide significantly greater lateral strength, stiffness, and energy dissipation than conventional CLT shear wall configurations.
This poster provides a comparison between the strong ground motions observed in the 22 February 2011 Mw6.3 Christchurch earthquake with those observed in Tokyo during the 11 March 2011 Mw9.0 Tohoku earthquake. The destuction resulting from both of these events has been well documented, although tsunami was the principal cause of damage in the latter event, and less attention has been devoted to the impact of earthquake-induced ground motions. Despite Tokyo being located over 100km from the nearest part of the causative rupture, the ground motions observed from the Tohoku earthquake were significant enough to cause structural damage and also significant liquefaction to loose reclaimed soils in Tokyo Bay. The author was fortunate enough (from the perspective of an earthquake engineer) to experience first-hand both of these events. Following the Tohoku event, the athor conducted various ground motion analyses and reconniassance of the Urayasu region in Tokyo Bay affected by liquefaction in collaboration with Prof. Kenji Ishihara. This conference is therefore a fitting opportunity in which to discuss some of authors insights obtained as a result of this first hand knowledge. Figure 1 illustrates the ground motions recorded in the Christchurch CBD in the 22 February 2011 and 4 September 2010 earthquakes, with that recorded in Tokyo Bay in the 11 March 2011 Tohoku earthquake. It is evident that these three ground motions vary widely in their amplitude and duration. The CBGS ground motion from the 22 February 2011 event has a very large amplitude (nearly 0.6g) and short duration (approx. 10s of intense shaking), as a result of the causal Mw6.3 rupture at short distance (Rrup=4km). The CBGS ground motion from the 4 September 2010 earthquake has a longer duration (approx. 30s of intense shaking), but reduced acceleration amplitude, as a result of the causal Mw7.1 rupture at a short-to-moderate distance (Rrup=14km). Finally, the Urayasu ground motion in Tokyo bay during the 11 March 2011 Tohoku earthquake exhibits an acceleration amplitude similar to the 4 September 2010 CBGS ground motion, but a significantly larger duration (approx 150s of intense shaking). Clearly, these three different ground motions will affect structures and soils in different ways depending on the vibration characteristics of the structures/soil, and the potential for strength and stiffness degradation due to cumulative effects. Figure 2 provides a comparison between the arias intensities of the several ground motion records from the three different events. It can be seen that the arias intensities of the ground motions in the Christchurch CBD from the 22 February 2011 earthquake (which is on average AI=2.5m/s) is approximately twice that from the 4 September 2010 earthquake (average AI≈1.25). This is consistent with a factor of approximately 1.6 obtained by Cubrinovski et al. (2011) using the stress-based (i.e.PGA-MSF) approach of liquefaction triggering. It can also be seen that the arias intensity of the ground motions recorded in Tokyo during the 2011 Tohoku earthquake are larger than ground motions in the Christchurch CBD from the 4 September 2011 earthquake, but smaller than those of the 22 February 2011 earthquake. Based on the arias intensity liquefaction triggering approach it can therefore be concluded that the ground motion severity, in terms of liquefaction potential, for the Tokyo ground motions is between those ground motions in Christchurch CBD from the 4 September 2010 and 22 February 2011 events.
Spatial variations in river facies exerted a strong influence on the distribution of liquefaction features observed in Christchurch during the 2010-11 Canterbury Earthquake Sequence (CES). Liquefaction and liquefaction-induced ground deformation was primarily concentrated near modern waterways and areas underlain by Holocene fluvial deposits with shallow water tables (< 1 to 2 m). In southern Christchurch, spatial variations of liquefaction and subsidence were documented in the suburbs within inner meander loops of the Heathcote River. Newly acquired geospatial data, geotechnical reports and eye-witness discussions are compiled to provide a detailed account of the surficial effects of CES liquefaction and ground deformation adjacent to the Heathcote River. LiDAR data and aerial photography are used to produce a new series of original figures which reveal the locations of recurrent liquefaction and subsidence. To investigate why variable liquefaction patterns occurred, the distribution of surface ejecta and associated ground damage is compared with near-surface sedimentologic, topographic, and geomorphic variability to seek relationships between the near-surface properties and observed ground damages. The most severe liquefaction was concentrated within a topographic low in the suburb of St Martins, an inner meander loop of the Heathcote River, with liquefaction only minor or absent in the surrounding areas. Subsurface investigations at two sites in St Martins enable documentation of fluvial stratigraphy, the expressions of liquefaction, and identification of pre-CES liquefaction features. Excavation to water table depths (~1.5 m below the surface) across sand boils reveals multiple generations of CES liquefaction dikes and sills that cross-cut Holocene fluvial and anthropogenic stratigraphy. Based on in situ geotechnical tests (CPT) indicating sediment with a factor of safety < 1, the majority of surface ejecta was sourced from well-sorted fine to medium sand at < 5 m depth, with the most damaging liquefaction corresponding with the location of a low-lying sandy paleochannel, a remnant river channel from the Holocene migration of the meander in St Martins. In the adjacent suburb of Beckenham, where migration of the Heathcote River has been laterally confined by topography associated with the volcanic lithologies of Banks Peninsula, severe liquefaction was absent with only minor sand boils occurring closest to the modern river channel. Auger sampling across the suburb revealed thick (>1 m) clay-rich overbank and back swamp sediments that produced a stratigraphy which likely confined the units susceptible to liquefaction and prevented widespread ejection of liquefied material. This analysis suggests river migration promotes the formation and preservation of fluvial deposits prone to liquefaction. Trenching revealed the strongest CES earthquakes with large vertical accelerations favoured sill formation and severe subsidence at highly susceptible locations corresponding with an abandoned channel. Less vulnerable sites containing deeper and thinner sand bodies only liquefied in the strongest and most proximal earthquakes forming minor localised liquefaction features. Liquefaction was less prominent and severe subsidence was absent where lateral confinement of a Heathcote meander has promoted the formation of fluvial stratum resistant to liquefaction. Correlating CES liquefaction with geomorphic interpretations of Christchurch’s Heathcote River highlights methods in which the performance of liquefaction susceptibility models can be improved. These include developing a reliable proxy for estimating soil conditions in meandering fluvial systems by interpreting the geology and geomorphology, derived from LiDAR data and modern river morphology, to improve the methods of accounting for the susceptibility of an area. Combining geomorphic interpretations with geotechnical data can be applied elsewhere to identify regional liquefaction susceptibilities, improve existing liquefaction susceptibility datasets, and predict future earthquake damage.
Sewerage systems convey sewage, or wastewater, from residential or commercial buildings through complex reticulation networks to treatment plants. During seismic events both transient ground motion and permanent ground deformation can induce physical damage to sewerage system components, limiting or impeding the operability of the whole system. The malfunction of municipal sewerage systems can result in the pollution of nearby waterways through discharge of untreated sewage, pose a public health threat by preventing the use of appropriate sanitation facilities, and cause serious inconvenience for rescuers and residents. Christchurch, the second largest city in New Zealand, was seriously affected by the Canterbury Earthquake Sequence (CES) in 2010-2011. The CES imposed widespread damage to the Christchurch sewerage system (CSS), causing a significant loss of functionality and serviceability to the system. The Christchurch City Council (CCC) relied heavily on temporary sewerage services for several months following the CES. The temporary services were supported by use of chemical and portable toilets to supplement the damaged wastewater system. The rebuild delivery agency -Stronger Christchurch Infrastructure Rebuild Team (SCIRT) was created to be responsible for repair of 85 % of the damaged horizontal infrastructure (i.e., water, wastewater, stormwater systems, and roads) in Christchurch. Numerous initiatives to create platforms/tools aiming to, on the one hand, support the understanding, management and mitigation of seismic risk for infrastructure prior to disasters, and on the other hand, to support the decision-making for post-disaster reconstruction and recovery, have been promoted worldwide. Despite this, the CES in New Zealand highlighted that none of the existing platforms/tools are either accessible and/or readable or usable by emergency managers and decision makers for restoring the CSS. Furthermore, the majority of existing tools have a sole focus on the engineering perspective, while the holistic process of formulating recovery decisions is based on system-wide approach, where a variety of factors in addition to technical considerations are involved. Lastly, there is a paucity of studies focused on the tools and frameworks for supporting decision-making specifically on sewerage system restoration after earthquakes. This thesis develops a decision support framework for sewerage pipe and system restoration after earthquakes, building on the experience and learning of the organisations involved in recovering the CSS following the CES in 2010-2011. The proposed decision support framework includes three modules: 1) Physical Damage Module (PDM); 2) Functional Impact Module (FIM); 3) Pipeline Restoration Module (PRM). The PDM provides seismic fragility matrices and functions for sewer gravity and pressure pipelines for predicting earthquake-induced physical damage, categorised by pipe materials and liquefaction zones. The FIM demonstrates a set of performance indicators that are categorised in five domains: structural, hydraulic, environmental, social and economic domains. These performance indicators are used to assess loss of wastewater system service and the induced functional impacts in three different phases: emergency response, short-term recovery and long-term restoration. Based on the knowledge of the physical and functional status-quo of the sewerage systems post-earthquake captured through the PDM and FIM, the PRM estimates restoration time of sewer networks by use of restoration models developed using a Random Forest technique and graphically represented in terms of restoration curves. The development of a decision support framework for sewer recovery after earthquakes enables decision makers to assess physical damage, evaluate functional impacts relating to hydraulic, environmental, structural, economic and social contexts, and to predict restoration time of sewerage systems. Furthermore, the decision support framework can be potentially employed to underpin system maintenance and upgrade by guiding system rehabilitation and to monitor system behaviours during business-as-usual time. In conjunction with expert judgement and best practices, this framework can be moreover applied to assist asset managers in targeting the inclusion of system resilience as part of asset maintenance programmes.
Christchurch City Council (Council) is undertaking the Land Drainage Recovery Programme in order to assess the effects of the earthquakes on flood risk to Christchurch. In the course of these investigations it has become better understood that floodplain management should be considered in a multi natural hazards context. Council have therefore engaged the Jacobs, Beca, University of Canterbury, and HR Wallingford project team to investigate the multihazards in eastern areas of Christchurch and develop flood management options which also consider other natural hazards in that context (i.e. how other hazards contribute to flooding both through temporal and spatial coincidence). The study has three stages: Stage 1 Gap Analysis – assessment of information known, identification of gaps and studies required to fill the gaps. Stage 2 Hazard Studies – a gap filling stage with the studies identified in Stage 1. Stage 3 Collating, Optioneering and Reporting – development of options to manage flood risk. This present report is to document findings of Stage 1 and recommends the studies that should be completed for Stage 2. It has also been important to consider how Stage 3 would be delivered and the gaps are prioritised to provide for this. The level of information available and hazards to consider is extensive; requiring this report to be made up of five parts each identifying individual gaps. A process of identifying information for individual hazards in Christchurch has been undertaken and documented (Part 1) followed by assessing the spatial co-location (Part 2) and probabilistic presence of multi hazards using available information. Part 3 considers multi hazard presence both as a temporal coincidence (e.g. an earthquake and flood occurring at one time) and as a cascade sequence (e.g. earthquake followed by a flood at some point in the future). Council have already undertaken a number of options studies for managing flood risk and these are documented in Part 4. Finally Part 5 provides the Gap Analysis Summary and Recommendations to Council. The key findings of Stage 1 gap analysis are: - The spatial analysis showed eastern Christchurch has a large number of hazards present with only 20% of the study area not being affected by any of the hazards mapped. Over 20% of the study area is exposed to four or more hazards at the frequencies and data available. - The majority of the Residential Red Zone is strongly exposed to multiple hazards, with 86% of the area being exposed to 4 or more hazards, and 24% being exposed to 6 or more hazards. - A wide number of gaps are present; however, prioritisation needs to consider the level of benefit and risks associated with not undertaking the studies. In light of this 10 studies ranging in scale are recommended to be done for the project team to complete the present scope of Stage 3. - Stage 3 will need to consider a number of engineering options to address hazards and compare with policy options; however, Council have not established a consistent policy on managed retreat that can be applied for equal comparison; without which substantial assumptions are required. We recommend Council undertake a study to define a managed retreat framework as an option for the city. - In undertaking Stage 1 with floodplain management as the focal point in a multi hazards context we have identified that Stage 3 requires consideration of options in the context of economics, implementation and residual risk. Presently the scope of work will provide a level of definition for floodplain options; however, this will not be at equal levels of detail for other hazard management options. Therefore, we recommend Council considers undertaking other studies with those key hazards (e.g. Coastal Hazards) as a focal point and identifies the engineering options to address such hazards. Doing so will provide equal levels of information for Council to make an informed and defendable decision on which options are progressed following Stage 3.
This report presents the simplified seismic assessment of a case study reinforced concrete (RC) building following the newly developed and refined NZSEE/MBIE guidelines on seismic assessment (NZSEE/MBIE, semi-final draft 26 October 2016). After an overview of the step-by-step ‘diagnostic’ process, including an holistic and qualitative description of the expected vulnerabilities and of the assessment strategy/methodology, focus is given, whilst not limited, to the implementation of a Detailed Seismic Assessment (DSA) (NZSEE/MBIE, 2016c). The DSA is intended to provide a more reliable and consistent outcome than what can be provided by an initial seismic assessment (ISA). In fact, while the Initial Seismic Assessment (ISA), of which the Initial Evaluation Procedure is only a part of, is the more natural and still recommended first step in the overall assessment process, it is mostly intended to be a coarse evaluation involving as few resources as reasonably possible. It is thus expected that an ISA will be followed by a Detailed Seismic Assessment (DSA) not only where the threshold of 33%NBS is not achieved but also where important decisions are intended that are reliant on the seismic status of the building. The use of %NBS (% New Building Standard) as a capacity/demand ratio to describe the result of the seismic assessment at all levels of assessment procedure (ISA through to DSA) is deliberate by the NZSEE/MBIE guidelines (Part A) (NZSEE/MBIE 2016a). The rating for the building needs only be based on the lowest level of assessment that is warranted for the particular circumstances. Discussion on how the %NBS rating is to be determined can be found in Section A3.3 (NZSEE/MBIE 2016a), and, more specifically, in Part B for the ISA (NZSEE/MBIE 2016b) and Part C for the DSA (NZSEE/MBIE 2016c). As per other international approaches, the DSA can be based on several analysis procedures to assess the structural behaviour (linear, nonlinear, static or dynamic, force or displacement-based). The significantly revamped NZSEE 2016 Seismic Assessment Guidelines strongly recommend the use of an analytical (basically ‘by hand’) method, referred to the Simple Lateral Mechanism Analysis (SLaMA) as a first phase of any other numerically-based analysis method. Significant effort has thus been dedicated to provide within the NZSEE 2016 guidelines (NZSEE/MBIE 2016c) a step-by-step description of the procedure, either in general terms (Chapter 2) or with specific reference to Reinforced Concrete Buildings (Chapter 5). More specifically, extract from the guidelines, NZSEE “recommend using the Simple Lateral Mechanism Analysis (SLaMA) procedure as a first step in any assessment. While SLaMA is essentially an analysis technique, it enables assessors to investigate (and present in a simple form) the potential contribution and interaction of a number of structural elements and their likely effect on the building’s global capacity. In some cases, the results of a SLaMA will only be indicative. However, it is expected that its use should help assessors achieve a more reliable outcome than if they only carried out a detailed analysis, especially if that analysis is limited to the elastic range For complex structural systems, a 3D dynamic analysis may be necessary to supplement the simplified nonlinear Simple Lateral Mechanism Analysis (SLaMA).” This report presents the development of a full design example for the the implementation of the SLaMA method on a case study buildings and a validation/comparison with a non-linear static (pushover) analysis. The step-by-step-procedure, summarized in Figure 1, will be herein demonstrated from a component level (beams, columns, wall elements) to a subassembly level (hierarchy of strength in a beam-column joint) and to a system level (frame, C-Wall) assuming initially a 2D behaviour of the key structural system, and then incorporating a by-hand 3D behaviour (torsional effects).
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.
This dissertation addresses a diverse range of topics in the area of physics-based ground motion simulation with particular focus on the Canterbury, New Zealand region. The objectives achieved provide the means to perform hybrid broadband ground motion simulation and subsequently validates the simulation methodology employed. In particu- lar, the following topics are addressed: the development of a 3D seismic velocity model of the Canterbury region for broadband ground motion simulation; the development of a 3D geologic model of the interbedded Quaternary formations to provide insight on observed ground motions; and the investigation of systematic effects through ground motion sim- ulation of small-to-moderate magnitude earthquakes. The paragraphs below outline each contribution in more detail. As a means to perform hybrid broadband ground motion simulation, a 3D model of the geologic structure and associated seismic velocities in the Canterbury region is devel- oped utilising data from depth-converted seismic reflection lines, petroleum and water well logs, cone penetration tests, and implicitly guided by existing contour maps and geologic cross sections in data sparse subregions. The model explicitly characterises five significant and regionally recognisable geologic surfaces that mark the boundaries between geologic units with distinct lithology and age, including the Banks Peninsula volcanics, which are noted to strongly influence seismic wave propagation. The Basement surface represents the base of the Canterbury sedimentary basin, where a large impedance contrast exists re- sulting in basin-generated waves. Seismic velocities for the lithological units between the geologic surfaces are derived from well logs, seismic reflection surveys, root mean square stacking velocities, empirical correlations, and benchmarked against a regional crustal model, thus providing the necessary information for a Canterbury velocity model for use in broadband seismic wave propagation. A 3D high-resolution model of the Quaternary geologic stratigraphic sequence in the Canterbury region is also developed utilising datasets of 527 high-quality water well logs, and 377 near-surface cone penetration test records. The model, developed using geostatistical Kriging, represents the complex interbedded regional Quaternary geology by characterising the boundaries between significant interbedded geologic formations as 3D surfaces including explicit modelling of the formation unconformities resulting from the Banks Peninsula volcanics. The stratigraphic layering present can result in complex wave propagation. The most prevalent trend observed in the surfaces was the downward dip from inland to the eastern coastline as a result of the dominant fluvial depositional environment of the terrestrial gravel formations. The developed model provides a benefi- cial contribution towards developing a comprehensive understanding of recorded ground motions in the region and also providing the necessary information for future site char- acterisation and site response analyses. To highlight the practicality of the model, an example illustrating the role of the model in constraining surface wave analysis-based shear wave velocity profiling is illustrated along with the calculation of transfer functions to quantify the effect of the interbedded geology on wave propagation. Lastly, an investigation of systematic biases in the (Graves and Pitarka, 2010, 2015) ground motion simulation methodology and the specific inputs used for the Canterbury region is presented considering 144 small-to-moderate magnitude earthquakes. In the simulation of these earthquakes, the 3D Canterbury Velocity Model, developed as a part of this dissertation, is used for the low-frequency simulation, and a regional 1D velocity model for the high-frequency simulation. Representative results for individual earthquake sources are first presented to highlight the characteristics of the small-to-moderate mag- nitude earthquake simulations through waveforms, intensity measure scaling with source- to-site distance, and spectral bias of the individual events. Subsequently, a residual de- composition is performed to examine the between- and within-event residuals between observed data, and simulated and empirical predictions. By decomposing the residuals into between- and within-event residuals, the biases in source, path and site effects, and their causes, can be inferred. The residuals are comprehensively examined considering their aggregated characteristics, dependence on predictor variables, spatial distribution, and site-specific effects. The results of the simulation are also benchmarked against empir- ical ground motion models, where their similarities manifest from common components in their prediction. Ultimately, suggestions to improve the predictive capability of the simulations are presented as a result of the analysis.
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.
This thesis is concerned with modelling rockfall parameters associated with cliff collapse debris and the resultant “ramp” that formed following the high peak ground acceleration (PGA) events of 22 February 2011 and 13 June 2011. The Christchurch suburb of Redcliffs, located at the base of the Port Hills on the northern side of Banks Peninsula, New Zealand, is comprised of Miocene-age volcanics with valley-floor infilling marine sediments. The area is dominated by basaltic lava flows of the Mt Pleasant Formation, which is a suite of rocks forming part of the Lyttelton Volcanic Group that were erupted 11.0-10.0Ma. Fresh exposure enabled the identification of a basaltic ignimbrite unit at the study site overlying an orange tuff unit that forms a marker horizon spanning the length of the field area. Prior to this thesis, basaltic ignimbrite on Banks Peninsula has not been recorded, so descriptions and interpretations of this unit are the first presented. Mapping of the cliff face by remote observation, and analysis of hand samples collected from the base of the debris slopes, has identified a very strong (>200MPa), columnar-jointed, welded unit, and a very weak (<5MPa), massive, so-called brecciated unit that together represent the end-member components of the basaltic ignimbrite. Geochemical analysis shows the welded unit is picrite basalt, and the brecciated unit is hawaiite, making both clearly distinguishable from the underlying trachyandesite tuff. RocFall™ 4.0 was used to model future rockfalls at Redcliffs. RocFall™ is a two-dimensional (2D), hybrid, probabilistic modelling programme for which topographical profile data is used to generate slope profiles. GNS Science collected the data used for slope profile input in March 2011. An initial sensitivity analysis proved the Terrestrial Laser Scan (TLS)-derived slope to be too detailed to show any results when the slope roughness parameter was tested. A simplified slope profile enabled slope roughness to be varied, however the resulting model did not correlate with field observations as well. By using slope profile data from March 2011, modelled rockfall behaviour has been calibrated with observed rockfall runout at Redcliffs in the 13 June 2011 event to create a more accurate rockfall model. The rockfall model was developed on a single slope profile (Section E), with the chosen model then applied to four other section lines (A-D) to test the accuracy of the model, and to assess future rockfall runout across a wider area. Results from Section Lines A, B, and E correlate very well with field observations, with <=5% runout exceeding the modelled slope, and maximum bounce height at the toe of the slope <=1m. This is considered to lie within observed limits given the expectation that talus slopes will act as a ramp on which modelled rocks travel further downslope. Section Lines C and D produced higher runout percentage values than the other three section lines (23% and 85% exceeding the base of the slope, respectively). Section D also has a much higher maximum bounce height at the toe of the slope (~8.0m above the slope compared to <=1.0m for the other four sections). Results from modelling of all sections shows the significance of the ratio between total cliff height (H) and horizontal slope distance (x), and of maximum drop height to the top of the talus (H*) and horizontal slope distance (x). H/x can be applied to the horizontal to vertical ratio (H:V) as used commonly to identify potential slope instability. Using the maximum value from modelling at Redcliffs, the future runout limit can be identified by applying a 1.4H:1V ratio to the remainder of the cliff face. Additionally, the H*/x parameter shows that when H*/x >=0.6, the percentage of rock runout passing the toe of the slope will exceed 5%. When H*/x >=0.75, the maximum bounce height at the toe of the slope can be far greater than when H*/x is below this threshold. Both of these parameters can be easily obtained, and can contribute valuable guideline data to inform future land-use planning decisions. This thesis project has demonstrated the applicability of a 2D probabilistic-based model (RocFall™ 4.0) to evaluate rockfall runout on the talus slope (or ramp) at the base of ~35-70m high cliff with a basaltic ignimbrite source. Limitations of the modelling programme have been identified, in particular difficulties with adjusting modelled roughness of the slope profile and the inability to consider fragmentation. The runout profile using RocFall™ has been successfully calibrated against actual profiles and some anomalous results have been identified.
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’.
Deconstruction, at the end of the useful life of a building, produces a considerable amount of materials which must be disposed of, or be recycled / reused. At present, in New Zealand, most timber construction and demolition (C&D) material, particularly treated timber, is simply waste and is placed in landfills. For both technical and economic reasons (and despite the increasing cost of landfills), this position is unlikely to change in the next 10 – 15 years unless legislation dictates otherwise. Careful deconstruction, as opposed to demolition, can provide some timber materials which can be immediately re-used (eg. doors and windows), or further processed into other components (eg. beams or walls) or recycled (‘cascaded’) into other timber or composite products (e.g. fibre-board). This reusing / recycling of materials is being driven slowly in NZ by legislation, the ‘greening’ of the construction industry and public pressure. However, the recovery of useful material can be expensive and uneconomic (as opposed to land-filling). In NZ, there are few facilities which are able to sort and separate timber materials from other waste, although the soon-to-be commissioned Burwood Resource Recovery Park in Christchurch will attempt to deal with significant quantities of demolition waste from the recent earthquakes. The success (or otherwise) of this operation should provide good information as to how future C&D waste will be managed in NZ. In NZ, there are only a few, small scale facilities which are able to burn waste wood for energy recovery (e.g. timber mills), and none are known to be able to handle large quantities of treated timber. Such facilities, with constantly improving technology, are being commissioned in Europe (often with Government subsidies) and this indicates that similar bio-energy (co)generation will be established in NZ in the future. However, at present, the NZ Government provides little assistance to the bio-energy industry and the emergence worldwide of shale-gas reserves is likely to push the economic viability of bio-energy further into the future. The behaviour of timber materials placed in landfills is complex and poorly understood. Degrading timber in landfills has the potential to generate methane, a potent greenhouse gas, which can escape to the atmosphere and cancel out the significant benefits of carbon sequestration during tree growth. Improving security of landfills and more effective and efficient collection and utilisation of methane from landfills in NZ will significantly reduce the potential for leakage of methane to the atmosphere, acting as an offset to the continuing use of underground fossil fuels. Life cycle assessment (LCA), an increasingly important methodology for quantifying the environmental impacts of building materials (particularly energy, and global warming potential (GWP)), will soon be incorporated into the NZ Green Building Council Greenstar rating tools. Such LCA studies must provide a level playing field for all building materials and consider the whole life cycle. Whilst the end-of-life treatment of timber by LCA may establish a present-day base scenario, any analysis must also present a realistic end-of-life scenario for the future deconstruction of any 6 new building, as any building built today will be deconstructed many years in the future, when very different technologies will be available to deal with construction waste. At present, LCA practitioners in NZ and Australia place much value on a single research document on the degradation of timber in landfills (Ximenes et al., 2008). This leads to an end-of-life base scenario for timber which many in the industry consider to be an overestimation of the potential negative effects of methane generation. In Europe, the base scenario for wood disposal is cascading timber products and then burning for energy recovery, which normally significantly reduces any negative effects of the end-of-life for timber. LCA studies in NZ should always provide a sensitivity analysis for the end-of-life of timber and strongly and confidently argue that alternative future scenarios are realistic disposal options for buildings deconstructed in the future. Data-sets for environmental impacts (such as GWP) of building materials in NZ are limited and based on few research studies. The compilation of comprehensive data-sets with country-specific information for all building materials is considered a priority, preferably accounting for end-of-life options. The NZ timber industry should continue to ‘champion’ the environmental credentials of timber, over and above those of the other major building materials (concrete and steel). End-of-life should not be considered the ‘Achilles heel’ of the timber story.
This thesis describes the management process of innovation through construction infrastructure projects. This research focuses on the innovation management process at the project level from four views. These are categorised into the separate yet related areas of: “innovation definition”, “Project time”, “project team motivation” and “Project temporary organisation”. A practical knowledge is developed for each of these research areas that enables project practitioners to make the best decision for the right type of innovation at the right phase of projects, through a capable project organisation. The research developed a holistic view on both innovation and the construction infrastructure project as two complex phenomena. An infrastructure project is a long-term capital investment, highly risky and an uncertain. Infrastructure projects can play a key role in innovation and performance improvement throughout the construction industry. The delivery of an infrastructure project is affected in most cases by critical issues of budget constraint, programme delays and safety Where the business climate is characterized by uncertainty, risk and a high level of technological change, construction infrastructure projects are unable to cope with the requirement to develop innovation. Innovation in infrastructure projects, as one of the key performance indicators (KPI) has been identified as a critical capability for performance improvement through the industry. However, in spite of the importance of infrastructure projects in improving innovation, there are a few research efforts that have developed a comprehensive view on the project context and its drivers and inhibitors for innovation in the construction industry. Two main reasons are given as the inhibitors through the process of comprehensive research on innovation management in construction. The first reason is the absence of an understanding of innovation itself. The second is a bias towards research at a firm and individual level, so a comprehensive assessment of project-related factors and their effects on innovation in infrastructure projects has not been undertaken. This study overcomes these issues by adopting as a case study approach of a successful infrastructure project. This research examines more than 500 construction innovations generated by a unique infrastructure alliance. SCIRT (Stronger Christchurch Infrastructure Rebuild Team) is a temporary alliancing organisation that was created to rebuild and recover the damaged infrastructure after the Christchurch 2011 earthquake. Researchers were given full access to the innovation project information and innovation systems under a contract with SCIRT Learning Legacy, provided the research with material which is critical for understanding innovations in large, complex alliancing infrastructure organisation. In this research, an innovation classification model was first constructed. Clear definitions have been developed for six types of construction innovation with a variety of level of novelties and benefits. The innovation classification model was applied on the SCIRT innovation database and the resultant trends and behaviours of different types of innovation are presented. The trends and behaviours through different types of SCIRT innovations developed a unique opportunity to research the projectrelated factors and their effect on the behaviour of different classified types of innovation throughout the project’s lifecycle. The result was the identification of specific characteristics of an infrastructure project that affect the innovation management process at the project level. These were categorised in four separate chapters. The first study presents the relationship between six classified types of innovation, the level of novelty and the benefit they come up with, by applying the innovation classification model on SCIRT innovation database. The second study focused on the innovation potential and limitations in different project lifecycle phases by using a logic relationship between the six classified types of innovation and the three classified phases of the SCIRT project. The third study result develops a holistic view of different elements of the SCIRT motivation system and results in a relationship between the maturity level of definition developed for innovation as one of the KPIs and a desire though the SCIRT innovation incentive system to motivate more important innovations throughout the project. The fourth study is about the role of the project’s temporary organisation that finally results in a multiple-view innovation model being developed for project organisation capability assessment in the construction industry. The result of this thesis provides practical and instrumental knowledge to be used by a project practitioner. Benefits of the current thesis could be categorized in four groups. The first group is the innovation classification model that provides a clear definition for six classified types of innovation with four levels of novelty and specifically defined outcomes and the relationship between the innovation types, novelty and benefit. The second is the ability that is provided for the project practitioner to make the best decision for the right type of innovation at the right phases of a project’s lifecycle. The third is an optimisation that is applied on the SCIRT innovation motivation system that enables the project practitioner to incentivize the right type of innovation with the right level of financial gain. This drives the project teams to develop a more important innovation instead of a simple problemsolving one. Finally, the last and probably more important benefit is the recommended multiple-view innovation model. This is a tool that could be used by a project practitioner in order to empower the project team to support innovation throughout the project.
One of the great challenges facing human systems today is how to prepare for, manage, and adapt successfully to the profound and rapid changes wreaked by disasters. Wellington, New Zealand, is a capital city at significant risk of devastating earthquake and tsunami, potentially requiring mass evacuations with little or short notice. Subsequent hardship and suffering due to widespread property damage and infrastructure failure could cause large areas of the Wellington Region to become uninhabitable for weeks to months. Previous research has shown that positive health and well-being are associated with disaster-resilient outcomes. Preventing adverse outcomes before disaster strikes, through developing strengths-based skill sets in health-protective attitudes and behaviours, is increasingly advocated in disaster research, practise, and management. This study hypothesised that well-being constructs involving an affective heuristic play vital roles in pathways to resilience as proximal determinants of health-protective behaviours. Specifically, this study examined the importance of health-related quality of life and subjective well-being in motivating evacuation preparedness, measured in a community sample (n=695) drawn from the general adult population of Wellington’s isolated eastern suburbs. Using a quantitative epidemiological approach, the study measured the prevalence of key quality of life indicators (physical and mental health, emotional well-being or “Sense of Coherence”, spiritual well-being, social well-being, and life satisfaction) using validated psychometric scales; analysed the strengths of association between these indicators and the level of evacuation preparedness at categorical and continuous levels of measurement; and tested the predictive power of the model to explain the variance in evacuation preparedness activity. This is the first study known to examine multi-dimensional positive health and global well-being as resilient processes for engaging in evacuation preparedness behaviour. A cross-sectional study design and quantitative survey were used to collect self-report data on the study variables; a postal questionnaire was fielded between November 2008 and March 2009 to a sampling frame developed through multi-stage cluster randomisation. The survey response rate was 28.5%, yielding a margin of error of +/- 3.8% with 95% confidence and 80% statistical power to detect a true correlation coefficient of 0.11 or greater. In addition to the primary study variables, data were collected on demographic and ancillary variables relating to contextual factors in the physical environment (risk perception of physical and personal vulnerability to disaster) and the social environment (through the construct of self-determination), and other measures of disaster preparedness. These data are reserved for future analyses. Results of correlational and regression analyses for the primary study variables show that Wellingtonians are highly individualistic in how their well-being influences their preparedness, and a majority are taking inadequate action to build their resilience to future disaster from earthquake- or tsunami-triggered evacuation. At a population level, the conceptual multi-dimensional model of health-related quality of life and global well-being tested in this study shows a positive association with evacuation preparedness at statistically significant levels. However, it must be emphasised that the strength of this relationship is weak, accounting for only 5-7% of the variability in evacuation preparedness. No single dimension of health-related quality of life or well-being stands out as a strong predictor of preparedness. The strongest associations for preparedness are in a positive direction for spiritual well-being, emotional well-being, and life satisfaction; all involve a sense of existential meaningfulness. Spiritual well-being is the only quality of life variable making a statistically significant unique contribution to explaining the variance observed in the regression models. Physical health status is weakly associated with preparedness in a negative direction at a continuous level of measurement. No association was found at statistically significant levels for mental health status and social well-being. These findings indicate that engaging in evacuation preparedness is a very complex, holistic, yet individualised decision-making process, and likely involves highly subjective considerations for what is personally relevant. Gender is not a factor. Those 18-24 years of age are least likely to prepare and evacuation preparedness increases with age. Multidimensional health and global well-being are important constructs to consider in disaster resilience for both pre-event and post-event timeframes. This work indicates a need for promoting self-management of risk and building resilience by incorporating a sense of personal meaning and importance into preparedness actions, and for future research into further understanding preparedness motivations.
Deformational properties of soil, in terms of modulus and damping, exert a great influence on seismic response of soil sites. However, these properties for sands containing some portion of fines particles have not been systematically addressed. In addition, simultaneous modelling of the modulus and damping behaviour of soils during cyclic loading is desirable. This study presents an experimental and computational investigation into the deformational properties of sands containing fines content in the context of site response analysis. The experimental investigation is carried on sandy soils sourced from Christchurch, New Zealand using a dynamic triaxial apparatus while the computational aspect is based on the framework of total-stress one-dimensional (1D) cyclic behaviour of soil. The experimental investigation focused on a systematic study on the deformational behaviour of sand with different amounts of fines content (particle diameter ≤ 75µm) under drained conditions. The silty sands were prepared by mixing clean sand with three different percentages of fines content. A series of bender element tests at small-strain range and stress-controlled dynamic triaxial tests at medium to high-strain ranges were conducted on samples of clean sand and silty sand. This allowed measurements of linear and nonlinear deformational properties of the same specimen for a wide strain range. The testing program was designed to quantify the effects of void ratio and fines content on the low-strain stiffness of the silty sand as well as on the nonlinear stress-strain relationship and corresponding shear modulus and damping properties as a function of cyclic shear strains. Shear wave velocity, Vs, and maximum shear modulus, Gmax, of silty sand was shown to be significantly smaller than the respective values for clean sands measured at the same void ratio, e, or same relative density, Dr. However, the test results showed that the difference in the level of nonlinearity between clean sand and silty sands was small. For loose samples prepared at an identical relative density, the behaviour of clean sand was slightly less nonlinear as compared to sandy soils with higher fines content. This difference in the nonlinear behaviour of clean sand and sandy soils was negligible for dense soils. Furthermore, no systematic influence of fines content on the material damping curve was observed for sands with fines content FC = 0 to 30%. In order to normalize the effects of fines on moduli of sands, equivalent granular void ratio, e*, was employed. This was done through quantifying the participation of fines content in the force transfer chain of the sand matrix. As such, a unified framework for modelling of the variability of shear wave velocity, Vs, (or shear modulus, Gmax) with void ratio was achieved for clean sands and sands with fines, irrespective of their fines content. Furthermore, modelling of the cyclic stress-strain behaviour based on this experimental program was investigated. The modelling effort focused on developing a simple constitutive model which simultaneously models the soil modulus and damping relationships with shear strains observed in laboratory tests. The backbone curve of the cyclic model was adopted based on a modified version of Kondner and Zelasko (MKZ) hyperbolic function, with a curvature coefficient, a. In order to simulate the hysteretic cycles, the conventional Masing rules (Pyke 1979) were revised. The parameter n, in the Masing’s criteria was assumed to be a function of material damping, h, measured in the laboratory. As such the modulus and damping produced by the numerical model could match the stress-strain behaviour observed in the laboratory over the course of this study. It was shown that the Masing parameter n, is strain-dependent and generally takes values of n ≤ 2. The model was then verified through element test simulations under different cyclic loadings. It was shown that the model could accurately simulate the modulus and the damping simultaneously. The model was then incorporated within the OpenSees computational platform and was used to scrutinize the effects of damping on one-dimensional seismic site response analysis. For this purpose, several strong motion stations which recorded the Canterbury earthquake sequence were selected. The soil profiles were modelled as semi-infinite horizontally layered deposits overlying a uniform half-space subjected to vertically propagating shear waves. The advantages and limitations of the nonlinear model in terms of simulating soil nonlinearity and associated material damping were further scrutinized. It was shown that generally, the conventional Masing criteria unconservatively may underestimate some response parameters such as spectral accelerations. This was shown to be due to larger hysteretic damping modelled by using conventional Masing criteria. In addition, maximum shear strains within the soil profiles were also computed smaller in comparison to the values calculated by the proposed model. Further analyses were performed to study the simulation of backbone curve beyond the strain ranges addressed in the experimental phase of this study. A key issue that was identified was that relying only on the modulus reduction curves to simulate the stress-strain behaviour of soil may not capture the actual soil strength at larger strains. Hence, strength properties of the soil layer should also be incorporated to accurately simulate the backbone curve.
Natural hazard disasters often have large area-wide impacts, which can cause adverse stress-related mental health outcomes in exposed populations. As a result, increased treatment-seeking may be observed, which puts a strain on the limited public health care resources particularly in the aftermath of a disaster. It is therefore important for public health care planners to know whom to target, but also where and when to initiate intervention programs that promote emotional wellbeing and prevent the development of mental disorders after catastrophic events. A large body of literature assesses factors that predict and mitigate disaster-related mental disorders at various time periods, but the spatial component has rarely been investigated in disaster mental health research. This thesis uses spatial and spatio-temporal analysis techniques to examine when and where higher and lower than expected mood and anxiety symptom treatments occurred in the severely affected Christchurch urban area (New Zealand) after the 2010/11 Canterbury earthquakes. High-risk groups are identified and a possible relationship between exposure to the earthquakes and their physical impacts and mood and anxiety symptom treatments is assessed. The main research aim is to test the hypothesis that more severely affected Christchurch residents were more likely to show mood and anxiety symptoms when seeking treatment than less affected ones, in essence, testing for a dose-response relationship. The data consisted of mood and anxiety symptom treatment information from the New Zealand Ministry of Health’s administrative databases and demographic information from the National Health Index (NHI) register, when combined built a unique and rich source for identifying publically funded stress-related treatments for mood and anxiety symptoms in almost the whole population of the study area. The Christchurch urban area within the Christchurch City Council (CCC) boundary was the area of interest in which spatial variations in these treatments were assessed. Spatial and spatio-temporal analyses were done by applying retrospective space-time and spatial variation in temporal trends analysis using SaTScan™ software, and Bayesian hierarchical modelling techniques for disease mapping using WinBUGS software. The thesis identified an overall earthquake-exposure effect on mood and anxiety symptom treatments among Christchurch residents in the context of the earthquakes as they experienced stronger increases in the risk of being treated especially shortly after the catastrophic 2011 Christchurch earthquake compared to the rest of New Zealand. High-risk groups included females, elderly, children and those with a pre-existing mental illness with elderly and children especially at-risk in the context of the earthquakes. Looking at the spatio-temporal distribution of mood and anxiety symptom treatments in the Christchurch urban area, a high rates cluster ranging from the severely affected central city to the southeast was found post-disaster. Analysing residential exposure to various earthquake impacts found that living in closer proximity to more affected areas was identified as a risk factor for mood and anxiety symptom treatments, which largely confirms a dose-response relationship between level of affectedness and mood and anxiety symptom treatments. However, little changes in the spatial distribution of mood and anxiety symptom treatments occurred in the Christchurch urban area over time indicating that these results may have been biased by pre-existing spatial disparities. Additionally, the post-disaster mobility activity from severely affected eastern to the generally less affected western and northern parts of the city seemed to have played an important role as the strongest increases in treatment rates occurred in less affected northern areas of the city, whereas the severely affected eastern areas tended to show the lowest increases. An investigation into the different effects of mobility confirmed that within-city movers and temporary relocatees were generally more likely to receive care or treatment for mood or anxiety symptoms, but moving within the city was identified as a protective factor over time. In contrast, moving out of the city from minor, moderately or severely damaged plain areas of the city, which are generally less affluent than Port Hills areas, was identified as a risk factor in the second year post-disaster. Moreover, residents from less damaged plain areas of the city showed a decrease in the likelihood of receiving care or treatment for mood or anxiety symptoms compared to those from undamaged plain areas over time, which also contradicts a possible dose-response relationship. Finally, the effects of the social and physical environment, as well as community resilience on mood and anxiety symptom treatments among long-term stayers from Christchurch communities indicate an exacerbation of pre-existing mood and anxiety symptom treatment disparities in the city, whereas exposure to ‘felt’ earthquake intensities did not show a statistically significant effect. The findings of this thesis highlight the complex relationship between different levels of exposure to a severe natural disaster and adverse mental health outcomes in a severely affected region. It is one of the few studies that have access to area-wide health and impact information, are able to do a pre-disaster / post-disaster comparison and track their sample population to apply spatial and spatio-temporal analysis techniques for exposure assessment. Thus, this thesis enhances knowledge about the spatio-temporal distribution of adverse mental health outcomes in the context of a severe natural disaster and informs public health care planners, not only about high-risk groups, but also where and when to target health interventions. The results indicate that such programs should broadly target residents living in more affected areas as they are likely to face daily hardship by living in a disrupted environment and may have already been the most vulnerable ones before the disaster. Special attention should be focussed on women, elderly, children and people with pre-existing mental illnesses as they are most likely to receive care or treatment for stress-related mental health symptoms. Moreover, permanent relocatees from affected areas and temporarily relocatees shortly after the disaster may need special attention as they face additional stressors due to the relocation that may lead to the development of adverse mental health outcomes needing treatment.
Mixed conifer, beech and hardwood forests are relatively common in Aotearoa/New Zealand, but are not well studied. This thesis investigates the coexistence, regeneration dynamics and disturbance history of a mixed species forest across an environmental gradient of drainage and soil development in north Westland. The aim was to investigate whether conifers, beech and non-beech hardwood species were able to coexist on surfaces that differed in their underlying edaphic conditions, and if so to understand the mechanisms that influenced their regeneration on both poorly drained and well drained soils. The site selected was an area of high tree species diversity on a lowland 0.8 km² post-glacial terrace at the base of Mount Harata in the Grey River Valley. My approach was to use forest stand history reconstruction at two spatial scales: an intensive within-plot study of stand dynamics (chapter 1) and a whole-landform approach (chapter 2) that examined whether the dynamics identified at the smaller within-plot scale reflected larger patterns across the terrace. In chapter 1, three large permanent plots (0.3-0.7 ha) were placed at different points along the drainage gradient, one plot situated in each of the mainly well-drained, poorly drained and very poorly drained areas along the terrace. Information was gathered on species age and size structures, spatial distributions of tree ages, species interactions, microsite establishment preferences, patterns of stand mortality, and disturbance history in each plot. There were differences in stand structure, composition and relative abundance of species found between the well drained plot and the two poorer drained plots. On the well drained site conifers were scarce, the beeches Nothofagus fusca and N. menziesii dominated the canopy, and in the subcanopy the hardwood species Weinmannia racemosa and Quintinia acutifolia were abundant. As drainage became progressively poorer, the conifers Dacrydium cupressinum and Dacrycarpus dacrydioides became more abundant and occupied the emergent tier over a beech canopy. The hardwoods W. racemosa and Q. acutifolia became gradually less abundant in the subcanopy, whereas the hardwood Elaeocarpus hookerianus became more so. In the well drained plot, gap partitioning for light between beeches and hardwoods enabled coexistence in response to a range of different sized openings resulting from disturbances of different extent. In the two more poorly drained plots, species also coexisted by partitioning microsite establishment sites according to drainage. There were several distinct periods where synchronous establishment of different species occurred in different plots, suggesting there were large disturbances: c. 100yrs, 190-200 yrs, 275-300 yrs and 375-425 yrs ago. Generally after the same disturbance, different species regenerated in different plots reflecting the underlying drainage gradient. However, at the same site after different disturbances, different sets of species regenerated, suggesting the type and extent of disturbances and the conditions left behind influenced species regeneration at some times but not others. The regeneration of some species (e.g., N. fusca in the well-drained plot, and Dacrydium in the poorer drained plots) was periodic and appeared to be closely linked to these events. In the intervals between these disturbances, less extensive disturbances resulted in the more frequent N. menziesii and especially hardwood regeneration. The type of tree death caused by different disturbances favoured different species, with dead standing tree death favouring the more shade-tolerant N. menziesii and hardwoods, whereas uprooting created a mosaic of microsite conditions and larger gap sizes that enabled Dacrycarpus, N. fusca and E. hookerianus to maintain themselves in the poorly drained areas. In chapter 2, 10 circular plots (c. 0.12 ha) were placed in well drained areas and 10 circular plots (c. 0.2 ha) in poorly drained plots to collect information on species population structures and microsite preferences. The aims were to reconstruct species' regeneration responses to a range of disturbances of different type and extent across the whole terrace, and to examine whether there were important differences in the effects of these disturbances. At this landform scale, the composition and relative abundances of species across the drainage gradient reflected those found in chapter 1. There were few scattered conifers in well drained areas, despite many potential regeneration opportunities created from a range of different stand destroying and smaller scale disturbances. Three of the four periods identified in chapter 1 reflected distinct terrace-wide periods of regeneration 75-100 yrs, 200-275 yrs and 350-450 yrs ago, providing strong evidence of periodic large, infrequent disturbances that occurred at intervals of 100-200 yrs. These large, infrequent disturbances have had a substantial influence in determining forest history, and have had long term effects on forest structure and successional processes. Different large, infrequent disturbances had different effects across the terrace, with the variability in conditions that resulted enabling different species to regenerate at different times. For example, the regeneration of distinct even-aged Dacrydium cohorts in poorly drained areas was linked to historical Alpine Fault earthquakes, but not to more recent storms. The variation in the intensity of different large, infrequent disturbances at different points along the environmental drainage gradient, was a key factor influencing the scale of impacts. In effect, the underlying edaphic conditions influenced species composition along the drainage gradient and disturbance history regulated the relative abundances of species. The results presented here further emphasise the importance of large scale disturbances as a mechanism that allows coexistence of different tree species in mixed forest, in particular for the conifers Dacrydium, Dacrycarpus and the beech N. fusca, by creating much of the environmental variation to which these species responded. This study adds to our understanding of the effects of historical earthquakes in the relatively complex forests of north Westland, and further illustrates their importance in the Westland forest landscape as the major influential disturbance on forest pattern and history. These results also further develop the 'two-component' model used to describe conifer/angiosperm dynamics, by identifying qualitative differences in the impacts of different large, infrequent disturbances across an environmental gradient that allowed for coexistence of different species. In poorer drained areas, these forests may even be thought of as 'three-component' systems with conifers, beeches and hardwoods exhibiting key differences in their regeneration patterns after disturbances of different type and extent, and in their microsite preferences.
Environmental stress and disturbance can affect the structure and functioning of marine ecosystems by altering their physical, chemical and biological features. In estuaries, benthic invertebrate communities play important roles in structuring sediments, influencing primary production and biogeochemical flux, and occupying key food web positions. Stress and disturbance can reduce species diversity, richness and abundance, with ecological theory predicting that biodiversity will be at its lowest soon after a disturbance with assemblages dominated by opportunistic species. The Avon-Heathcote Estuary in Christchurch New Zealand has provided a novel opportunity to examine the effects of stress, in the form of eutrophication, and disturbance, in the form of cataclysmic earthquake events, on the structure and functioning of an estuarine ecosystem. For more than 50 years, large quantities (up to 500,000m3/day) of treated wastewater were released into this estuary but in March 2010 this was diverted to an ocean outfall, thereby reducing the nutrient loading by around 90% to the estuary. This study was therefore initially focussed on the reversal of eutrophication and consequent effects on food web structure in the estuary as it responded to lower nutrients. In 2011, however, Christchurch was struck with a series of large earthquakes that greatly changed the estuary. Massive amounts of liquefied sediments, covering up to 65% of the estuary floor, were forced up from deep below the estuary, the estuary was tilted by up to a 50cm rise on one side and a corresponding drop on the other, and large quantities of raw sewage from broken wastewater infrastructure entered the estuary for up to nine months. This study was therefore a test of the potentially synergistic effects of nutrient reduction and earthquake disturbance on invertebrate communities, associated habitats and food web dynamics. Because there was considerable site-to-site heterogeneity in the estuary, the sites in this study were selected to represent a eutrophication gradient from relatively “clean” (where the influence of tidal flows was high) to highly impacted (near the historical discharge site). The study was structured around these sites, with components before the wastewater diversion, after the diversion but before the earthquakes, and after the earthquakes. The eutrophication gradient was reflected in the composition and isotopic chemistry of primary producer and invertebrate communities and the characteristics of sediments across the sample sites. Sites closest to the former wastewater discharge pipe were the most eutrophic and had cohesive organic -rich, fine sediments and relatively depauperate communities dominated by the opportunistic taxa Capitellidae. The less-impacted sites had coarser, sandier sediments with fewer pollutants and far less organic matter than at the eutrophic sites, relatively high diversity and lower abundances of micro- and macro-algae. Sewage-derived nitrogen had became incorporated into the estuarine food web at the eutrophic sites, starting at the base of the food chain with benthic microalgae (BMA), which were found to use mostly sediment-derived nitrogen. Stable isotopic analysis showed that δ13C and δ15N values of most food sources and consumers varied spatially, temporally and in relation to the diversion of wastewater, whereas the earthquakes did not appear to affect the overall estuarine food web structure. This was seen particularly at the most eutrophic site, where isotopic signatures became more similar to the cleaner sites over two-and-a-half years after the diversion. New sediments (liquefaction) produced by the earthquakes were found to be coarser, have lower concentrations of heavy metals and less organic matter than old (existing) sediments. They also had fewer macroinvertebrate inhabitants initially after the earthquakes but most areas recovered to pre-earthquake abundance and diversity within two years. Field experiments showed that there were higher amounts of primary production and lower amounts of nutrient efflux from new sediments at the eutrophic sites after the earthquakes. Primary production was highest in new sediments due to the increased photosynthetic efficiency of BMA resulting from the increased permeability of new sediments allowing increased light penetration, enhanced vertical migration of BMA and the enhanced transport of oxygen and nutrients. The reduced efflux of NH4-N in new sediments indicated that the capping of a large portion of eutrophic old sediments with new sediments had reduced the release of legacy nutrients (originating from the historical discharge) from the sediments to the overlying water. Laboratory experiments using an array of species and old and new sediments showed that invertebrates altered levels of primary production and nutrient flux but effects varied among species. The mud snail Amphibola crenata and mud crab Austrohelice crassa were found to reduce primary production and BMA biomass through the consumption of BMA (both species) and its burial from bioturbation and the construction of burrows (Austrohelice). In contrast, the cockle Austrovenus stutchburyi did not significantly affect primary production and BMA biomass. These results show that changes in the structure of invertebrate communities resulting from disturbances can also have consequences for the functioning of the system. The major conclusions of this study were that the wastewater diversion had a major effect on food web dynamics and that the large quantities of clean and unpolluted new sediments introduced to the estuary during the earthquakes altered the recovery trajectory of the estuary, accelerating it at least throughout the duration of this study. This was largely through the ‘capping’ effect of the new liquefied, coarser-grained sediments as they dissipated across the estuary and covered much of the old organic-rich eutrophic sediments. For all aspects of this study, the largest changes occurred at the most eutrophic sites; however, the surrounding habitats were important as they provided the context for recovery of the estuary, particularly because of the very strong influence of sediments, their biogeochemistry, microalgal and macroalgal dynamics. There have been few studies documenting system level responses to eutrophication amelioration and to the best on my knowledge there are no other published studies examining the impacts of large earthquakes on benthic communities in an estuarine ecosystem. This research gives valuable insight and advancements in the scientific understanding of the effects that eutrophication recovery and large-scale disturbances can have on the ecology of a soft-sediment ecosystem.