scholarly journals Estimation of effective porosity in large-scale groundwater models by combining particle tracking, auto-calibration and <sup>14</sup>C dating

2018 ◽  
Author(s):  
Rena Meyer ◽  
Peter Engesgaard ◽  
Klaus Hinsby ◽  
Jan A. Piotrowski ◽  
Torben O. Sonnenborg

Abstract. Effective porosity plays an important role in contaminant management. However, the effective porosity is often assumed constant in space and hence either neglected or simplified in transport model calibration. Based on a calibrated highly parametrized flow model, a three-dimensional advective transport model (MODPATH) of a 1300 km2-large coastal area of southern Denmark and northern Germany is presented. A detailed voxel model represents the highly heterogeneous geological composition of the area. Inverse modelling of advective transport is used to estimate seven, spatially distributed, effective porosity units based on apparent groundwater ages inferred from 11 14C measurements in Pleistocene and Miocene aquifers, corrected for the effects of diffusion and geochemical reactions. The match between the observed and simulated ages is improved significantly by the calibration of seven porosity units by a reduction of ME of 99 % and RMS of 82 % compared to a uniform porosity approach. Groundwater ages range from a few hundred years in the Pleistocene to several thousand years in Miocene aquifers. The advective age distributions derived from particle tracking at each sampling well show unimodal (for younger ages) to multimodal (for older ages) shapes and thus reflect the heterogeneity that particles encounter along their travel path. The estimated effective porosity field, with values ranging between 4.3 % in clay and 45 % in sand formations, is used in a direct simulation of distributed mean groundwater ages. Although the absolute ages are affected by various uncertainties, a unique insight into the complex three-dimensional age distribution pattern and potential advance of young contaminated groundwater in the investigated regional aquifer system is provided, highlighting the importance of estimating effective porosity in groundwater transport modelling and the implications for groundwater quantity and quality assessment and management.

2018 ◽  
Vol 22 (9) ◽  
pp. 4843-4865 ◽  
Author(s):  
Rena Meyer ◽  
Peter Engesgaard ◽  
Klaus Hinsby ◽  
Jan A. Piotrowski ◽  
Torben O. Sonnenborg

Abstract. Effective porosity plays an important role in contaminant management. However, the effective porosity is often assumed to be constant in space and hence heterogeneity is either neglected or simplified in transport model calibration. Based on a calibrated highly parametrized flow model, a three-dimensional advective transport model (MODPATH) of a 1300 km2 coastal area of southern Denmark and northern Germany is presented. A detailed voxel model represents the highly heterogeneous geological composition of the area. Inverse modelling of advective transport is used to estimate the effective porosity of 7 spatially distributed units based on apparent groundwater ages inferred from 11 14C measurements in Pleistocene and Miocene aquifers, corrected for the effects of diffusion and geochemical reactions. By calibration of the seven effective porosity units, the match between the observed and simulated ages is improved significantly, resulting in a reduction of ME of 99 % and RMS of 82 % compared to a uniform porosity approach. Groundwater ages range from a few hundred years in the Pleistocene to several thousand years in Miocene aquifers. The advective age distributions derived from particle tracking at each sampling well show unimodal (for younger ages) to multimodal (for older ages) shapes and thus reflect the heterogeneity that particles encounter along their travel path. The estimated effective porosity field, with values ranging between 4.3 % in clay and 45 % in sand formations, is used in a direct simulation of distributed mean groundwater ages. Although the absolute ages are affected by various uncertainties, a unique insight into the complex three-dimensional age distribution pattern and potential advance of young contaminated groundwater in the investigated regional aquifer system is provided, highlighting the importance of estimating effective porosity in groundwater transport modelling and the implications for groundwater quantity and quality assessment and management.


2011 ◽  
Vol 11 (1) ◽  
pp. 363-373 ◽  
Author(s):  
H. Bencherif ◽  
L. El Amraoui ◽  
G. Kirgis ◽  
J. Leclair De Bellevue ◽  
A. Hauchecorne ◽  
...  

Abstract. This paper reports on an increase of ozone event observed over Kerguelen (49.4° S, 70.3° E) in relationship with large-scale isentropic transport. This is evidenced by ground-based observations (co-localised radiosonde and SAOZ experiments) together with satellite global observations (Aura/MLS) assimilated into MOCAGE, a Méteo-France model. The study is based on the analyses of the first ozonesonde experiment never recorded at the Kerguelen site within the framework of a French campaign called ROCK that took place from April to August 2008. Comparisons and interpretations of the observed event are supported by co-localised SAOZ observations, by global mapping of tracers (O3, N2O and columns of O3) from Aura/MLS and Aura/OMI experiments, and by model simulations of Ertel Potential Vorticity initialised by the ECMWF (European Centre for Medium-Range Weather Forecasts) data reanalyses. Satellite and ground-based observational data revealed a consistent increase of ozone in the local stratosphere by mid-April 2008. Additionally, Ozone (O3) and nitrous oxide (N2O) profiles obtained during January–May 2008 using the Microwave Limb Sounder (MLS) aboard the Aura satellite are assimilated into MOCAGE (MOdèle de Chimie Atmosphérique à Grande Echelle), a global three-dimensional chemistry transport model of Météo-France. The assimilated total O3 values are consistent with SAOZ ground observations (within ±5%), and isentropic distributions of O3 match well with maps of advected potential vorticity (APV) derived from the MIMOSA model, a high-resolution advection transport model, and from the ECMWF reanalysis. The event studied seems to be related to the isentropic transport of air masses that took place simultaneously in the lower- and middle-stratosphere, respectively from the polar region and from the tropics to the mid-latitudes. In fact, the ozone increase observed by mid April 2008 resulted simultaneously: (1) from an equator-ward departure of polar air masses characterised with a high-ozone layer in the lower stratosphere (near the 475 K isentropic level), and (2) from a reverse isentropic transport from the tropics to mid- and high-latitudes in the upper stratosphere (nearby the 700 K level). The increase of ozone observed over Kerguelen from the 16-April ozonesonde profile is thus attributed to a concomitant isentropic transport of ozone in two stratospheric layers: the tropical air moving southward and reaching over Kerguelen in the upper stratosphere, and the polar air passing over the same area but in the lower stratosphere.


2009 ◽  
Vol 9 (5) ◽  
pp. 18511-18543 ◽  
Author(s):  
J. Aschmann ◽  
B. M. Sinnhuber ◽  
E. L. Atlas ◽  
S. M. Schauffler

Abstract. The transport of very short-lived substances into the tropical upper troposphere and lower stratosphere is investigated by a three-dimensional chemical transport model using archived convective updraft mass fluxes (or detrainment rates) from the European Centre for Medium-Range Weather Forecast's ERA-Interim reanalysis. Large-scale vertical velocities are calculated from diabatic heating rates. With this approach we explicitly model the large scale subsidence in the tropical troposphere with convection taking place in fast and isolated updraft events. The model calculations agree generally well with observations of bromoform and methyl iodide from aircraft campaigns and with ozone and water vapor from sonde and satellite observations. Using a simplified treatment of dehydration and bromine product gas washout we give a range of 1.6 to 3 ppt for the contribution of bromoform to stratospheric bromine, assuming a uniform source in the boundary layer of 1 ppt. We show that the most effective region for VSLS transport into the stratosphere is the West Pacific, accounting for about 55% of the bromine from bromoform transported into the stratosphere under the supposition of a uniformly distributed source.


2009 ◽  
Vol 9 (3) ◽  
pp. 13889-13916 ◽  
Author(s):  
A. Voulgarakis ◽  
O. Wild ◽  
N. H. Savage ◽  
G. D. Carver ◽  
J. A. Pyle

Abstract. We use a three-dimensional chemical transport model to examine the shortwave radiative effects of clouds on the tropospheric ozone budget. In addition to looking at changes in global concentrations as previous studies have done, we examine changes in ozone chemical production and loss caused by clouds and how these vary in different parts of the troposphere. On a global scale, we find that clouds have a modest effect on ozone chemistry, but on a regional scale their role is much more significant, with the size of the response dependent on the region. The largest averaged changes in chemical budgets (±10–14%) are found in the marine troposphere, where cloud optical depths are high. We demonstrate that cloud effects are small on average in the middle troposphere because this is a transition region between reduction and enhancement in photolysis rates. We show that increases in boundary layer ozone due to clouds are driven by large-scale changes in downward ozone transport from higher in the troposphere rather than by decreases in in-situ ozone chemical loss rates. Increases in upper tropospheric ozone are caused by higher production rates due to backscattering of radiation and consequent increases in photolysis rates, mainly J(NO2). The global radiative effect of clouds on isoprene is stronger than on ozone. Tropospheric isoprene lifetime increases by 7% when taking clouds into account. We compare the importance of clouds in contributing to uncertainties in the global ozone budget with the role of other radiatively-important factors. The budget is most sensitive to the overhead ozone column, while surface albedo and clouds have smaller effects. However, uncertainty in representing the spatial distribution of clouds may lead to a large sensitivity on regional scales.


Author(s):  
Cheng He ◽  
Padala Chittibabu ◽  
Matthew Graham

a three-dimensional hydrodynamic and mud transport modelling system MIKE 3-FM is used to study the flow circulation and the sediment plume fate in the Randle Reef (RR) area of Hamilton Harbour.It was found that the simulated plume extents can be very sensitive to transport model mesh resolution because of the additional numerical advection induced in the calculations. Detailed discussion of this particular issue has not been observed in other publications.The results of the simulation provide useful information for remediation project managers to plan and guide the environmental monitoring of dredging operations. In addition, the methodology used in this study can be adapted to other dredging projects.


2010 ◽  
Vol 10 (7) ◽  
pp. 17727-17751
Author(s):  
H. Bencherif ◽  
L. El Amraoui ◽  
G. Kirgis ◽  
J. Leclair De Bellevue ◽  
A. Hauchecorne ◽  
...  

Abstract. This paper reports on an increase of ozone event observed over Kerguelen (49.4° S, 70.3° E) in relationship with large-scale isentropic transport. It is evidenced from ground-based observations, together with satellite global observations and assimilated fields. The study is based on the analyses of the first ozonesonde experiment never recorded at the Kerguelen site in the framework of a French campaign called ROCK that took place from April to August 2008. Comparisons and interpretations of the observed event are supported by co-localised SAOZ observations, by global mapping of tracers (O3, N2O and columns of O3) from Aura/MLS and Aura/OMI experiments, and by model simulations of Ertel Potential Vorticity initialised by ECMWF (European Centre for Medium-Range Weather Forecasts) data reanalyses. Satellite and ground-based observational data revealed a consistent increase of ozone in the local stratosphere by mid-April 2008. Additionally, Ozone (O3) and nitrous oxide (N2O) profiles obtained during January–May 2008 by the Microwave Lamb Sounder (MLS) aboard the Aura satellite are assimilated into MOCAGE (MOdèle de Chimie Atmosphérique à Grande Echelle), a global three-dimensional chemistry transport model of Météo-France. The assimilated total O3 values are consistent with SAOZ ground observations (within ±5%), and isentropic distributions of O3 are matching well with maps of advected potential vorticity (APV) derived from the MIMOSA model, a high-resolution advection transport model, and from ECMWF reanalysis. The studied event seems to be related to isentropic transport of air masses that took place simultaneously in the lower- and middle-stratosphere, respectively from the polar region and from tropics to the mid-latitudes. In fact, the studied ozone increase by mid April 2008 results simultaneously: (1) from an equator-ward departure of polar air masses characterised with a high-ozone layer in the lower stratosphere (nearby the 475 K isentropic level), and (2) from a reverse isentropic transport from tropics to mid- and high-latitudes in the upper stratosphere (nearby the 700 K level). The increase of ozone observed over Kerguelen from the 16-April ozonesonde profile is then attributed to a concomitant isentropic transport of ozone in two stratospheric layers: the tropical air moving southward and reaches over Kerguelen in the upper stratosphere, and the polar air passing over the same area but in the lower stratosphere.


2009 ◽  
Vol 9 (21) ◽  
pp. 8235-8246 ◽  
Author(s):  
A. Voulgarakis ◽  
O. Wild ◽  
N. H. Savage ◽  
G. D. Carver ◽  
J. A. Pyle

Abstract. We use a three-dimensional chemical transport model to examine the shortwave radiative effects of clouds on the tropospheric ozone budget. In addition to looking at changes in global concentrations as previous studies have done, we examine changes in ozone chemical production and loss caused by clouds and how these vary in different parts of the troposphere. On a global scale, we find that clouds have a modest effect on ozone chemistry, but on a regional scale their role is much more significant, with the size of the response dependent on the region. The largest averaged changes in chemical budgets (±10–14%) are found in the marine troposphere, where cloud optical depths are high. We demonstrate that cloud effects are small on average in the middle troposphere because this is a transition region between reduction and enhancement in photolysis rates. We show that increases in boundary layer ozone due to clouds are driven by large-scale changes in downward ozone transport from higher in the troposphere rather than by decreases in in-situ ozone chemical loss rates. Increases in upper tropospheric ozone are caused by higher production rates due to backscattering of radiation and consequent increases in photolysis rates, mainly J(NO2). The global radiative effect of clouds on isoprene, through decreases of OH in the lower troposphere, is stronger than on ozone. Tropospheric isoprene lifetime increases by 7% when taking clouds into account. We compare the importance of clouds in contributing to uncertainties in the global ozone budget with the role of other radiatively-important factors. The budget is most sensitive to the overhead ozone column, while surface albedo and clouds have smaller effects. However, uncertainty in representing the spatial distribution of clouds may lead to a large sensitivity of the ozone budget components on regional scales.


2019 ◽  
Vol 8 (2) ◽  
Author(s):  
Valentina Vincenzi ◽  
Michele Sapigni

In the framework of the final design of a new hydroelectric plant close to the existing dam on the Adige River in the municipality of San Giovanni Lupatoto (VR, Italy), the design of deep excavations in loose sediments below groundwater level, and the connected risks of piping and excavation bottom heave required a detailed study of hydrogeological characterization, aimed to an accurate aquifer parameterization and so a confident design of excavation supports and dewatering systems. Particularly a three-dimensional analysis of the problem was necessary, because the hydrogeological system was complicated by many “artificial” boundary conditions (dam on the river, deep foundations of the existing and new engineering works) causing vertical not negligible flow paths. The conceptual model has been developed applying, in parallel, different prospecting tools: a preliminary Step Drawdown Test, a long duration Constant Rate Test, a multiple tracing test and the implementation of a three-dimensional numerical model of groundwater flow (MODFLOW-2005 code) and advective transport (MODPATH code). The modeling has allowed a good estimate of the aquifer parameters such as hydraulic conductivity and effective porosity through calibration of the pumping/tracing tests under transient conditions. Moreover, the simulation of post operam conditions and the comparison with ante operam, in terms of flow paths length and filtration velocity, showed that the construction of the new works foundations would have determined a sensitive reduction of groundwater velocities closed to the existing dam.


Water ◽  
2020 ◽  
Vol 12 (5) ◽  
pp. 1480
Author(s):  
Xingwei Liu ◽  
Qiulan Zhang ◽  
Tangpei Cheng

To overcome the large time and memory consumption problems in large-scale high-resolution contaminant transport simulations, an efficient approach was presented to parallelize the modular three-dimensional transport model for multi-species (MT3DMS) (University of Alabama, Tuscaloosa, AL, USA) program on J adaptive structured meshes applications infrastructures (JASMIN). In this approach, a domain decomposition method and a stencil-based method were used to accomplish parallel implementation, while a ghost cell strategy was used for communication. The MODFLOW-MT3DMS coupling mode was optimized to achieve the parallel coupling of flow and contaminant transport. Five types of models were used to verify the correctness and test the parallel performance of the method. The developed parallel program JMT3D (China University of Geosciences (Beijing), Beijing, China) can increase the speed by up to 31.7 times, save memory consumption by 96% with 46 processors, and ensure that the solution accuracy and convergence do not decrease as the number of domains increases. The BiCGSTAB (Bi-conjugate gradient variant algorithm) method required the least amount of time and achieved high speedup in most cases. Coupling the flow and contaminant transport further improved the efficiency of the simulations, with a 33.45 times higher speedup achieved on 46 processors. The AMG (algebraic multigrid) method achieved a good scalability, with an efficiency above 100% on hundreds of processors for the simulation of tens of millions of cells.


Author(s):  
Khodyar Javadi ◽  
Mohammad Taeibi-Rahni ◽  
Masoud Darbandi

This work is conducted with evaluation of different turbulence models capabilities in predicting three dimensional jet-into-crossflow (JICF) interactions. For this purpose, first of all, comprehensive discussions on the near wall flow complexities due to discharge of a jet into a crossflow are presented. In this regards, large scale coherent structures such as: counter rotating vortex pairs (CRVP’s), near wall secondary motions, horseshoe vortices, and wall jets like are discussed. Secondly, the abilities of different turbulence models in predicting such flows (JICF) are evaluated. Our evaluation is based on three points of view including: 1) JICF characteristics, 2) computed location, and 3) sensitivity to different flow variables. In this regard, the turbulence models such as k-ε, k-ω, shear stress transport model (SST), and Reynolds stress model (RSM) are employed. Their related results are compared to credential available experimental/numerical data as well themselves. Since the same basic code with the same grid density as well as numerical discretization scheme is used, it is save to conclude that, any differences in the results are due to the abilities of turbulence models. The flow field computation was governed by Reynolds Averaged Navier-Stokes (RANS) equations performing finite volume method with SIMPLE algorithm over a non-uniform structured grid.


Sign in / Sign up

Export Citation Format

Share Document