T.J. Heimovaara
Please Note
25 records found
1
Although numerous projects worldwide have shown encouraging results, landfill stabilization remains challenged by the inherent heterogeneity of the waste body. This complexity limits uniform treatment and leaves uncertainties about the physical, chemical, and biological interactions at play. Addressing these knowledge gaps, this thesis investigates the effectiveness of aeration and water recirculation in three Dutch pilot landfills: Braambergen and Wieringermeer (aerated), and Kragge (water recirculation).
The pilots revealed that the effects of aeration are highly variable in space and time. At Braambergen, variability in aeration performance revealed the strong influence of site heterogeneity. Differences in water levels in aeration wells affected gas composition and flow, yet high water columns alone could not explain the observed contrast between compartments. Other factors, such as spatial variability in gas permeability within the waste body, also played a role. Where aeration was more effective, higher gas extraction, elevated temperatures, and greater settlement indicated enhanced microbial activity and carbon mineralization.
Beyond gas monitoring, stabilization was assessed by comparing the carbon generation of waste samples under aerobic and anaerobic conditions with model predictions and with carbon actually recovered on-site. The heterogeneity of the waste samples was reflected in the carbon potential and decay rate constants (k-values). Aerated pilots showed reduced aerobic carbon potential, reflecting advanced stabilization, while the recirculated pilot retained substantial degradable organic matter. These results highlight both the large potential of aeration to accelerate stabilization and the persistence of heterogeneity that complicates prediction and management.
A further focus was placed on building a comprehensive carbon and nitrogen balance across the solid, aqueous, and gas phases at field scale. Over seven years, aerated pilots exhibited higher organic matter degradation than the anaerobic pilot with a significant share of carbon and nitrogen released through the gas phase. In contrast, the recirculated pilot retained larger amounts of degradable carbon and poorly mobilizable nitrogen. Importantly, the analysis revealed that a substantial fraction of nitrogen remains fixed in solid or microbial pools, potentially delaying compliance with leachate emission targets.
Taken together, these findings advance understanding of how aeration and water recirculation influence landfill stabilization. They demonstrate the benefits of aeration for accelerating degradation while also underlining the challenges posed by spatial variability and persistent nitrogen pools. Such insights are crucial for improving the design and implementation of in-situ stabilization strategies and for reducing the long-term aftercare needs of landfills. ...
Although numerous projects worldwide have shown encouraging results, landfill stabilization remains challenged by the inherent heterogeneity of the waste body. This complexity limits uniform treatment and leaves uncertainties about the physical, chemical, and biological interactions at play. Addressing these knowledge gaps, this thesis investigates the effectiveness of aeration and water recirculation in three Dutch pilot landfills: Braambergen and Wieringermeer (aerated), and Kragge (water recirculation).
The pilots revealed that the effects of aeration are highly variable in space and time. At Braambergen, variability in aeration performance revealed the strong influence of site heterogeneity. Differences in water levels in aeration wells affected gas composition and flow, yet high water columns alone could not explain the observed contrast between compartments. Other factors, such as spatial variability in gas permeability within the waste body, also played a role. Where aeration was more effective, higher gas extraction, elevated temperatures, and greater settlement indicated enhanced microbial activity and carbon mineralization.
Beyond gas monitoring, stabilization was assessed by comparing the carbon generation of waste samples under aerobic and anaerobic conditions with model predictions and with carbon actually recovered on-site. The heterogeneity of the waste samples was reflected in the carbon potential and decay rate constants (k-values). Aerated pilots showed reduced aerobic carbon potential, reflecting advanced stabilization, while the recirculated pilot retained substantial degradable organic matter. These results highlight both the large potential of aeration to accelerate stabilization and the persistence of heterogeneity that complicates prediction and management.
A further focus was placed on building a comprehensive carbon and nitrogen balance across the solid, aqueous, and gas phases at field scale. Over seven years, aerated pilots exhibited higher organic matter degradation than the anaerobic pilot with a significant share of carbon and nitrogen released through the gas phase. In contrast, the recirculated pilot retained larger amounts of degradable carbon and poorly mobilizable nitrogen. Importantly, the analysis revealed that a substantial fraction of nitrogen remains fixed in solid or microbial pools, potentially delaying compliance with leachate emission targets.
Taken together, these findings advance understanding of how aeration and water recirculation influence landfill stabilization. They demonstrate the benefits of aeration for accelerating degradation while also underlining the challenges posed by spatial variability and persistent nitrogen pools. Such insights are crucial for improving the design and implementation of in-situ stabilization strategies and for reducing the long-term aftercare needs of landfills.
In this dissertation, I present an integrated framework to estimate and predict landfill emission potentials by combining stochastic modeling, Bayesian uncertainty quantification, data assimilation, and hydrogeophysical measurements. The research introduces a stochastic Lagrangian-based travel time modeling approach to simulate the heterogeneous water flow and solute transport within landfill bodies. This method, unlike traditional grid-based models, captures preferential flow phenomena and accommodates the spatial variability inherent in landfill waste structures.
The model calibration is performed using Bayesian inference, employing long-term observational data of leachate production and quality from the Braambergen landfill in the Netherlands. This probabilistic calibration explicitly quantifies uncertainties in model parameters and outputs, providing more credible risk assessments and long-term predictions of leachate emissions.
Recognizing the risk of error accumulation in history-matching methods, I further implement data assimilation techniques, notably the Weakly Coupled Particle Filter (WCPF) and a hybrid Particle Filter–Markov Chain Monte Carlo (PF-MCMC) method. These approaches enable sequential updating of model parameters and system states as new data become available, improving the predictive performance and reducing uncertainty over time. The PF-MCMC method, in particular, can estimate parameters and hidden processes, which is very helpful for understanding the dynamics in the landfill.
To further enhance the accuracy of emission potential estimations, the framework integrates hydrogeophysical data obtained through Electrical Resistivity Tomography (ERT). Using a Bayesian evidential learning approach, resistivity measurements are directly mapped into probabilistic water storage estimates within landfill waste bodies. This additional constraint strengthens the characterization of subsurface hydrological conditions, distinguishing between leachable and isolated water fractions.
The dissertation is structured across six chapters, beginning with an overview of the landfill aftercare problem, followed by the development of the stochastic modeling framework, the application of particle filtering and PF-MCMC, the incorporation of ERT data through Bayesian evidential learning, and concluding with a synthesis of findings and recommendations for future research.
Overall, this work advances the scientific understanding of landfill emission dynamics by offering a unified methodological framework that integrates stochastic modeling, data assimilation, and hydrogeophysical surveying. The contributions herein support the development of more robust, data-driven, and cost-effective strategies for landfill aftercare, ensuring long-term environmental protection and sustainability. ...
In this dissertation, I present an integrated framework to estimate and predict landfill emission potentials by combining stochastic modeling, Bayesian uncertainty quantification, data assimilation, and hydrogeophysical measurements. The research introduces a stochastic Lagrangian-based travel time modeling approach to simulate the heterogeneous water flow and solute transport within landfill bodies. This method, unlike traditional grid-based models, captures preferential flow phenomena and accommodates the spatial variability inherent in landfill waste structures.
The model calibration is performed using Bayesian inference, employing long-term observational data of leachate production and quality from the Braambergen landfill in the Netherlands. This probabilistic calibration explicitly quantifies uncertainties in model parameters and outputs, providing more credible risk assessments and long-term predictions of leachate emissions.
Recognizing the risk of error accumulation in history-matching methods, I further implement data assimilation techniques, notably the Weakly Coupled Particle Filter (WCPF) and a hybrid Particle Filter–Markov Chain Monte Carlo (PF-MCMC) method. These approaches enable sequential updating of model parameters and system states as new data become available, improving the predictive performance and reducing uncertainty over time. The PF-MCMC method, in particular, can estimate parameters and hidden processes, which is very helpful for understanding the dynamics in the landfill.
To further enhance the accuracy of emission potential estimations, the framework integrates hydrogeophysical data obtained through Electrical Resistivity Tomography (ERT). Using a Bayesian evidential learning approach, resistivity measurements are directly mapped into probabilistic water storage estimates within landfill waste bodies. This additional constraint strengthens the characterization of subsurface hydrological conditions, distinguishing between leachable and isolated water fractions.
The dissertation is structured across six chapters, beginning with an overview of the landfill aftercare problem, followed by the development of the stochastic modeling framework, the application of particle filtering and PF-MCMC, the incorporation of ERT data through Bayesian evidential learning, and concluding with a synthesis of findings and recommendations for future research.
Overall, this work advances the scientific understanding of landfill emission dynamics by offering a unified methodological framework that integrates stochastic modeling, data assimilation, and hydrogeophysical surveying. The contributions herein support the development of more robust, data-driven, and cost-effective strategies for landfill aftercare, ensuring long-term environmental protection and sustainability.
Coupling ORCHESTRA to Python
A modelling framework for simulating waste degradation in a bioreactor
mol/L) affect the hydrolysis rate by toxicity inhibition. This translates to a model that is sensitive to a variety of input parameters such as the initial C/N ratio, the Kd (degree of adsorption) value and the initial organic fraction. A comparison with literature studies, however, implies that the model lacks fundamental processes such as gas and fluid transport or changes in volume and temperature. In addition to this, we discuss that the Freundlich isotherms may be inappropriate to capture adsorption in a coupled framework as no ion exchange is taken into account. Models that do so, such as the NICA-Donnan model, could improve the model results significantly. ...
mol/L) affect the hydrolysis rate by toxicity inhibition. This translates to a model that is sensitive to a variety of input parameters such as the initial C/N ratio, the Kd (degree of adsorption) value and the initial organic fraction. A comparison with literature studies, however, implies that the model lacks fundamental processes such as gas and fluid transport or changes in volume and temperature. In addition to this, we discuss that the Freundlich isotherms may be inappropriate to capture adsorption in a coupled framework as no ion exchange is taken into account. Models that do so, such as the NICA-Donnan model, could improve the model results significantly.
The leachate levels indicated perched leachate zones: above the basal drainage system and below the injection drains at the top. Leachate injected through the infiltration drains cannot efficiently infiltrate the waste body, and the effects of the infiltration events were not picked up in the wells and piezometers throughout the landfill, implying little hydraulic connectivity.
The results of the dilution tests indicated horizontal and vertical flow within the landfill. Vertical velocities measured in the wells were estimated to be considerably higher (77 - 225 m/d) than the average horizontal velocities (0.02 - 1.0 m/d). The wells provide a path for vertical flow. Apart from the highest horizontal velocities measured in deeper sections of the landfill (15-18 m below ground level), velocities varied without a clear relation to landfill depth, indicating preferential flow paths. Uniform dilution tests performed with the infiltration drains turned off suggested that leachate infiltration does not increase the horizontal velocities. This research suggests that waste stabilization through recirculation is not optimal at De Kragge.
\noindent Due to an overall lack of understanding about the well and filter pack installation and the high spatial heterogeneity of the waste, the calculated velocities are uncertain. Further research into possible error sources (e.g., the borehole correction factor) is recommended. Additional tracer tests, including tests on the neighboring compartment without stabilization measures, are recommended to further assess the effectiveness of the recirculation system. ...
The leachate levels indicated perched leachate zones: above the basal drainage system and below the injection drains at the top. Leachate injected through the infiltration drains cannot efficiently infiltrate the waste body, and the effects of the infiltration events were not picked up in the wells and piezometers throughout the landfill, implying little hydraulic connectivity.
The results of the dilution tests indicated horizontal and vertical flow within the landfill. Vertical velocities measured in the wells were estimated to be considerably higher (77 - 225 m/d) than the average horizontal velocities (0.02 - 1.0 m/d). The wells provide a path for vertical flow. Apart from the highest horizontal velocities measured in deeper sections of the landfill (15-18 m below ground level), velocities varied without a clear relation to landfill depth, indicating preferential flow paths. Uniform dilution tests performed with the infiltration drains turned off suggested that leachate infiltration does not increase the horizontal velocities. This research suggests that waste stabilization through recirculation is not optimal at De Kragge.
\noindent Due to an overall lack of understanding about the well and filter pack installation and the high spatial heterogeneity of the waste, the calculated velocities are uncertain. Further research into possible error sources (e.g., the borehole correction factor) is recommended. Additional tracer tests, including tests on the neighboring compartment without stabilization measures, are recommended to further assess the effectiveness of the recirculation system.
Urban Soils as Hotspots for Carbon Storage
A case study of the urban greenspaces of The Hague, the Netherlands
Although urban soil can be highly disturbed or altered by anthropogenic activities, the high C densities in The Hague suggested that its potential to store C appeared unaffected. Along the transect, a mean C density of 88 t/ha, of which 82 t/ha was considered organic C, was detected, which was higher than the values currently assigned to urban soils in national C inventories. The urban soil C storage was dependent on the type of vegetation, urbanization extent and land ownership. The hypothesized links between land use and soil type were not apparent in this case study, suggesting that processes driving soil C storage are controlled by different factors.
The total soil C storage of the upper 30 cm of the greenspaces in The Hague was estimated at 18.8 kt of C. The use of high spatial resolution GIS data with a scale of 10 x 10 m enabled the inclusion of small patches in the total soil C storage of The Hague, which proved to be significant as the smaller urban greenspaces, which are typical for dense urban centres, contained similar soil C density as the larger urban greenspaces, such as urban forests.
Soil C storage in urban ecosystems is highly variable. How generalizable these results are across other Dutch cities requires further research. Moreover, to translate current soil C stocks to annual C fluxes further research is required. This study found that urban soil C stocks are underestimated, which potentially also is the case for urban soil sequestration rates that are currently applied in C modelling studies. ...
Although urban soil can be highly disturbed or altered by anthropogenic activities, the high C densities in The Hague suggested that its potential to store C appeared unaffected. Along the transect, a mean C density of 88 t/ha, of which 82 t/ha was considered organic C, was detected, which was higher than the values currently assigned to urban soils in national C inventories. The urban soil C storage was dependent on the type of vegetation, urbanization extent and land ownership. The hypothesized links between land use and soil type were not apparent in this case study, suggesting that processes driving soil C storage are controlled by different factors.
The total soil C storage of the upper 30 cm of the greenspaces in The Hague was estimated at 18.8 kt of C. The use of high spatial resolution GIS data with a scale of 10 x 10 m enabled the inclusion of small patches in the total soil C storage of The Hague, which proved to be significant as the smaller urban greenspaces, which are typical for dense urban centres, contained similar soil C density as the larger urban greenspaces, such as urban forests.
Soil C storage in urban ecosystems is highly variable. How generalizable these results are across other Dutch cities requires further research. Moreover, to translate current soil C stocks to annual C fluxes further research is required. This study found that urban soil C stocks are underestimated, which potentially also is the case for urban soil sequestration rates that are currently applied in C modelling studies.
The control of infiltration and seepage of water is one ofthe most challenging tasks in water management and civil-engineering and, inan attempt to control this, methods for forming a waterimpermeable layer inthe soil have been widely practised in soil engineering (Laumann et al., 2018;Proto et al., 2016). The use of natural processes to modify the engineeringproperties of the subsurface could help to develop cost-effective, robust andsustainable engineering technologies and is attracting increasing attentionfrom the industry (Zhou, 2020). This research aims to reduce the permeabilityby using aluminium (Al) and organic matter (OM) precipitates mixed ex-situwith porous media to create a horizontal barrier. The Al-OM precipitates wereexpected to clog the pore space with a reduction in permeability and hydraulicconductivity as result. To find out if it is feasible to use Al-OMprecipitates mixed ex-situ with porous media for a permeability reductivelayer, the Al-OM precipitates and the permeability of the medium wereresearched. To characterise the flocculation reaction, experiments onthe flocculation were performed. The yield of the reaction was obtained byadding certain amounts of Al and OM solutions to form particular amounts of drymass of flocs. This experiment confirmed the hypothesis that 85% of the addedmass of OM will result in dry mass of flocs. The concept of a critical metal tocarbon ratio (M/C-ratio), indicating flocculation regardless of the inputconcentrations, was tested by measuring the pH over an increasing M/C-ratio.From this titration curve, the found critical M/C ratio is between 0.0230.031,and the pH stabilises at a level lower than pH 4. This result proves theconcept behind the numerical scenarios describing the titration of OM solutionswith Al3+, with a critical molar M/C ratio independent of the inputconcentrations (Veerkamp, 2018; Zhou, 2020). After determining the yield of thereaction and concluding that the concentration of Al and OM was not ofinfluence, the by-products of the Al-OM reaction were quantified. The ionicstrength of the supernatant of an increasing density of flocs in solution wasdetermined by measuring the electrical conductivity (EC). The results showedthat the ionic strength increased linearly with an increasing density of flocs.The relationship between the ionic strength and the density of flocs was coupledto the linear relationship between concentrations potassium chloride (KCl) andits EC. From the results, the measured EC can be used as a tracer since the K+and Cl are nonreactive. The hydraulic conductivity measurements were conducted by afalling head test to be able to make an indication of the change inpermeability when adding the Al-OM precipitates to the sand. To find theoptimal method to mix the Al, OM and porous media, the influence of differentmethods of producing, adding and mixing the materials on the permeabilityreduction is explored. In the first mixing method, the Al and OM were added insolution, the solution containing Al-OM precipitates was centrifuged until thereduced ionic strength was at an EC value less than 700 µm/cm. The hydraulicconductivity measurements were used to obtain the relationship between thehydraulic conductivity reduction over an increasing concentration of flocsretained by one kilogram of sand. The hydraulic conductivity was reducedexponentially over an increasing concentration of flocs up to a magnitude of 3.The results imply a large variability in the achieved reduction dominated bythe amount of retained flocs. The second method is using Al and OM in powderformat and adding them to the dry sand and adding 500 ml of water to thismixture. This mixing method resulted in a completely different floc structure.The flocs produced by mixing in solution have a shear dependency feature, whiledry mixing created particles that have a constant size. For this method, theincrease of the reduction is linear over an increase of concentration of flocsretained by the soil. The highest reduction for this method was found to be ofa magnitude of two, measured at 50 grams of flocs retained by one kilogram ofsand. Finally, this research gave proof of principle of using Al-OMprecipitates mixed directly with sand could reduce the permeability up to amagnitude of 3. These results present a new road to research on this Al-OM-sandmixture’s strength parameters and compaction over time over an increasing flocdensity, since these parameters are critical for using the layer in practice. ...
The control of infiltration and seepage of water is one ofthe most challenging tasks in water management and civil-engineering and, inan attempt to control this, methods for forming a waterimpermeable layer inthe soil have been widely practised in soil engineering (Laumann et al., 2018;Proto et al., 2016). The use of natural processes to modify the engineeringproperties of the subsurface could help to develop cost-effective, robust andsustainable engineering technologies and is attracting increasing attentionfrom the industry (Zhou, 2020). This research aims to reduce the permeabilityby using aluminium (Al) and organic matter (OM) precipitates mixed ex-situwith porous media to create a horizontal barrier. The Al-OM precipitates wereexpected to clog the pore space with a reduction in permeability and hydraulicconductivity as result. To find out if it is feasible to use Al-OMprecipitates mixed ex-situ with porous media for a permeability reductivelayer, the Al-OM precipitates and the permeability of the medium wereresearched. To characterise the flocculation reaction, experiments onthe flocculation were performed. The yield of the reaction was obtained byadding certain amounts of Al and OM solutions to form particular amounts of drymass of flocs. This experiment confirmed the hypothesis that 85% of the addedmass of OM will result in dry mass of flocs. The concept of a critical metal tocarbon ratio (M/C-ratio), indicating flocculation regardless of the inputconcentrations, was tested by measuring the pH over an increasing M/C-ratio.From this titration curve, the found critical M/C ratio is between 0.0230.031,and the pH stabilises at a level lower than pH 4. This result proves theconcept behind the numerical scenarios describing the titration of OM solutionswith Al3+, with a critical molar M/C ratio independent of the inputconcentrations (Veerkamp, 2018; Zhou, 2020). After determining the yield of thereaction and concluding that the concentration of Al and OM was not ofinfluence, the by-products of the Al-OM reaction were quantified. The ionicstrength of the supernatant of an increasing density of flocs in solution wasdetermined by measuring the electrical conductivity (EC). The results showedthat the ionic strength increased linearly with an increasing density of flocs.The relationship between the ionic strength and the density of flocs was coupledto the linear relationship between concentrations potassium chloride (KCl) andits EC. From the results, the measured EC can be used as a tracer since the K+and Cl are nonreactive. The hydraulic conductivity measurements were conducted by afalling head test to be able to make an indication of the change inpermeability when adding the Al-OM precipitates to the sand. To find theoptimal method to mix the Al, OM and porous media, the influence of differentmethods of producing, adding and mixing the materials on the permeabilityreduction is explored. In the first mixing method, the Al and OM were added insolution, the solution containing Al-OM precipitates was centrifuged until thereduced ionic strength was at an EC value less than 700 µm/cm. The hydraulicconductivity measurements were used to obtain the relationship between thehydraulic conductivity reduction over an increasing concentration of flocsretained by one kilogram of sand. The hydraulic conductivity was reducedexponentially over an increasing concentration of flocs up to a magnitude of 3.The results imply a large variability in the achieved reduction dominated bythe amount of retained flocs. The second method is using Al and OM in powderformat and adding them to the dry sand and adding 500 ml of water to thismixture. This mixing method resulted in a completely different floc structure.The flocs produced by mixing in solution have a shear dependency feature, whiledry mixing created particles that have a constant size. For this method, theincrease of the reduction is linear over an increase of concentration of flocsretained by the soil. The highest reduction for this method was found to be ofa magnitude of two, measured at 50 grams of flocs retained by one kilogram ofsand. Finally, this research gave proof of principle of using Al-OMprecipitates mixed directly with sand could reduce the permeability up to amagnitude of 3. These results present a new road to research on this Al-OM-sandmixture’s strength parameters and compaction over time over an increasing flocdensity, since these parameters are critical for using the layer in practice.
The project 'introductie Duurzaam Stortbeheer (iDS) aims to create aftercare with an ending, or finite aftercare. It is looking for possibilities to stabilise the landfill mass, i.e. to eliminate the threat of pollutants by treating the wastebodies.
The goal of the CURE project is to develop fundamental insight into landfill processes in order to research the feasibility of wastebody stabilisation. To predict behaviour of landfills, and to monitor the processes, as much information as possible needs to be gathered through measurements. These consist of measurement of gas concentration and production rates together with the variation in leachate quality and volumes, as well as many more. This research, as part of the CURE project, presents the applicability of settlement as and addition to these measurements.
The main goal of waste body stabilisation is to actively reduce the amount of organic matter in a landfill. Uncontrolled landfills produce considerable methane emissions as well as high concentrations of nitrogen and heavy metals in the leachate, leading to groundwater hazards in the environment.
McDougall (2007) introduced the fundamental conceptual model upon which this research is based. This model is not publicly available, so one of the goals of this research is to provide a fundamental conceptual model. The relation between degradation and settlement has been studied and implemented in a 1D model which allows for hydrological systems, oxygen penetration, degradation and cell strain. The outcome is a relation between mechanics, biochemistry and hydrology. A few assumptions and simplifications were made to make the model versatile and adjustible, but also easy to read. For instance, the model uses oxygen as absolute limiting factor in degradation, neglecting all processes concerning anaerobic degradation. The model also assumes that at the beginning of the simulation, each cell has the same composition, because there is no conclusive data about the distribution of waste types. The model has been created in one dimension, thereby neglecting all multi-dimensional processes and limitations.
In the final scenario, the main limitation to degradation is the availability of dissolved oxygen. As irrigation provides infiltration of 5 mm water per day with 10 mg dissolved oxygen per liter, the daily reduction of oxygen demand is 50 mg. The total modelled oxygen demand of the landfill is over 9600 kg, which concludes that a different tactic needs to happen in order to stabilise the wastebody. To increase the degradation rate, more oxygen needs to be applied. An approximation of the effect of applying a partial vacuum above the water table in the landfill to attract air from the environment is modelled, giving the oxygen more transfer area into unsaturated water. An increase by a factor 10^{4} is applied to the effect of oxygen on organic matter to model this enhanced irrigation system. The results are modeled over 100 years, after which approximately 26.17% of the degradation has happened. This shows that aerobic degradation over a wastebody takes too long with conventional degradation, and the unsaturated voids might help the process.
The model shows credible results for an unspecified landfill with deterministic parameters. Settlement has been brought in relation to a simplified form of biodegradation. The need for further research, with a spatial fluctuation of these parameters and detailed multidimensional water flow is needed to predict landfill behaviour in more detail. \pagebreak ...
The project 'introductie Duurzaam Stortbeheer (iDS) aims to create aftercare with an ending, or finite aftercare. It is looking for possibilities to stabilise the landfill mass, i.e. to eliminate the threat of pollutants by treating the wastebodies.
The goal of the CURE project is to develop fundamental insight into landfill processes in order to research the feasibility of wastebody stabilisation. To predict behaviour of landfills, and to monitor the processes, as much information as possible needs to be gathered through measurements. These consist of measurement of gas concentration and production rates together with the variation in leachate quality and volumes, as well as many more. This research, as part of the CURE project, presents the applicability of settlement as and addition to these measurements.
The main goal of waste body stabilisation is to actively reduce the amount of organic matter in a landfill. Uncontrolled landfills produce considerable methane emissions as well as high concentrations of nitrogen and heavy metals in the leachate, leading to groundwater hazards in the environment.
McDougall (2007) introduced the fundamental conceptual model upon which this research is based. This model is not publicly available, so one of the goals of this research is to provide a fundamental conceptual model. The relation between degradation and settlement has been studied and implemented in a 1D model which allows for hydrological systems, oxygen penetration, degradation and cell strain. The outcome is a relation between mechanics, biochemistry and hydrology. A few assumptions and simplifications were made to make the model versatile and adjustible, but also easy to read. For instance, the model uses oxygen as absolute limiting factor in degradation, neglecting all processes concerning anaerobic degradation. The model also assumes that at the beginning of the simulation, each cell has the same composition, because there is no conclusive data about the distribution of waste types. The model has been created in one dimension, thereby neglecting all multi-dimensional processes and limitations.
In the final scenario, the main limitation to degradation is the availability of dissolved oxygen. As irrigation provides infiltration of 5 mm water per day with 10 mg dissolved oxygen per liter, the daily reduction of oxygen demand is 50 mg. The total modelled oxygen demand of the landfill is over 9600 kg, which concludes that a different tactic needs to happen in order to stabilise the wastebody. To increase the degradation rate, more oxygen needs to be applied. An approximation of the effect of applying a partial vacuum above the water table in the landfill to attract air from the environment is modelled, giving the oxygen more transfer area into unsaturated water. An increase by a factor 10^{4} is applied to the effect of oxygen on organic matter to model this enhanced irrigation system. The results are modeled over 100 years, after which approximately 26.17% of the degradation has happened. This shows that aerobic degradation over a wastebody takes too long with conventional degradation, and the unsaturated voids might help the process.
The model shows credible results for an unspecified landfill with deterministic parameters. Settlement has been brought in relation to a simplified form of biodegradation. The need for further research, with a spatial fluctuation of these parameters and detailed multidimensional water flow is needed to predict landfill behaviour in more detail. \pagebreak
At the start of this project, we have a number of water level measurements obtained from various wells in the landfill. Straightforward spatial interpolation of this data leads to unexpected results. Most likely this is caused by the highly complex heterogeneity in this porous system. For this reason, this research aims to apply Electrical Resistivity Tomography (ERT) technology to explain the water distribution variations between wells. The apparent resistivity along several lines are measured over depth using different arrays. Some scripts written in Python with 'pyBERT' and 'pyGIMLi' packages are used to get electrical resistivity inversion results from the apparent resistivity. It is known that the decrease in the water content leads to a significant increase in the resistivity. Therefore, the possible existence of saturated and unsaturated blocks in the waste body can be visualized from the inversion maps.
Initially, the interface between the saturated and unsaturated zones is expected to be identified from Laplacian edge detection, while the results indicate that this technique fails to represent the area boundaries under highly-heterogeneous situations. Subsequently, Archie's law and van Genuchten equation are coupled to give a relation between the resistivity and water pressure head. Archie's law is used to compute the resistivity from water content and van Genuchten equation is used to compute the water content from the water pressure head. There are two hypotheses during this analysis: (a) where the resistivity is 20(ohm-m) gives the interface of dry and wet zones; and (b) the landfill leachate is under hydrostatic condition. Then the water pressure head is the distance from the interface, which can be read from the inversion maps. By selecting a certain range of empirical parameters, the computed resistivity-pressure head curves provide relatively good fits to the measured results. ...
At the start of this project, we have a number of water level measurements obtained from various wells in the landfill. Straightforward spatial interpolation of this data leads to unexpected results. Most likely this is caused by the highly complex heterogeneity in this porous system. For this reason, this research aims to apply Electrical Resistivity Tomography (ERT) technology to explain the water distribution variations between wells. The apparent resistivity along several lines are measured over depth using different arrays. Some scripts written in Python with 'pyBERT' and 'pyGIMLi' packages are used to get electrical resistivity inversion results from the apparent resistivity. It is known that the decrease in the water content leads to a significant increase in the resistivity. Therefore, the possible existence of saturated and unsaturated blocks in the waste body can be visualized from the inversion maps.
Initially, the interface between the saturated and unsaturated zones is expected to be identified from Laplacian edge detection, while the results indicate that this technique fails to represent the area boundaries under highly-heterogeneous situations. Subsequently, Archie's law and van Genuchten equation are coupled to give a relation between the resistivity and water pressure head. Archie's law is used to compute the resistivity from water content and van Genuchten equation is used to compute the water content from the water pressure head. There are two hypotheses during this analysis: (a) where the resistivity is 20(ohm-m) gives the interface of dry and wet zones; and (b) the landfill leachate is under hydrostatic condition. Then the water pressure head is the distance from the interface, which can be read from the inversion maps. By selecting a certain range of empirical parameters, the computed resistivity-pressure head curves provide relatively good fits to the measured results.
measurements of chloride and ammonium. By comparing the SARIMAX interpolated data and the data with reduced size, the results indicate that directly dropping half of the measurements can be regarded as an acceptable way to reduce the measurement frequency, as the data properties are well preserved and the errors in estimating the mass of substances leaching out are in the acceptable range. However, interpolating using SARIMAX model doesn’t have significant improvements in preserving the data properties. Further quartering the data can lead to large deviations in data properties. ...
measurements of chloride and ammonium. By comparing the SARIMAX interpolated data and the data with reduced size, the results indicate that directly dropping half of the measurements can be regarded as an acceptable way to reduce the measurement frequency, as the data properties are well preserved and the errors in estimating the mass of substances leaching out are in the acceptable range. However, interpolating using SARIMAX model doesn’t have significant improvements in preserving the data properties. Further quartering the data can lead to large deviations in data properties.
Gas flow through methane oxidation systems
A laboratory and numerical study for optimising system design
Landfills have been indicated as a major methane source. Methane oxidation systems are `low technology' systems that can treat these methane emissions. Yet, methane oxidation systems are herein still sub-optimal and leave room for improvement. A numerical model was established to research the effective gas permeability ratio between the gas distribution layer and the methane oxidation layer, and the centre-to-centre distance of the gas inlet points, necessary to achieve a spatial homogeneous methane load. In order to relate the permeability ratio to the design choice for the materials, laboratory experiments were performed to asses the influences of compaction level, hydraulic conditions and physical properties of a soil on the effective permeability of that soil. Overall, it is concluded that the effective permeability is predominantly influenced by the compaction level and soil texture. The water saturation only has a significant influence at near saturated levels. This means that the choice of suitable material and adequate construction practice has more effect on the effective permeability than seasonal changes in saturation levels in moderate climates. Furthermore, it is concluded that there are two parameters that govern the spatial homogeneity of the methane fluxes from the gas distribution layer into the methane oxidation layers: the permeability ratio between these layers, and the centre-to-centre distance between the inlet points. The required permeability ratio increases quadratically with an increasing centre-to-centre distance.
Nederlands:
Stortplaatsen zijn aangeduid als een belangrijke bron van methaanemissies. Methaanoxidatiesystemen zijn technisch simpele systemen die het stortgas kunnen saneren. Op dit moment zijn de bestaande methaanoxidatiesystemen nog weinig efficiënt en is het noodzakelijk om de homogeniteit van de laterale distributie van het stortgas te optimaliseren. Er is een numeriek model gegenereerd om inzicht te geven in de benodigde ratio tussen de effectieve permeabiliteit voor gas van de gasdistributielaag en van de methaanoxidatielaag, en in de maximale hart-op-hartafstand tussen de gasinlaatpunten om deze lateraal homogene methaan distributie te bereiken. Om de effectieve permeabiliteitsratio te relateren aan de materiaalselectie zijn er laboratoriumexperimenten uitgevoerd, die de invloeden van het compactieniveau, de hydraulische condities en de fysieke grondeigenschappen op de effectieve permeabiliteit voor gas vaststellen. Al met al kan worden geconcludeerd dat de effectieve permeabiliteit hoofdzakelijk wordt beïnvloed door het compactieniveau en de grondtextuur. Het watergehalte blijkt alleen significante invloed te hebben onder bijna verzadigde omstandigheden. Dit betekent dat een geschikte materiaalkeuze en adequate constructie de effectieve permeabiliteit voor gas meer beïnvloeden dan de seizoensgerelateerde veranderingen in het watergehalte in gematigde klimaten. Daarnaast kan worden geconcludeerd dat twee parameters bepalend zijn voor de laterale homogeniteit van de methaanstroom van de gasdistributielaag naar de methaanoxidatielaag: de ratio tussen de effectieve permeabiliteit voor gas tussen deze twee lagen, en de hart-op-hartafstand van de gasinlaatpunten. De benodigde permeabiliteitsratio neemt kwadratisch toe bij een toenemende hart-op-hartafstand. ...
Landfills have been indicated as a major methane source. Methane oxidation systems are `low technology' systems that can treat these methane emissions. Yet, methane oxidation systems are herein still sub-optimal and leave room for improvement. A numerical model was established to research the effective gas permeability ratio between the gas distribution layer and the methane oxidation layer, and the centre-to-centre distance of the gas inlet points, necessary to achieve a spatial homogeneous methane load. In order to relate the permeability ratio to the design choice for the materials, laboratory experiments were performed to asses the influences of compaction level, hydraulic conditions and physical properties of a soil on the effective permeability of that soil. Overall, it is concluded that the effective permeability is predominantly influenced by the compaction level and soil texture. The water saturation only has a significant influence at near saturated levels. This means that the choice of suitable material and adequate construction practice has more effect on the effective permeability than seasonal changes in saturation levels in moderate climates. Furthermore, it is concluded that there are two parameters that govern the spatial homogeneity of the methane fluxes from the gas distribution layer into the methane oxidation layers: the permeability ratio between these layers, and the centre-to-centre distance between the inlet points. The required permeability ratio increases quadratically with an increasing centre-to-centre distance.
Nederlands:
Stortplaatsen zijn aangeduid als een belangrijke bron van methaanemissies. Methaanoxidatiesystemen zijn technisch simpele systemen die het stortgas kunnen saneren. Op dit moment zijn de bestaande methaanoxidatiesystemen nog weinig efficiënt en is het noodzakelijk om de homogeniteit van de laterale distributie van het stortgas te optimaliseren. Er is een numeriek model gegenereerd om inzicht te geven in de benodigde ratio tussen de effectieve permeabiliteit voor gas van de gasdistributielaag en van de methaanoxidatielaag, en in de maximale hart-op-hartafstand tussen de gasinlaatpunten om deze lateraal homogene methaan distributie te bereiken. Om de effectieve permeabiliteitsratio te relateren aan de materiaalselectie zijn er laboratoriumexperimenten uitgevoerd, die de invloeden van het compactieniveau, de hydraulische condities en de fysieke grondeigenschappen op de effectieve permeabiliteit voor gas vaststellen. Al met al kan worden geconcludeerd dat de effectieve permeabiliteit hoofdzakelijk wordt beïnvloed door het compactieniveau en de grondtextuur. Het watergehalte blijkt alleen significante invloed te hebben onder bijna verzadigde omstandigheden. Dit betekent dat een geschikte materiaalkeuze en adequate constructie de effectieve permeabiliteit voor gas meer beïnvloeden dan de seizoensgerelateerde veranderingen in het watergehalte in gematigde klimaten. Daarnaast kan worden geconcludeerd dat twee parameters bepalend zijn voor de laterale homogeniteit van de methaanstroom van de gasdistributielaag naar de methaanoxidatielaag: de ratio tussen de effectieve permeabiliteit voor gas tussen deze twee lagen, en de hart-op-hartafstand van de gasinlaatpunten. De benodigde permeabiliteitsratio neemt kwadratisch toe bij een toenemende hart-op-hartafstand.
The suitability of gibbsite as aluminum source to complex with dissolved organic matter
In the interest of soil permeability reduction
A literature study was conducted to investigate the interplay between gibbsite and dissolved organic matter (DOM). Four processes are primarily studied in this research, they are: the release of free Al via dissolution of gibbsite; the complexation between free Al and DOM; the protonation and subsequently precipitation of DOM and the adsorption of H+/DOM onto the gibbsite surface. Thereafter the following hypotheses were formulated:
•The extent of gibbsite dissolution is larger at low pH. However, due to its low solubility and slow dissolution kinetics, a low concentration of free Al is expected.
•The organic matter (OM) source used in this study is a humic acid, therefore it is expected to become completely insoluble when the pH is below two.
•The contradicting pH dependency of gibbsite dissolution and OM solubility leads a narrow transitional pH range, that is favourable for the occurrence of complexation between free Al and DOM.
These hypotheses are tested by performing experiments in the laboratory, where the pH varies from 2 to 8 and the aluminum/carbon (mol/mol) ratio is 0; 0.05; 0.1; 0.3 and 1. A synthesized crystalline gibbsite powder and a potassium humate are used.
\bigskip
The results are subdivided into three pH ranges in which similar behaviour was found. In the low pH range (pH = 2 - 2.5), the protonation and subsequently precipitation of DOM is the dominant process in short term. However, when considering the long term behaviour, gibbsite dissolution is the dominant process as the DOM has precipitated. Due to the difference in time scale, there is very little complexation between free Al and DOM. In the transitional pH range (pH = 3 - 3.5), the competition for H+ between DOM and gibbsite is profound. This increases the solubility of DOM and slows the dissolution of gibbsite down. As a consequence, this favours the gradual complexation between free Al and DOM. Nevertheless, due to the ambiguity of free Al measurement method used in this research, it is difficult to quantify the complexation. In neutral to high pH range (pH = 4 - 8) DOM undergoes protonation over time but stays soluble. The release of free Al via the dissolution of gibbsite is negligible in this pH range. Therefore complexation between free Al and DOM should not be expected.
The results presented in this research indicate that gibbsite may not be suitable as the aluminum source to form precipitates with organic matter in engineering solutions that require fast soil permeability reduction. However, the slow release of free Al from gibbsite in combination with the gradual formation of precipitates with organic matter might be interesting for robust problems that require a self-healing ability. For efficiency concern, another OM source that is rich in fulvic acids (higher solubility at low pH than humic acid) and a natural gibbsite source (higher solubility than synthesized gibbsite) should be investigated.
...
A literature study was conducted to investigate the interplay between gibbsite and dissolved organic matter (DOM). Four processes are primarily studied in this research, they are: the release of free Al via dissolution of gibbsite; the complexation between free Al and DOM; the protonation and subsequently precipitation of DOM and the adsorption of H+/DOM onto the gibbsite surface. Thereafter the following hypotheses were formulated:
•The extent of gibbsite dissolution is larger at low pH. However, due to its low solubility and slow dissolution kinetics, a low concentration of free Al is expected.
•The organic matter (OM) source used in this study is a humic acid, therefore it is expected to become completely insoluble when the pH is below two.
•The contradicting pH dependency of gibbsite dissolution and OM solubility leads a narrow transitional pH range, that is favourable for the occurrence of complexation between free Al and DOM.
These hypotheses are tested by performing experiments in the laboratory, where the pH varies from 2 to 8 and the aluminum/carbon (mol/mol) ratio is 0; 0.05; 0.1; 0.3 and 1. A synthesized crystalline gibbsite powder and a potassium humate are used.
\bigskip
The results are subdivided into three pH ranges in which similar behaviour was found. In the low pH range (pH = 2 - 2.5), the protonation and subsequently precipitation of DOM is the dominant process in short term. However, when considering the long term behaviour, gibbsite dissolution is the dominant process as the DOM has precipitated. Due to the difference in time scale, there is very little complexation between free Al and DOM. In the transitional pH range (pH = 3 - 3.5), the competition for H+ between DOM and gibbsite is profound. This increases the solubility of DOM and slows the dissolution of gibbsite down. As a consequence, this favours the gradual complexation between free Al and DOM. Nevertheless, due to the ambiguity of free Al measurement method used in this research, it is difficult to quantify the complexation. In neutral to high pH range (pH = 4 - 8) DOM undergoes protonation over time but stays soluble. The release of free Al via the dissolution of gibbsite is negligible in this pH range. Therefore complexation between free Al and DOM should not be expected.
The results presented in this research indicate that gibbsite may not be suitable as the aluminum source to form precipitates with organic matter in engineering solutions that require fast soil permeability reduction. However, the slow release of free Al from gibbsite in combination with the gradual formation of precipitates with organic matter might be interesting for robust problems that require a self-healing ability. For efficiency concern, another OM source that is rich in fulvic acids (higher solubility at low pH than humic acid) and a natural gibbsite source (higher solubility than synthesized gibbsite) should be investigated.
A water balance model for landfill De Kragge II
Supporting the development of a sustainable aftercare approach
This research focuses on adapting the prediction models, created for two demonstration projects, to describe the water balance for the third demonstration field at De Kragge II in Bergen op Zoom. Compared with the other two projects, the water balance for this field is a bit more complicated, since it has a different layout of the drainage system and lateral flow to and from adjacent fields is possible. The aim of this research is to model the water dynamics in the landfill with minimal uncertainty.
Available input and output data are: rainfall, evaporation potential, leachate levels and leachate outflow, available for the pilot field and the adjacent compartment. The leachate outflow is controlled by valves, level meters and pumps in the flow system, also the operator has influence on which compartment is drained. Another complication with the outflow data is that the data from weighed trucks transporting the leachate indicate that sometimes leachate was directly pumped from the landfill, herewith bypassing the flowmeter.
The model consists of three layers, a recultivation-, waste- and drainage layer. In the recultivation layer infiltration into the waste layer is calculated by balancing rainfall, evapotranspiration and storage. The water volume infiltrating the waste layer is distributed stochastically, according to a travel time distribution, that discretizes the infiltrated water to faster and slower moving regimes. In the drainage layer model, the balance of water inflow, leachate outflow, sideflow and storage is calculated. Resulting
in a volume of water ex-filtrating the landfill. To evaluate model uncertainty, visual and quantative criteria are used. The fits of modelled on measured data is used as a visual check. The quantitative analysis consists of evaluating the Kullback-Leibler divergence and the marginalized likelihood. The Kullback-leibler divergence estimates how much information is gained from the parameters, while the marginalized likelihood determines the balance between information content and complexity.
In order to find a model that describes the water dynamics with minimal uncertainty, three different model implementations were evaluated. In the first approach both leachate level and outflow were used for calibration with measured data. Evaluation of the model performance showed that the outflow could not be determined with acceptable error. This was indicated by large standard deviations of the model and measurement error with respect to outflow measurements. Likely the reason for this is the gap in the water balance and the erratic patterns in the outflow data. A second approach was therefore modelled in which only the leachate levels were fitted with measured data and the outflow data was given as input. This increased the leachate level fits slightly, also the quantitative criteria showed that approach 2 is better than approach 1. Some of the parameters of approach 2 had large uncertainty and the model is quite complex given the available measured data. Therefore a third, simpler and faster model was implemented. In this model the waste layer calculations were simplified using the circular convolution function of MATLAB. This function calculates the travel time distribution continuously instead of discretizing the function over given retention times, which was done in the first two approaches. This model approach gave the best leachate level fits. The Kullback-leibler divergence indicated higher information content compared to approach 2.
In addition, each approach was evaluated with different model scenarios in which the waste compartments where either coupled or uncoupled. For each approach the uncoupled models performed visually and quantitatively better than the coupled models.
Approach 2 gave the best insight in the water dynamics given its complexity, shown by the higher values for the marginalized likelihood. Approach 3 gave the best fits, of the leachate levels, the highest 퐷KL values and is the fastest model. From these results it would be advised to use approach 3 to analyse the water dynamics of the landfill. Given the results of the different scenarios, it would be advised to use the uncoupled models for the analysis of water dynamics inside the landfill, since these models showed the best fits, highest 퐷KL and marginalized likelihood values.
Based on the obtained results, the following insights could be drawn about the landfill dynamics. The sideflow between the two compartments is about 5 to 25 m3/day. The model showed that water in the landfill, flows fast from the cover layer to the drainage layer. The infiltration flux of the recultivation layer model seemed to be dominated by rainfall and evaporation. Therefore this model could be simplified by omitting water storage and flow through the layer. ...
This research focuses on adapting the prediction models, created for two demonstration projects, to describe the water balance for the third demonstration field at De Kragge II in Bergen op Zoom. Compared with the other two projects, the water balance for this field is a bit more complicated, since it has a different layout of the drainage system and lateral flow to and from adjacent fields is possible. The aim of this research is to model the water dynamics in the landfill with minimal uncertainty.
Available input and output data are: rainfall, evaporation potential, leachate levels and leachate outflow, available for the pilot field and the adjacent compartment. The leachate outflow is controlled by valves, level meters and pumps in the flow system, also the operator has influence on which compartment is drained. Another complication with the outflow data is that the data from weighed trucks transporting the leachate indicate that sometimes leachate was directly pumped from the landfill, herewith bypassing the flowmeter.
The model consists of three layers, a recultivation-, waste- and drainage layer. In the recultivation layer infiltration into the waste layer is calculated by balancing rainfall, evapotranspiration and storage. The water volume infiltrating the waste layer is distributed stochastically, according to a travel time distribution, that discretizes the infiltrated water to faster and slower moving regimes. In the drainage layer model, the balance of water inflow, leachate outflow, sideflow and storage is calculated. Resulting
in a volume of water ex-filtrating the landfill. To evaluate model uncertainty, visual and quantative criteria are used. The fits of modelled on measured data is used as a visual check. The quantitative analysis consists of evaluating the Kullback-Leibler divergence and the marginalized likelihood. The Kullback-leibler divergence estimates how much information is gained from the parameters, while the marginalized likelihood determines the balance between information content and complexity.
In order to find a model that describes the water dynamics with minimal uncertainty, three different model implementations were evaluated. In the first approach both leachate level and outflow were used for calibration with measured data. Evaluation of the model performance showed that the outflow could not be determined with acceptable error. This was indicated by large standard deviations of the model and measurement error with respect to outflow measurements. Likely the reason for this is the gap in the water balance and the erratic patterns in the outflow data. A second approach was therefore modelled in which only the leachate levels were fitted with measured data and the outflow data was given as input. This increased the leachate level fits slightly, also the quantitative criteria showed that approach 2 is better than approach 1. Some of the parameters of approach 2 had large uncertainty and the model is quite complex given the available measured data. Therefore a third, simpler and faster model was implemented. In this model the waste layer calculations were simplified using the circular convolution function of MATLAB. This function calculates the travel time distribution continuously instead of discretizing the function over given retention times, which was done in the first two approaches. This model approach gave the best leachate level fits. The Kullback-leibler divergence indicated higher information content compared to approach 2.
In addition, each approach was evaluated with different model scenarios in which the waste compartments where either coupled or uncoupled. For each approach the uncoupled models performed visually and quantitatively better than the coupled models.
Approach 2 gave the best insight in the water dynamics given its complexity, shown by the higher values for the marginalized likelihood. Approach 3 gave the best fits, of the leachate levels, the highest 퐷KL values and is the fastest model. From these results it would be advised to use approach 3 to analyse the water dynamics of the landfill. Given the results of the different scenarios, it would be advised to use the uncoupled models for the analysis of water dynamics inside the landfill, since these models showed the best fits, highest 퐷KL and marginalized likelihood values.
Based on the obtained results, the following insights could be drawn about the landfill dynamics. The sideflow between the two compartments is about 5 to 25 m3/day. The model showed that water in the landfill, flows fast from the cover layer to the drainage layer. The infiltration flux of the recultivation layer model seemed to be dominated by rainfall and evaporation. Therefore this model could be simplified by omitting water storage and flow through the layer.
Groundwater and solute transport modelling study Vosdonk Noord at Etten-Leur
Examining the effect of two implementation methodologies for highly heterogenic shallow subsurface characteristics
Causes and prevention of irrigation blockages in the bioleaching process
Case study: Terrafame mine