Estuarine ecosystems are under enormous stress due to rapid coastal developments and climate change. Proper management of these important ecosystems requires a good understanding of their key processes. In this thesis, riparian groundwater-surface water interaction is explored for an aquifer-estuary system primarily by a series of numerical experiments. The work focuses on riparian-scale groundwater flow and salinization. The overall aim of the study was to extend our understanding of aquifer-estuary exchange, which is currently centered on the lower marine estuarine reach, to middle estuaries (i.e., the estuary reach that has variable salinity). The numerical experiments were guided by previous studies and observations made from an exploratory field investigation conducted in and next to Sandy Creek, a macro-tidal estuary incised in the alluvial aquifer of the Pioneer Valley, North-eastern Australia (Longitude 49.11°, Latitude -21.27°). The following observations were made from the field investigation: Sandy Creek estuary experiences a variable salinity regime in its mid reaches that consists of periods of 1) freshwater flushing due to up catchment-derived flooding, 2) persistent freshwater conditions for at least 2 months following the flooding, 3) tidal salinity fluctuations and 4) constant near-seawater salinity; laterally extensive and disconnected aquitards were found to occur at the field site; Sandy Creek had an essentially ‘vertical’ bank slope. Numerical simulations were conducted using the finite element modeling code FEFLOW for saturated unsaturated, variable-density groundwater flow and solute transport, to examine the influence of the following factors on aquifer-estuary exchange: a tidally varying estuarine salinity and hydraulic head, a seasonal freshwater flush (i.e., estuary with freshwater and an elevated stage due to an up catchment sourced flood), near estuary aquitard layers, lateral asymmetry (about the estuary centerline) in hydraulic conductivity and regional hydraulic gradients. The simulations neglected seepage face development after numerical experiments showed that for a vertical bank estuary interacting with a sandy loam aquifer, seepage face effects on groundwater flow and associated salinity distribution were minimal. The following observations were drawn from the range of numerical experiments considered. Tidal salinity fluctuations in the estuary (varying between 0 and 1 - i.e., using a relative salinity scale where a salinity of 1 is seawater) produced flow paths and residence times that were distinctly different to the constant seawater salinity case. While the constant average 0.5 salinity case and the corresponding tidally-varying salinity case (i.e., salinity varying between 0 and 1) produced somewhat comparable results in terms of RUC and RLC (RUC represents groundwater discharge to the estuary that originated from recharge to the estuary bank and RLC groundwater discharge to the estuary that originated from recharge through the estuary bed), whereas flow paths and the total salt mass in the aquifer differed. Freshwater flushing simulations indicated that the near-estuary aquifer responds rapidly to a 2-day ‘wet season’ flushing event with a short-lived freshwater lens created through freshening of the hyporheic zone. Annual cycling of the seasonal flushing led to significant disruption of the estuary water circulation in the aquifer thereby impacting on residence times, transport pathways, and RUC and RLC, and acting to potentially remobilize groundwater and contaminants previously trapped in continuous and semi-continuous re-circulation cells. Although groundwater flow paths determined using tide-averaged velocity vectors were representative of flow paths from transient tidally driven flow vector field, residence times calculated from the two flow fields were markedly different. The influence of riparian scale aquitards and lateral asymmetry (about the estuary centreline) in hydraulic gradients and hydraulic conductivity on groundwater flow and associated salinity distribution was also found to be sensitive to estuarine salinity conditions. The results indicate that observations made about aquifer-estuary interaction in the lower estuary may not be directly applicable to the middle estuary. According to the simulations, tidal salinity variations in the estuary are important factors that affect hyporheic-riparian salt transport processes and that the use of a time averaged estuarine salinity as an approximation to variable salinity conditions is unsuitable for the accurate prediction of the near-estuary dynamics in middle estuaries. This study was based on a two dimensional representation of the riparian scale interaction and it is clear that future research needs to focus on the three-dimensionality of the aquifer-estuary system, incorporating spatially and temporally varying flow and transport characteristics. That is, many estuaries are tortuous and the aquifer geology spatially complex such that assumptions required for the two-dimensional section will most likely restrict application to the field. The tidal dynamics in the middle estuary is also expected to generate three dimensional aspects to the aquifer-estuary interaction. Thus further investigation that explicitly models the hydrodynamics and salt transport in the estuary and estuarine morphology is required to refine the insight provided by the simple conceptual model adopted in this study.

Estuarine ecosystems are under enormous stress due to rapid coastal developments and climate change. Proper management of these important ecosystems requires a good understanding of their key processes. In this thesis, riparian groundwater-surface water interaction is explored for an aquifer-estuary system primarily by a series of numerical experiments. The work focuses on riparian-scale groundwater flow and salinization. The overall aim of the study was to extend our understanding of aquifer-estuary exchange, which is currently centered on the lower marine estuarine reach, to middle estuaries (i.e., the estuary reach that has variable salinity). The numerical experiments were guided by previous studies and observations made from an exploratory field investigation conducted in and next to Sandy Creek, a macro-tidal estuary incised in the alluvial aquifer of the Pioneer Valley, North-eastern Australia (Longitude 49.11°, Latitude -21.27°). The following observations were made from the field investigation: Sandy Creek estuary experiences a variable salinity regime in its mid reaches that consists of periods of 1) freshwater flushing due to up catchment-derived flooding, 2) persistent freshwater conditions for at least 2 months following the flooding, 3) tidal salinity fluctuations and 4) constant near-seawater salinity; laterally extensive and disconnected aquitards were found to occur at the field site; Sandy Creek had an essentially ‘vertical’ bank slope. Numerical simulations were conducted using the finite element modeling code FEFLOW for saturated unsaturated, variable-density groundwater flow and solute transport, to examine the influence of the following factors on aquifer-estuary exchange: a tidally varying estuarine salinity and hydraulic head, a seasonal freshwater flush (i.e., estuary with freshwater and an elevated stage due to an up catchment sourced flood), near estuary aquitard layers, lateral asymmetry (about the estuary centerline) in hydraulic conductivity and regional hydraulic gradients. The simulations neglected seepage face development after numerical experiments showed that for a vertical bank estuary interacting with a sandy loam aquifer, seepage face effects on groundwater flow and associated salinity distribution were minimal. The following observations were drawn from the range of numerical experiments considered. Tidal salinity fluctuations in the estuary (varying between 0 and 1 - i.e., using a relative salinity scale where a salinity of 1 is seawater) produced flow paths and residence times that were distinctly different to the constant seawater salinity case. While the constant average 0.5 salinity case and the corresponding tidally-varying salinity case (i.e., salinity varying between 0 and 1) produced somewhat comparable results in terms of RUC and RLC (RUC represents groundwater discharge to the estuary that originated from recharge to the estuary bank and RLC groundwater discharge to the estuary that originated from recharge through the estuary bed), whereas flow paths and the total salt mass in the aquifer differed. Freshwater flushing simulations indicated that the near-estuary aquifer responds rapidly to a 2-day ‘wet season’ flushing event with a short-lived freshwater lens created through freshening of the hyporheic zone. Annual cycling of the seasonal flushing led to significant disruption of the estuary water circulation in the aquifer thereby impacting on residence times, transport pathways, and RUC and RLC, and acting to potentially remobilize groundwater and contaminants previously trapped in continuous and semi-continuous re-circulation cells. Although groundwater flow paths determined using tide-averaged velocity vectors were representative of flow paths from transient tidally driven flow vector field, residence times calculated from the two flow fields were markedly different. The influence of riparian scale aquitards and lateral asymmetry (about the estuary centreline) in hydraulic gradients and hydraulic conductivity on groundwater flow and associated salinity distribution was also found to be sensitive to estuarine salinity conditions. The results indicate that observations made about aquifer-estuary interaction in the lower estuary may not be directly applicable to the middle estuary. According to the simulations, tidal salinity variations in the estuary are important factors that affect hyporheic-riparian salt transport processes and that the use of a time averaged estuarine salinity as an approximation to variable salinity conditions is unsuitable for the accurate prediction of the near-estuary dynamics in middle estuaries. This study was based on a two dimensional representation of the riparian scale interaction and it is clear that future research needs to focus on the three-dimensionality of the aquifer-estuary system, incorporating spatially and temporally varying flow and transport characteristics. That is, many estuaries are tortuous and the aquifer geology spatially complex such that assumptions required for the two-dimensional section will most likely restrict application to the field. The tidal dynamics in the middle estuary is also expected to generate three dimensional aspects to the aquifer-estuary interaction. Thus further investigation that explicitly models the hydrodynamics and salt transport in the estuary and estuarine morphology is required to refine the insight provided by the simple conceptual model adopted in this study.

Estuarine ecosystems are under enormous stress due to rapid coastal developments and climate change. Proper management of these important ecosystems requires a good understanding of their key processes. In this thesis, riparian groundwater-surface water interaction is explored for an aquifer-estuary system primarily by a series of numerical experiments. The work focuses on riparian-scale groundwater flow and salinization. The overall aim of the study was to extend our understanding of aquifer-estuary exchange, which is currently centered on the lower marine estuarine reach, to middle estuaries (i.e., the estuary reach that has variable salinity). The numerical experiments were guided by previous studies and observations made from an exploratory field investigation conducted in and next to Sandy Creek, a macro-tidal estuary incised in the alluvial aquifer of the Pioneer Valley, North-eastern Australia (Longitude 49.11°, Latitude -21.27°). The following observations were made from the field investigation: Sandy Creek estuary experiences a variable salinity regime in its mid reaches that consists of periods of 1) freshwater flushing due to up catchment-derived flooding, 2) persistent freshwater conditions for at least 2 months following the flooding, 3) tidal salinity fluctuations and 4) constant near-seawater salinity; laterally extensive and disconnected aquitards were found to occur at the field site; Sandy Creek had an essentially ‘vertical’ bank slope. Numerical simulations were conducted using the finite element modeling code FEFLOW for saturated unsaturated, variable-density groundwater flow and solute transport, to examine the influence of the following factors on aquifer-estuary exchange: a tidally varying estuarine salinity and hydraulic head, a seasonal freshwater flush (i.e., estuary with freshwater and an elevated stage due to an up catchment sourced flood), near estuary aquitard layers, lateral asymmetry (about the estuary centerline) in hydraulic conductivity and regional hydraulic gradients. The simulations neglected seepage face development after numerical experiments showed that for a vertical bank estuary interacting with a sandy loam aquifer, seepage face effects on groundwater flow and associated salinity distribution were minimal. The following observations were drawn from the range of numerical experiments considered. Tidal salinity fluctuations in the estuary (varying between 0 and 1 - i.e., using a relative salinity scale where a salinity of 1 is seawater) produced flow paths and residence times that were distinctly different to the constant seawater salinity case. While the constant average 0.5 salinity case and the corresponding tidally-varying salinity case (i.e., salinity varying between 0 and 1) produced somewhat comparable results in terms of RUC and RLC (RUC represents groundwater discharge to the estuary that originated from recharge to the estuary bank and RLC groundwater discharge to the estuary that originated from recharge through the estuary bed), whereas flow paths and the total salt mass in the aquifer differed. Freshwater flushing simulations indicated that the near-estuary aquifer responds rapidly to a 2-day ‘wet season’ flushing event with a short-lived freshwater lens created through freshening of the hyporheic zone. Annual cycling of the seasonal flushing led to significant disruption of the estuary water circulation in the aquifer thereby impacting on residence times, transport pathways, and RUC and RLC, and acting to potentially remobilize groundwater and contaminants previously trapped in continuous and semi-continuous re-circulation cells. Although groundwater flow paths determined using tide-averaged velocity vectors were representative of flow paths from transient tidally driven flow vector field, residence times calculated from the two flow fields were markedly different. The influence of riparian scale aquitards and lateral asymmetry (about the estuary centreline) in hydraulic gradients and hydraulic conductivity on groundwater flow and associated salinity distribution was also found to be sensitive to estuarine salinity conditions. The results indicate that observations made about aquifer-estuary interaction in the lower estuary may not be directly applicable to the middle estuary. According to the simulations, tidal salinity variations in the estuary are important factors that affect hyporheic-riparian salt transport processes and that the use of a time averaged estuarine salinity as an approximation to variable salinity conditions is unsuitable for the accurate prediction of the near-estuary dynamics in middle estuaries. This study was based on a two dimensional representation of the riparian scale interaction and it is clear that future research needs to focus on the three-dimensionality of the aquifer-estuary system, incorporating spatially and temporally varying flow and transport characteristics. That is, many estuaries are tortuous and the aquifer geology spatially complex such that assumptions required for the two-dimensional section will most likely restrict application to the field. The tidal dynamics in the middle estuary is also expected to generate three dimensional aspects to the aquifer-estuary interaction. Thus further investigation that explicitly models the hydrodynamics and salt transport in the estuary and estuarine morphology is required to refine the insight provided by the simple conceptual model adopted in this study.

In this study, The Beach Dewatering System, a relatively recent technology to combat beach erosion, which is proposed as a practical alternative to more traditional shoreline stabilization methods, is investigated and an informative overview on the genesis, development and recent use of this technique is provided. On the basis of the link existing between the elevation of beach groundwater and erosional or accretionary trends at the beach face, a numerical model that simulates groundwater flow in a coastal aquifer under beach drainage is presented. In this model, the seaward boundary of the domain is considered to be tidally fluctuating in a large scale to represent the occurrence of seepage face significantly. The unsteady groundwater flow equation is solved numerically using the method of finite differences. The results clearly showed that the water table being lowered caused the reduction of the seepage face which is the main aim of Beach Dewatering projects. The positional design parameters, i.e. horizontal and vertical location of the drain, are also investigated by utilizing an efficiency index. It is observed that the system efficiency decreased as the drain is shifted landward. The results also indicated that, the efficiency slightly increased with the vertical drain elevation.

Non-overflow flow-through rockfill structures are river engineering elements used to attenuate and delay inflow hydrographs. They represent expedient places to deposit rather enormous quantities of waste rock at mountainous mine sites. Their application has become so common that matters of safety regarding their design have been laid out in Section 8.5 of the Canadian Dam Safety Guidelines (CDA 2007). The research described herein was directed at investigating the different aspects of the hydraulics of these flow-through rockfill structures. In order to assess the potential for an unraveling failure of flow-through rockfill dams, a systematic study of the hydraulic design of these structures was conducted and the non-linear nature of flow through these structures was dealt with using a p-LaPlacian-like partial differential equation. Subsequently, factors of safety against this type of failure are presented for a range of downstream slopes, thus showing the unsafe combinations of embankment slope and particle diameter. Three different index gradients within the toe of such structures were investigated. In this regard, the gradient most suitable for independently computing the height of the point of first flow emergence on the downstream face is examined and a method for independently computing the variation in hydraulic head within that vertical (which allows for the toe of the structure to be isolated) is presented. An additional gradient that allows for the independent estimation of the default tailwater depth is proposed. In order to provide better tools to assess the behavior of these embankments at the toe, laboratory and analytical studies were undertaken. In this regard, the hydraulics associated with the zone of the downstream toe were studied. The depth variation of the seepage-face was computationally modeled, and two approaches for solving the spatially varied flow (SVF) condition problem within the toe region undertaken. The results show that a dual linear variation in depth can be used to good accuracy, without inducing any unrealistic exit gradients in the zone of primary concern with respect to unraveling. It is hoped that these techniques and computational tools provided herein will aid in facilitating the design and assessment of these flow-through rockfill structures.