Modelling solute transport in structured soils: performance evaluation of the ADR and TRM models

The movement of chemicals through the soil to the groundwater or discharged to surface waters represents a degradation of these resources. In many cases, serious human and stock health implications are associated with this form of pollution. The chemicals of interest include nutrients, pesticides, salts, and industrial wastes. Recent studies have shown that current models and methods do not adequately describe the leaching of nutrients through soil, often underestimating the risk of groundwater contamination by surface-applied chemicals and overestimating the concentration of resident solutes. This inaccuracy results primarily from ignoring soil structure and nonequilibrium between soil constituents, water, and solutes. A multiple sample percolation system (MSPS), consisting of 25 individual collection wells, was constructed to study the effects of localized soil heterogeneities on the transport of nutrients (NO−3, Cl−, PO3−4) in the vadose zone of an agricultural soil predominantly dominated by clay. Very significant variations in drainage patterns across a small spatial scale were observed (one-way ANOVA, p < 0.001 indicating considerable heterogeneity in water flow patterns and nutrient leaching. Using data collected from the multiple sample percolation experiments, this paper compares the performance of two mathematical models for predicting solute transport, the advective-dispersion model with a reaction term (ADR), and a two-region preferential flow model (TRM) suitable for modelling nonequilibrium transport. These results have implications for modelling solute transport and predicting nutrient loading on a larger scale.

Quantifying effects of soil heterogeneity on groundwater pollution at four sites in the USA

Four sites located in the north-eastern region of the United States of America have been chosen to investigate the impacts of soil heterogeneity in the transport of solutes (bromide and chloride) through the vadose zone (the zone in the soil that lies below the root zone and above the permanent saturated groundwater). A recently proposed mathematical model based on the cumulative beta distribution has been deployed to compare and contrast the regions' heterogeneity from multiple sample percolation experiments. Significant differences in patterns of solute leaching were observed even over a small spatial scale, indicating that traditional sampling methods for solute transport, for example the gravity pan or suction Iysimeters, or more recent inventions such as the multiple sample percolation systems may not be effective in estimating solute fluxes in soils when a significant degree of soil heterogeneity is present. Consequently, ignoring soil heterogeneity in solute transport studies will likely result in under- or overprediction of leached fluxes and potentially lead to serious pollution of soils and/or groundwater. The cumulative beta distribution technique is found to be a versatile and simple technique of gaining valuable information regarding soil heterogeneity effects on solute transport. It is also an excellent tool for guiding future decisions of experimental designs particularly in regard to the number of samples within one site and the number of sampling locations between sites required to obtain a representative estimate of field solute or drainage flux.

Concepts and dimensionality in modeling unsaturated water flow and solute transport

Many environmental studies require accurate simulation of water and solute fluxes in the unsaturated zone. This paper evaluates one- and multi-dimensional approaches for soil water flow as well as different spreading mechanisms to model solute behavior at different scales. For quantification of soil water fluxes,Richards equation has become the standard. Although current numerical codes show perfect water balances, the calculated soil water fluxes in case of head boundary conditions may depend largely on the method used for spatial averaging of the hydraulic conductivity. Atmospheric boundary conditions, especially in the case of phreatic groundwater levels fluctuating above and below a soil surface, require sophisticated solutions to ensure convergence. Concepts for flow in soils with macro pores and unstable wetting fronts are still in development. One-dimensional flow models are formulated to work with lumped parameters in order to account for the soil heterogeneity and preferential flow. They can be used at temporal and spatial scales that are of interest to water managers and policymakers. Multi-dimensional flow models are hampered by data and computation requirements.Their main strength is detailed analysis of typical multi-dimensional flow problems, including soil heterogeneity and preferential flow. Three physically based solute-transport concepts have been proposed to describe solute spreading during unsaturated flow: The stochastic-convective model (SCM), the convection-dispersion equation (CDE), and the fraction aladvection-dispersion equation (FADE). A less physical concept is the continuous-time random-walk process (CTRW). Of these, the SCM and the CDE are well established, and their strengths and weaknesses are identified. The FADE and the CTRW are more recent,and only a tentative strength weakness opportunity threat (SWOT)analysis can be presented at this time. We discuss the effect of the number of dimensions in a numerical model and the spacing between model nodes on solute spreading and the values of the solute-spreading parameters. In order to meet the increasing complexity of environmental problems, two approaches of model combination are used: Model integration and model coupling. Amain drawback of model integration is the complexity of there sulting code. Model coupling requires a systematic physical domain and model communication analysis. The setup and maintenance of a hydrologic framework for model coupling requires substantial resources, but on the other hand, contributions can be made by many research groups.

Process-based modelling of mass transport in the subsurface

Modern societies rely on natural water pathways that include subsurface flow of water and dissolved chemicals. The thesis presents a range of numerical and analytical models for simulating physical, chemical and biological processes in the subsurface, including coastal aquifers, the near-surface vadose zone, and solute transport in biogeochemically active aquifers.

Investigation of seismic effects on disturbed rock slope stability assessments

For general stability analysis of rock slopes, rock mass strength and rock mass disturbance are definitely should be considered. In addition, the impact of earthquakes must be taken into account. In fact, the rock mass strength is very difficult to be assessed which causes the difficulty of analysing rock slope stability. Therefore, an empirical failure criterion, the Hoek-Brown failure criterion, has been proposed. It is one of the most widely accepted approaches to estimate rock mass strength. The rock mass disturbance is important and was found having significant influence on evaluating rock slope stability, especially for rock slope with poor quality rock mass. In the Hoek-Brown failure criterion, the disturbance factor can represent the level of the rock mass disturbance which would provide a reasonable basis for estimating rock mass strength. This research will not only discuss the slope factor of safety, but also consider the influence of the seismic force on rock slope stability assessment using pseudo-static method. In practice, only horizontal seismic coefficient is used. Various magnitudes of the disturbance factor and recommended blasting damage zone thickness are also taken into account. The blasting damage zone thickness considered ranges from 0.5 to 2.5 times of slope height. The research results have potential to be extended and then sets of comprehensive stability charts can be provided for the rock slope stability evaluations. They will be convenient tools for practising engineers. In this study, finite element upper bound and lower bound limit analysis methods are employed. Their applicability has been investigated in some previous studies. The differences between upper bound and lower bound solutions are less than ±10% which would provide reasonable and acceptable range for rock slope stability safety factor estimation.

Retinal tissue thickness in type 1 and type 2 diabetes

BACKGROUND: The objective was to investigate full retinal and inner retinal thickness in individuals with type 1 and type 2 diabetes. METHODS: Eighty-four individuals with type 1 diabetes (T1DM), 67 individuals with type 2 diabetes (T2DM) and 42 non-diabetic individuals (control group) were enrolled. Participants underwent full retinal thickness evaluation in the central retinal, parafoveal and perifoveal zones and in the retinal nerve fibre layer (RNFL) and ganglion cell complex (GCC), using spectral domain optical coherence tomography. As a preliminary step, the key variables of interest - age, sex, diabetic retinopathy (DR), duration of diabetes and HbA1c levels - were analysed and compared between the three groups. Full retinal thickness, RNFL and GCC thicknesses were also compared between the groups. The relationship between the type of diabetes and retinal tissue thickness was explored, adjusting for the five potential confounders. RESULTS: Compared to individuals with T1DM, individuals with T2DM had significantly reduced full retinal thickness in the parafovea and perifovea and reduced RNFL and GCC thickness. The mean differences were six (p = 0.020), seven (p = 0.008), six (p = 0.021) and four micrometres (p = 0.013) for the parafovea, perifovea, RNFL and GCC thicknesses, respectively. Thicknesses within the central zone (p = 0.018) and at the parafovea (p = 0.007) were significantly reduced in T2DM when compared to the control group. After adjusting for age, sex, diabetic retinopathy, duration of diabetes and HbA1c levels, the relationship between type of diabetes and retinal tissue thickness was not statistically significant (p > 0.056). CONCLUSION: Retinal tissue thickness is not significantly different between type 1 and type 2 diabetes, when adjusted for age, sex, diabetic retinopathy, duration of diabetes and HbA1c levels.

Retinal thickness profile of individuals with diabetes

PURPOSE: To examine the retinal thickness profiles of individuals with and without diabetic retinopathy (DR).

METHODS: Full retinal thickness in the central zone, overall and hemisphere thicknesses of the parafovea and perifovea, ganglion cell complex (GCC) thickness and retinal nerve fibre layer (RNFL) thickness were assessed in 185 individuals using spectral domain optical coherence tomography (88 individuals with diabetes but no DR, 55 with DR, and 42 non-diabetic controls). The DR group comprised of 60% of participants with very mild non-proliferative diabetic retinopathy (NPDR) (representing microaneurysms only) and 40% with mild NPDR (hard exudates, cotton-wool spots, and/or mild retinal haemorrhages). Regression analysis was performed to determine the factors associated with retinal tissue thickness, taking into account, age, sex, presence of DR, duration of diabetes, HbA1c levels and type of diabetes.

RESULTS: The mean (S.D.) of the overall parafoveal thickness was 306 (16) in the DR group and 314 (14) in the control group (p = 0.02). The mean (S.D.) of the superior hemisphere parafoveal thickness was 309 (16) in the DR group and 318 (14) in the control group (p = 0.02). The mean (S.D.) of the inferior hemisphere parafoveal thickness was 303 (17) in the DR group and 311 (15) in the control group (p = 0.02). There were no significant differences in retinal thickness between groups in the central zone (p = 0.27) or perifovea (p > 0.41). Neither the overall nor the hemisphere RNFL (p > 0.75) and GCC thickness (p > 0.37) were significantly different between the groups. Regression analysis revealed that parafoveal thickness in diabetic individuals was reduced in association with presence of DR (B = -5.9 μm, p = 0.02) and with advancing age (B = -4.5 μm, p = 0.004, for every 10 year increase in age) when adjusted for sex, duration of diabetes, HbA1c levels and type of diabetes.

CONCLUSION: The inner macula is thinner in the presence of clinical signs of diabetic retinopathy and is compounded by advancing age. The influence of any macular oedema or that by cotton-wool spots could not be ruled out and may still confound these results.