Recent experiences from the Darfield and Canterbury, New Zealand earthquakes have shown that the soft soil condition of saturated liquefiable sand has a profound effect on seismic response of buildings, bridges and other lifeline infrastructure. For detailed evaluation of seismic response three dimensional integrated analysis comprising structure, foundation and soil is required; such an integrated analysis is referred to as Soil Foundation Structure Interaction (SFSI) in literatures. SFSI is a three-dimensional problem because of three primary reasons: first, foundation systems are three-dimensional in form and geometry; second, ground motions are three-dimensional, producing complex multiaxial stresses in soils, foundations and structure; and third, soils in particular are sensitive to complex stress because of heterogeneity of soils leading to a highly anisotropic constitutive behaviour. In literatures the majority of seismic response analyses are limited to plane strain configuration because of lack of adequate constitutive models both for soils and structures, and computational limitation. Such two-dimensional analyses do not represent a complete view of the problem for the three reasons noted above. In this context, the present research aims to develop a three-dimensional mathematical formulation of an existing plane-strain elasto-plastic constitutive model of sand developed by Cubrinovski and Ishihara (1998b). This model has been specially formulated to simulate liquefaction behaviour of sand under ground motion induced earthquake loading, and has been well-validated and widely implemented in verifcation of shake table and centrifuge tests, as well as conventional ground response analysis and evaluation of case histories. The approach adopted herein is based entirely on the mathematical theory of plasticity and utilises some unique features of the bounding surface plasticity formalised by Dafalias (1986). The principal constitutive parameters, equations, assumptions and empiricism of the existing plane-strain model are adopted in their exact form in the three-dimensional version. Therefore, the original two-dimensional model can be considered as a true subset of the three-dimensional form; the original model can be retrieved when the tensorial quantities of the three dimensional version are reduced to that of the plane-strain configuration. Anisotropic Drucker-Prager type failure surface has been adopted for the three-dimensional version to accommodate triaxial stress path. Accordingly, a new mixed hardening rule based on Mroz’s approach of homogeneous surfaces (Mroz, 1967) has been introduced for the virgin loading surface. The three-dimensional version is validated against experimental data for cyclic torsional and triaxial stress paths.
Background and methodology The Mw 7.8, 14th November 2016 earthquake centred (item b, figure 1) in the Hurunui District of the South Island, New Zealand, damaged critical infrastructure across North Canterbury and Marlborough. We investigate the impacts to infrastructure and adaptations to the resulting service disruption in four small rural towns (figure 1): Culverden (a), Waiau (c), Ward (d) and Seddon (e). This is accomplished though literary research, interviews and geospatial analysis. Illustrating our methods, we have displayed here a Hurunui District hazard map (figure 2b) and select infrastructure inventories (figures 2a, 3).
This study examines the performance of nonlinear total-stress wave-propagation site response analysis for modelling site effects in physics-based ground motion simulations of the 2010-2011 Canterbury, New Zealand earthquake sequence. This approach allows for explicit modeling of 3-dimensional ground motion phenomena at the regional scale, as well as detailed site effects and soil nonlinearity at the local scale. The approach is compared to a more commonly used empirical VS30 (30 m time-averaged shear wave velocity)-based method for computing site amplification as proposed by Graves and Pitarka (2010, 2015).
Background This study examines the performance of site response analysis via nonlinear total-stress 1D wave-propagation for modelling site effects in physics-based ground motion simulations of the 2010-2011 Canterbury, New Zealand earthquake sequence. This approach allows for explicit modeling of 3D ground motion phenomena at the regional scale, as well as detailed nonlinear site effects at the local scale. The approach is compared to a more commonly used empirical VS30 (30 m time-averaged shear wave velocity)-based method for computing site amplification as proposed by Graves and Pitarka (2010, 2015), and to empirical ground motion prediction via a ground motion model (GMM).
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.
Liquefaction-induced lateral spreading in large seismic events often results in pervasive and costly damage to engineering structures and lifelines, making it a critical component of engineering design. However, the complex nature of this phenomenon leads to designing for such a hazard extremely challenging and there is a clear for an improved understanding and predicting liquefaction-induced lateral spreading. The 2010-2011 Canterbury (New Zealand) Earthquakes triggered severe liquefaction-induced lateral spreading along the streams and rivers of the Christchurch region, causing extensive damage to roads, bridges, lifelines, and structures in the vicinity. The unfortunate devastation induced from lateral spreading in these events also rendered the rare opportunity to gain an improved understanding of lateral spreading displacements specific to the Christchurch region. As part of this thesis, the method of ground surveying was employed following the 4 September 2010 Darfield (Mw 7.1) and 22 February 2011 Christchurch (Mw 6.2) earthquakes at 126 locations (19 repeated) throughout Christchurch and surrounding suburbs. The method involved measurements and then summation of crack widths along a specific alignment (transect) running approximately perpendicular to the waterway to indicate typically a maximum lateral displacement at the bank and reduction of the magnitude of displacements with distance from the river. Rigorous data processing and comparisons with alternative measurements of lateral spreading were performed to verify results from field observations and validate the method of ground surveying employed, as well as highlight the complex nature of lateral spreading displacements. The welldocumented field data was scrutinized to gain an understanding of typical magnitudes and distribution patterns (distribution of displacement with distance) of lateral spreading observed in the Christchurch area. Maximum displacements ranging from less than 10 cm to over 3.5 m were encountered at the sites surveyed and the area affected by spreading ranged from less than 20 m to over 200 m from the river. Despite the highly non-uniform displacements, four characteristic distribution patterns including large, distributed ground displacements, block-type movements, large and localized ground displacements, and areas of little to no displacements were identified. Available geotechnical, seismic, and topographic data were collated at the ground surveying sites for subsequent analysis of field measurements. Two widely-used empirical models (Zhang et al. (2004), Youd et al. (2002)) were scrutinized and applied to locations in the vicinity of field measurements for comparison with model predictions. The results indicated generally poor correlation (outside a factor of two) with empirical predictions at most locations and further validated the need for an improved, analysis- based method of predicting lateral displacements that considers the many factors involved on a site-specific basis. In addition, the development of appropriate model input parameters for the Youd et al. (2002) model led to a site-specific correlation of soil behavior type index, Ic, and fines content, FC, for sites along the Avon River in Christchurch that matched up well with existing Ic – FC relationships commonly used in current practice. Lastly, a rigorous analysis was performed for 25 selected locations of ground surveying measurements along the Avon River where ground slope conditions are mild (-1 to 2%) and channel heights range from about 2 – 4.5 m. The field data was divided into categories based on the observed distribution pattern of ground displacements including: large and distributed, moderate and distributed, small to negligible, and large and localized. A systematic approach was applied to determine potential critical layers contributing to the observed displacement patterns which led to the development of characteristic profiles for each category considered. The results of these analyses outline an alternative approach to the evaluation of lateral spreading in which a detailed geotechnical analysis is used to identify the potential for large spreading displacements and likely spatial distribution patterns of spreading. Key factors affecting the observed magnitude and distribution of spreading included the thickness of the critical layer, relative density, soil type and layer continuity. It was found that the large and distributed ground displacements were associated with a thick (1.5 – 2.5 m) deposit of loose, fine to silty sand (qc1 ~4-7 MPa, Ic 1.9-2.1, qc1n_cs ~50-70) that was continuous along the bank and with distance from the river. In contrast, small to negligible displacements were characterized by an absence of or relatively thin (< 1 m), discontinuous critical layer. Characteristic features of the moderate and distributed displacements were found to be somewhere between these two extremes. The localized and large displacements showed a characteristic critical layer similar to that observed in the large and distributed sites but that was not continuous and hence leading to the localized zone of displacement. The findings presented in this thesis illustrate the highly complex nature of lateral displacements that cannot be captured in simplified models but require a robust geotechnical analysis similar to that performed for this research.
This research briefing reports on the key findings of a computer-assisted text analysis of records from The Press newspaper related to the Earthquake Commission (EQC) from 2010 to 2019. The briefing has been prepared as a submission to the Public Inquiry into the Earthquake Commission. The aim of producing this research briefing is to provide the Public Inquiry with preliminary findings of a large-scale overview of media coverage on EQC and to identify and quantify key features and trends in public discourse about EQC over time. This research, which aggregates many stories and voices over time, offers a unique lens to view how EQC has been collectively represented, understood and experienced by the people of Canterbury.
This paper summarizes the development of a high-resolution surficial shear wave velocity model based on the combination of the large high-spatial-density database of cone penetration test (CPT) logs in and around Christchurch, New Zealand and a recently-developed Christchurch-specific empirical correlation between soil shear wave velocity and CPT. This near-surface shear wave velocity model has applications for site characterization efforts via the development of maps of time-averaged shear wave velocities over specific depths, as well as use in site response analysis and ground motion simulation.
This poster discusses several possible approaches by which the nonlinear response of surficial soils can be explicitly modelled in physics-based ground motion simulations, focusing on the relative advantages and limitations of the various methodologies. These methods include fully-coupled 3D simulation models that directly allow soil nonlinearity in surficial soils, the domain reduction method for decomposing the physical domain into multiple subdomains for separate simulation, conventional site response analysis uncoupled from the simulations, and finally, the use of simple empirically based site amplification factors We provide the methodology for an ongoing study to explicitly incorporate soil nonlinearity into hybrid broadband simulations of the 2010-2011 Canterbury, New Zealand earthquakes.
The Sendai Framework for Disaster Risk Reduction 2015-2030 finds that, despite progress in disaster risk reduction over the last decade “evidence indicates that exposure of persons and assets in all countries has increased faster than vulnerability has decreased, thus generating new risk and a steady rise in disaster losses” (p.4, UNISDR 2015). Fostering cooperation among relevant stakeholders and policy makers to “facilitate a science-policy interface for effective decisionmaking in disaster risk management” is required to achieve two priority areas for action, understanding disaster risk and enhancing disaster preparedness (p. 13, p. 23, UNISDR 2015). In other topic areas, the term science-policy interface is used interchangeably with the term boundary organisation. Both terms are usually used refer to systematic collaborative arrangements used to manage the intersection, or boundary, between science and policy domains, with the aim of facilitating the joint construction of knowledge to inform decision-making. Informed by complexity theory, and a constructivist focus on the functions and processes that minimize inevitable tensions between domains, this conceptual framework has become well established in fields where large complex issues have significant economic and political consequences, including environmental management, biodiversity, sustainable development, climate change and public health. To date, however, there has been little application of this framework in the disaster risk reduction field. In this doctoral project the boundary management framework informs an analysis of the research response to the 2010-2011 Canterbury Earthquake Sequence, focusing on the coordination role of New Zealand’s national Natural Hazards Research Platform. The project has two aims. It uses this framework to tell the nuanced story of the way this research coordination role evolved in response to both the complexity of the unfolding post-disaster environment, and to national policy and research developments. Lessons are drawn from this analysis for those planning and implementing arrangements across the science-policy boundary to manage research support for disaster risk reduction decision-making, particularly after disasters. The second aim is to use this case study to test the utility of the boundary management framework in the disaster risk reduction context. This requires that terminology and concepts are explained and translated in terms that make this analysis as accessible as possible across the disciplines, domains and sectors involved in disaster risk reduction. Key findings are that the focus on balance, both within organisations, and between organisations and domains, and the emphasis on systemic effects, patterns and trends, offer an effective and productive alternative to the more traditional focus on individual or organisational performance. Lessons are drawn concerning the application of this framework when planning and implementing boundary organisations in the hazard and disaster risk management context.
The purpose of this research is to investigate men’s experiences of the 2016 7.8 magnitude Kaikōura earthquake and Tsunami. While, research into the impacts of the earthquake has been conducted, few studies have examined how gender shaped people’s experiences of this natural hazard event. Analysing disasters through a gender lens has significantly contributed to disaster scholarship in identifying the resilience and vulnerabilities of individuals and communities pre- and post-disaster (Fordham, 2012; Bradshaw, 2013). This research employs understandings of masculinities (Connell, 2005), to examine men’s strengths and challenges in responding, recovering, and coping following the earthquake. Qualitative inquiry was carried out in Northern Canterbury and Marlborough involving 18 face-to-face interviews with men who were impacted by the Kaikōura earthquake and its aftermath. Interview material is being analysed using thematic and narrative analysis. Some of the preliminary findings have shown that men took on voluntary roles in addition to their fulltime paid work resulting in long hours, poor sleep and little time spent with family. Some men assisted wives and children to high ground then drove into the tsunami zone to check on relatives or to help evacuate people. Although analysis of the findings is currently ongoing, preliminary findings have identified that the men who participated in the study have been negatively impacted by the 2016 Kaikōura earthquake. A theme identified amongst participants was an avoidance to seek support with the challenges they were experiencing due to the earthquake. The research findings align with key characteristics of masculinity, including demonstrating risky behaviours and neglecting self or professional care. This study suggests that these behaviours affect men’s overall resilience, and thus the resilience of the wider community.
After a high-intensity seismic event, inspections of structural damages need to be carried out as soon as possible in order to optimize the emergency management, as well as improving the recovery time. In the current practice, damage inspections are performed by an experienced engineer, who physically inspect the structures. This way of doing not only requires a significant amount of time and high skilled human resources, but also raises the concern about the inspector’s safety. A promising alternative is represented using new technologies, such as drones and artificial intelligence, which can perform part of the damage classification task. In fact, drones can safely access high hazard components of the structures: for instance, bridge piers or abutments, and perform the reconnaissance by using highresolution cameras. Furthermore, images can be automatically processed by machine learning algorithms, and damages detected. In this paper, the possibility of applying such technologies for inspecting New Zealand bridges is explored. Firstly, a machine-learning model for damage detection by performing image analysis is presented. Specifically, the algorithm was trained to recognize cracks in concrete members. A sensitivity analysis was carried out to evaluate the algorithm accuracy by using database images. Depending on the confidence level desired,i.e. by allowing a manual classification where the alghortim confidence is below a specific tolerance, the accuracy was found reaching up to 84.7%. In the second part, the model is applied to detect the damage observed on the Anzac Bridge (GPS coordinates -43.500865, 172.701138) in Christchurch by performing a drone reconnaissance. Reults show that the accuracy of the damage detection was equal to 88% and 63% for cracking and spalling, respectively.
This thesis is concerned with springs that appeared in the Hillsborough, Christchurch during the 2010-2011 Canterbury Earthquake Sequence, and which have continued to discharge groundwater to the surface to the present time. Investigations have evolved, measurements of discharge at selected sites, limited chemical data on anions and isotope analysis. The springs are associated with earthquake generated fissures (extensional) and compression zones, mostly in loess-colluvium soils of the valley floor and lower slopes. Extensive peat swamps are present in the Hillsborough valley, with a groundwater table at ~1m below ground. The first appearance of the ‘new’ springs took place following the Mw 7.1 Darfield Earthquake on 4 September 2010, and discharges increased both in volume and extent of the Christchurch Mw 6.3 Earthquake of 22 February 2011. Five monitored sites show flow rates in the range of 4.2-14.4L/min, which have remained effectively constant for the duration of the study (2014-2015). Water chemistry analysis shows that the groundwater discharges are sourced primarily from volcanic bedrocks which underlies the valley at depths ≤50m below ground level. Isotope values confirm similarities with bedrock-sourced groundwater, and the short term (hours-days) influence of extreme rainfall events. Cyclone Lusi (2013-2014) affects were monitored and showed recovery of the bedrock derived water signature within 72 hours. Close to the mouth of the valley sediments interfinger with Waimakiriri River derived alluvium bearing a distinct and different isotope signature. Some mixing is evident at certain locations, but it is not clear if there is any influence from the Huntsbury reservoir which failed in the Port Hills Earthquake (22 February 2011) and stored groundwater from the Christchurch artesian aquifer system (Riccarton Gravel).
In recent work on commons and commoning, scholars have argued that we might delink the practice of commoning from property ownership, while paying attention to modes of governance that enable long-term commons to emerge and be sustained. Yet commoning can also occur as a temporary practice, in between and around other forms of use. In this article we reflect on the transitional commoning practices and projects enabled by the Christchurch post-earthquake organisation Life in Vacant Spaces, which emerged to connect and mediate between landowners of vacant inner city demolition sites and temporary creative or entrepreneurial users. While these commons are often framed as transitional or temporary, we argue they have ongoing reverberations changing how people and local government in Christchurch approach common use. Using the cases of the physical space of the Victoria Street site “The Commons” and the virtual space of the Life in Vacant Spaces website, we show how temporary commoning projects can create and sustain the conditions of possibility required for nurturing commoner subjectivities. Thus despite their impermanence, temporary commoning projects provide a useful counter to more dominant forms of urban development and planning premised on property ownership and “permanent” timeframes, in that just as the physical space of the city being opened to commoning possibilities, so too are the expectations and dispositions of the city’s inhabitants, planners, and developers.
In major seismic events, a number of plan-asymmetric buildings which experienced element failure or structural collapse had twisted significantly about their vertical axis during the earthquake shaking. This twist, known as “building torsion”, results in greater demands on one side of a structure than on the other side. The Canterbury Earthquakes Royal Commission’s reports describe the response of a number of buildings in the February 2011 Christchurch earthquakes. As a result of the catastrophic collapse of one multi-storey building with significant torsional irregularity, and significant torsional effects also in other buildings, the Royal Commission recommended that further studies be undertaken to develop improved simple and effective guides to consider torsional effects in buildings which respond inelastically during earthquake shaking. Separately from this, as building owners, the government, and other stakeholders, are planning for possible earthquake scenarios, they need good estimates of the likely performance of both new and existing buildings. These estimates, often made using performance based earthquake engineering considerations and loss estimation techniques, inform decision making. Since all buildings may experience torsion to some extent, and torsional effects can influence demands on building structural and non-structural elements, it is crucial that demand estimates consider torsion. Building seismic response considering torsion can be evaluated with nonlinear time history analysis. However, such analysis involves significant computational effort, expertise and cost. Therefore, from an engineers’ point of view, simpler analysis methods, with reasonable accuracy, are beneficial. The consideration of torsion in simple analysis methods has been investigated by many researchers. However, many studies are theoretical without direct relevance to structural design/assessment. Some existing methods also have limited applicability, or they are difficult to use in routine design office practice. In addition, there has been no consensus about which method is best. As a result, there is a notable lack of recommendations in current building design codes for torsion of buildings that respond inelastically. There is a need for building torsion to be considered in yielding structures, and for simple guidance to be developed and adopted into building design standards. This study aims to undertaken to address this need for plan-asymmetric structures which are regular over their height. Time history analyses are first conducted to quantify the effects of building plan irregularity, that lead to torsional response, on the seismic response of building structures. Effects of some key structural and ground motion characteristics (e.g. hysteretic model, ground motion duration, etc.) are considered. Mass eccentricity is found to result in rather smaller torsional response compared to stiffness/strength eccentricity. Mass rotational inertia generally decreases the torsional response; however, the trend is not clearly defined for torsionally restrained systems (i.e. large λty). Systems with EPP and bilinear models have close displacements and systems with Takeda, SINA, and flag-shaped models yield almost the same displacements. Damping has no specific effect on the torsional response for the single-storey systems with the unidirectional eccentricity and excitation. Displacements of the single-storey systems subject to long duration ground motion records are smaller than those for short duration records. A method to consider torsional response of ductile building structures under earthquake shaking is then developed based on structural dynamics for a wide range of structural systems and configurations, including those with low and high torsional restraint. The method is then simplified for use in engineering practice. A novel method is also proposed to simply account for the effects of strength eccentricity on response of highly inelastic systems. A comparison of the accuracy of some existing methods (including code-base equivalent static method and model response spectrum analysis method), and the proposed method, is conducted for single-storey structures. It is shown that the proposed method generally provides better accuracy over a wide range of parameters. In general, the equivalent static method is not adequate in capturing the torsional effects and the elastic modal response spectrum analysis method is generally adequate for some common parameters. Record-to-record variation in maximum displacement demand on the structures with different degrees of torsional response is considered in a simple way. Bidirectional torsional response is then considered. Bidirectional eccentricity and excitation has varying effects on the torsional response; however, it generally increases the weak and strong edges displacements. The proposed method is then generalized to consider the bidirectional torsion due to bidirectional stiffness/strength eccentricity and bidirectional seismic excitation. The method is shown to predict displacements conservatively; however, the conservatism decreases slightly for cases with bidirectional excitation compared to those subject to unidirectional excitation. In is shown that the roof displacement of multi-storey structures with torsional response can be predicted by considering the first mode of vibration. The method is then further generalized to estimate torsional effects on multi-storey structure displacement demands. The proposed procedure is tested multi-storey structures and shown to predict the displacements with a good accuracy and conservatively. For buildings which twist in plan during earthquake shaking, the effect of P-Δλ action is evaluated and recommendations for design are made. P-Δλ has more significant effects on systems with small post- yield stiffness. Therefore, system stability coefficient is shown not to be the best indicator of the importance of P-Δλ and it is recommended to use post-yield stiffness of system computed with allowance for P-Δλ effects. For systems with torsional response, the global system stability coefficient and post- yield stiffness ration do not reflect the significance of P-Δλ effects properly. Therefore, for torsional systems individual seismic force resisting systems should be considered. Accuracy of MRSA is investigated and it is found that the MRSA is not always conservative for estimating the centre of mass and strong edge displacements as well as displacements of ductile systems with strength eccentricity larger than stiffness eccentricity. Some modifications are proposed to get the MRSA yields a conservative estimation of displacement demands for all cases.
As a global phenomenon, many cities are undergoing urban renewal to accommodate rapid growth in urban population. However, urban renewal can struggle to balance social, economic, and environmental outcomes, whereby economic outcomes are often primarily considered by developers. This has important implications for urban forests, which have previously been shown to be negatively affected by development activities. Urban forests serve the purpose of providing ecosystem services and thus are beneficial to human wellbeing. Better understanding the effect of urban renewal on city trees may help improve urban forest outcomes via effective management and policy strategies, thereby maximising ecosystem service provision and human wellbeing. Though the relationship between certain aspects of development and urban forests has received consideration in previous literature, little research has focused on how the complete property redevelopment cycle affects urban forest dynamics over time. This research provides an opportunity to gain a comprehensive understanding of the effect of residential property redevelopment on urban forest dynamics, at a range of spatial scales, in Christchurch, New Zealand following a series of major earthquakes which occurred in 2010 – 2011. One consequence of the earthquakes is the redevelopment of thousands of properties over a relatively short time-frame. The research quantifies changes in canopy cover city-wide, as well as, tree removal, retention, and planting on individual residential properties. Moreover, the research identifies the underlying reasons for these dynamics, by exploring the roles of socio-economic and demographic factors, the spatial relationships between trees and other infrastructure, and finally, the attitudes of residential property owners. To quantify the effect of property redevelopment on canopy cover change in Christchurch, this research delineated tree canopy cover city-wide in 2011 and again in 2015. An object-based image analysis (OBIA) technique was applied to aerial imagery and LiDAR data acquired at both time steps, in order to estimate city-wide canopy cover for 2011 and 2015. Changes in tree canopy cover between 2011 and 2015 were then spatially quantified. Tree canopy cover change was also calculated for all meshblocks (a relatively fine-scale geographic boundary) in Christchurch. The results show a relatively small magnitude of tree canopy cover loss, city-wide, from 10.8% to 10.3% between 2011 and 2015, but a statistically significant change in mean tree canopy cover across all the meshblocks. Tree canopy cover losses were more likely to occur in meshblocks containing properties that underwent a complete redevelopment cycle, but the loss was insensitive to the density of redevelopment within meshblocks. To explore property-scale individual tree dynamics, a mixed-methods approach was used, combining questionnaire data and remote sensing analysis. A mail-based questionnaire was delivered to residential properties to collect resident and household data; 450 residential properties (321 redeveloped, 129 non- redeveloped) returned valid questionnaires and were identified as analysis subjects. Subsequently, 2,422 tree removals and 4,544 tree retentions were identified within the 450 properties; this was done by manually delineating individual tree crowns, based on aerial imagery and LiDAR data, and visually comparing the presence or absence of these trees between 2011 and 2015. The tree removal rate on redeveloped properties (44.0%) was over three times greater than on non-redeveloped properties (13.5%) and the average canopy cover loss on redeveloped properties (52.2%) was significantly greater than on non-redeveloped properties (18.8%). A classification tree (CT) analysis was used to model individual tree dynamics (i.e. tree removal, tree retention) and candidate explanatory variables (i.e. resident and household, economic, land cover, and spatial variables). The results indicate that the model including land cover, spatial, and economic variables had the best predicting ability for individual tree dynamics (accuracy = 73.4%). Relatively small trees were more likely to be removed, while trees with large crowns were more likely to be retained. Trees were most likely to be removed from redeveloped properties with capital values lower than NZ$1,060,000 if they were within 1.4 m of the boundary of a redeveloped building. Conversely, trees were most likely to be retained if they were on a property that was not redeveloped. The analysis suggested that the resident and household factors included as potential explanatory variables did not influence tree removal or retention. To conduct a further exploration of the relationship between resident attitudes and actions towards trees on redeveloped versus non-redeveloped properties, this research also asked the landowners from the 450 properties that returned mail questionnaires to indicate their attitudes towards tree management (i.e. tree removal, tree retention, and tree planting) on their properties. The results show that residents from redeveloped properties were more likely to remove and/or plant trees, while residents from non- redeveloped properties were more likely to retain existing trees. A principal component analysis (PCA) was used to explore resident attitudes towards tree management. The results of the PCA show that residents identified ecosystem disservices (e.g. leaf litter, root damage to infrastructure) as common reasons for tree removal; however, they also noted ecosystem services as important reasons for both tree planting and tree retention on their properties. Moreover, the reasons for tree removal and tree planting varied based on whether residents’ property had been redeveloped. Most tree removal occurred on redeveloped properties because trees were in conflict with redevelopment, but occurred on non- redeveloped properties because of perceived poor tree health. Residents from redeveloped properties were more likely to plant trees due to being aesthetically pleasing or to replace trees removed during redevelopment. Overall, this research adds to, and complements, the existing literature on the effects of residential property redevelopment on urban forest dynamics. The findings of this research provide empirical support for developing specific legislation or policies about urban forest management during residential property redevelopment. The results also imply that urban foresters should enhance public education on the ecosystem services provided by urban forests and thus minimise the potential for tree removal when undertaking property redevelopment.
The previously unknown Greendale Fault was buried beneath the Canterbury Plains and ruptured in the September 4th 2010 moment magnitude (Mw) 7.1 Darfield Earthquake. The Darfield Earthquake and subsequent Mw 6 or greater events that caused damage to Christchurch highlight the importance of unmapped faults near urban areas. This thesis examines the morphology, age and origin of the Canterbury Plains together with the paleoseismology and surface-rupture displacement distributions of the Greendale Fault. It offers new insights into the surface-rupture characteristics, paleoseismology and recurrence interval of the Greendale Fault and related structures involved in the 2010 Darfield Earthquake. To help constrain the timing of the penultimate event on the Greendale Fault the origin and age of the faulted glacial outwash deposits have been examined using sedimentological analysis of gravels and optically stimulated luminescence (OSL) dating combined with analysis of GPS and LiDAR survey data. OSL ages from this and other studies, and the analysis of surface paleochannel morphology and subsurface gravel deposits indicate distinct episodes of glacial outwash activity across the Canterbury Plains, at ~20 to 24 and ~28 to 33 kyr separated by a hiatus in sedimentation possibly indicating an interstadial period. These data suggest multiple glacial periods between ~18 and 35 kyr which may have occurred throughout the Canterbury region and wider New Zealand. A new model for the Waimakariri Fan is proposed where aggradation is mainly achieved during episodic sheet flooding with the primary river channel location remaining approximately fixed. The timing, recurrence interval and displacements of the penultimate surface-rupturing earthquake on the Greendale Fault have been constrained by trenching the scarp produced in 2010 at two locations. These excavations reveal a doubling of the magnitude of surface displacement at depths of 2-4 m. Aided by OSL ages of sand lenses in the gravel deposits, this factor-of-two increase is interpreted to indicate that in the central section of the Greendale Fault the penultimate surface-rupturing event occurred between ca. 20 and 30 kyr ago. The Greendale Fault remained undetected prior to the Darfield earthquake because the penultimate fault scarp was eroded and buried during Late Pleistocene alluvial activity. The Darfield earthquake rupture terminated against the Hororata Anticline Fault (HAF) in the west and resulted in up to 400 mm of uplift on the Hororata Anticline immediately above the HAF. Folding in 2010 is compared to Quaternary and younger deformation across the anticline recorded by a seismic reflection line, GPS-measured topographic profiles along fluvial surfaces, and river channel sinuosity and morphology. It is concluded that the HAF can rupture during earthquakes dissimilar to the 2010 event that may not be triggered by slip on the Greendale Fault. Like the Greendale Fault geomorphic analyses provide no evidence for rupture of the HAF in the last 18 kyr, with the average recurrence interval for the late Quaternary inferred to be at least ~10 kyr. Surface rupture of the Greendale Fault during the Darfield Earthquake produced one of the most accessible and best documented active fault displacement and geometry datasets in the world. Surface rupture fracture patterns and displacements along the fault were measured with high precision using real time kinematic (RTK) GPS, tape and compass, airborne light detection and ranging (LiDAR), and aerial photos. This allowed for detailed analysis of the cumulative strike-slip displacement across the fault zone, displacement gradient (ground shear strain) and the type of displacement (i.e. faulting or folding). These strain profiles confirm that the rupture zone is generally wide (~30 to ~300 metres) with >50% of displacement (often 70-80%) accommodated by ground flexure rather than discrete fault slip and ground cracking. The greatest fault-zone widths and highest proportions of folding are observed at fault stepovers.
We present the initial findings from a study of adaptive resilience of lifelines organisations providing essential infrastructure services, in Christchurch, New Zealand following the earthquakes of 2010-2011. Qualitative empirical data was collected from 200 individuals in 11 organisations. Analysis using a grounded theory method identified four major factors that aid organisational response, recovery and renewal following major disruptive events. Our data suggest that quality of top and middle-level leadership, quality of external linkages, level of internal collaboration, ability to learn from experience, and staff well-being and engagement influence adaptive resilience. Our data also suggest that adaptive resilience is a process or capacity, not an outcome and that it is contextual. Post-disaster capacity/resources and post-disaster environment influence the nature of adaptive resilience.
Disasters are rare events with major consequences; yet comparatively little is known about managing employee needs in disaster situations. Based on case studies of four organisations following the devastating earthquakes of 2010 - 2011 in Christchurch, New Zealand, this paper presents a framework using redefined notions of employee needs and expectations, and charting the ways in which these influence organisational recovery and performance. Analysis of in-depth interview data from 47 respondents in four organisations highlighted the evolving nature of employee needs and the crucial role of middle management leadership in mitigating the effects of disasters. The findings have counterintuitive implications for human resource functions in a disaster, suggesting that organisational justice forms a central framework for managing organisational responses to support and engage employees for promoting business recovery.
The Canterbury earthquake and aftershock sequence in New Zealand during 2010-2011 subjected the city’s structures to a significant accumulated cyclic demand and raised significant questions regarding the low-cycle fatigue demands imposed upon the structures. There is a significant challenge to quantify the level of cumulative demand imposed on structures and to assess the percentage of a structure's fatigue life that has been consumed as a result of this earthquake sequence. It is important to be able to quantify the cumulative demand to determine how a building will perform in a subsequent large earthquake and inform repair and re-occupancy decisions. This paper investigates the cumulative fatigue demand for a structure located within the Christchurch Central Business District (CBD). Time history analysis and equivalent cycle counting methods are applied across the Canterbury earthquake sequence, using key events from September 4th 2010 and February 22nd , 2011 main shocks. The estimate of the cumulative fatigue demand is then compared to the expected capacity of a case study reinforced concrete bridge pier, to undertake a structure-specific fatigue assessment. The analysis is undertaken to approximate the portion of the structural fatigue capacity that has been consumed, and how much residual capacity remains. Results are assessed for recordings at the four Christchurch central city strong motion recording sites installed by the GeoNet programme, to provide an estimate of variation in results. The computed cyclic demand results are compared to code-based design methods and as assessment of the inelastic displacement demand of the reinforcing steel. Results are also presented in a fragility context where a de minimis (inconsequential), irreparable damage and full fatigue fracture are defined to provide a probabilistic assessment of the fatigue damage incurred. This methodology can provide input into the overall assessment of fatigue demands and residual capacity.
The initial goal of this research was to explore how SME business models change in response to a crisis. Keeping this in mind, the business model canvas (Osterwalder & Pigneur, 2010) was used as a tool to analyse SME business models in the Canterbury region of New Zealand. The purpose was to evaluate the changes SMEs instituted in their business models after being hit by a series of earthquakes in 2010 and 2011. The idea was to conduct interviews with business owners and analyse them using grounded theory methods. As this method is iterative and requires simultaneous data collection and analysis, a tentative model was proposed after first phase of the data collection and analysis. However, as a result of this process, it became apparent that owner-specific characteristics, action orientation and networks were more prominent in the data than business model elements. Although the SMEs in this study experienced several operational changes in their business models, such as a change of location, modifications to their payment terms or expanded/restricted target markets, the suggested framework highlights how owner-specific attributes ensured the recovery of their businesses. After the initial framework was suggested, subsequent interviews were conducted to test, verify, and modify the tentative model. Three aspects of business recovery emerged: (a) cognitive coping – the business owner’s mind-set and motive; (b) adaptive coping – the ability of business owner to take corrective actions; and (c) social capital – the social network of a business owner, including formal and informal connections and their significance. Three distinct groups were identified; self-sufficient SMEs, socially-based SMEs and surviving SMEs. This thesis proposes a grounded theory of business recovery for SMEs following a disaster. Cognitive coping and social capital enabled the owners to take actions, which eventually led to the desired outcomes for the businesses.
The Covid-19 pandemic has brought to the foreground the importance of social connectedness for wellbeing, at the individual, community and societal level. Within the context of the local community, pro-connection facilities are fundamental to foster community development, resilience and public health. Through identifying the gap in social connectedness literature for Māori, this has created space for new opportunities and to reflect on what is already occurring in Ōtautahi. It is well documented that Māori experience unequal societal impacts across all health outcomes. Therefore, narrowing the inequities between Māori and non-Māori across a spectrum of dimensions is a priority. Evaluating the #WellconnectedNZ project, which explores the intersections between social connection and wellbeing is one way to trigger these conversations. This was achieved by curating a dissimilar set of community pro-connection facilities and organizing them into a Geographic Information System (GIS). Which firstly involved, the collecting and processing of raw data, followed by spatial analysis through creating maps, this highlighted the alignment between the distribution of places, population and social data. Secondly, statistical analysis focusing on the relationship between deprivation and accessibility. Finally, semi-structured interviews providing perceptions of community experience. This study describes findings following a kaupapa Māori research approach. Results demonstrated that, in general some meshblocks in Ōtautahi benefit from a high level of accessibility to pro-connection facilities; but with an urban-rural gradient (as is expected, further from the central business district (CBD) are less facilities). Additionally, more-deprived meshblocks in the Southern and Eastern suburbs of Christchurch have poorer accessibility, suggesting underlying social and spatial inequalities, likely exacerbated by Covid-19 and the Christchurch earthquakes. In this context, it is timely to (re)consider pro-connection places and their role in the development of social infrastructure for connected communities, in the community facility planning space. ‘We are all interwoven, we just need to make better connections’.
Probabilistic Structural Fire Engineering (PSFE) has been introduced to overcome the limitations of current conventional approaches used for the design of fire-exposed structures. Current structural fire design investigates worst-case fire scenarios and include multiple thermal and structural analyses. PSFE permits buildings to be designed to a level of life safety or economic loss that may occur in future fire events with the help of a probabilistic approach. This thesis presents modifications to the adoption of a Performance-Based Earthquake Engineering (PBEE) framework in Probabilistic Structural Fire Engineering (PSFE). The probabilistic approach runs through a series of interrelationships between different variables, and successive convolution integrals of these interrelationships result in probabilities of different measures. The process starts with the definition of a fire severity measure (FSM), which best relates fire hazard intensity with structural response. It is identified by satisfying efficiency and sufficiency criteria as described by the PBEE framework. The relationship between a fire hazard and corresponding structural response is established by analysis methods. One method that has been used to quantify this relationship in PSFE is Incremental Fire Analysis (IFA). The existing IFA approach produces unrealistic fire scenarios, as fire profiles may be scaled to wide ranges of fire severity levels, which may not physically represent any real fires. Two new techniques are introduced in this thesis to limit extensive scaling. In order to obtain an annual rate of exceedance of fire hazard and structural response for an office building, an occurrence model and an attenuation model for office fires are generated for both Christchurch city and New Zealand. The results show that Christchurch city is 15% less likely to experience fires that have the potential to cause structural failures in comparison to all of New Zealand. In establishing better predictive relationships between fires and structural response, cumulative incident radiation (a fire hazard property) is found to be the most appropriate fire severity measure. This research brings together existing research on various sources of uncertainty in probabilistic structural fire engineering, such as elements affecting post-flashover fire development factors (fuel load, ventilation, surface lining and compartment geometry), fire models, analysis methods and structural reliability. Epistemic uncertainty and aleatory uncertainty are investigated in the thesis by examining the uncertainty associated with modelling and the factors that influence post-flashover development of fires. A survey of 12 buildings in Christchurch in combination with recent surveys in New Zealand produced new statistical data on post-flashover development factors in office buildings in New Zealand. The effects of these parameters on temperature-time profiles are evaluated. The effects of epistemic uncertainty due to fire models in the estimation of structural response is also calculated. Parametric fires are found to have large uncertainty in the prediction of post-flashover fires, while the BFD curves have large uncertainties in prediction of structural response. These uncertainties need to be incorporated into failure probability calculations. Uncertainty in structural modelling shows that the choices that are made during modelling have a large influence on realistic predictions of structural response.
This dissertation addresses a diverse range of applied aspects in ground motion simulation validation via the response of complex structures. In particular, the following topics are addressed: (i) the investigation of similarity between recorded and simulated ground motions using code-based 3D irregular structural response analysis, (ii) the development of a framework for ground motion simulations validation to identify the cause of differences between paired observed and simulated dataset, and (iii) the illustration of the process of using simulations for seismic performance-based assessment. The application of simulated ground motions is evaluated for utilisation in engineering practice by considering responses of 3D irregular structures. Validation is performed in a code-based context when the NZS1170.5 (NZS1170.5:2004, 2004) provisions are followed for response history analysis. Two real buildings designed by engineers and physically constructed in Christchurch before the 2010-2011 Canterbury earthquake sequence are considered. The responses are compared when the buildings are subjected to 40 scaled recorded and their subsequent simulated ground motions selected from 22 February 2011 Christchurch. The similarity of recorded and simulated responses is examined using statistical methods such as bootstrapping and hypothesis testing to determine whether the differences are statistically significant. The findings demonstrate the applicability of simulated ground motion when the code-based approach is followed in response history analysis. A conceptual framework is developed to link the differences between the structural response subjected to simulated and recorded ground motions to the differences in their corresponding intensity measures. This framework allows the variability to be partitioned into the proportion that can be “explained” by the differences in ground motion intensity measures and the remaining “unexplained” variability that can be attributed to different complexities such as dynamic phasing of multi-mode response, nonlinearity, and torsion. The application of this framework is examined through a hierarchy of structures reflecting a range of complexity from single-degree-of-freedom to 3D multi-degree-of-freedom systems with different materials, dynamic properties, and structural systems. The study results suggest the areas that ground motion simulation should focus on to improve simulations by prioritising the ground motion intensity measures that most clearly account for the discrepancies in simple to complex structural responses. Three approaches are presented to consider recorded or simulated ground motions within the seismic performance-based assessment framework. Considering the applications of ground motions in hazard and response history analyses, different pathways in utilising ground motions in both areas are explored. Recorded ground motions are drawn from a global database (i.e., NGA-West2 Ancheta et al., 2014). The NZ CyberShake dataset is used to obtain simulations. Advanced ground motion selection techniques (i.e., generalized conditional intensity measure, GCIM) are used for ground motion selection at a few intensity levels. The comparison is performed by investigating the response of an example structure (i.e., 12-storey reinforced concrete special moment frame) located in South Island, NZ. Results are compared and contrasted in terms of hazard, groundmotion selection, structural responses, demand hazard, and collapse risk, then, the probable reasons for differences are discussed. The findings from this study highlight the present opportunities and shortcomings in using simulations in risk assessment. i
This paper develops representative ground motion ensembles for several major earthquake scenarios in New Zealand. Cases considered include representative ground motions for the occurrence of Alpine, Hope, and Porters Pass earthquakes in Christchurch, and the occurrence of Wellington, Wairarapa, and Ohariu, fault ruptures in Wellington. Challenges in the development of ground motion ensembles for subduction zone earthquakes are also highlighted. The ground motions are selected based on the generalized conditional intensity measure (GCIM) approach, ensuring that the ground motion ensembles represent both the mean, and distribution of ground motion intensity which such scenarios could impose. These scenario-based ground motion sets can be used to complement ground motions which are often selected in conjunction with probabilistic seismic hazard analysis, in order to understand the performance of structures for the question “what if this fault ruptures?”
SeisFinder is an open-source web service developed by QuakeCoRE and the University of Canterbury, focused on enabling the extraction of output data from computationally intensive earthquake resilience calculations. Currently, SeisFinder allows users to select historical or future events and retrieve ground motion simulation outputs for requested geographical locations. This data can be used as input for other resilience calculations, such as dynamic response history analysis. SeisFinder was developed using Django, a high-level python web framework, and uses a postgreSQL database. Because our large-scale computationally-intensive numerical ground motion simulations produce big data, the actual data is stored in file systems, while the metadata is stored in the database.
Overview of SeisFinder SeisFinder is an open-source web service developed by QuakeCoRE and the University of Canterbury, focused on enabling the extraction of output data from computationally intensive earthquake resilience calculations. Currently, SeisFinder allows users to select historical or future events and retrieve ground motion simulation outputs for requested geographical locations. This data can be used as input for other resilience calculations, such as dynamic response history analysis. SeisFinder was developed using Django, a high-level python web framework, and uses a postgreSQL database. Because our large-scale computationally-intensive numerical ground motion simulations produce big data, the actual data is stored in file systems, while the metadata is stored in the database. The basic SeisFinder architecture is shown in Figure 1.
The purpose of this thesis is to evaluate the seismic response of the UC Physics Building based on recorded ground motions during the Canterbury earthquakes, and to use the recorded response to evaluate the efficacy of various conventional structural analysis modelling assumptions. The recorded instrument data is examined and analysed to determine how the UC Physics Building performed during the earthquake-induced ground motions. Ten of the largest earthquake events from the 2010-11 Canterbury earthquake sequence are selected in order to understand the seismic response under various levels of demand. Peak response amplitude values are found which characterise the demand from each event. Spectral analysis techniques are utilised to find the natural periods of the structure in each orthogonal direction. Significant torsional and rocking responses are also identified from the recorded ground motions. In addition, the observed building response is used to scrutinise the adequacy of NZ design code prescriptions for fundamental period, response spectra, floor acceleration and effective member stiffness. The efficacy of conventional numerical modelling assumptions for representing the UC Physics Building are examined using the observed building response. The numerical models comprise of the following: a one dimensional multi degree of freedom model, a two dimensional model along each axis of the building and a three dimensional model. Both moderate and strong ground motion records are used to examine the response and subsequently clarify the importance of linear and non-linear responses and the inclusion of base flexibility. The effects of soil-structure interaction are found to be significant in the transverse direction but not the longitudinal direction. Non-linear models predict minor in-elastic behaviour in both directions during the 4 September 2010 Mw 7.1 Darfield earthquake. The observed torsional response is found to be accurately captured by the three dimensional model by considering the interaction between the UC Physics Building and the adjacent structure. With the inclusion of adequate numerical modelling assumptions, the structural response is able to be predicted to within 10% for the majority of the earthquake events considered.
The 2010-2011 Canterbury Earthquakes brought devastation to the city of Christchurch and has irrevocably affected the lives of the city’s residents. Years after the conclusion of these earthquakes, Christchurch and its residents are well on the path to recovery. Crime has proven an ongoing topic of discussion throughout this period, with news reports of increased burglary and arson in areas left largely abandoned by earthquake damage, and a rise in violent crime in suburban areas of Christchurch. Following the body of research that has considered the reaction of crime to natural disasters, this research has sought to comprehensively examine and understand the effects that the Canterbury Earthquakes had on crime. Examining Christchurch-wide offending, crime rates fell over the study period (July 2008 to June 2013), with the exception of domestic violence. Aside from a momentary increase in burglary in the days immediately following the Christchurch Earthquake, crime rates (as of 2013) have remained largely below pre-earthquake levels. Using Dual Kernel Density Estimation Analysis, a distinct spatial change in pre-earthquake crime hotspots was observed. These changes included an enormous decrease in central city offences, a rise in burglary in the eastern suburbs, and an increase in assault in areas outside of the central city. Logistic regression analysis, using a time-compensated dependent variable, identified a number of statistically-significant relationships between per CAU crime rate change and factors measuring socio-demographic characteristics, community cohesion, and the severity of disaster effects. The significance of these findings was discussed using elements of Social Disorganisation Theory, Routine Activity Theory, and Strain Theory. Consistent with past findings, social order was largely maintained following the Canterbury Earthquakes, with suggestion that increased collective efficacy and therapeutic communities had a negative influence on crime in the post-earthquake period. Areas of increased burglary and assault were associated with large population decreases, suggesting a link with the dissolution of communities and the removal of their inherent informal guardianship. Though observed, the increase in domestic violence was not associated with most neighbourhood-level variables. Trends in crime after the Canterbury Earthquakes were largely consistent with past research, and the media’s portrayal.
Disasters are rare events with major consequences; yet comparatively little is known about managing employee needs in disaster situations. Based on case studies of four organisations following the devastating earthquakes of 2010 - 2011 in Christchurch, New Zealand, this paper presents a framework using redefined notions of employee needs and expectations, and charting the ways in which these influence organisational recovery and performance. Analysis of in-depth interview data from 47 respondents in four organisations highlighted the evolving nature of employee needs and the crucial role of middle management leadership in mitigating the effects of disasters. The findings have counterintuitive implications for human resource functions in a disaster, suggesting that organisational justice forms a central framework for managing organisational responses to support and engage employees for promoting business recovery.
