This research has been undertaken to further our knowledge of decade-to-century timescale braided, gravel-bed river delta growth dynamics. The study included: a review of available literature; field studies; the development of microscale models for two study deltas; and the development of a simple numerical model incorporating movement of braided river channels across a delta topset (varying the location of sediment delivery to the delta). Results from the microscale modelling showed that successful physical modelling requires well-defined fixed boundaries and, ideally, good historical aerial photography for the estimation of the model time scale. A complex braided gravel-bed river delta system composed of two merging deltas entering a deep, low-energy receiving basins was able to be successfully modelled to provide valuable information on delta growth dynamics. However, a microscale model of a delta prograding into shallow receiving basins, with a large supply of fine sediment, was more difficult to calibrate and assess (partly due to limited field data), and was considered less reliable. The simple rule-based numerical model ‘DELGROW’, developed to simulate a braided river system entering a deep, low-energy body of water, requires a known sediment supply rate, as well as information on the braided river topography, submerged delta foreset, and lakebed bathymetry. Unlike simple 1-d width-averaged geometric models, DELGROW takes into consideration barriers (e.g. islands) as well as relatively complex converging braided river delta configurations. By changing the sediment supply, or modifying the river system, the response of the river system to various scenarios can also be assessed. Microscale models and DELGROW appear to realistically simulate decade-to-century timescale growth of braided gravel-bed river deltas entering a deep, low-energy, receiving basin. Both of these modelling methods initially use the supplied sediment to try and eliminate any riverbed irregularities (e.g. low areas), before continuing to advance and deposit sediment in a more evenly-distributed manner, whilst taking into consideration irregularities due to barriers, and asymmetric sediment sources such as merging deltas. Neither model can reliably predict locations of bank erosion, or channel avulsions that divert flow and sediment outside of the fixed model boundaries.

A thermal conductive heating (TCH) pilot test was conducted at the Naval Air Warfare Center (NAWC) in West Trenton, New Jersey in 2009 in collaboration with TerraTherm, Inc., the Naval Facilities Engineering Services Center and the United States Geological Survey. The NAWC site was historically used as a jet engine testing facility from the mid-1950s to the late 1990s. During this time, the subsurface was contaminated with trichloroethylene (TCE) which was a common solvent used at the facility. The pilot test consisted of 15 heater/extraction wells installed to a depth of 16.8 m in weathered mudstone and operated for 102 days. Rock core samples were taken pre- and post-remediation to measure the initial TCE concentrations and evaluate the effect the TCH pilot test had. The data collected during the pilot test was used to create a two-dimensional (2D) finite difference model using TMVOC. TMVOC is part of the TOUGH 2 family of codes and is a numerical model that is capable of simulating multiphase flow, heat transfer and transport of volatile organic compounds in three-dimensional heterogenous porous media or fractured rock. The 2D model was used as a screening model to investigate TCE removal from the rock matrix when heating for 100 days with a similar heating pattern to what was employed at the NAWC site. The numerical domain incorporated three primary fractures with competent bedrock in between. As the test pilot was conducted in the weathered bedrock zone, a sensitivity analysis was first completed on the matrix permeability to help to match the TCE removal from the pilot test. The pilot test had a 63.5% removal of TCE from the study area compared to 67% from the baseline model. A limited sensitivity analysis was completed which investigated how the matrix porosity and rate of energy application would have on the success of TCE removal from the rock matrix. It revealed that the TCE removal increases with increased matrix porosity and increased rate of energy application. / Thesis (Master, Civil Engineering) -- Queen's University, 2013-01-02 22:32:05.274

In this thesis, I focused on numerical modelling of groundwater flow in the vicinity of groundwater circulation well (GCW). The thesis can be devided to two parts. In the first one, I created a simple numerical model of GCW. I used a cylindrical flow domain because of the radial symmetry of the problem. The task was to find out which parameters have the biggest influence on groundwater flow. It was found that the hydrogeological settings and hydraulic properties have more considerable effect on GW flow than the well construction parameters. Distance between the well screens is the only important construction parameter, other parameters can be neglected. However, we cannot neglect the influence of hydraulic properties of the porous medium. Considering the radius of influence of GCW, the presence of aquitards and anisotropy of layers affect the groundwater flow around the well seriously. These two parameters cause the extension of the radius of influence, which leads to a longer path of the particle of water. The travel time of particles increases too this slows down the remediation. A big influence has also an intensity of natural groundwater flow in the area because it decreases the radius of influence of the well. In the second part I created a numerical model of pilot site Velká Hleďsebe to investigate...

Mountain watersheds recently burned by wildfire often experience greater amounts of runoff and increased rates of sediment transport relative to similar unburned areas. Given the sedimentation and debris flow threats caused by increases in erosion, more work is needed to better understand the physical mechanisms responsible for the observed increase in sediment transport in burned environments and the time scale over which a heightened geomorphic response can be expected. In this study, we quantified the relative importance of different hillslope erosion mechanisms during two postwildfire rainstorms at a drainage basin in Southern California by combining terrestrial laser scanner-derived maps of topographic change, field measurements, and numerical modeling of overland flow and sediment transport. Numerous debris flows were initiated by runoff at our study area during a long-duration storm of relatively modest intensity. Despite the presence of a well-developed rill network, numerical model results suggest that the majority of eroded hillslope sediment during this long-duration rainstorm was transported by raindrop-induced sediment transport processes, highlighting the importance of raindrop-driven processes in supplying channels with potential debris flow material. We also used the numerical model to explore relationships between postwildfire storm characteristics, vegetation cover, soil infiltration capacity, and the total volume of eroded sediment from a synthetic hillslope for different end-member erosion regimes. This study adds to our understanding of sediment transport in steep, postwildfire landscapes and shows how data from field monitoring can be combined with numerical modeling of sediment transport to isolate the processes leading to increased erosion in burned areas.

Numerische Experimente wurden durchgeführt, um den Einfluß von Gittermaschenweite und subskaliger Heterogenität auf die Vorhersage der am Wasserkreislauf beteiligten Größen zu untersuchen. Die Modellergebnisse zeigen, daß die Evapotranspiration, Bewölkung und der Niederschlag von der Gittermaschenweite und der Heterogenität beeinflußt werden. Es zeigte sich, daß bei Verwendung gröberer Maschenweiten unter Einbezug der verschiedenen Landnutzungstypen innerhalb der Gittermasche die Obertlächenprozesse und Phänomene (z.B. Wärmeinseleffekt) realistischer beschrieben werden, als wenn nur ein Landnutzungstyp für das gesamte Gitterelement angenommen wird. / Numerical experiments were performed to investigate the influence of grid resolution and subgrid heterogeneity on the prediction of the quantities of the water cycle. The results were compared with each other and with those provided by a simulation using the same surface parameterization scheme but taking subgrid scale surface heterogeneity into account. The model results substantiate that the evapotranspiration, cloudiness and precipitation are affected by the grid resolution and the heterogeneity. lt was found that increasing the grid size but including the heterogeneity describes more realistically the surface processes and phenomena (e.g„ heat island effect) than assuming one land use type for the whole grid element.

No description available.

The hypoxic area in the Gulf of Mexico is the second largest in the world, which has received extensive scientific study and management interest. Previous modeling studies have concluded that the increased hypoxic area in the Gulf of Mexico was caused by the increased anthropogenic nitrogen loading of the Mississippi River; however, the nitrogen-area relationship is complicated by many other factors, such as wind, river discharge, and the ratio of Mississippi to Atchafalaya River flow. These factors are related to large-scale climate variability, and thus will not be affected by regional nitrogen reduction efforts. In the research presented here, both statistical (regression) and numerical models are used to study the influence of anthropogenic and climate factors on the hypoxic area variability in the Gulf of Mexico. The numerical model is a three-dimensional, coupled hydrological-biogeochemical model (ROMS-Fennel). Results include: (1) the west wind duration during the summer explain 55% of the hypoxic area variability since 1993. Combined wind duration and nitrogen loading explain over 70% of the variability, and combined wind duration and river discharge explain over 85% of the variability. (2) The numerical model captures the temporal variability, but overestimates the bottom oxygen concentrations. The model shows that the simulated hypoxic area is in agreement with the observations from the year 1991, as long as hypoxia is defined as oxygen concentrations below 3 mg/L rather than below 2 mg/L. (3) The first three modes from an Empirical Orthogonal Function (EOF) analysis of the numerical model output results explain 62%, 8.1% and 4.9% of the variability of the hypoxic area. The Principle Component time series is cross-correlated with wind, dissolved inorganic nitrogen concentration and river discharge. (4) Scenario experiments with the same nitrogen loading, but different duration of upwelling favorable wind, indicate that the upwelling favorable wind is important for hypoxic area development. However, a long duration of upwelling wind decreases the area. (5) Scenario experiments with the same nitrogen loading, but different discharges, indicate that increasing river discharge by 50% increases the area by 42%. Additionally, scenario experiments with the same river discharge, but different nitrogen concentrations, indicate that reducing the nitrogen concentration by 50% decreases the area by 75%. (6) Scenario experiments with the same nitrogen loading, but different flow diver- sions, indicate that if the Atchafalaya River discharges increased to 66.7%, the total hypoxic area increases the hypoxic area by 30%, and most of the hypoxic area moved from east to west Louisiana shelf. Additionally, if the Atchafalaya River discharge decreased to zero, the total hypoxic area increases by 13%. (7) Scenario experiments with the same nitrogen loading, but different nitrogen forms, indicate that if all the nitrogen was in the inorganic forms, the hypoxic area increases by 15%. These results have multiple implications for understanding the mechanisms that control the oxygen dynamics, reevaluating management strategies, and improving the observational methods.

In view of the foreign mariculture is gradually diversified, and even has a trend that the fish cage aquaculture combines with algae culture at the same facility. A submersible single-point-mooring (SPM) cage system was modified and installed in-situ to investigate the feasibility of the cage system. A numercial model was established to simulate the cage dynamic motion as well as the mooring line tension. A detailed cage construction process was described in this study. Tension meter was used to keep track of mooring line tension, while the ADCP( Acoustic Doppler Current Profiler) was utilized to record the sea state during the test period. Those data were used to validate the numerical model. The field experiements were carried out at a location north to Xiaoliuqiu island. Since the testing period was in winter, the wave height was relativly calm and found to be between 0.5 and 1.2 m, wave period 4~7 seconds, and wave current about 0.2~0.6 m/s. The numerical results indicate that the maximum mooring line tension has good agreement with the meauresments of the tension meter. These comparisons verify that this numerical model is sufficient to simulate this kind of alage cage systems.

Corpus Christi Bay (TX, USA) is a shallow wind-driven bay which is designated as a National Estuary due to its impact on the economy. But this bay experiences periodic hypoxia (dissolved oxygen <2 mg/l) which threatens aerobic aquatic organisms. Development of the Coastal Margin Observation and Assessment System (CMOAS) through integration of real-time observations with numerical modeling helps to understand the processes causing hypoxia in this energetic bay. CMOAS also serves as a template for the implementation of observational systems in other dynamic ecosystems for characterizing and predicting other episodic events such as harmful algal blooms, accidental oil spills, sediment resuspension events, etc. State-of-the-art sensor technologies are involved in real-time monitoring of hydrodynamic, meteorological and water quality parameters in the bay. Three different platform types used for the installation of sensor systems are: 1) Fixed Robotic, 2) Mobile, and 3) Remote. An automated profiler system, installed on the fixed robotic platform, vertically moves a suite of in-situ sensors within the water column for continuous measurements. An Integrated Data Acquisition, Communication and Control system has been configured on our mobile platform (research vessel) for the synchronized measurements and real-time visualization of hydrodynamic and water quality parameters at greater spatial resolution. In addition, a high frequency (HF) radar system has been installed on remote platforms to generate surface current maps for Corpus Christi (CC) Bay and its offshore area. This data is made available to stakeholders in real-time through the development of cyberinfrastructure which includes establishment of communication network, software development, web services, database development, etc. Real-time availability of measured datasets assists in implementing an integrated sampling scheme for our monitoring systems installed at different platforms. With our integrated system, we were able to capture evidence of an hypoxic event in Summer 2007. Data collected from our monitoring systems are used to drive and validate numerical models developed in this study. The analysis of observational datasets and developed 2-D hydrodynamic model output suggests that a depth-integrated model is not able to capture the water current structure of CC Bay. Also, the development of a threedimensional mechanistic dissolved oxygen model and a particle aggregation transport model (PAT) helps to clarify the critical processes causing hypoxia in the bay. The various numerical models and monitoring systems developed in this study can serve as valuable tools for the understanding and prediction of various episodic events dominant in other dynamic ecosystems.

A 3D numerical model ¡V COHERENS has been applied to construct a coupled hydrodynamic and ecological model for studying Tapeng Bay, which is a coastal lagoon situated in southwest of Taiwan. The simulations have been carried out to study the influences and their interacting mechanisms among the tidal currents, nutrients and micro planktons in the Lagoon. Model results have been compiled for calculating the nutrient fluxes and the primary productions in the Tapeng Bay. Tapeng Bay is a semi-enclosed coastal lagoon, which has only one tidal inlet for exchanging lagoon water with the coastal currents along the Kaoping coast on the narrow shelf in southwest of Taiwan. The study area is situated in the tropical climate zone where has sunshine through out the year except the rainy days concentrated in the summer season, which is influenced by the southwest monsoon. There are several drainage channels that collect the untreated domestic sewerage and wastewater discharged from the fish farms surround the lagoon. The discharges in these channels are usually low during the dry season. The solid contained in the water are mostly settled on the channel beds. During the raining season, high discharges due to the storm rainfalls re-suspend the sediments and carry into the lagoon. These sediments, which contain high concentrations of suspended solids and nutrients, cause the Bay water highly eutrophied. Therefore, the Bay is fully influenced by the seasonal variations. There are a lot of aquaculture, i.e. oyster farming and fish cage, in the Bay area since the water is calm and rich. But the balance between the nature and the anthropogenic disturbance is breaking. Besides the water level variation generated from the tidal inlet, the fresh water inflow from 3 major channels are included in the model to simulate their influences to the hydrodynamics and the density driven circulation due to changing salinities and temperatures from these inlets. Plankton, detritus, dissolved nutrients and dissolved oxygen is taking into account as the model variables for this marine eco-system. The plankton growth is mainly generated due to temperature, light intensity and nutrient level. Only the nitrogen cycle has been considered in the model by assuming there are enough supply of phosphate and silicate. Model runs have been carried out according to different seasonal situations of the boundary conditions. Furthermore, climates (heats, lights, winds, etc) are also included in the model to distinct seasonal characteristics. It is shown, from the model results, that the currents mainly dominate the distribution of nutrients in the Tapeng Bay. The nutrient level controls plankton growth. The nutrient sources are mainly coming from the coastal currents (through tidal inlet) in the wintertime, whereas the summer source was from the drainage channels due to the wash out by the high discharge rates. Beside these, dissolved oxygen concentrations in the Bay water are strongly influenced by the plankton growth rate, faster the photosynthesis higher the DO concentrations. The eutrophication levels of the Tapeng Bay water have been compiled using the plankton carbon level modeled at various situations. According to the Nixon standard (1995), Tapeng Bay has eutrophication through out the year. Mesotrophic condition can be observed during the wintertime, whereas the hypereutrophic level can be concluded during the raining season.