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.
Following the 22nd February 2011, Mw 6.2 earthquake located along a previously unknown fault beneath the Port Hills of Christchurch, surface cracking was identified in contour parallel locations within fill material at Quarry Road on the lower slopes of Mount Pleasant. GNS Science, in the role of advisor to the Christchurch City Council, concluded that these cracks were a part of a potential rotational mass movement (named zone 11A) within the fill and airfall loess material present. However, a lack of field evidence for slope instability and an absence of laboratory geotechnical data on which slope stability analysis was based, suggested this conclusion is potentially incorrect. It was hypothesised that ground cracking was in fact due to earthquake shaking, and not mass movement within the slope, thus forming the basis of this study. Three soil units were identified during surface and subsurface investigations at Quarry Road: fill derived from quarry operations in the adjacent St. Andrews Quarry (between 1893 and 1913), a buried topsoil, and underlying in-situ airfall loess. The fill material was identified by the presence of organic-rich topsoil “clods” that were irregular in both size (∼10 – 200 mm) and shape, with variable thicknesses of 1 – 10 m. Maximum thickness, as indicated by drill holes and geophysical survey lines, was identified below 6 Quarry Road and 7 The Brae where it is thought to infill a pre-existing gully formed in the underlying airfall loess. Bearing strength of the fill consistently exceeded 300 kPa ultimate below ∼500 mm depth. The buried topsoil was 200 – 300 mm thick, and normally displayed a lower bearing strength when encountered, but not below 300 kPa ultimate (3 – 11 blows per 100mm or ≥100 kPa allowable). In-situ airfall loess stood vertically in outcrop due to its characteristic high dry strength and also showed Scala penetrometer values of 6 – 20+ blows per 100 mm (450 – ≥1000 kPa ultimate). All soils were described as being moist to dry during subsurface investigations, with no groundwater table identified during any investigation into volcanic bedrock. In-situ moisture contents were established using bulk disturbed samples from hand augers and test pitting. Average moisture contents were low at 9% within the fill, 11 % within the buried topsoil, and 8% within the airfall loess: all were below the associated average plastic limit of 17, 15, and 16, respectively, determined during Atterberg limit analysis. Particle size distributions, identified using the sieve and pipette method, were similar between the three soil units with 11 – 20 % clay, 62 – 78 % silt, and 11 – 20 % fine sand. Using these results and the NZGS soil classification, the loess derived fill and in-situ airfall loess are termed SILT with some clay and sand, and the buried topsoil is SILT with minor clay and sand. Dispersivity of the units was found using the Emerson crumb test, which established that the fill can be non- to completely dispersive (score 0 – 4). The buried topsoil was always non-dispersive (score 0), and airfall loess completely dispersive (score 4). Values for cohesion (c) and internal friction angle (φ) of the three soil units were established using the direct shear box at field moisture contents. Results showed all soil units had high shear strengths at the moisture contents tested (c = 18 – 24 kPa and φ = 42 – 50°), with samples behaving in a brittle fashion. Moisture content was artificially increased to 16% within the buried topsoil, which reduced the shear strength (c = 10 kPa, φ = 18°) and allowed it to behave plastically. Observational information indicating stability at Quarry Road included: shallow, discontinuous, cracks that do not display vertical offset; no scarp features or compressional zones typical of landsliding; no tilted or deformed structures; no movement in inclinometers; no basal shear zone identified in logged core to 20 m depth; low field moisture contents; no groundwater table; and high soil strength using Scala penetrometers. Limit equilibrium analysis of the slope was conducted using Rocscience software Slide 5.0 to verify the slope stability identified by observational methods. Friction, cohesion, and density values determined during laboratory were input into the two slope models investigated. Results gave minimum static factor of safety values for translational (along buried topsoil) and rotational (in the fill) slides of 2.4 – 4.2. Sensitivity of the slope to reduced shear strength parameters was analysed using c = 10 kPa and φ = 18° for the translational buried topsoil plane, and a cohesion of 0 kPa within the fill for the rotational plane. The only situation that gave a factor of safety <1.0 was in nonengineered fill at 0.5 m depth. Pseudostatic analysis based on previous peak ground acceleration (PGA) values for the Canterbury Earthquake Sequence, and predicted PGAs for future Alpine Fault and Hope Fault earthquakes established minimum factor of safety values between 1.2 and 3.3. Yield acceleration PGAs were computed to be between 0.8g and 1.6g. Based on all information gathered, the cracking at Quarry Road is considered to be shallow deformation in response to earthquake shaking, and not due to deep-seated landsliding. It is recommended that the currently bare site be managed by smoothing the land, installing contour drainage, and bioremediation of the surface soils to reduce surface water infiltration and runoff. Extensive earthworks, including removal of the fill, are considered unnecessary. Any future replacement of housing would be subject to site-specific investigations, and careful foundation design based on those results.
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.
Ongoing climate change triggers increasing temperature and more frequent extreme events which could limit optimal performance of haliotids, affect their physiology and biochemistry as well as influencing their population structure. Haliotids are a valuable nearshore fishery in a number of countries and many are showing a collapse of stocks because of overexploitation, environmental changes, loss of habitat, and disease. The haliotid in New Zealand commonly referred to as the blackfoot pāua (Haliotis iris) contribute a large and critical cultural, recreational and economic resource. Little was known about pāua responses to increasing temperature and acute environmental factors, as well as information about population size structure in Kaikoura after the earthquake 2016 and in Banks Peninsula. The aims of this study were to investigate the effects of temperature on scope for growth (SfG); physiological and biochemical responses of pāua subjected to different combined stressors including acute temperature, acute salinity and progressive hypoxia; and describe population size structure and shell morphology in different environments in Kaikoura and Banks Peninsula. The main findings of the present study found that population size structures of pāua were site-specific, and the shell length and shell height ratio of 3.25 could distinguish between stunted and non-stunted populations. The study found that high water temperature resulted in a reduction in absorbed energy from food, an increase in respiration energy, and ammonia excretion energy. Surveys were conducted at six study sites around the Canterbury Region over three years in order to better understand the population size structure and shell morphology of pāua. The findings found that the population size structure at 6 sites differed. Both juveniles and adults were found in intertidal areas at five sites. However, at Cape Three Points, pāua were found only in subtidal zones. One of the sites, Little Port Cooper, had a stunted population where only two pāua reached 125 mm in length over three years. In addition, most pāua in Little Port Cooper and Cape Three Points were adults, while Seal Reef had mostly juveniles. Wakatu Quay and Omihi had a full size range of pāua. Oaro population was dominated with juveniles and sub-adults. Recruitment and growth of pāua were successful after the earthquake in 2016. Research into pāua shell morphologies also determined that shell dimensions differed between sites. The relationships of shell length to shell width were linear and the relationship of shell length to shell height was curvilinear. Interestingly, SL:SH ratio of 3.25 is able to be used to identify stunted and non-stunted populations for pāua larger than 90 mm in length. Little Port Cooper was a stunted population with mean SL:SH ratio being 3.16. In the laboratory, scope for growth of pāua was investigated at four different temperatures of 12oC, 15oC, 18oC and 21oC over four weeks’ acclimation. The current study has found that SfG of pāua highly depended on temperature. Absorbed energy and respiration energy accounted for the highest proportion of the SfG of pāua. The respiration energy of pāua accounted for approximately 36%, 40%, 49% and 69% of the absorbed energy at 12°C, 15°C, 18°C and 21°C, respectively. The pāua at all acclimation temperatures had a positive scope for growth. The study suggested that the SfG was highest at 15°C, while the value at 21°C was the lowest. However, SfG at 18°C and 21°C decreased after 14 days of acclimation. Because of maintaining almost unchanged oxygen consumption over four weeks’ acclimation, pāua showed their poor abilities to acclimate to an increase in temperature. Therefore, they may be more vulnerable in future warming scenarios. The physiological and biochemical responses of pāua toward different combined stressors included three experiments. In terms of the acute temperature experiment, pāua were acclimated at 12oC, 15oC, 18oC or 21oC for two weeks before stepwise exposure to four temperatures of 12oC, 15oC, 18oC and 21oC every 4 hours. The acute salinity change, pāua were acclimated at 12oC, 15oC or 18oC over two weeks. Pāua were then exposed to a stepwise decrease of salinity of 2‰ every two hours from 34 – 22‰. Regarding the declining oxygen level, pāua were acclimated at 15 oC or 18oC for two weeks before exposure to one of four temperatures at 12oC, 15oC, 18oC or 21oC in one hour. After that acute progressive hypoxia was studied in closed respirometers for around six hours. The findings showed that there were interactions between combined stressors, affecting physiology of pāua (metabolism and heart rate). This suggests that environmental factors do not have a separate effect, but they also have interactions that enhance negative effects on pāua. Also, both oxygen uptake and heart rate responded quickly to temperature change and increased with rising temperature. On the other hand, oxygen uptake and heart rate decreased with reducing salinity and progressive hypoxia (before critical oxygen tension - Pcrit). Pcrit over four acute temperature exposures, ranged between 30.2 and 80.0 mmHg, depending on the exposure temperature. Acclimation temperature, combined with acute temperature, salinity or hypoxia stress affected the biochemistry of pāua. Pāua are osmoconformers so decreased salinity resulted in reducing haemolymph ionic concentration and increasing body volume. They were hypo-ionic with respect to sodium and potassium over the salinity ranges of 34 - 22‰. Haemocyanin accounts for a large pecentage of haemolymph protein, so trends of protein followed haemocyanin. Pāua tended to store oxygen in haemocyanin under extreme salinity stress at 22‰ and extreme hypoxia around 10 mmHg, rather than in oxygen transport. In conclusion, pāua at different sites had different population structures and morphologies. Pāua are sensitive to environmental stressors. They consumed more oxygen at high temperatures because they do not have thermal acclimation capacity. They are also osmoconformers with haemolymph sodium and potassium decreasing with salinity medium. Under progressive hypoxia, pāua could regulate oxygen and heart rate until Pcrit depending on temperature. Acute environmental changes also disturbed haemolyph parameters. 12°C and 15°C could be in the range of optimal temperature with higher SfG and less stress when exposed to acute environmental changes. Meanwhile long term exposure to 21°C is likely to be outside of the optimal range for the pāua. With ongoing climate change, pāua populations are more vulnerable so conservation is necessary. The research contributes to improving fishery management, providing insights into different environmental stressors affecting the energy demand and physiological and biochemical responses of pāua. It also allow to predicting the growth patterns and responses of pāua to adapt to climate change.