@article{ErdalBaroniSanchezLeonetal.2019,
  author    = {Erdal, Daniel and Baroni, Gabriele and S{\´a}nchez Le{\´o}n, Eduardo Emilio and Cirpka, Olaf A.},
  title     = {The value of simplified models for spin up of complex models with an application to subsurface hydrology},
  series = {Computers \& geosciences : an international journal devoted to the publication of papers on all aspects of geocomputation and to the distribution of computer programs and test data sets ; an official journal of the International Association for Mathematical Geology},
  volume    = {126},
  journal   = {Computers \& geosciences : an international journal devoted to the publication of papers on all aspects of geocomputation and to the distribution of computer programs and test data sets ; an official journal of the International Association for Mathematical Geology},
  publisher = {Elsevier},
  address   = {Oxford},
  issn      = {0098-3004},
  doi       = {10.1016/j.cageo.2019.01.014},
  pages     = {62 -- 72},
  year      = {2019},
  abstract  = {Spinning up large-scale coupled surface-subsurface numerical models can be a time and resource consuming task. If an uninformed initial condition is chosen, the spin-up can easily require 20 years of repeated simulations on high-performance computing machines. In this paper we compare the classical approach of starting from a fixed shallow depth to groundwater (here 3 m) with three more informed approaches for the definition of initial conditions in the spin up. In the first of these three approaches, we start from a known-steady state groundwater table, calculated with a 2-D groundwater model and the yearly net recharge, and combine it with an unsaturated zone that assumes hydrostatic conditions. In the second approach, we start from the same groundwater table combined with vertical profiles in the unsaturated zone with uniform vertical flow identical to the groundwater recharge. In the third approach we calculate a dynamic steady state from a simplified subsurface model combining a transient 2-D groundwater model with a limited number of 1-D transient unsaturated zone columns on top. Results for spinning-up a 3-D Parflow-CLM model using the different initial conditions show that large gains can be made by considering states in groundwater and the vadose zone that are consistent, i.e. where groundwater recharge and the vertical flux in the vadose zone agree. By this, the spin-up time was reduced from about 10 years to about 3 years of simulated time. In the light of seasonal fluctuations of net recharge, using the transient approach showed more stable results.},
  language  = {en}
}
@article{DeBiaseCarminatiOswaldetal.2013,
  author    = {De Biase, Cecilia and Carminati, Andrea and Oswald, Sascha Eric and Thullner, Martin},
  title     = {Numerical modeling analysis of VOC removal processes in different aerobic vertical flow systems for groundwater remediation},
  series = {Journal of contaminant hydrology},
  volume    = {154},
  journal   = {Journal of contaminant hydrology},
  number    = {11},
  publisher = {Elsevier},
  address   = {Amsterdam},
  issn      = {0169-7722},
  doi       = {10.1016/j.jconhyd.2013.07.007},
  pages     = {53 -- 69},
  year      = {2013},
  abstract  = {Vertical flow systems filled with porous medium have been shown to efficiently remove volatile organic contaminants (VOCs) from contaminated groundwater. To apply this semi-natural remediation strategy it is however necessary to distinguish between removal due to biodegradation and due to volatile losses to the atmosphere. Especially for (potentially) toxic VOCs, the latter needs to be minimized to limit atmospheric emissions. In this study, numerical simulation was used to investigate quantitatively the removal of volatile organic compounds in two pilot-scale water treatment systems: an unplanted vertical flow filter and a planted one, which could also be called a vertical flow constructed wetland, both used for the treatment of contaminated groundwater. These systems were intermittently loaded with contaminated water containing benzene and MTBE as main VOCs. The highly dynamic but permanently unsaturated conditions in the porous medium facilitated aerobic biodegradation but could lead to volatile emissions of the contaminants. Experimental data from porous material analyses, flow rate measurements, solute tracer and gas tracer test, as well as contaminant concentration measurements at the boundaries of the systems were used to constrain a numerical reactive transport modeling approach. Numerical simulations considered unsaturated water flow, transport of species in the aqueous and the gas phase as well as aerobic degradation processes, which made it possible to quantify the rates of biodegradation and volatile emissions and calculating their contribution to total contaminant removal. A range of degradation rates was determined using experimental results of both systems under two operation modes and validated by field data obtained at different operation modes applied to the filters. For both filters, simulations and experimental data point to high biodegradation rates, if the flow filters have had time to build up their removal capacity. For this case volatile emissions are negligible and total removal can be attributed to biodegradation, only. The simulation study thus supports the use of both of these vertical flow systems for the treatment of groundwater contaminated with VOCs and the use of reactive transport modeling for the assessment of VOCs removal and operation modes in these high performance treatment systems.},
  language  = {en}
}