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.

Runoff, sediment loss and water quality from forest roads in a southeast Queensland coastal plain Pinus plantation

Runoff and sediment loss from forest roads were monitored for a two-year period in a Pinus plantation in southeast Queensland. Two classes of road were investigated: a gravelled road, which is used as a primary daily haulage route for the logging area, and an ungravelled road, which provides the main access route for individual logging compartments and is intensively used as a haulage route only during the harvest of these areas (approximately every 30 years). Both roads were subjected to routine traffic loads and maintenance during the study. Surface runoff in response to natural rainfall was measured and samples taken for the determination of sediment and nutrient (total nitrogen, total phosphorus, dissolved organic carbon and total iron) loads from each road. Results revealed that the mean runoff coefficient (runoff depth/rainfall depth) was consistently higher from the gravelled road plot with 0.57, as compared to the ungravelled road with 0.38. Total sediment loss over the two-year period was greatest from the gravelled road plot at 5.7 t km−1 compared to the ungravelled road plot with 3.9 t km−1. Suspended solids contributed 86% of the total sediment loss from the gravelled road, and 72% from the ungravelled road over the two years. Nitrogen loads from the two roads were both relatively constant throughout the study, and averaged 5.2 and 2.9 kg km−1 from the gravelled and ungravelled road, respectively. Mean annual phosphorus loads were 0.6 kg km−1 from the gravelled road and 0.2 kg km−1 from the ungravelled road. Organic carbon and total iron loads increased in the second year of the study, which was a much wetter year, and are thought to reflect the breakdown of organic matter in roadside drains and increased sediment generation, respectively. When road and drain maintenance (grading) was performed runoff and sediment loss were increased from both road types. Additionally, the breakdown of the gravel road base due to high traffic intensity during wet conditions resulted in the formation of deep (10 cm) ruts which increased erosion. The Water Erosion Prediction Project (WEPP):Road model was used to compare predicted to observed runoff and sediment loss from the two road classes investigated. For individual rainfall events, WEPP:Road predicted output showed strong agreement with observed values of runoff and sediment loss. WEPP:Road predictions for annual sediment loss from the entire forestry road network in the study area also showed reasonable agreement with the extrapolated observed values.

Optimal crop planning and conjunctive use of water resources in a coastal river basin

Due to increasing trend of intensive rice cultivation in a coastal river basin, crop planning and groundwater management are imperative for the sustainable agriculture. For effective management, two models have been developed viz. groundwater balance model and optimum cropping and groundwater management model to determine optimum cropping pattern and groundwater allocation from private and government tubewells according to different soil types (saline and non-saline), type of agriculture (rainfed and irrigated) and seasons (monsoon and winter). A groundwater balance model has been developed considering mass balance approach. The components of the groundwater balance considered are recharge from rainfall, irrigated rice and non-rice fields, base flow from rivers and seepage flow from surface drains. In the second phase, a linear programming optimization model is developed for optimal cropping and groundwater management for maximizing the economic returns. The models developed were applied to a portion of coastal river basin in Orissa State, India and optimal cropping pattern for various scenarios of river flow and groundwater availability was obtained.

Prediction of Ground Water Levels in the Uplands of a Tropical Coastal Riparian Wetland using Artificial Neural Networks

Artificial Neural Networks (ANNs) have been found to be a robust tool to model many non-linear hydrological processes. The present study aims at evaluating the performance of ANN in simulating and predicting ground water levels in the uplands of a tropical coastal riparian wetland. The study involves comparison of two network architectures, Feed Forward Neural Network (FFNN) and Recurrent Neural Network (RNN) trained under five algorithms namely Levenberg Marquardt algorithm, Resilient Back propagation algorithm, BFGS Quasi Newton algorithm, Scaled Conjugate Gradient algorithm, and Fletcher Reeves Conjugate Gradient algorithm by simulating the water levels in a well in the study area. The study is analyzed in two cases-one with four inputs to the networks and two with eight inputs to the networks. The two networks-five algorithms in both the cases are compared to determine the best performing combination that could simulate and predict the process satisfactorily. Ad Hoc (Trial and Error) method is followed in optimizing network structure in all cases. On the whole, it is noticed from the results that the Artificial Neural Networks have simulated and predicted the water levels in the well with fair accuracy. This is evident from low values of Normalized Root Mean Square Error and Relative Root Mean Square Error and high values of Nash-Sutcliffe Efficiency Index and Correlation Coefficient (which are taken as the performance measures to calibrate the networks) calculated after the analysis. On comparison of ground water levels predicted with those at the observation well, FFNN trained with Fletcher Reeves Conjugate Gradient algorithm taken four inputs has outperformed all other combinations.

Hydrographical characteristics and oxygen isotopic signatures of water in a coastal environment (Mangalore) along the southeastern Arabian Sea

Coastal marine environments are important links between the continents and the open ocean. The coast off Mangalore forms part of the upwelling zone along the southeastern Arabian Sea. The temperature, salinity, density, dissolved oxygen and stable oxygen isotope ratio (delta O-18) of surface waters as well as those of bottom waters off coastal Mangalore were studied every month from October 2010 to May 2011. The coastal waters were stratified in October and November due to precipitation and runoff. The region was characterised by upwelled bottom waters in October, whereas the region exhibited a temperature inversion in November. The surface and bottom waters presented almost uniform properties from December until April. The coastal waters were observed to be most dense in January and May. Comparatively cold and poorly oxygenated bottom waters during the May sampling indicated the onset of upwelling along the region. delta O-18 of the coastal waters successfully documented the observed variations in the hydrographical characteristics of the Mangalore coast during the monthly sampling period. We also noted that the monthly variability in the properties of the coastal waters of Mangalore was related to the hydrographical characteristics of the adjacent open ocean inferred from satellite-derived surface winds, sea surface height anomaly data and sea surface temperatures.

Water quality and harmful algae in Southeastern coastal stormwater ponds

Several long-term monitoring studies describing the water quality and biological condition of Southeastern estuaries (National Estuarine Eutrophication Assessment Project, South Carolina Estuarine and Coastal Assessment Program (SCECAP), Environmental Monitoring and Assessment Program (EMAP), South Carolina Harmful Algal Bloom Program (SCHAB), South Carolina Tidal Creek Project, and others) have been developed. Many of the same water quality issues determined for open estuaries are also found in coastal stormwater ponds, and there are important interactions between the man-made ponds and the natural systems. Researchers have highlighted problems such as nutrient eutrophication, bacterial and chemical contamination, hypoxia, and harmful algal blooms (HABs). This technical memorandum summarizes the state-of-the-knowledge of water quality indicators (dissolved oxygen, nutrients, and chlorophyll a), and harmful algae in Southeastern coastal stormwater ponds. (PDF contains 31 pages)

Assessing the sustainability of brackish water aquaculture systems in the Philippines. Working Paper 4: Aquaculture and poverty - a case study of five coastal communities in the Philippines

After reviewing the rather thin literature on the subject, we investigate the relationship between aquaculture and poverty based on a case study of five coastal communities in the Philippines. The analysis relies on a data set collated through a questionnaire survey of 148 households randomly selected in these five communities. The methodological approach combines the qualitative analysis of how this relationship is perceived by the surveyed households and a quantitative analysis of the levels and determinants of poverty and inequality in these communities. There is overwhelming evidence that aquaculture benefits the poor in important ways and that it is perceived very positively by the poor and non-poor alike. In particular, the poor derive a relatively larger share of their income from aquaculture than the rich, and a lowering of the poverty line only reinforces this result. Further, a Gini decomposition exercise shows unambiguously that aquaculture represents an inequality-reducing source of income. We believe that the pro-poor character of brackish water aquaculture in the study areas is explained by the fact that the sector provides employment to a large number of unskilled workers in communities characterized by large surpluses of labour. Our results also suggest that the analysis of the relationship between aquaculture and poverty should not focus exclusively on the socio-economic status of the farm operator/owner, as has often been the case in the past. [PDF contains 51 pages]

Towed vehicles: undulating platforms as tools for maping coastal processes and water quality assessment, Seaside, California, February 5-7,2007. Workshop proceedings

The Alliance for Coastal Technologies (ACT) Workshop on Towed Vehicles: Undulating Platforms As Tools for Mapping Coastal Processes and Water Quality Assessment was convened February 5-7,2007 at The Embassy Suites Hotel, Seaside, California and sponsored by the ACT-Pacific Coast partnership at the Moss Landing Marine Laboratories (MLML). The TUV workshop was co-chaired by Richard Burt (Chelsea Technology Group) and Stewart Lamerdin (MLML Marine Operations). Invited participants were selected to provide a uniform representation of the academic researchers, private sector product developers, and existing and potential data product users from the resource management community to enable development of broad consensus opinions on the application of TUV platforms in coastal resource assessment and management. The workshop was organized to address recognized limitations of point-based monitoring programs, which, while providing valuable data, are incapable of describing the spatial heterogeneity and the extent of features distributed in the bulk solution. This is particularly true as surveys approach the coastal zone where tidal and estuarine influences result in spatially and temporally heterogeneous water masses and entrained biological components. Aerial or satellite based remote sensing can provide an assessment of the aerial extent of plumes and blooms, yet provide no information regarding the third dimension of these features. Towed vehicles offer a cost-effective solution to this problem by providing platforms, which can sample in the horizontal, vertical, and time-based domains. Towed undulating vehicles (henceforth TUVs) represent useful platforms for event-response characterization. This workshop reviewed the current status of towed vehicle technology focusing on limitations of depth, data telemetry, instrument power demands, and ship requirements in an attempt to identify means to incorporate such technology more routinely in monitoring and event-response programs. Specifically, the participants were charged to address the following: (1) Summarize the state of the art in TUV technologies; (2) Identify how TUV platforms are used and how they can assist coastal managers in fulfilling their regulatory and management responsibilities; (3) Identify barriers and challenges to the application of TUV technologies in management and research activities, and (4) Recommend a series of community actions to overcome identified barriers and challenges. A series of plenary presentation were provided to enhance subsequent breakout discussions by the participants. Dave Nelson (University of Rhode Island) provided extensive summaries and real-world assessment of the operational features of a variety of TUV platforms available in the UNOLs scientific fleet. Dr. Burke Hales (Oregon State University) described the modification of TUV to provide a novel sampling platform for high resolution mapping of chemical distributions in near real time. Dr. Sonia Batten (Sir Alister Hardy Foundation for Ocean Sciences) provided an overview on the deployment of specialized towed vehicles equipped with rugged continuous plankton recorders on ships of opportunity to obtain long-term, basin wide surveys of zooplankton community structure, enhancing our understanding of trends in secondary production in the upper ocean. [PDF contains 32 pages]

From the individual to the community and beyond: water quality, stress indicators and key species in coastal waters

This review examines water quality and stress indicators at levels of organisation from the individual to the community and beyond by means of three case studies concentrating on rocky shores within the north-east Atlantic. Responses of dogwhelks (Nucella) to tributyltin pollution from antifouling paints is examined as the main case study. There are effects at the individual level (development of male sexual characteristics in the female leading to effective sterility) and population level (reduction in juveniles, few females and eventual population disappearance of dogwhelks in badly contaminated areas) but information on community level effects of dogwhelk demise is sparse. Such effects were simulated by dogwhelk removal experiments on well studied, moderately exposed ledges on shores on the Isle of Man. The removal of dogwhelks reduced the size and longevity of newly established Fucus clumps that had escaped grazing. Removal of dogwhelks also increased the likelihood of algal escapes. In a factorial experiment dogwhelks were shown to be less important than limpets \{Patella) in structuring communities but still had a significant modifying effect by increasing the probability of algal escapes. Community level responses to stress on rocky shores are then explored by reference to catastrophic impacts such as oil spills, using the Torrey Canyon as a case study. Recovery of the system in response to this major perturbation took between 10-15 years through a series of damped oscillations. The final case study is that of indicators of ecosystem level change in response to climate fluctuations, using ratios of northern \{Semibalanus balanoides) and southern (Chthamalus spp.) barnacles. Indices derived from counts on the shore show good correlations with inshore sea-water temperatures after a 2-year lag phase. The use of barnacles to measure offshore changes is reviewed. The discussion considers the use of bioindicators at various levels of organisation.

Coastal States and Distant-water Fishing Nation Relations: An Economist's Perspective

The widespread implementation of Extended Fisheries Jurisdiction (EFJ) has confronted coastal states with several resource management problems. One of these consists of the economic relations, if any, that the coastal state should establish with distant-water fishing nations (DWFN's) seeking access to the coastal state's 200-mile zone. Several of the other papers presented here deal with specific aspects of the issue. This paper, on the other hand, will concern itself with the question of the analytical framework to be used by economists in studying this issue. It will offer some suggestions with respect to possible components of the framework. In doing so, the paper will restrict itself to the coastal state's perspective of EFJ and the management issues arising therefrom. It goes without saying, of course, that an enlightened coastal state will attempt to acquaint itself with the DWFN view of the world.

Using wildlife activity and antibiotic resistance analysis to model bacterial water quality in coastal ponds

Models that help predict fecal coliform bacteria (FCB) levels in environmental waters can be important tools for resource managers. In this study, we used animal activity along with antibiotic resistance analysis (ARA), land cover, and other variables to build models that predict bacteria levels in coastal ponds that discharge into an estuary. Photographic wildlife monitoring was used to estimate terrestrial and aquatic wildlife activity prior to sampling. Increased duck activity was an important predictor of increased FCB in coastal ponds. Terrestrial animals like deer and raccoon, although abundant, were not significant in our model. Various land cover types, rainfall, tide, solar irradiation, air temperature, and season parameters, in combination with duck activity, were significant predictors of increased FCB. It appears that tidal ponds allow for settling of bacteria under most conditions. We propose that these models can be used to test different development styles and wildlife management techniques to reduce bacterial loading into downstream shellfish harvesting and contact recreation areas.

Salt and water balance in the Mundel Lake: a strongly choked coastal lagoon

The Mundel Lake is an extremely shallow lagoon on the west coast of Sri Lanka. It is connected to the Puttalam Lagoon through 15 km long Dutch Canal. Salinity measurements and daily sea level data were obtained fortnightly from January 1993 to March 1994 and they were used to quantify the salt and water budget along with precipitation, evaporation and freshwater runoff. Extreme fluctuations of salinity and sea level are striking features of the system. Salinity of the Mundel Lake and Dutch Canal varied from 5-46.5 and 6 61 ppt respectively while the sea level ranged from -0.25 to +1.2 m. Tidal variations were not seen in the lagoon due to its long narrow canal system. Salt budget showed that the deposition of salt on the lagoon bottom during periods of decreasing water level. During increasing water level, salt is dissolved again. Flow of water through the Dutch Canal between the Puttalam Lagoon and Mundel Lake is driven by the changes in sea level. These changes are mainly due to seasonal changes of net freshwater supply and, to a lesser degree, to seasonal changes in sea surface height. As the flow rates are small due to the long and narrow canal, the residence time ranges between two months and several months in the Mundel Lake, except during season of high freshwater supply. As the water exchange is weak, the Mundel Lake becomes hyper saline with strong fluctuations in salinity. This implies a stress to all lagoon dwelling aquatic organisms and also to aquaculture practices in the area.