Assessment of groundwater quality recovery strategies using nitrate transport modelling. Application to the Saône alluvial formations (Tournus, Saône-et-Loire)

To assess the efficiency of different agro-environmental strategies used to reduce groundwater pollution by nitrates, transport modelling in soils and groundwater has been carried out on two withdrawal areas in an alluvial plain. In a first time, the agro-environmental model AgriFlux allowed the simulation of water and nitrates fluxes flowing to groundwater. This model was calibrated for each agro-pedological unit of the studied territory. In a second time, the application of the hydrogeological model MODFLOW-MT3D allowed the simulation of nitrate transport in groundwater for the 1980-2004 period. This soil-groundwater coupled modelling has shown that soil nature is the first factor that conditions the vulnerability to nitrates. Thus, nitrate leaching occurs preferentially under sandy soils. Efficiency of different agro-environmental operations for groundwater quality recovery was quantified. The best results are obtained by combination of (1) grassland re-installation on sandy agricultural lots located in near well protection perimeter and (2) fertilization reduction on sandy agricultural lots located in the well alimentation area upstream the near protection perimeter. On other soils, the effect of grassland on groundwater quality improvement is more limited. Nevertheless, the control of nitrate fertilisation remains essential and is justified in both near and far well protection perimeters. Modelling thus allows optimising and priorizing agro-environmental actions in alluvial agricultural zones. [Comte J.-C., Banton O., Kockmann F., Villard A., Creuzot G. (2006), Assessment of groundwater quality recovery strategies using nitrate transport modelling. Application to the Saône alluvial formations (Tournus, Saône-et-Loire), Ingénieries Eau-Agriculture-Territoires, 45, 15-28]

Comparison and visualisation of high-resolution transport modelling with aircraft measurements

This article presents a case study of a comparison of an Eulerian chemical transport model (CTM) and Lagrangian chemical model with measurements taken by aircraft. High-resolution Eulerian integrations produce improved point-by-point comparisons between model results and measurements compared to low resolution. The Lagrangian model requires mixing to be introduced in order to model the measurements.

Calibration Strategies Designed For Long Term Sediment Transport Modelling

The presented work deals with the calibration of a 2D numerical model for the simulation of long term bed load transport. A settled basin along an alpine stream was used as a case study. The focus is to parameterise the used multi fractional transport model such that a dynamically balanced behavior regarding erosion and deposition is reached. The used 2D hydrodynamic model utilizes a multi-fraction multi-layer approach to simulate morphological changes and bed load transport. The mass balancing is performed between three layers: a top mixing layer, an intermediate subsurface layer and a bottom layer. Using this approach bears computational limitations in calibration. Due to the high computational demands, the type of calibration strategy is not only crucial for the result, but as well for the time required for calibration. Brute force methods such as Monte Carlo type methods may require a too large number of model runs. All here tested calibration strategies used multiple model runs utilising the parameterization and/or results from previous run. One concept was to reset to initial bed elevations after each run, allowing the resorting process to convert to stable conditions. As an alternative or in combination, the roughness was adapted, based on resulting nodal grading curves, from the previous run. Since the adaptations are a spatial process, the whole model domain is subdivided in homogeneous sections regarding hydraulics and morphological behaviour. For a faster optimization, the adaptation of the parameters is made section wise. Additionally, a systematic variation was done, considering results from previous runs and the interaction between sections. The used approach can be considered as similar to evolutionary type calibration approaches, but using analytical links instead of random parameter changes.

Transport modelling: macro and micro simulation for the studied case of Funchal

The work done in this thesis attempts to demonstrate the importance of using models that can predict and represent the mobility of our society. To answer the proposed challenges two models were examined, the first corresponds to macro simulation with the intention of finding a solution to the frequency of the bus company Horários do Funchal, responsible for transport in the city of Funchal, and some surrounding areas. Where based on a simplified model of the city it was possible to increase the frequency of journeys getting an overall reduction in costs. The second model concerns the micro simulation of Avenida do Mar, where currently is being built a new roundabout (Praça da Autonomia), which connects with this avenue. Therefore it was proposed to study the impact on local traffic, and the implementation of new traffic lights for this purpose. Four possible situations in which was seen the possibility of increasing the number of lanes on the roundabout or the insertion of a bus lane were created. The results showed that having a roundabout with three lanes running is the best option because the waiting queues are minimal, and at environmental level this model will project fewer pollutants. Thus, this thesis presents two possible methods of urban planning. Transport modelling is an area that is under constant development, the global goal is to encourage more and more the use of these models, and as such it is important to have more people to devote themselves to studying new ways of addressing current problems, so that we can have more accurate models and increasing their credibility.

Biogeochemical processes in a clay formation in situ experiment: Part F - Reactive transport modelling

Reactive transport modelling was used to simulate solute transport, thermodynamic reactions, ion exchange and biodegradation in the Porewater Chemistry (PC) experiment at the Mont Terri Rock Laboratory. Simulations show that the most important chemical processes controlling the fluid composition within the borehole and the surrounding formation during the experiment are ion exchange, biodegradation and dissolution/precipitation reactions involving pyrite and carbonate minerals. In contrast, thermodynamic mineral dissolution/precipitation reactions involving alumo-silicate minerals have little impact on the fluid composition on the time-scale of the experiment. With the accurate description of the initial chemical condition in the formation in combination with kinetic formulations describing the different stages of bacterial activities, it has been possible to reproduce the evolution of important system parameters, such as the pH, redox potential, total organic C. dissolved inorganic C and SO(4) concentration. Leaching of glycerol from the pH-electrode may be the primary source of organic material that initiated bacterial growth, which caused the chemical perturbation in the borehole. Results from these simulations are consistent with data from the over-coring and demonstrate that the Opalinus Clay has a high buffering capacity in terms of chemical perturbations caused by bacterial activity. This buffering capacity can be attributed to the carbonate system as well as to the reactivity of clay surfaces.

Reactive mass transport modelling of a three-dimensional vertical fault zone with a finger-like convective flow regime

For a three-dimensional vertically-oriented fault zone, we consider the coupled effects of fluid flow, heat transfer and reactive mass transport, to investigate the patterns of fluid flow, temperature distribution, mineral alteration and chemically induced porosity changes. We show, analytically and numerically, that finger-like convection patterns can arise in a vertically-oriented fault zone. The onset and patterns of convective fluid flow are controlled by the Rayleigh number which is a function of the thermal properties of the fluid and the rock, the vertical temperature gradient, and the height and the permeability of the fault zone. Vigorous fluid flow causes low temperature gradients over a large region of the fault zone. In such a case, flow across lithological interfaces becomes the most important mechanism for the formation of sharp chemical reaction fronts. The degree of rock buffering, the extent and intensity of alteration, the alteration mineralogy and in some cases the formation of ore deposits are controlled by the magnitude of the flow velocity across these compositional interfaces in the rock. This indicates that alteration patterns along compositional boundaries in the rock may provide some insights into the convection pattern. The advective mass and heat exchanges between the fault zone and the wallrock depend on the permeability contrast between the fault zone and the wallrock. A high permeability contrast promotes focussed convective flow within the fault zone and diffusive exchange of heat and chemical reactants between the fault zone and the wallrock. However, a more gradual permeability change may lead to a regional-scale convective flow system where the flow pattern in the fault affects large-scale fluid flow, mass transport and chemical alteration in the wallrocks