Predicting flow and transport in highly heterogeneous alluvial aquifers

Successful prediction of groundwater flow and solute transport through highly heterogeneous aquifers has remained elusive due to the limitations of methods to characterize hydraulic conductivity (K) and generate realistic stochastic fields from such data. As a result, many studies have suggested that the classical advective-dispersive equation (ADE) cannot reproduce such transport behavior. Here we demonstrate that when high-resolution K data are used with a fractal stochastic method that produces K fields with adequate connectivity, the classical ADE can accurately predict solute transport at the macrodispersion experiment site in Mississippi. This development provides great promise to accurately predict contaminant plume migration, design more effective remediation schemes, and reduce environmental risks. Key Points Non-Gaussian transport behavior at the MADE site is unraveledADE can reproduce tracer transport in heterogeneous aquifers with no calibrationNew fractal method generates heterogeneous K fields with adequate connectivity

Co-optimisation of indoor environmental quality and energy consumption within urban office buildings

This study aimed to develop a multi-component model that can be used to maximise indoor environmental quality inside mechanically ventilated office buildings, while minimising energy usage. The integrated model, which was developed and validated from fieldwork data, was employed to assess the potential improvement of indoor air quality and energy saving under different ventilation conditions in typical air-conditioned office buildings in the subtropical city of Brisbane, Australia. When operating the ventilation system under predicted optimal conditions of indoor environmental quality and energy conservation and using outdoor air filtration, average indoor particle number (PN) concentration decreased by as much as 77%, while indoor CO2 concentration and energy consumption were not significantly different compared to the normal summer time operating conditions. Benefits of operating the system with this algorithm were most pronounced during the Brisbane’s mild winter. In terms of indoor air quality, average indoor PN and CO2 concentrations decreased by 48% and 24%, respectively, while potential energy savings due to free cooling went as high as 108% of the normal winter time operating conditions. The application of such a model to the operation of ventilation systems can help to significantly improve indoor air quality and energy conservation in air-conditioned office buildings.

Applying HYDRUS to flow in a sodic clay soil with solution composition–dependent hydraulic conductivity

Management of sodic soils under irrigation often requires application of chemical ameliorants to improve permeability combined with leaching of excess salts. Modeling irrigation, soil treatments, and leaching in these sodic soils requires a model that can adequately represent the physical and chemical changes in the soil associated with the amelioration process. While there are a number of models that simulate reactive solute transport, UNSATCHEM and HYDRUS-1D are currently the only models that also include an ability to simulate the impacts of soil chemistry on hydraulic conductivity. Previous researchers have successfully applied these models to simulate amelioration experiments on a sodic loam soil. To further gauge their applicability, we extended the previous work by comparing HYDRUS simulations of sodic soil amelioration with the results from recently published laboratory experiments on a more reactive, repacked sodic clay soil. The general trends observed in the laboratory experiments were able to be simulated using HYDRUS. Differences between measured and simulated results were attributed to the limited flexibility of the function that represents chemistry-dependent hydraulic conductivity in HYDRUS. While improvements in the function could be made, the present work indicates that HYDRUS-UNSATCHEM captures the key changes in soil hydraulic properties that occur during sodic clay soil amelioration and thus extends the findings of previous researchers studying sodic loams.

The risk of pollen-mediated gene flow from exotic Corymbia plantations into native Corymbia populations in Australia

Sectors of the forest plantation industry in Australia are set to expand in the near future using species or hybrids of the spotted gums (Corymbia, Section Politaria). Plantations of these taxa have already been introduced across temperate and subtropical Australia, representing locally exotic introductions from native stands in Queensland and New South Wales. A literature review was undertaken to provide insights into the potential for pollen-mediated gene flow from these plantations into native populations. Three factors suggest that such gene flow is likely; (1) interspecific hybridisation within the genus has frequently been recorded, including between distantly related species from different sections, (2) apparent high levels of vertebrate pollinator activity may result in plantation pollen being moved over hundreds of kilometres, (3) much of the plantation estate is being established among closely related taxa and therefore few barriers to gene flow are expected. Across Australia, 20 of the 100 native Corymbia taxa were found to have regional level co-occurrence with plantations. These were located most notably within regions of north-east New South Wales and south-east Queensland, however, co-occurrence was also found in south-west Western Australia and eastern Victoria. The native species found to have co-occurrence were then assessed for the presence of reproductive barriers at each step in the process of gene flow that may reduce the number of species at risk even further. The available data suggest three risk categories exist for Corymbia. The highest risk was for gene flow from plantations of spotted gums to native populations of spotted gums. This was based on the expected limited existence of pre- and post-zygotic barriers, substantial long-distance pollen dispersal and an apparent broad period of flowering in Corymbia citriodora subsp. variegata plantations. The following risk category focussed on gene flow from Corymbia torelliana × C. c. variegata hybrid plantations into native C. c. variegata, as the barriers associated with the production and establishment of F1 hybrids have been circumvented. For the lowest risk category, Corymbia plantations may present a risk to other non-spotted gum species, however, further investigation of the particular cross-combinations is required. A list of research directions is provided to better quantify these risks. Empirical data will need to be combined within a risk assessment framework that will not only estimate the likelihood of exotic gene flow, but also consider the conservation status/value of the native populations. In addition, the potential impacts of pollen flow from plantations will need to be weighed up against their various economic and environmental benefits.

The effect of habitat and environmental history on otolith chemistry of barramundi Lates calcarifer in estuarine populations of a regulated tropical river

We examine the microchemistry of otoliths of cohorts of a fished shed population of the large catadromous fish, barramundi Lates calcarifer from the estuary of a large tropical river. Barramundi from the estuary of the large, heavily regulated Fitzroy River, north eastern Australia were analysed by making transects of 87Sr/86Sr isotope and trace metal/Ca ratios from the core to the outer edge. Firstly, we examined the Sr/Ca, Ba/Ca, Mg/Ca and Mn/Ca and 87Sr/86Sr isotope ratios in otoliths of barramundi tagged in either freshwater or estuarine habitats that were caught by the commercial fishery in the estuary. We used 87Sr/86Sr isotope ratios to identify periods of freshwater residency and assess whether trace metal/Ca ratios varied between habitats. Only Sr/Ca consistently varied between known periods of estuarine or freshwater residency. The relationships between trace metal/Ca and river flow, salinity, temperature were examined in fish tagged and recaptured in the estuary. We found weak and inconsistent patterns in relationships between these variables in the majority of fish. These results suggest that both individual movement history within the estuary and the scale of environmental monitoring were reducing our ability to detect any patterns. Finally, we examined fish in the estuary from two dominant age cohorts (4 and 7 yr old) before and after a large flood in 2003 to ascertain if the flood had enabled fish from freshwater habitats to migrate to the estuary. There was no difference in the proportion of fish in the estuary that had accessed freshwater after the flood. Instead, we found that larger individuals with each age cohort were more likely to have spent a period in freshwater. This highlights the need to maintain freshwater flows in rivers. About half the fish examined had accessed freshwater habitats before capture. Of these, all had spent at least their first two months in marine salinity waters before entering freshwater and some did not enter freshwater until four years of age. This contrasts with the results of several previous studies in other parts of the range that found that access to freshwater swamps by larval barramundi was important for enhanced population productivity and recruitment.

The role of pollen in the changing environmental conditions of Scots pine

Rapid change in climate is challenge for the adaptation of forest trees in the future. In wind pollinated tree species pollen mediated long distance gene flow may provide alleles that are (pre)adapted to a future climate. In order to examine the long distance pollen flow in Scots pine (Pinus sylvestris L.), we measured the amount and viability of airborne pollen and flowering phenology in central, northern, and northernmost Finland during four years. Viable airborne pollen grains were detected during female flowering and before local pollen shedding in all study sites. The situation when there was nonlocal pollen in the air lasted from one to four days depending on the year and study site. The amount of nonlocal airborne pollen varied also between years and study sites, the total amount of nonlocal viable pollen in the air was 2.3% from all detected viable pollen grains. The effect of pollen origin on seeds siring ability was studied with artificial pollination experiments. Pollen genotypes originating from southern Finland sired 76% and 48 % of the analysed seeds in competition studies where both pollen origin were introduced simultaneously into the female strobili. We examined the importance of arrival order of pollen grains in to the strobili in a study where pollen genotypes of different origin were introduced in two hours interval. Northern genotypes sired 76% of the analysed seeds when it was injected first, but in the "southern first" experiment both pollen types sired equal amount of seeds. The first pollen grain in the pollen chamber do not always fertilizes the ovum, instead there likely is more complex way of competition between pollen grains. To examine chemically mediated pollen-pollen interactions we conducted in vitro germination experiment where different pollen genotypes had chemical but not physical contact. Both positive and negative effects of interactions were found. We found highly negative effects in germinability of northern pollen grains when they were germinating with southern pollen, and increase in the germinability of southern pollen. There were no variation in the size of the dry pollen grains between pollen origins, and minor variation between different genotypes. After hydration and germination northern pollen grains were larger than southern pollen. Pollen genotypes having high hydration rates had low germinability and tube growth rate, however, germinated pollen grains were larger in size than nongerminated. This supports the suggestion that the early germination and growth of pollen tube is dependent on pollen storage materialsand less dependent on water intake and hydration. Long distance pollen movements and good competition ability of southern pollen makes gene flow possible, although rising temperature and timing of pollen movements may affect pollen competition and the amount of gene flow.

Spatial and temporal variability in midge (Nematocera) assemblages in shallow Finnish lakes (60−70 °N) : community-based modelling of past environmental change

Multi- and intralake datasets of fossil midge assemblages in surface sediments of small shallow lakes in Finland were studied to determine the most important environmental factors explaining trends in midge distribution and abundance. The aim was to develop palaeoenvironmental calibration models for the most important environmental variables for the purpose of reconstructing past environmental conditions. The developed models were applied to three high-resolution fossil midge stratigraphies from southern and eastern Finland to interpret environmental variability over the past 2000 years, with special focus on the Medieval Climate Anomaly (MCA), the Little Ice Age (LIA) and recent anthropogenic changes. The midge-based results were compared with physical properties of the sediment, historical evidence and environmental reconstructions based on diatoms (Bacillariophyta), cladocerans (Crustacea: Cladocera) and tree rings. The results showed that the most important environmental factor controlling midge distribution and abundance along a latitudinal gradient in Finland was the mean July air temperature (TJul). However, when the dataset was environmentally screened to include only pristine lakes, water depth at the sampling site became more important. Furthermore, when the dataset was geographically scaled to southern Finland, hypolimnetic oxygen conditions became the dominant environmental factor. The results from an intralake dataset from eastern Finland showed that the most important environmental factors controlling midge distribution within a lake basin were river contribution, water depth and submerged vegetation patterns. In addition, the results of the intralake dataset showed that the fossil midge assemblages represent fauna that lived in close proximity to the sampling sites, thus enabling the exploration of within-lake gradients in midge assemblages. Importantly, this within-lake heterogeneity in midge assemblages may have effects on midge-based temperature estimations, because samples taken from the deepest point of a lake basin may infer considerably colder temperatures than expected, as shown by the present test results. Therefore, it is suggested here that the samples in fossil midge studies involving shallow boreal lakes should be taken from the sublittoral, where the assemblages are most representative of the whole lake fauna. Transfer functions between midge assemblages and the environmental forcing factors that were significantly related with the assemblages, including mean air TJul, water depth, hypolimnetic oxygen, stream flow and distance to littoral vegetation, were developed using weighted averaging (WA) and weighted averaging-partial least squares (WA-PLS) techniques, which outperformed all the other tested numerical approaches. Application of the models in downcore studies showed mostly consistent trends. Based on the present results, which agreed with previous studies and historical evidence, the Medieval Climate Anomaly between ca. 800 and 1300 AD in eastern Finland was characterized by warm temperature conditions and dry summers, but probably humid winters. The Little Ice Age (LIA) prevailed in southern Finland from ca. 1550 to 1850 AD, with the coldest conditions occurring at ca. 1700 AD, whereas in eastern Finland the cold conditions prevailed over a longer time period, from ca. 1300 until 1900 AD. The recent climatic warming was clearly represented in all of the temperature reconstructions. In the terms of long-term climatology, the present results provide support for the concept that the North Atlantic Oscillation (NAO) index has a positive correlation with winter precipitation and annual temperature and a negative correlation with summer precipitation in eastern Finland. In general, the results indicate a relatively warm climate with dry summers but snowy winters during the MCA and a cool climate with rainy summers and dry winters during the LIA. The results of the present reconstructions and the forthcoming applications of the models can be used in assessments of long-term environmental dynamics to refine the understanding of past environmental reference conditions and natural variability required by environmental scientists, ecologists and policy makers to make decisions concerning the presently occurring global, regional and local changes. The developed midge-based models for temperature, hypolimnetic oxygen, water depth, littoral vegetation shift and stream flow, presented in this thesis, are open for scientific use on request.

Topographical, structural and geophysical characterization of fracture zones: implications for groundwater flow and vulnerability

The main objective of this study is to evaluate selected geophysical, structural and topographic methods on regional, local, and tunnel and borehole scales, as indicators of the properties of fracture zones or fractures relevant to groundwater flow. Such information serves, for example, groundwater exploration and prediction of the risk of groundwater inflow in underground construction. This study aims to address how the features detected by these methods link to groundwater flow in qualitative and semi-quantitative terms and how well the methods reveal properties of fracturing affecting groundwater flow in the studied sites. The investigated areas are: (1) the Päijänne Tunnel for water-conveyance whose study serves as a verification of structures identified on regional and local scales; (2) the Oitti fuel spill site, to telescope across scales and compare geometries of structural assessment; and (3) Leppävirta, where fracturing and hydrogeological environment have been studied on the scale of a drilled well. The methods applied in this study include: the interpretation of lineaments from topographic data and their comparison with aeromagnetic data; the analysis of geological structures mapped in the Päijänne Tunnel; borehole video surveying; groundwater inflow measurements; groundwater level observations; and information on the tunnel s deterioration as demonstrated by block falls. The study combined geological and geotechnical information on relevant factors governing groundwater inflow into a tunnel and indicators of fracturing, as well as environmental datasets as overlays for spatial analysis using GIS. Geophysical borehole logging and fluid logging were used in Leppävirta to compare the responses of different methods to fracturing and other geological features on the scale of a drilled well. Results from some of the geophysical measurements of boreholes were affected by the large diameter (gamma radiation) or uneven surface (caliper) of these structures. However, different anomalies indicating more fractured upper part of the bedrock traversed by well HN4 in Leppävirta suggest that several methods can be used for detecting fracturing. Fracture trends appear to align similarly on different scales in the zone of the Päijänne Tunnel. For example, similarities of patterns were found between the regional magnetic trends, correlating with orientations of topographic lineaments interpreted as expressions of fracture zones. The same structural orientations as those of the larger structures on local or regional scales were observed in the tunnel, even though a match could not be made in every case. The size and orientation of the observation space (patch of terrain at the surface, tunnel section, or borehole), the characterization method, with its typical sensitivity, and the characteristics of the location, influence the identification of the fracture pattern. Through due consideration of the influence of the sampling geometry and by utilizing complementary fracture characterization methods in tandem, some of the complexities of the relationship between fracturing and groundwater flow can be addressed. The flow connections demonstrated by the response of the groundwater level in monitoring wells to pressure decrease in the tunnel and the transport of MTBE through fractures in bedrock in Oitti, highlight the importance of protecting the tunnel water from a risk of contamination. In general, the largest values of drawdown occurred in monitoring wells closest to the tunnel and/or close to the topographically interpreted fracture zones. It seems that, to some degree, the rate of inflow shows a positive correlation with the level of reinforcement, as both are connected with the fracturing in the bedrock. The following geological features increased the vulnerability of tunnel sections to pollution, especially when several factors affected the same locations: (1) fractured bedrock, particularly with associated groundwater inflow; (2) thin or permeable overburden above fractured rock; (3) a hydraulically conductive layer underneath the surface soil; and (4) a relatively thin bedrock roof above the tunnel. The observed anisotropy of the geological media should ideally be taken into account in the assessment of vulnerability of tunnel sections and eventually for directing protective measures.

Gene flow from transgenic plant populations: models and applications for risk assessment

The future use of genetically modified (GM) plants in food, feed and biomass production requires a careful consideration of possible risks related to the unintended spread of trangenes into new habitats. This may occur via introgression of the transgene to conventional genotypes, due to cross-pollination, and via the invasion of GM plants to new habitats. Assessment of possible environmental impacts of GM plants requires estimation of the level of gene flow from a GM population. Furthermore, management measures for reducing gene flow from GM populations are needed in order to prevent possible unwanted effects of transgenes on ecosystems. This work develops modeling tools for estimating gene flow from GM plant populations in boreal environments and for investigating the mechanisms of the gene flow process. To describe spatial dimensions of the gene flow, dispersal models are developed for the local and regional scale spread of pollen grains and seeds, with special emphasis on wind dispersal. This study provides tools for describing cross-pollination between GM and conventional populations and for estimating the levels of transgenic contamination of the conventional crops. For perennial populations, a modeling framework describing the dynamics of plants and genotypes is developed, in order to estimate the gene flow process over a sequence of years. The dispersal of airborne pollen and seeds cannot be easily controlled, and small amounts of these particles are likely to disperse over long distances. Wind dispersal processes are highly stochastic due to variation in atmospheric conditions, so that there may be considerable variation between individual dispersal patterns. This, in turn, is reflected to the large amount of variation in annual levels of cross-pollination between GM and conventional populations. Even though land-use practices have effects on the average levels of cross-pollination between GM and conventional fields, the level of transgenic contamination of a conventional crop remains highly stochastic. The demographic effects of a transgene have impacts on the establishment of trangenic plants amongst conventional genotypes of the same species. If the transgene gives a plant a considerable fitness advantage in comparison to conventional genotypes, the spread of transgenes to conventional population can be strongly increased. In such cases, dominance of the transgene considerably increases gene flow from GM to conventional populations, due to the enhanced fitness of heterozygous hybrids. The fitness of GM plants in conventional populations can be reduced by linking the selectively favoured primary transgene to a disfavoured mitigation transgene. Recombination between these transgenes is a major risk related to this technique, especially because it tends to take place amongst the conventional genotypes and thus promotes the establishment of invasive transgenic plants in conventional populations.

Substance flow analysis in Finland : Four case studies on N and P flows

Nitrogen (N) and phosphorus (P) are essential elements for all living organisms. However, in excess, they contribute to several environmental problems such as aquatic and terrestrial eutrophication. Globally, human action has multiplied the volume of N and P cycling since the onset of industrialization. The multiplication is a result of intensified agriculture, increased energy consumption and population growth. Industrial ecology (IE) is a discipline, in which human interaction with the ecosystems is investigated using a systems analytical approach. The main idea behind IE is that industrial systems resemble ecosystems, and, like them, industrial systems can then be described using material, energy and information flows and stocks. Industrial systems are dependent on the resources provided by the biosphere, and these two cannot be separated from each other. When studying substance flows, the aims of the research from the viewpoint of IE can be, for instance, to elucidate the ways how the cycles of a certain substance could be more closed and how the flows of a certain substance could be decreased per unit of production (= dematerialization). In Finland, N and P are studied widely in different ecosystems and environmental emissions. A holistic picture comparing different societal systems is, however, lacking. In this thesis, flows of N and P were examined in Finland using substance flow analysis (SFA) in the following four subsystems: I) forest industry and use of wood fuels, II) food production and consumption, III) energy, and IV) municipal waste. A detailed analysis at the end of the 1990s was performed. Furthermore, historical development of the N and P flows was investigated in the energy system (III) and the municipal waste system (IV). The main research sources were official statistics, literature, monitoring data, and expert knowledge. The aim was to identify and quantify the main flows of N and P in Finland in the four subsystems studied. Furthermore, the aim was to elucidate whether the nutrient systems are cyclic or linear, and to identify how these systems could be more efficient in the use and cycling of N and P. A final aim was to discuss how this type of an analysis can be used to support decision-making on environmental problems and solutions. Of the four subsystems, the food production and consumption system and the energy system created the largest N flows in Finland. For the creation of P flows, the food production and consumption system (Paper II) was clearly the largest, followed by the forest industry and use of wood fuels and the energy system. The contribution of Finland to N and P flows on a global scale is low, but when compared on a per capita basis, we are one of the largest producers of these flows, with relatively high energy and meat consumption being the main reasons. Analysis revealed the openness of all four systems. The openness is due to the high degree of internationality of the Finnish markets, the large-scale use of synthetic fertilizers and energy resources and the low recycling rate of many waste fractions. Reduction in the use of fuels and synthetic fertilizers, reorganization of the structure of energy production, reduced human intake of nutrients and technological development are crucial in diminishing the N and P flows. To enhance nutrient recycling and replace inorganic fertilizers, recycling of such wastes as wood ash and sludge could be promoted. SFA is not usually sufficiently detailed to allow specific recommendations for decision-making to be made, but it does yield useful information about the relative magnitude of the flows and may reveal unexpected losses. Sustainable development is a widely accepted target for all human action. SFA is one method that can help to analyse how effective different efforts are in leading to a more sustainable society. SFA's strength is that it allows a holistic picture of different natural and societal systems to be drawn. Furthermore, when the environmental impact of a certain flow is known, the method can be used to prioritize environmental policy efforts.

Environmental Mechanics Research in China

By recalling mankind's path during past 50 years in the present article, we mainly highlight the significance of environmental issues today. In particular, two major factors leading to environment deterioration in China such as water resources and coal burning are stressed on. Present-day environmental issues are obviously interdisciplinary, of multiple scales and multi-composition in nature. Therefore, a process-based approach for environment research is absolutely necessarily. A series of sub-processes, either physical, chemical or biological, are subsequently analyzed in order to established reasonable parameterization scheme and credible comprehensive model. And we are now in a position to answer questions still open to us, improve existing somewhat empirical engineering approaches and enhance quantitative accuracy in prediction. To illustrate this process-based research approach, three typical examples associated with the Yangtze River Estuary, Loess Plateau and Tenggeli Desert environments have been dealt with respectively. A theoretical model of vertical flow field accounting for runoff and tide interaction has been established to delineate salinity and sediment motion which are responsible for the formation of mouth bar at the outlet and the ecological evolution there. A kinematic wave theory combined with the revised Green-Ampt infiltration formula is applied to the prediction of runoff generation and erosion in three types of erosion region on the Loess Plateau. Three approaches describing water motion in SPAC system in arid areas at different levels have been improved by introducing vegetation sub-models. However, we have found that the formation of a dry sandy layer and biological crust skin are additional primary causes leading to deterioration of water supply and succession of ecological system.

Thermocapillary flow in an annular liquid layer painted on a moving fiber

In the present paper, a theoretical model is studied on the flow in the liquid annular film, which is ejected from a vessel with relatively higher temperature and painted on the moving solid fiber. A temperature gradient, driving a thermocapillary flow, is formed on the free surface because of the heat transfer from the liquid with relatively higher temperature to the environmental gas with relatively lower temperature. The thermocapillary flow may change the radii profile of the liquid film. This process analyzed is based on the approximations of lubrication theory and perturbation theory, and the equation of the liquid layer radii and the process of thermal hydrodynamics in the liquid layer are solved for a temperature distribution on the solid fiber.