Resource management and Māori attitudes to water in southern New Zealand

Abstract:  Pre-contact Māori regarded land and water as a single entity, with a common regime of resource management practices. Underpinning these was a world-view that involved unique spiritual concepts, the most important of which was mauri: the notion that a body of water had its own life-force. Waters were classified according to the state of their mauri. The paper outlines traditional approaches and how they are applied today.     

Bar-top hollows: A new element in the architecture of sandy braided rivers

Discrete hollows in the bar tops of the South Saskatchewan River are described that form a newly-recognized morphological element of sandy braided rivers. These bar-top hollows, which are up to 1.7 m deep and may extend for 10–30 m down and across flow, have a circular to ovoid planform and are shown, through use of ground penetrating radar, to be filled by a series of distinct, often angle-of-repose, foresets. The hollows form by both erosion and bar-top deposition and may be generated by bar-tail accretion, cross-bar channel cutoff and subsequent fill or lateral accretion at the bar-head. Bar-top hollows occur in the upper part of the bar depositional sequence and may thus prove useful indicators for braid bar reconstruction in ancient sediments, and should not be confused with channel scour.     

Assessing catchment-scale spatial and temporal patterns of groundwater and stream salinity

Understanding catchment-scale patterns of groundwater and stream salinity are important in land- and water-salinity management. A large-scale assessment of groundwater and stream data was undertaken in the eastern Mt Lofty Ranges of South Australia using geographical information systems (GIS), regional scale hydrologic data, hydrograph separation and hydrochemical techniques. Results of the study show: (1) salts were mostly of marine origin (75%), while sulfate and bicarbonate from mineral weathering comprised most of the remainder, (2) elevated groundwater salinities and stable water isotopic compositions similar to mean rainfall indicated that plant transpiration was the primary salt accumulation mechanism, (3) key factors explaining groundwater salinity were geology and rainfall, with overall catchment salinity inversely proportional to average annual rainfall, and groundwater salinity ‘hotspots’ (EC >8 mS/cm) associated with geological formations comprising sulfidic marine siltstones and shales, (4) shallow groundwater correlated with elevated stream salinity, implying that baseflow contributed to stream salt loads, with most of the annual salt load (estimated to be 24,500 tonnes) occurring in winter when baseflow volume was highest. Salt-load analysis using stream data could be a practical, low-cost technique to rapidly target the investigation of problem areas within a catchment.     

Climatic Change, Wars and Dynastic Cycles in China Over the Last Millennium

In recent years, the phenomenon of global warming and its implications for the future of the human race have been intensively studied. In contrast, few quantitative studies have been attempted on the notable effects of past climatic changes upon human societies. This study explored the relationship between climatic change and war in China by comparing high-resolution paleo-climatic reconstructions with known war incidences in China in the last millennium. War frequencies showed a cyclic pattern that closely followed the global paleo-temperature changes. Strong and significant correlations were found between climatic change, war occurrence, harvest level, population size and dynastic transition. During cold phases, China suffered more often from frequent wars, population decline and dynastic changes. The quantitative analyses suggested that the reduction of thermal energy input during a cold phase would lower the land carrying capacity in the traditional agrarian society, and the population size, with significant accretions accrued in the previous warm phase, could not be sustained by the shrinking resource base. The stressed human-nature relationship generated a ‘push force’, leading to more frequent wars between states, regions and tribes, which could lead to the collapse of dynasties and collapses of human population size. War frequencies varied according to geographical locations (North, Central and South China) due to spatial variations in the physical environment and hence differential response to climatic change. Moreover, war occurrences demonstrated an obvious time lag after an episode of temperature fall, and the three geographical regions experienced different length of time lags. This research also shows that human population increases and collapses were correlated with the climatic phases and the social instabilities that were induced by climate changes during the last millennium. The findings proposed a new interpretation of human-nature relationship in the past, with implications for the impacts of anomalous global warming on future human conflicts.     

Study of the Probability Density Functions From a Large-Eddy Simulation for a Stably Stratified Boundary Layer

Turbulence in a non-strongly stably stratified large-eddy simulation (LES) case is studied through probability density functions (PDFs) to obtain additional information than that provided by classical LES averages. The PDFs are computed for one hour within the steady-state regime at three different levels: near the surface, in the middle and at the top of the boundary layer, for the wind components and the temperature. The physical significance of these PDFs from LES is discussed and they are compared to those obtained from observations. The analysis of the eddy structures within the stably stratified boundary layer is made through the combined study of the fields, the spectra and the statistical moments obtained from the PDFs and joint PDFs. The homogeneity of the fields is inspected through a comparison of the ensemble to the temporal and the spatial PDFs, showing that the ergodic theorem is not fulfilled. To this end, the sensitivity of the PDF moments to the LES resolution is explored.     

Binary upscaling—the role of connectivity and a new formula

Although recognized as important, measures of connectivity (i.e. the existence of high-conductivity paths that increase flow and allow for early solute arrival) have not yet been incorporated into methods for upscaling hydraulic conductivities of porous media. We present and evaluate a binary upscaling formula that utilizes connectivity information. The upscaled hydraulic conductivity (K) of binary media is determined as a function of the proportions and conductivities of the two materials, the geometry of the inclusions, and the mean distance between them. The use of a phase interchange theorem renders the formula equally applicable to two-dimensional media with inclusions of low K and high K as compared with the matrix. The new upscaling formula is tested on two-dimensional binary random fields spanning a broad range of spatial correlation structures and conductivity contrasts. The computed effective conductivities are compared to what is obtained using self-consistent effective medium theory, the coated ellipsoids approximation, and to a streamline approach. It is shown that, although simple, the proposed formula performs better than available methods for binary upscaling. The use of connectivity information leads to significantly improved behavior close to the percolation threshold. The proposed upscaling formula depends exclusively on parameters that are obtainable from field investigations.     

Emplacement of subaerial pahoehoe lava sheet flows into water: 1990 Kūpaianaha flow of Kilauea volcano at Kaimū Bay, Hawai`i

Episode 48 of the ongoing eruption of Kilauea, Hawai`i, began in July 1986 and continuously extruded lava for the next 5.5 years from a low shield, Kūpaianaha. The flows in March 1990 headed for Kalapana and inundated the entire town under 15–25 m of lava by the end of August. As the flows advanced eastward, they entered into Kaimū Bay, replacing it with a plain of lava that extends 300 m beyond the original shoreline. The focus of our study is the period from August 1 to October 31, 1990, when the lava buried almost 406,820 m² of the 5-m deep bay. When lava encountered the sea, it flowed along the shoreline as a narrow primary lobe up to 400 m long and 100 m wide, which in turn inflated to a thickness of 5–6 m. The flow direction of the primary lobes was controlled by the submerged delta below the lavas and by damming up lavas fed at low extrusion rates. Breakout flows through circumferential and axial inflation cracks on the inflating primary lobes formed smaller secondary lobes, burying the lows between the primary lobes and hiding their original outlines. Inflated flow lobes eventually ruptured at proximal and/or distal ends as well as mid-points between the two ends, feeding new primary lobes which were emplaced along and on the shore side of the previously inflated lobes. The flow lobes mapped with the aid of aerial photographs were correlated with daily observations of the growing flow field, and 30 primary flow lobes were dated. Excluding the two repose periods that intervened while the bay was filled, enlargement of the flow field took place at a rate of 2,440–22,640 square meters per day in the bay. Lobe thickness was estimated to be up to 11 m on the basis of cross sections of selected lobes measured using optical measurement tools, measuring tape and hand level. The total flow-lobe volume added in the bay during August 1–October 31 was approximately 3.95 million m³, giving an average supply rate of 0.86 m³/s.     

Quantifying Damage, Saturation and Anisotropy in Cracked Rocks by Inverting Elastic Wave Velocities

Crack damage results in a decrease of elastic wave velocities and in the development of anisotropy. Using non-interactive crack effective medium theory as a fundamental tool, we calculate dry and wet elastic properties of cracked rocks in terms of a crack density tensor, average crack aspect ratio and mean crack fabric orientation from the solid grains and fluid elastic properties. Using this same tool, we show that both the anisotropy and shear-wave splitting of elastic waves can be derived. Two simple crack distributions are considered for which the predicted anisotropy depends strongly on the saturation, reaching up to 60% in the dry case. Comparison with experimental data on two granites, a basalt and a marble, shows that the range of validity of the non-interactive effective medium theory model extends to a total crack density of approximately 0.5, considering symmetries up to orthorhombic. In the isotropic case, Kachanov's (1994) non-interactive effective medium model was used in order to invert elastic wave velocities and infer both crack density and aspect ratio evolutions. Inversions are stable and give coherent results in terms of crack density and aperture evolution. Crack density variations can be interpreted in terms of crack growth and/or changes of the crack surface contact areas as cracks are being closed or opened respectively. More importantly, the recovered evolution of aspect ratio shows an exponentially decreasing aspect ratio (and therefore aperture) with pressure, which has broader geophysical implications, in particular on fluid flow. The recovered evolution of aspect ratio is also consistent with current mechanical theories of crack closure. In the anisotropic cases—both transverse isotropic and orthorhombic symmetries were considered—anisotropy and saturation patterns were well reproduced by the modelling, and mean crack fabric orientations we recovered are consistent with in situ geophysical imaging. Our results point out that: (1) It is possible to predict damage, anisotropy and saturation in terms of a crack density tensor and mean crack aspect ratio and orientation; (2) using well constrained wave velocity data, it is possible to extrapolate the contemporaneous evolution of crack density, anisotropy and saturation using wave velocity inversion as a tool; 3) using such an inversion tool opens the door in linking elastic properties, variations to permeability.     

A parallelizable method for two‐phase flows in naturally‐fractured reservoirs

A parallelizable, semi‐implicit numerical method is proposed for the study of naturally‐fractured reservoir systems. It has proved to be computationally efficient in producing accurate numerical solutions for the dual‐porosity model for immiscible, two‐phase flow in such reservoirs. The method combines hybridized mixed finite elements, a new version of the modified method of characteristics, a sophisticated operator‐splitting procedure for separating the pressure calculation in the fractures from that of the saturation, another operator splitting to handle the interaction of the matrix blocks and the fractures, and domain decomposition iterative procedures for both the pressure and the saturation. It permits moderately long time steps for the pressure and the saturation in the fractures and matrix blocks by using short, inexpensive microsteps to treat the transport portion of the saturation equation in the fractures. This paper is devoted to the formulation of the method and a discussion of numerical results for five‐spot and vertical cross‐section examples.