February floods in the lower North Island, 2004: Catastrophe – causes and consequences

Abstract: The February 2004 Manawatu floods in New Zealand were the result of a naturally occurring, although unusual, storm. Up to 300 mm of rain fell on the already saturated ground of the lower North Island over two days, generating substantial and rapid runoff from catchment slopes. Rivers rose quickly, inundating unprotected farmland and properties and in places breaching stopbanks. There was widespread slope failure in the hill country of the lower North Island, affecting an area of ca. 7500 km². Slopes under scrub, plantation forest and native bush were not as badly affected as those under pasture, where slopes typically failed by shallow translational landsliding. Flooding caused catastrophic channel change in a number of small to medium sized channel systems in the upland fringes. Whilst the occurrence of landsliding and channel changes during an extreme event such as this is natural, the intensity of both landsliding and channel erosion was exacerbated by human activity within the catchments.     

Flexural and gravity modelling of the Mérida Andes and Barinas–Apure Basin, Western Venezuela

The kinematic evolution of the Barinas–Apure Basin and the southern Mérida Andes from Lower Miocene to the Present is numerically modelled using flexural isostatic theory and geophysical and geological data. Two published regional transects are used to build up a reference section, which is then used to constrain important parameters (e.g. shortenings and sedimentary thicknesses) for the flexural modelling. To control the location of the main fault system in the flexural model earthquake information is also used. The estimated flexural elastic thickness of the South American lithosphere beneath the Barinas–Apure Basin and the Mérida Andes Range is 25 km. The value for the final total shortening is 60 km. The flexural isostatic model shows that the Andean uplift has caused the South American lithosphere subsidence and the development of the Barinas–Apure Basin.In addition, gravity modelling was used to understand deep crustal features that could not be predicted by flexural theory. Consequently, the best-fit flexural model is used to build a gravity model across the Mérida Andes and the Barinas–Apure Basin preserving the best-controlled structural features from the flexural modelling (e.g. basin wavelength and depth) and slightly changing the main bodies density values and deep crustal structures. The final gravity model is intended to be representative of the major features affecting the gravity field in the study area. The predicted morphology in the lower crustal level of the final gravity model favours the hypothesis of a present delamination or megathrust of the Maracaibo crust over the South American Shield. This process would use the Conrad discontinuity as a main detachment surface within an incipient NW dipping continental subduction.     

Soil liquefaction during the Arequipa Mw 8.4, June 23, 2001 earthquake, southern coastal Peru

The Arequipa June 23, 2001, earthquake with a moment magnitude of Mw 8.4 struck southern Peru, northern Chile and western Bolivia. This shallow (29 km deep) interplate event, occurring in the coupled zone of the Nazca subduction next to the southeast of the subducting Nazca ridge, triggered very localized but widely outspread soil liquefaction. Although sand blows and lateral spreading of river banks and road bridge abutments were observed 390 km away from the epicenter in the southeast direction (nearing the town of Tacna, close to the Chile border), liquefaction features were only observed in major river valleys and delta and coastal plains in the meizoseismal area. This was strongly controlled by the aridity along the coastal strip of Southern Peru. From the sand blow distribution along the coastal area, a first relationship of isolated sand blow diameter versus epicentral distance for a single event is ever proposed. The most significant outcome from this liquefaction field reconnaissance is that energy propagation during the main June 23, 2001, event is further supported by the distribution and size of the isolated sand blows in the meizoseismal area. The sand blows are larger to the southeast of the epicenter than its northwestern equivalents. This can be stated in other words as well. The area affected by liquefaction to the northwest is less spread out than to the southeast. Implications of these results in future paleoliquefaction investigations for earthquake magnitude and epicentral determinations are extremely important. In cases of highly asymmetrical distribution of liquefaction features such as this one, where rupture propagation tends to be mono-directional, it can be reliably determined an epicentral distance (between earthquake and liquefaction evidence) and an earthquake magnitude only if the largest sand blow is found. Therefore, magnitude estimation using this uneven liquefaction occurrence will surely lead to underrating if only the shortest side of the meizoseismal area is unluckily studied, which can eventually be the only part exhibiting liquefaction evidence, depending on the earthquake location and the distribution of liquefaction-prone environments.     

Kamafugitic diatremes: their textures and field relationships with examples from the Goiás alkaline province, Brazil

Kamafugitic rocks intruded the Precambrian basement and Phanerozoic sediments at the northeast border of the Paraná basin as part of the Late Cretaceous Goiás alkaline province (GAP). Plutonic complexes dominate the north of the province, whereas lavas and pyroclastic rocks prevail in the south. The central GAP is characterized by kamafugitic diatremes, which may crop out continuously for up to 850 m and consist of a central breccia body, surrounded and overlain by lava flows and crosscut by dykes. The breccias contain some special spheroidal juvenile fragments—namely, accretionary and armored lapilli, frozen droplets, spinning droplets, and wrapped fragments—whose textural and mineralogical aspects are described in detail. Irregularly shaped tuff pockets that occur within the breccias contain textures and structures similar to those of subaerial surge deposits and formed in confined, high gas to solid+liquid ratio domains in the conduit. Diatreme emplacement affected the country rock through thermal metamorphism, development of columnar jointing, and formation of peperite-like mixtures. There is no evidence of phreatomagmatic activity in the diatremes, and CO₂, rather than H₂O, seems to have been the major volatile component of the kamafugitic magmas. This finding implies that features such as accretionary lapilli and peperites are not exclusively associated with H₂O-dominated processes.     

Application of Distributed Temperature Sensing for coupled mapping of sedimentation processes and spatio‐temporal variability of groundwater discharge in soft‐bedded streams

The delineation of groundwater discharge areas based on Distributed Temperature Sensing (DTS) data of the streambed can be difficult in soft‐bedded streams where sedimentation and scouring processes constantly change the position of the fibre optic cable relative to the streambed. Deposition‐induced temperature anomalies resemble the signal of groundwater discharge while scouring will cause the cable to float in the water column and measure stream water temperatures. DTS applied in a looped layout with nine fibre optic cable rows in a 70 × 5 m section of a soft‐bedded stream made it possible to detect variability in streambed temperatures between October 2011 and January 2012. Detailed monthly streambed elevation surveys were carried out to monitor the position of the fibre optic cable relative to the streambed and to quantify the effect of sedimentation processes on streambed temperatures. Based on the simultaneous interpretation of streambed temperature and elevation data, a method is proposed to delineate potential high‐groundwater discharge areas and identify deposition‐induced temperature anomalies in soft‐bedded streams. Potential high‐discharge sites were detected using as metrics the daily minimum, maximum and mean streambed temperatures as well as the daily amplitude and standard deviation of temperatures. The identified potential high‐discharge areas were mostly located near the channel banks, also showing temporal variability because of the scouring and redistribution of streambed sediments, leading to the relocation of pool‐riffle sequences. This study also shows that sediment deposits of 0.1 m thickness already resulted in an increase in daily minimum streambed temperatures and decrease in daily amplitude and standard deviation. Scouring sites showed lower daily minimum streambed temperatures and higher daily amplitude and standard deviation compared with areas without sedimentation and scouring. As a limitation of the approach, groundwater discharge occurring at depositional and scouring areas cannot be identified by the metrics applied. Copyright © 2015 John Wiley & Sons, Ltd.     

Application of a confluence‐based sediment‐fingerprinting approach to a dynamic sedimentary catchment,New Zealand

Fine sediment is a dynamic component of the fluvial system, contributing to the physical form, chemistry and ecological health of a river. It is important to understand rates and patterns of sediment delivery, transport and deposition. Sediment fingerprinting is a means of directly determining sediment sources via their geochemical properties, but it faces challenges in discriminating sources within larger catchments. In this research, sediment fingerprinting was applied to major river confluences in the Manawatu catchment as a broad‐scale application to characterizing sub‐catchment sediment contributions for a sedimentary catchment dominated by agriculture. Stepwise discriminant function analysis and principal component analysis of bulk geochemical concentrations and geochemical indicators were used to investigate sub‐catchment geochemical signatures. Each confluence displayed a unique array of geochemical variables suited for discrimination. Geochemical variation in upstream sediment samples was likely a result of the varying geological source compositions. The Tiraumea sub‐catchment provided the dominant signature at the major confluence with the Upper Manawatu and Mangatainoka sub‐catchments. Subsequent downstream confluences are dominated by the upstream geochemical signatures from the main stem of Manawatu River. Variability in the downstream geochemical signature is likely due to incomplete mixing caused in part by channel configuration. Results from this exploratory investigation indicate that numerous geochemical elements have the ability to differentiate fine sediment sources using a broad‐scale confluence‐based approach and suggest there is enough geochemical variation throughout a large sedimentary catchment for a full sediment fingerprint model. Combining powerful statistical procedures with other geochemical analyses is critical to understanding the processes or spatial patterns responsible for sediment signature variation within this type of catchment. Copyright © 2015 John Wiley & Sons, Ltd.     

Assessment of dissolved nutrients dispersal derived from offshore fish-farm using nitrogen stable isotope ratios (δ15N) in macroalgal bioassays

In this study, the dispersal of wastes from offshore fish farms was evaluated by analyzing nitrogen stable isotope ratios (δ¹⁵N) in macroalgae incubated in the water column at sites located at an increasing distance from the fish cages. Bioassays were performed at three fish farms situated in separate localities with different nutritional conditions (Canary Islands, Murcia and Catalonia) and varying in size, species of fish reared and annual production. Macroalgal bioassays were carried out in two different directions (DI and DII) and they were replicated at each distance in order to evaluate the effect of small-scale variability on the spatial extent of fish farm wastes. The results obtained with δ¹⁵N contribute to a better understanding of the application of nitrogen stable isotopes ratios in macroalgae as an effective bioindicator for tracing the dispersion of offshore fish farm wastes, and demonstrate that fish farm wastes can be traced even over distances of some km from the pollution source. In the Canary Islands, the maximum distance obtained for detection of fish farm wastes was between 450 and 700 m. Of the three installations studied, Murcia presented the greatest distance for detection of fish farm waste influence, ranging from between 1550 and 2450 m, whilst in Catalonia this distance was less than 120 m. In Catalonia, the results were masked by the influence of other sources of nitrogen, and thus fish farm wastes were detected at more reduced distances than expected. These results confirm that fish farm wastes can be traced using the nitrogen stable isotope ratios of macroalgae and that this method can also be useful for identifying areas of potential risk to some sensitive ecosystems, and as an early signal that changes in the community structure might occur.     

Improved algorithms for generalized finite element simulations of three-dimensional hydraulic fracture propagation

This paper reports improvements to algorithms for the simulation of 3-D hydraulic fracturing with the Generalized Finite Element Method (GFEM). Three optimizations are presented and analyzed. First, an improved initial guess based on solving a 3-D elastic problem with the pressure from the previous step is shown to decrease the number of Newton iterations and increase robustness. Second, an improved methodology to find the time step that leads to fracture propagation is proposed and shown to decrease significantly the number of iterations. Third, reduced computational cost is observed by properly recycling the linear part of the coupled stiffness matrix. Two representative examples are used to analyze these improvements. Additionally, a methodology to include the leak-off term is presented and verified against asymptotic analytical solutions. Conservation of mass is shown to be well satisfied in all examples.