El sistema cárstico de Sant Martí de Llémena

Descripció del sistema càrstic de Sant Martí de Llémena

Espeleotemas de ópalo y pseudocarst granítico

The amorphous silica (opal-A) speleothems associated to the open structural system of granitic rocks where the slow circulation of runoff is produced are mainly the result of the biological degradation of the rock. These speleothems are found in many different geographic, climatic and geological environments though always associated to granitoids. They show two different morphologies: cylindrical or long bodies and laminar or layered forms. They are internally formed by a mass of clasts and spheres of opal-Awith a porous texture that evolves to compact and massive due to the reiteration of the re-dissolution/re-precipitation of the amorphous opal by the water that circulates through it after each rainfall period. A final characteristic of each growth stage (end of rainy period) is the development of whiskers of minerals, normally gypsum, on the water output points of the speleothems. The dimensions of the speleothems are millimetric (length and/or thickness). In this paper their morphology and mineralogy are described based on their study by SEM, XRD and XRF, and there is established a new and more complete genetic hypothesis than the one that exists up to now to understand their formation

La zona lacustre de Banyoles (Girona, España)

Descripció de la zona lacustre de Banyoles des del punt de vista geològic i hidrogeologic

Lattice Boltzmann modeling of fluid flow and solute transport in karst aquifers

A novel modeling approach is applied to karst hydrology. Long-standing problems in karst hydrology and solute transport are addressed using Lattice Boltzmann methods (LBMs). These methods contrast with other modeling approaches that have been applied to karst hydrology. The motivation of this dissertation is to develop new computational models for solving ground water hydraulics and transport problems in karst aquifers, which are widespread around the globe. This research tests the viability of the LBM as a robust alternative numerical technique for solving large-scale hydrological problems. The LB models applied in this research are briefly reviewed and there is a discussion of implementation issues. The dissertation focuses on testing the LB models. The LBM is tested for two different types of inlet boundary conditions for solute transport in finite and effectively semi-infinite domains. The LBM solutions are verified against analytical solutions. Zero-diffusion transport and Taylor dispersion in slits are also simulated and compared against analytical solutions. These results demonstrate the LBM’s flexibility as a solute transport solver. The LBM is applied to simulate solute transport and fluid flow in porous media traversed by larger conduits. A LBM-based macroscopic flow solver (Darcy’s law-based) is linked with an anisotropic dispersion solver. Spatial breakthrough curves in one and two dimensions are fitted against the available analytical solutions. This provides a steady flow model with capabilities routinely found in ground water flow and transport models (e.g., the combination of MODFLOW and MT3D). However the new LBM-based model retains the ability to solve inertial flows that are characteristic of karst aquifer conduits. Transient flows in a confined aquifer are solved using two different LBM approaches. The analogy between Fick’s second law (diffusion equation) and the transient ground water flow equation is used to solve the transient head distribution. An altered-velocity flow solver with source/sink term is applied to simulate a drawdown curve. Hydraulic parameters like transmissivity and storage coefficient are linked with LB parameters. These capabilities complete the LBM’s effective treatment of the types of processes that are simulated by standard ground water models. The LB model is verified against field data for drawdown in a confined aquifer.

Lattice Boltzmann Modeling of Fluid Flow and Solute Transport in Karst Aquifers

A novel modeling approach is applied to karst hydrology. Long-standing problems in karst hydrology and solute transport are addressed using Lattice Boltzmann methods (LBMs). These methods contrast with other modeling approaches that have been applied to karst hydrology. The motivation of this dissertation is to develop new computational models for solving ground water hydraulics and transport problems in karst aquifers, which are widespread around the globe. This research tests the viability of the LBM as a robust alternative numerical technique for solving large-scale hydrological problems. The LB models applied in this research are briefly reviewed and there is a discussion of implementation issues. The dissertation focuses on testing the LB models. The LBM is tested for two different types of inlet boundary conditions for solute transport in finite and effectively semi-infinite domains. The LBM solutions are verified against analytical solutions. Zero-diffusion transport and Taylor dispersion in slits are also simulated and compared against analytical solutions. These results demonstrate the LBM’s flexibility as a solute transport solver. The LBM is applied to simulate solute transport and fluid flow in porous media traversed by larger conduits. A LBM-based macroscopic flow solver (Darcy’s law-based) is linked with an anisotropic dispersion solver. Spatial breakthrough curves in one and two dimensions are fitted against the available analytical solutions. This provides a steady flow model with capabilities routinely found in ground water flow and transport models (e.g., the combination of MODFLOW and MT3D). However the new LBM-based model retains the ability to solve inertial flows that are characteristic of karst aquifer conduits. Transient flows in a confined aquifer are solved using two different LBM approaches. The analogy between Fick’s second law (diffusion equation) and the transient ground water flow equation is used to solve the transient head distribution. An altered-velocity flow solver with source/sink term is applied to simulate a drawdown curve. Hydraulic parameters like transmissivity and storage coefficient are linked with LB parameters. These capabilities complete the LBM’s effective treatment of the types of processes that are simulated by standard ground water models. The LB model is verified against field data for drawdown in a confined aquifer.

Lattice Boltzmann Modeling and Specialized Laboratory Techniques to Determine the Permeability of Megaporous Karst Rock

The Pleistocene carbonate rock Biscayne Aquifer of south Florida contains laterally-extensive bioturbated ooltic zones characterized by interconnected touching-vug megapores that channelize most flow and make the aquifer extremely permeable. Standard petrophysical laboratory techniques may not be capable of accurately measuring such high permeabilities. Instead, innovative procedures that can measure high permeabilities were applied. These fragile rocks cannot easily be cored or cut to shapes convenient for conducting permeability measurements. For the laboratory measurement, a 3D epoxy-resin printed rock core was produced from computed tomography data obtained from an outcrop sample. Permeability measurements were conducted using a viscous fluid to permit easily observable head gradients (~2 cm over 1 m) simultaneously with low Reynolds number flow. For a second permeability measurement, Lattice Boltzmann Method flow simulations were computed on the 3D core renderings. Agreement between the two estimates indicates an accurate permeability was obtained that can be applied to future studies.

The hydrogeology, hydrology and hydrochemistry of the Condamine River Valley alluvium

The Cainozoic alluvium of the Condamine River valley is interpreted to consist of sediments deposited as floodplain and sheetwash deposits in bedrock valleys eroded into Mesozoic sedimentary rocks and tertiary volcanics. A maximum recorded sediment accumulation of 134 m is centred just south of Dalby. The lower section ofboth the flood plain and sheetwash alluvium is composed of variegated sandy and clayey sediments and the upper section of brown and grey sandy and clayey sediments.

Mapping of shallow coastal groundwaters, their hydrology and environmental geochemistry : Pumicestone catchment, Southeast Queensland

Blooms of the toxic cyanobacterium majuscula Lyngbya in the coastal waters of southeast Queensland have caused adverse impacts on both environmental health and human health, and on local economies such as fishing and tourism. A number of studies have confirmed that the main limiting nutrients (“nutrients of concern”) that contribute to these blooms area Fe, DOC, N, P and also pH. This study is conducted to establish the distribution of these parameters in a typical southeast Queensland coastal setting. The study maps the geochemistry of shallow groundwater in the mainland Pumicestone catchment with an emphasis on the nutrients of concern to understand how these nutrients relate to aquifer materials, landuse and anthropogenic activities. The results of the study form a GIS information layer which will be incorporated into a larger GIS model being produced by Queensland Department of Environment and Resource Management (DERM) to support landuse management to avoid/minimize blooms of Lyngbya in Moreton Bay, southeast Queensland, and other similar settings. A total of 38 boreholes were established in the mainland Pumicestone region and four sampling rounds of groundwater carried out in both dry and wet conditions. These groundwater samples were measured in the field for physico-chemical parameters, and in the laboratory analyses for the nutrients of concern, and other major and minor ions. Aquifer materials were confirmed using the Geological Survey of Queensland digital geology map, and geomaterials were assigned to seven categories which are A (sands), B (silts, sandy silts), C (estuarine mud, silts), D (humid soils), E (alluvium), F (sandstone) and G (other bedrock). The results of the water chemistry were examined by use of the software package AquaChem/AqQA, and divided into six groundwater groups, based on groundwater chemical types and location of boreholes. The type of aquifer material and location, and proximity to waterways was found to be important because they affected physico-chemical properties and concentrations of nutrients of concern and dissolved ions. The analytical results showed that iron concentrations of shallow groundwaters were high due to acid sulfate soils, and also mud and silt, but were lower in sand materials. DOC concentrations of these shallow groundwaters in the sand material were high probably due to rapid infiltration. In addition, DOC concentrations in some boreholes were high because they were installed in organic rich wetlands. The pH values of boreholes were from acidic to near neutral; some boreholes with pH values were low (< 4), showing acid sulfate soils in these boreholes. Concentrations of total nitrogen and total phosphorus of groundwaters were generally low, and the main causes of elevated concentrations of total nitrogen and total phosphorus are largely due to animal and human wastes and tend to be found in localized source areas. Comparison of the relative percentage of nitrogen species (NH3/NH4< Org-N, NO3-N and NO2-N) demonstrated that they could be related to sources such as animal waste, residential and agricultural fertilizers, forest and vegetation, mixed residents and farms, and variable setting and vegetation covers. Total concentrations of dissolved ions in sampling round 3 (dry period) were higher than those in sampling round 2 (wet period) due to both evaporation of groundwater in the dry period and the dilution of rainfall in the wet period. This showed that the highest concentrations of nutrients of concern were due to acid sulfate soils, aquifer materials, landuse and anthropogenic activities and were typically in aquifer materials of E (alluvium) and C (estuarine muds) and locations of Burpengary, Caboolture, and Glass Mountain catchments.

Mechanism based emulation of dynamic simulation models : concept and application in hydrology

Many model-based investigation techniques, such as sensitivity analysis, optimization, and statistical inference, require a large number of model evaluations to be performed at different input and/or parameter values. This limits the application of these techniques to models that can be implemented in computationally efficient computer codes. Emulators, by providing efficient interpolation between outputs of deterministic simulation models, can considerably extend the field of applicability of such computationally demanding techniques. So far, the dominant techniques for developing emulators have been priors in the form of Gaussian stochastic processes (GASP) that were conditioned with a design data set of inputs and corresponding model outputs. In the context of dynamic models, this approach has two essential disadvantages: (i) these emulators do not consider our knowledge of the structure of the model, and (ii) they run into numerical difficulties if there are a large number of closely spaced input points as is often the case in the time dimension of dynamic models. To address both of these problems, a new concept of developing emulators for dynamic models is proposed. This concept is based on a prior that combines a simplified linear state space model of the temporal evolution of the dynamic model with Gaussian stochastic processes for the innovation terms as functions of model parameters and/or inputs. These innovation terms are intended to correct the error of the linear model at each output step. Conditioning this prior to the design data set is done by Kalman smoothing. This leads to an efficient emulator that, due to the consideration of our knowledge about dominant mechanisms built into the simulation model, can be expected to outperform purely statistical emulators at least in cases in which the design data set is small. The feasibility and potential difficulties of the proposed approach are demonstrated by the application to a simple hydrological model.

The structure of Norway spruce (Picea abies [L.] Karst.) stems in relation to wood properties of sawn timber

An important challenge in forest industry is to get the appropriate raw material out from the forests to the wood processing industry. Growth and stem reconstruction simulators are therefore increasingly integrated in industrial conversion simulators, for linking the properties of wooden products to the three-dimensional structure of stems and their growing conditions. Static simulators predict the wood properties from stem dimensions at the end of a growth simulation period, whereas in dynamic approaches, the structural components, e.g. branches, are incremented along with the growth processes. The dynamic approach can be applied to stem reconstruction by predicting the three-dimensional stem structure from external tree variables (i.e. age, height) as a result of growth to the current state. In this study, a dynamic growth simulator, PipeQual, and a stem reconstruction simulator, RetroSTEM, are adapted to Norway spruce (Picea abies [L.] Karst.) to predict the three-dimensional structure of stems (tapers, branchiness, wood basic density) over time such that both simulators can be integrated in a sawing simulator. The parameterisation of the PipeQual and RetroSTEM simulators for Norway spruce relied on the theoretically based description of tree structure developing in the growth process and following certain conservative structural regularities while allowing for plasticity in the crown development. The crown expressed both regularity and plasticity in its development, as the vertical foliage density peaked regularly at about 5 m from the stem apex, varying below that with tree age and dominance position (Study I). Conservative stem structure was characterized in terms of (1) the pipe ratios between foliage mass and branch and stem cross-sectional areas at crown base, (2) the allometric relationship between foliage mass and crown length, (3) mean branch length relative to crown length and (4) form coefficients in branches and stem (Study II). The pipe ratio between branch and stem cross-sectional area at crown base, and mean branch length relative to the crown length may differ in trees before and after canopy closure, but the variation should be further analysed in stands of different ages and densities with varying site fertilities and climates. The predictions of the PipeQual and RetroSTEM simulators were evaluated by comparing the simulated values to measured ones (Study III, IV). Both simulators predicted stem taper and branch diameter at the individual tree level with a small bias. RetroSTEM predictions of wood density were accurate. For focusing on even more accurate predictions of stem diameters and branchiness along the stem, both simulators should be further improved by revising the following aspects in the simulators: the relationship between foliage and stem sapwood area in the upper stem, the error source in branch sizes, the crown base development and the height growth models in RetroSTEM. In Study V, the RetroSTEM simulator was integrated in the InnoSIM sawing simulator, and according to the pilot simulations, this turned out to be an efficient tool for readily producing stand scale information about stem sizes and structure when approximating the available assortments of wood products.

Soil hydrological properties and conditions, site preparation, and the long-term performance of planted Scots pine (Pinus sylvestris L.) on upland forest sites in Finnish Lapland

Scots pine (Pinus sylvestris L.) and Norway spruce (Picea abies (L.) Karst.) forests dominate in Finnish Lapland. The need to study the effect of both soil factors and site preparation on the performance of planted Scots pine has increased due to the problems encountered in reforestation, especially on mesic and moist, formerly spruce-dominated sites. The present thesis examines soil hydrological properties and conditions, and effect of site preparation on them on 10 pine- and 10 spruce-dominated upland forest sites. Finally, the effects of both the site preparation and reforestation methods, and soil hydrology on the long-term performance of planted Scots pine are summarized. The results showed that pine and spruce sites differ significantly in their soil physical properties. Under field capacity or wetter soil moisture conditions, planted pines presumably suffer from excessive soil water and poor soil aeration on most of the originally spruce sites, but not on the pine sites. The results also suggested that site preparation affects the soil-water regime and thus prerequisites for forest growth over two decades after site preparation. High variation in the survival and mean height of planted pine was found. The study suggested that on spruce sites, pine survival is the lowest on sites that dry out slowly after rainfall events, and that height growth is the fastest on soils that reach favourable aeration conditions for root growth soon after saturation, and/or where the average air-filled porosity near field capacity is large enough for good root growth. Survival, but not mean height can be enhanced by employing intensive site preparation methods on spruce sites. On coarser-textured pine sites, site preparation methods don t affect survival, but methods affecting soil fertility, such as prescribed burning and ploughing, seem to enhance the height growth of planted Scots pines over several decades. The use of soil water content in situ as the sole criterion for sites suitable for pine reforestation was tested and found to be a relatively uncertain parameter. The thesis identified new potential soil variables, which should be tested using other data in the future.

Carbon dynamics in peatlands under changing hydrology : Effects of water level drawdown on litter quality, microbial enzyme activities and litter decomposition rates

Pristine peatlands are carbon (C) accumulating wetland ecosystems sustained by a high water level (WL) and consequent anoxia that slows down decomposition. Persistent WL drawdown as a response to climate and/or land-use change directly affects decomposition: increased oxygenation stimulates decomposition of the old C (peat) sequestered under prior anoxic conditions. Responses of the new C (plant litter) in terms of quality, production and decomposability, and the consequences for the whole C cycle of peatlands are not fully understood. WL drawdown induces changes in plant community resulting in shift in dominance from Sphagnum and graminoids to shrubs and trees. There is increasing evidence that the indirect effects of WL drawdown via the changes in plant communities will have more impact on the ecosystem C cycling than any direct effects. The aim of this study is to disentangle the direct and indirect effects of WL drawdown on the new C by measuring the relative importance of 1) environmental parameters (WL depth, temperature, soil chemistry) and 2) plant community composition on litter production, microbial activity, litter decomposition rates and, consequently, on the C accumulation. This information is crucial for modelling C cycle under changing climate and/or land-use. The effects of WL drawdown were tested in a large-scale experiment with manipulated WL at two time scales and three nutrient regimes. Furthermore, the effect of climate on litter decomposability was tested along a north-south gradient. Additionally, a novel method for estimating litter chemical quality and decomposability was explored by combining Near infrared spectroscopy with multivariate modelling. WL drawdown had direct effects on litter quality, microbial community composition and activity and litter decomposition rates. However, the direct effects of WL drawdown were overruled by the indirect effects via changes in litter type composition and production. Short-term (years) responses to WL drawdown were small. In long-term (decades), dramatically increased litter inputs resulted in large accumulation of organic matter in spite of increased decomposition rates. Further, the quality of the accumulated matter greatly changed from that accumulated in pristine conditions. The response of a peatland ecosystem to persistent WL drawdown was more pronounced at sites with more nutrients. The study demonstrates that the shift in vegetation composition as a response to climate and/or land-use change is the main factor affecting peatland ecosystem C cycle and thus dynamic vegetation is a necessity in any models applied for estimating responses of C fluxes to changes in the environment. The time scale for vegetation changes caused by hydrological changes needs to extend to decades. This study provides grouping of litter types (plant species and part) into functional types based on their chemical quality and/or decomposability that the models could utilize. Further, the results clearly show a drop in soil temperature as a response to WL drawdown when an initially open peatland converts into a forest ecosystem, which has not yet been considered in the existing models.

Turvemaan männyn (Pinus sylvestris L.) ja kuusen (Picea abies [L.] Karst.) kanto- ja juuripuun ominaisuudet ja asetoniliukoisten uuteaineiden määrä

Tutkimus on osa Metsäklusteri Oy:n Future Biorefinery –tutkimusohjelmaa, jossa kartoitetaan mahdollisuuksia hyödyntää metsäteollisuuden raaka-aineita aiempaa tarkemmin ja uusissa tuotteissa. Tutkimuksen tavoitteena on selvittää männyn (Pinus sylvestris L.) ja kuusen (Picea abies [L.] Karst.) juurten ja kantopuun rakenne ja ominaisuudet. Tutkimuksessa selvitetään löytyykö männyn ja kuusen juurista reaktiopuuta ja määritetään asetoniliukoisten uuteaineiden osuus kantoja juuripuussa. Tutkimusaineistona oli viisi eri-ikäistä mäntyä ja viisi eri-ikäistä kuusta. Juuri- ja kantoaineisto kerättiin Metsäntutkimuslaitoksen toimesta Parkanon seudulta (62.017°N, 23.017°E) hakkuun jälkeen. Maanalaisista juurista otettiin näytteet kolmelta eri etäisyydeltä juurenniskaan nähden. Kummankaan lajin juurista ei löytynyt varsinaista reaktiopuuta, mutta joissakin näytteissä havaittiin lievää reaktiopuuta. Lievää reaktiopuuta löytyi enemmän männyn kuin kuusen juurista, eikä sitä löytynyt lainkaan kaikkein ohuimmista, noin 2 cm paksuisista juurenosista. Männyn kannoissa uuteaineprosentti oli korkeampi kuin kuusen. Männyn juurissa uuteaineprosentti kasvoi edettäessä kohti juuren kärkeä. Kuusella uuteaineprosentti laski aluksi, mutta lähellä juuren kärkeä taas kasvoi. Kuoren uuteainepitoisuus oli molemmilla puulajeilla korkeampi kuin puuaineen. Tutkimusaineisto oli suppea, eikä tutkimuksessa pyritty tilastolliseen yleistettävyyteen. Laajemmasta aineistosta tehdylle tutkimukselle on tarvetta, sillä turvekankailta saatavan puun tarjonta on Suomessa kasvussa, mutta juurten uuteaine- ja reaktiopuututkimuksia on tehty vain kivennäismailta kerätyistä aineistoista.