Modelling sea water intrusion in coastal aquifers using heterogeneous computing

The objective of this PhD research program is to investigate numerical methods for simulating variably-saturated flow and sea water intrusion in coastal aquifers in a high-performance computing environment. The work is divided into three overlapping tasks: to develop an accurate and stable finite volume discretisation and numerical solution strategy for the variably-saturated flow and salt transport equations; to implement the chosen approach in a high performance computing environment that may have multiple GPUs or CPU cores; and to verify and test the implementation. The geological description of aquifers is often complex, with porous materials possessing highly variable properties, that are best described using unstructured meshes. The finite volume method is a popular method for the solution of the conservation laws that describe sea water intrusion, and is well-suited to unstructured meshes. In this work we apply a control volume-finite element (CV-FE) method to an extension of a recently proposed formulation (Kees and Miller, 2002) for variably saturated groundwater flow. The CV-FE method evaluates fluxes at points where material properties and gradients in pressure and concentration are consistently defined, making it both suitable for heterogeneous media and mass conservative. Using the method of lines, the CV-FE discretisation gives a set of differential algebraic equations (DAEs) amenable to solution using higher-order implicit solvers. Heterogeneous computer systems that use a combination of computational hardware such as CPUs and GPUs, are attractive for scientific computing due to the potential advantages offered by GPUs for accelerating data-parallel operations. We present a C++ library that implements data-parallel methods on both CPU and GPUs. The finite volume discretisation is expressed in terms of these data-parallel operations, which gives an efficient implementation of the nonlinear residual function. This makes the implicit solution of the DAE system possible on the GPU, because the inexact Newton-Krylov method used by the implicit time stepping scheme can approximate the action of a matrix on a vector using residual evaluations. We also propose preconditioning strategies that are amenable to GPU implementation, so that all computationally-intensive aspects of the implicit time stepping scheme are implemented on the GPU. Results are presented that demonstrate the efficiency and accuracy of the proposed numeric methods and formulation. The formulation offers excellent conservation of mass, and higher-order temporal integration increases both numeric efficiency and accuracy of the solutions. Flux limiting produces accurate, oscillation-free solutions on coarse meshes, where much finer meshes are required to obtain solutions with equivalent accuracy using upstream weighting. The computational efficiency of the software is investigated using CPUs and GPUs on a high-performance workstation. The GPU version offers considerable speedup over the CPU version, with one GPU giving speedup factor of 3 over the eight-core CPU implementation.

Sun, sea, sex, sand and soap

A recent issue of Young People Now (November 1995) mentioned the new (UK) Channel 4 television soap opera Hollyoaks by Phil Redmond, which raises the issue of the role of ‘soaps’ in the daily lives of young people. Australian soaps are especially popular in Britain and of interest to those who work with young people, because they have a high proportion of youthful looking actors and actresses and frequently depict scenes involving young people and apparent ‘real’ teenage dilemmas.

Sun, sea, sex, sand and conservatism in the Australian television soap opera

A recent issue of Young People Now (November 1995) mentioned the new (UK) television soap opera Hollyoaks by Phil Redmond, which raises the issue of the role of ‘soap operas’ (hereafter referred to as soaps) in the daily lives of young people. The term ‘soap’ originates with the sponsorship of radio and television programmes by companies such as Proctor and Gamble who in America in 1932 used a daytime radio domestic comedy, The Puddle Family to advertise Oxydol, a washing powder. The first British television soap was The Grove Family (BBC 1954-7) was followed by Emergency Ward Ten (ATV 1957-67), Coronation Street (Granada Television 1960-present) and Eastenders (BBC 1985-present). Australian soaps are especially popular in Britain and of potential interest to those who work with young people, because they have a high proportion of youthful looking actors and actresses and frequently depict scenes involving young people and apparent ‘real’ teenage dilemmas. On one level it may be commendable that actors who are young(ish) somewhere between the ages of 14 and 25 play roles that are ostensibly about young people and their alleged problems. However, the casting of young, largely unknown, actors reflects more the political economy of soaps in their relative cheapness and dispensability, rather than any genuine attempt to create an oppositional text for, about and by young people (Paterson 1986).

Past 140-year environmental record in the northern South China Sea : evidence from coral skeletal trace metal variations

About 140-year changes in the trace metals in Porites coral samples from two locations in the northern South China Sea were investigated. Results of PCA analyses suggest that near the coast, terrestrial input impacted behavior of trace metals by 28.4%, impact of Sea Surface Temperature (SST) was 19.0%, contribution of war and infrastructure were 14.4% and 15.6% respectively. But for a location in the open sea, contribution of War and SST reached 33.2% and 16.5%, while activities of infrastructure and guano exploration reached 13.2% and 14.7%. While the spatiotemporal change model of Cu, Cd and Pb in seawater of the north area of South China Sea during 1986–1997 were reconstructed. It was found that in the sea area Cu and Cd contaminations were distributed near the coast while areas around Sanya, Hainan had high Pb levels because of the well-developed tourism related activities.

Litigating law of the sea disputes using the UNCLOS Dispute Settlement System

The interpretation and application of the United Nations Convention on the Law of the Sea (UNCLOS) may be the source of many disputes. UNCLOS introduced an à la carte menu for dispute settlement with a number of options for international dispute resolution, including a compulsory procedure entailing binding decisions. While drafting this ambitious and complex system of dispute settlement, the drafters had to negotiate many delicate compromises to secure a system for the uniform interpretation of the Convention. The aim of this paper r is to explore why litigation using the UNCLOS dispute settlement system is, or is not, a preferred mode of settlement for law of the sea disputes.

Mid-Holocene sea-level and coral reef demise : U-Th dating of subfossil corals in Moreton Bay, Australia

It is increasingly apparent that sea-level data (e.g. microfossil transfer functions, dated coral microatolls and direct observations from satellite and tidal gauges) vary temporally and spatially at regional to local scales, thus limiting our ability to model future sea-level rise for many regions. Understanding sealevel response at ‘far-field’ locations at regional scales is fundamental for formulating more relevant sea-level rise susceptibility models within these regions under future global change projections. Fossil corals and reefs in particular are valuable tools for reconstructing past sea levels and possible environmental phase shifts beyond the temporal constraints of instrumental records. This study used abundant surface geochronological data based on in situ subfossil corals and precise elevation surveys to determine previous sea level in Moreton Bay, eastern Australia, a far-field site. A total of 64 U-Th dates show that relative sea level was at least 1.1 m above modern lowest astronomical tide (LAT) from at least ˜6600 cal. yr BP. Furthermore, a rapid synchronous demise in coral reef growth occurred in Moreton Bay ˜5800 cal. yr BP, coinciding with reported reef hiatus periods in other areas around the Indo-Pacific region. Evaluating past reef growth patterns and phases allows for a better interpretation of anthropogenic forcing versus natural environmental/climatic cycles that effect reef formation and demise at all scales and may allow better prediction of reef response to future global change.

Estimating present day extreme water level exceedance probabilities around the coastline of Australia : tides, extra-tropical storm surges and mean sea level

The occurrence of extreme water levels along low-lying, highly populated and/or developed coastlines can lead to considerable loss of life and billions of dollars of damage to coastal infrastructure. Therefore it is vitally important that the exceedance probabilities of extreme water levels are accurately evaluated to inform risk-based flood management, engineering and future land-use planning. This ensures the risk of catastrophic structural failures due to under-design or expensive wastes due to over-design are minimised. This paper estimates for the first time present day extreme water level exceedence probabilities around the whole coastline of Australia. A high-resolution depth averaged hydrodynamic model has been configured for the Australian continental shelf region and has been forced with tidal levels from a global tidal model and meteorological fields from a global reanalysis to generate a 61-year hindcast of water levels. Output from this model has been successfully validated against measurements from 30 tide gauge sites. At each numeric coastal grid point, extreme value distributions have been fitted to the derived time series of annual maxima and the several largest water levels each year to estimate exceedence probabilities. This provides a reliable estimate of water level probabilities around southern Australia; a region mainly impacted by extra-tropical cyclones. However, as the meteorological forcing used only weakly includes the effects of tropical cyclones, extreme water level probabilities are underestimated around the western, northern and north-eastern Australian coastline. In a companion paper we build on the work presented here and more accurately include tropical cyclone-induced surges in the estimation of extreme water level. The multi-decadal hindcast generated here has been used primarily to estimate extreme water level exceedance probabilities but could be used more widely in the future for a variety of other research and practical applications.

Groundwater table response to sea level rise and its impact on pavement structure

It has been predicted that sea level will rise about 0.8 m by 2100. Consequently, seawater can intrude into the coastal aquifers and change the level of groundwater table. A raise in groundwater table due to seawater intrusion threats the coastal infrastructure such as road pavements. The mechanical properties of subgrade materials will change due to elevated rise of groundwater table, leading to pavement weakening and decreasing the subgrade strength and stiffness. This paper presents an assessment of the vulnerability of subgrade in coastal areas to change in groundwater table due to sea-level rise. A simple bathtub approach is applied for estimating the groundwater level changes according to sea-level rise. Then the effect of groundwater level changes on the soil water content (SWC) of a single column of fine-sand soil is simulated using MIKE SHE. The impact of an increase in moisture content on subgrade strength/stiffness is assessed for a number of scenarios.

The SEA-CALIPSO volcano imaging experiment at Montserrat : plans, campaigns at sea and on land, scientific results, and lessons learned

Since 1995 the eruption of the andesitic Soufrière Hills Volcano (SHV), Montserrat, has been studied in substantial detail. As an important contribution to this effort, the Seismic Experiment with Airgunsource-Caribbean Andesitic Lava Island Precision Seismo-geodetic Observatory (SEA-CALIPSO) experiment was devised to image the arc crust underlying Montserrat, and, if possible, the magma system at SHV using tomography and reflection seismology. Field operations were carried out in October–December 2007, with deployment of 238 seismometers on land supplementing seven volcano observatory stations, and with an array of 10 ocean-bottom seismometers deployed offshore. The RRS James Cook on NERC cruise JC19 towed a tuned airgun array plus a digital 48-channel streamer on encircling and radial tracks for 77 h about Montserrat during December 2007, firing 4414 airgun shots and yielding about 47 Gb of data. The main objecctives of the experiment were achieved. Preliminary analyses of these data published in 2010 generated images of heterogeneous high-velocity bodies representing the cores of volcanoes and subjacent intrusions, and shallow areas of low velocity on the flanks of the island that reflect volcaniclastic deposits and hydrothermal alteration. The resolution of this preliminary work did not extend beyond 5 km depth. An improved three-dimensional (3D) seismic velocity model was then obtained by inversion of 181 665 first-arrival travel times from a more-complete sampling of the dataset, yielding clear images to 7.5 km depth of a low-velocity volume that was interpreted as the magma chamber which feeds the current eruption, with an estimated volume 13 km3. Coupled thermal and seismic modelling revealed properties of the partly crystallized magma. Seismic reflection analyses aimed at imaging structures under southern Montserrat had limited success, and suggest subhorizontal layering interpreted as sills at a depth of between 6 and 19 km. Seismic reflection profiles collected offshore reveal deep fans of volcaniclastic debris and fault offsets, leading to new tectonic interpretations. This chapter presents the project goals and planning concepts, describes in detail the campaigns at sea and on land, summarizes the major results, and identifies the key lessons learned.

Road infrastructure vulnerability to groundwater table variation due to sea level rise

This study was conducted to assess the vulnerability of coastal road infrastructures due to climate change induced sea level rise and extreme weather conditions through the estimation of road subgrade strength reduction as a result of changes in soil moisture content. The study area located in the Gold Coast, Australia highlighted that the risk is significant. In wet seasons or areas with wet condition, the groundwater table is already high, so even a small change in the groundwater table can raise the risk of inundation; particularly, in areas with existing shallow groundwater. The predicted risk of a high groundwater table on road infrastructure is a long-term hazard. Therefore, there is time to undertake some management plans to decrease the possible risks, for instance, some deep root plants could be planted along the roads with a high level of risk, to decrease the groundwater table elevation.

The relationship between N2O, NO, and N2 fluxes from fertilized and irrigated dryland soils of the Aral Sea Basin, Uzbekistan

Microbial respiratory reduction of nitrous oxide (N2O) to dinitrogen (N2) via denitrification plays a key role within the global N-cycle since it is the most important process for converting reactive nitrogen back into inert molecular N2. However, due to methodological constraints, we still lack a comprehensive, quantitative understanding of denitrification rates and controlling factors across various ecosystems. We investigated N2, N2O and NO emissions from irrigated cotton fields within the Aral Sera Basin using the He/O2 atmosphere gas flow soil core technique and an incubation assay. NH4NO3 fertilizer, equivalent to 75 kg ha−1 and irrigation water, adjusting the water holding capacity to 70, 100 and 130% were applied to the incubation vessels to assess its influence on gaseous N emissions. Under soil conditions as they are naturally found after concomitant irrigation and fertilization, denitrification was the dominant process and N2 the main end product of denitrification. The mean ratios of N2/N2O emissions increased with increasing soil moisture content. N2 emissions exceeded N2O emissions by a factor of 5 ± 2 at 70% soil water holding capacity (WHC) and a factor of 55 ± 27 at 130% WHC. The mean ratios of N2O/NO emissions varied between 1.5 ± 0.4 (70% WHC) and 644 ± 108 (130% WHC). The magnitude of N2 emissions for irrigated cotton was estimated to be in the range of 24 ± 9 to 175 ± 65 kg-N ha−1season−1, while emissions of NO were only of minor importance (between 0.1 to 0.7 kg-N ha−1 season−1). The findings demonstrate that for irrigated dryland soils in the Aral Sera Basin, denitrification is a major pathway of N-loss and that substantial amounts of N-fertilizer are lost as N2 to the atmosphere for irrigated dryland soils.

Nitrous oxide emissions from fertilized, irrigated cotton (Gossypium hirsutum L.) in the Aral Sea Basin, Uzbekistan : influence of nitrogen applications and irrigation practices

Nitrous oxide emissions were monitored at three sites over a 2-year period in irrigated cotton fields in Khorezm, Uzbekistan, a region located in the arid deserts of the Aral Sea Basin. The fields were managed using different fertilizer management strategies and irrigation water regimes. N2O emissions varied widely between years, within 1 year throughout the vegetation season, and between the sites. The amount of irrigation water applied, the amount and type of N fertilizer used, and topsoil temperature had the greatest effect on these emissions. Very high N2O emissions of up to 3000 μg N2O-N m−2 h−1 were measured in periods following N-fertilizer application in combination with irrigation events. These “emission pulses” accounted for 80–95% of the total N2O emissions between April and September and varied from 0.9 to 6.5 kg N2O-N ha−1.. Emission factors (EF), uncorrected for background emission, ranged from 0.4% to 2.6% of total N applied, corresponding to an average EF of 1.48% of applied N fertilizer lost as N2O-N. This is in line with the default global average value of 1.25% of applied N used in calculations of N2O emissions by the Intergovernmental Panel on Climate Change. During the emission pulses, which were triggered by high soil moisture and high availability of mineral N, a clear diurnal pattern of N2O emissions was observed, driven by daily changes in topsoil temperature. For these periods, air sampling from 8:00 to 10:00 and from 18:00 to 20:00 was found to best represent the mean daily N2O flux rates. The wet topsoil conditions caused by irrigation favored the production of N2O from NO3− fertilizers, but not from NH4+ fertilizers, thus indicating that denitrification was the main process causing N2O emissions. It is therefore argued that there is scope for reducing N2O emission from irrigated cotton production; i.e. through the exclusive use of NH4+ fertilizers. Advanced application and irrigation techniques such as subsurface fertilizer application, drip irrigation and fertigation may also minimize N2O emission from this regionally dominant agro-ecosystem.

Methane and nitrous oxide fluxes in annual and perennial land-use systems of the irrigated areas in the Aral Sea Basin

Land use and agricultural practices can result in important contributions to the global source strength of atmospheric nitrous oxide (N2O) and methane (CH4). However, knowledge of gas flux from irrigated agriculture is very limited. From April 2005 to October 2006, a study was conducted in the Aral Sea Basin, Uzbekistan, to quantify and compare emissions of N2O and CH4 in various annual and perennial land-use systems: irrigated cotton, winter wheat and rice crops, a poplar plantation and a natural Tugai (floodplain) forest. In the annual systems, average N2O emissions ranged from 10 to 150 μg N2O-N m−2 h−1 with highest N2O emissions in the cotton fields, covering a similar range of previous studies from irrigated cropping systems. Emission factors (uncorrected for background emission), used to determine the fertilizer-induced N2O emission as a percentage of N fertilizer applied, ranged from 0.2% to 2.6%. Seasonal variations in N2O emissions were principally controlled by fertilization and irrigation management. Pulses of N2O emissions occurred after concomitant N-fertilizer application and irrigation. The unfertilized poplar plantation showed high N2O emissions over the entire study period (30 μg N2O-N m−2 h−1), whereas only negligible fluxes of N2O (<2 μg N2O-N m−2 h−1) occurred in the Tugai. Significant CH4 fluxes only were determined from the flooded rice field: Fluxes were low with mean flux rates of 32 mg CH4 m−2 day−1 and a low seasonal total of 35.2 kg CH4 ha−1. The global warming potential (GWP) of the N2O and CH4 fluxes was highest under rice and cotton, with seasonal changes between 500 and 3000 kg CO2 eq. ha−1. The biennial cotton–wheat–rice crop rotation commonly practiced in the region would average a GWP of 2500 kg CO2 eq. ha−1 yr−1. The analyses point out opportunities for reducing the GWP of these irrigated agricultural systems by (i) optimization of fertilization and irrigation practices and (ii) conversion of annual cropping systems into perennial forest plantations, especially on less profitable, marginal lands.