Autoecología de los hayedos catalanes

Este estudio está basado en el muestreo de campo y posterior análisis de 24 parcelas de hayedo seleccionadas mediante una estratificación de su área de distribución basada en la clasificación CLATERES de la Ecorregión Catalano-Aragonesa. En cada parcela se han evaluado 3 parámetros fisiográficos, 15 climáticos y 18 edáficos, a partir de los cuales se han establecido sus valores paramétricos centrales y marginales que permiten definir los hábitats fisiográfico, climático y edáfico de las masas de Fagus sylvatica L. en Cataluña. Los hayedos catalanes se presentan sobre substratos litológicos muy diversos (plutonitas, vulcanitas, metamorfitas y sedimentitas, tanto ácidas como básicas), con texturas predominantes francas, franco-arenosas o franco-limosas. Los suelos, según FAO, son mayoritariamente cambisoles. A pesar de que la capacidad de retención de agua de sus suelos es escasa, la sequía fisiológica es reducida. Los humus predominantemente pertenecen a los tipos mull forestal y mull cálcico. Además, se presentan una serie de parámetros selvícolas ( Densidad de pies y densidad de chirpiales, Area basimétrica, Altura Total dominante, Índices de Hart-Becking, Índice de Calidad de Estación y Edad de la masa) que al correlacionarlos con los ecológicos nos ha permitido comprobar que los mejores hayedos se encuentran en las localizaciones más térmicas, en las que incluso se podría producir sequía fisiológica si no fuera por que existen suficientes precipitaciones estivales.

The Hercynian intrusive rocks of the Catalonian Coastal Ranges (NE Spain).

The hercynian outcrops of the Catalonian Coastal Ranges (NE Spain) consist mainly of Lower Permian-Upper Carboniferous, post-tectonic, epizona1 granitoid intrusions which form a typical applutonic calc-alkaline suite ranging from mafic hornblende gabbros and ultramafic olivine homblendites throught on alites and granodioritcs to leucogranites. This suite displays major andtrace-element characteristics and Sr isotope ratios similar to volcanic arc and post-collision magmatism oceanic lithosphere and to have been modified by contamination and is therefore believed to have formed above subducted with melts from the crust

Hercynian Metamorphism in the Catalonian Coastal Ranges

Paleozoic rocks in the Catalonian Coastal Ranges are in their largestpart affectedby alow-tovery-low grade Hercynian metamorphism. Amphibolite facies conditions are only found in restricted areas such as the southwestern part of the Guilleries massif where upper amphibolite facies conditions are reached. Metamorphic grade increases from top to bottom of the Paleozoic stratigraphic sequence and the metamorphic peak is diachronous, being progressively older in the lower grade metamorphic zones. The isograd pattern, mineral assemblages, mineral chemistry and preserved reaction textures are consistent with a low pressure metamorphism possibly evolving from a previous Barrovian type event. The metamorphic climax in the high grade zone was reached after the seconddeformational phase. Calculatedpeak P-Tconditions are 620-640 OC and around 3.5 Kb . A latter episode of decompression from the maximum conditions to 1-2 Kb, with an associated temperature decrease to 530-550 OC, is recognized. The intrusion of late Hercynian granitoids produced contact metamorphic aureoles where the pyroxene-hornfels facies is locally reached.

The Sottunga-Jurmo shear zone : structure and deformation history of a crustal-scale ductile shear zone in SW Finland

The aim of this study is to gain a better understanding of the structure and the deformation history of a NW-SE trending regional, crustal-scale shear structure in the Åland archipelago, SW Finland, called the Sottunga-Jurmo shear zone (SJSZ). Approaches involving e.g. structural geology, geochronology, geochemistry and metamorphic petrology were utilised in order to reconstruct the overall deformation history of the study area. The study therefore describes several features of the shear zone including structures, kinematics and lithologies within the study area, the ages of the different deformation phases (ductile to brittle) within the shear zone, as well as some geothermobarometric results. The results indicate that the SJSZ outlines a major crustal discontinuity between the extensively migmatized rocks NE of the shear zone and the unmigmatised, amphibolite facies rocks SW of the zone. The main SJSZ shows overall dextral lateral kinematics with a SW-side up vertical component and deformation partitioning into pure shear and simple shear dominated deformation styles that was intensified toward later stages of the deformation history. The deformation partitioning resulted in complex folding and refolding against the SW margin of the SJSZ, including conical and sheath folds, and in a formation of several minor strike-slip shear zones both parallel and conjugate to the main SJSZ in order to accommodate the regional transpressive stresses. Different deformation phases within the study area were dated by SIMS (zircon U-Pb), ID-TIMS (titanite U-Pb) and 40Ar/39Ar (pseudotachylyte wholerock) methods. The first deformation phase within the ca. 1.88 Ga rocks of the study area is dated at ca. 1.85 Ga, and the shear zone was reactivated twice within the ductile regime (at ca. 1.83 Ga and 1.79 Ga), during which the strain was successively increasingly partitioned into the main SJSZ and the minor shear zones. The age determinations suggest that the orogenic processes within the study area did not occur in a temporal continuum; instead, the metamorphic zircon rims and titanites show distinct, 10-20 Ma long breaks in deformation between phases of active deformation. The results of this study further imply slow cooling of the rocks through 600-700ºC so that at 1.79 Ga, 2 the temperature was still at least 600ºC. The highest recorded metamorphic pressures are 6.4-7.1 kbar. At the late stages or soon after the last ductile phase (ca. 1.79 Ga), relatively high-T mylonites and ultramylonites were formed, witnessing extreme deformation partitioning and high strain rates. After the rocks reached lower amphibolite facies to amphibolite-greenschist facies transitional conditions (ca. 500-550ºC), they cooled rapidly, probably due to crustal uplift and exhumation. The shear zone was reactivated at least once within the semi-brittle to brittle regime between ca. 1.79 Ga and 1.58 Ga, as evidenced by cataclasites and pseudotachylytes. In summary, the results of this study suggest that the Sottunga-Jurmo shear zone (and the South Finland shear zone) defines a major crustal discontinuity, and played a central role in accommodating the regional stresses during and after the Svecofennian orogeny.

Improved predictive mapping of indoor radon concentrations using ensemble regression trees based on automatic clustering of geological units.

PURPOSE: According to estimations around 230 people die as a result of radon exposure in Switzerland. This public health concern makes reliable indoor radon prediction and mapping methods necessary in order to improve risk communication to the public. The aim of this study was to develop an automated method to classify lithological units according to their radon characteristics and to develop mapping and predictive tools in order to improve local radon prediction. METHOD: About 240 000 indoor radon concentration (IRC) measurements in about 150 000 buildings were available for our analysis. The automated classification of lithological units was based on k-medoids clustering via pair-wise Kolmogorov distances between IRC distributions of lithological units. For IRC mapping and prediction we used random forests and Bayesian additive regression trees (BART). RESULTS: The automated classification groups lithological units well in terms of their IRC characteristics. Especially the IRC differences in metamorphic rocks like gneiss are well revealed by this method. The maps produced by random forests soundly represent the regional difference of IRCs in Switzerland and improve the spatial detail compared to existing approaches. We could explain 33% of the variations in IRC data with random forests. Additionally, the influence of a variable evaluated by random forests shows that building characteristics are less important predictors for IRCs than spatial/geological influences. BART could explain 29% of IRC variability and produced maps that indicate the prediction uncertainty. CONCLUSION: Ensemble regression trees are a powerful tool to model and understand the multidimensional influences on IRCs. Automatic clustering of lithological units complements this method by facilitating the interpretation of radon properties of rock types. This study provides an important element for radon risk communication. Future approaches should consider taking into account further variables like soil gas radon measurements as well as more detailed geological information.

Precisiones acerca de la significación petrológica y estructural de las rocas gneísicas y cataclásticas del Maresme (prov. de Barcelona)

Se estudian pequeños afloramientos de rocas metamórficas y cataclásticas ubicadas en el granito herciniano de la Cadena Costera Catalana al NE de Barcelona. Tras una breve íntesis de las diversas ideas emitidas sobre la génesis de estas rocas se realiza un estudio comparativo de las mismas. De él se infiere la presencia de dos grupos de rocas de significación petrogenética y tectónica distinta: 1) Los gneises de Mataró, que resultan del metamorfismo polifísico progresivo herciniano de sedimentos del Paleozoico inferior y 2) las rocas cataclásticas de Caldetes, ligadas a deformaciones internas en el batolito del granito herciniano postectónico, acompañadas de transformaciones hidrotermales en mayor o menor grado.

Seismic characterization of heterogeneous sedimentary and fractured rocks based on Biot's theory of poroelasticity

Understanding and quantifying seismic energy dissipation, which manifests itself in terms of velocity dispersion and attenuation, in fluid-saturated porous rocks is of considerable interest, since it offers the perspective of extracting information with regard to the elastic and hydraulic rock properties. There is increasing evidence to suggest that wave-induced fluid flow, or simply WIFF, is the dominant underlying physical mechanism governing these phenomena throughout the seismic, sonic, and ultrasonic frequency ranges. This mechanism, which can prevail at the microscopic, mesoscopic, and macroscopic scale ranges, operates through viscous energy dissipation in response to fluid pressure gradients and inertial effects induced by the passing wavefield. In the first part of this thesis, we present an analysis of broad-band multi-frequency sonic log data from a borehole penetrating water-saturated unconsolidated glacio-fluvial sediments. An inherent complication arising in the interpretation of the observed P-wave attenuation and velocity dispersion is, however, that the relative importance of WIFF at the various scales is unknown and difficult to unravel. An important generic result of our work is that the levels of attenuation and velocity dispersion due to the presence of mesoscopic heterogeneities in water-saturated unconsolidated clastic sediments are expected to be largely negligible. Conversely, WIFF at the macroscopic scale allows for explaining most of the considered data while refinements provided by including WIFF at the microscopic scale in the analysis are locally meaningful. Using a Monte-Carlo-type inversion approach, we compare the capability of the different models describing WIFF at the macroscopic and microscopic scales with regard to their ability to constrain the dry frame elastic moduli and the permeability as well as their local probability distribution. In the second part of this thesis, we explore the issue of determining the size of a representative elementary volume (REV) arising in the numerical upscaling procedures of effective seismic velocity dispersion and attenuation of heterogeneous media. To this end, we focus on a set of idealized synthetic rock samples characterized by the presence of layers, fractures or patchy saturation in the mesocopic scale range. These scenarios are highly pertinent because they tend to be associated with very high levels of velocity dispersion and attenuation caused by WIFF in the mesoscopic scale range. The problem of determining the REV size for generic heterogeneous rocks is extremely complex and entirely unexplored in the given context. In this pilot study, we have therefore focused on periodic media, which assures the inherent self- similarity of the considered samples regardless of their size and thus simplifies the problem to a systematic analysis of the dependence of the REV size on the applied boundary conditions in the numerical simulations. Our results demonstrate that boundary condition effects are absent for layered media and negligible in the presence of patchy saturation, thus resulting in minimum REV sizes. Conversely, strong boundary condition effects arise in the presence of a periodic distribution of finite-length fractures, thus leading to large REV sizes. In the third part of the thesis, we propose a novel effective poroelastic model for periodic media characterized by mesoscopic layering, which accounts for WIFF at both the macroscopic and mesoscopic scales as well as for the anisotropy associated with the layering. Correspondingly, this model correctly predicts the existence of the fast and slow P-waves as well as quasi and pure S-waves for any direction of wave propagation as long as the corresponding wavelengths are much larger than the layer thicknesses. The primary motivation for this work is that, for formations of intermediate to high permeability, such as, for example, unconsolidated sediments, clean sandstones, or fractured rocks, these two WIFF mechanisms may prevail at similar frequencies. This scenario, which can be expected rather common, cannot be accounted for by existing models for layered porous media. Comparisons of analytical solutions of the P- and S-wave phase velocities and inverse quality factors for wave propagation perpendicular to the layering with those obtained from numerical simulations based on a ID finite-element solution of the poroelastic equations of motion show very good agreement as long as the assumption of long wavelengths remains valid. A limitation of the proposed model is its inability to account for inertial effects in mesoscopic WIFF when both WIFF mechanisms prevail at similar frequencies. Our results do, however, also indicate that the associated error is likely to be relatively small, as, even at frequencies at which both inertial and scattering effects are expected to be at play, the proposed model provides a solution that is remarkably close to its numerical benchmark. -- Comprendre et pouvoir quantifier la dissipation d'énergie sismique qui se traduit par la dispersion et l'atténuation des vitesses dans les roches poreuses et saturées en fluide est un intérêt primordial pour obtenir des informations à propos des propriétés élastique et hydraulique des roches en question. De plus en plus d'études montrent que le déplacement relatif du fluide par rapport au solide induit par le passage de l'onde (wave induced fluid flow en anglais, dont on gardera ici l'abréviation largement utilisée, WIFF), représente le principal mécanisme physique qui régit ces phénomènes, pour la gamme des fréquences sismiques, sonique et jusqu'à l'ultrasonique. Ce mécanisme, qui prédomine aux échelles microscopique, mésoscopique et macroscopique, est lié à la dissipation d'énergie visqueuse résultant des gradients de pression de fluide et des effets inertiels induits par le passage du champ d'onde. Dans la première partie de cette thèse, nous présentons une analyse de données de diagraphie acoustique à large bande et multifréquences, issues d'un forage réalisé dans des sédiments glaciaux-fluviaux, non-consolidés et saturés en eau. La difficulté inhérente à l'interprétation de l'atténuation et de la dispersion des vitesses des ondes P observées, est que l'importance des WIFF aux différentes échelles est inconnue et difficile à quantifier. Notre étude montre que l'on peut négliger le taux d'atténuation et de dispersion des vitesses dû à la présence d'hétérogénéités à l'échelle mésoscopique dans des sédiments clastiques, non- consolidés et saturés en eau. A l'inverse, les WIFF à l'échelle macroscopique expliquent la plupart des données, tandis que les précisions apportées par les WIFF à l'échelle microscopique sont localement significatives. En utilisant une méthode d'inversion du type Monte-Carlo, nous avons comparé, pour les deux modèles WIFF aux échelles macroscopique et microscopique, leur capacité à contraindre les modules élastiques de la matrice sèche et la perméabilité ainsi que leur distribution de probabilité locale. Dans une seconde partie de cette thèse, nous cherchons une solution pour déterminer la dimension d'un volume élémentaire représentatif (noté VER). Cette problématique se pose dans les procédures numériques de changement d'échelle pour déterminer l'atténuation effective et la dispersion effective de la vitesse sismique dans un milieu hétérogène. Pour ce faire, nous nous concentrons sur un ensemble d'échantillons de roches synthétiques idéalisés incluant des strates, des fissures, ou une saturation partielle à l'échelle mésoscopique. Ces scénarios sont hautement pertinents, car ils sont associés à un taux très élevé d'atténuation et de dispersion des vitesses causé par les WIFF à l'échelle mésoscopique. L'enjeu de déterminer la dimension d'un VER pour une roche hétérogène est très complexe et encore inexploré dans le contexte actuel. Dans cette étude-pilote, nous nous focalisons sur des milieux périodiques, qui assurent l'autosimilarité des échantillons considérés indépendamment de leur taille. Ainsi, nous simplifions le problème à une analyse systématique de la dépendance de la dimension des VER aux conditions aux limites appliquées. Nos résultats indiquent que les effets des conditions aux limites sont absents pour un milieu stratifié, et négligeables pour un milieu à saturation partielle : cela résultant à des dimensions petites des VER. Au contraire, de forts effets des conditions aux limites apparaissent dans les milieux présentant une distribution périodique de fissures de taille finie : cela conduisant à de grandes dimensions des VER. Dans la troisième partie de cette thèse, nous proposons un nouveau modèle poro- élastique effectif, pour les milieux périodiques caractérisés par une stratification mésoscopique, qui prendra en compte les WIFF à la fois aux échelles mésoscopique et macroscopique, ainsi que l'anisotropie associée à ces strates. Ce modèle prédit alors avec exactitude l'existence des ondes P rapides et lentes ainsi que les quasis et pures ondes S, pour toutes les directions de propagation de l'onde, tant que la longueur d'onde correspondante est bien plus grande que l'épaisseur de la strate. L'intérêt principal de ce travail est que, pour les formations à perméabilité moyenne à élevée, comme, par exemple, les sédiments non- consolidés, les grès ou encore les roches fissurées, ces deux mécanismes d'WIFF peuvent avoir lieu à des fréquences similaires. Or, ce scénario, qui est assez commun, n'est pas décrit par les modèles existants pour les milieux poreux stratifiés. Les comparaisons des solutions analytiques des vitesses des ondes P et S et de l'atténuation de la propagation des ondes perpendiculaires à la stratification, avec les solutions obtenues à partir de simulations numériques en éléments finis, fondées sur une solution obtenue en 1D des équations poro- élastiques, montrent un très bon accord, tant que l'hypothèse des grandes longueurs d'onde reste valable. Il y a cependant une limitation de ce modèle qui est liée à son incapacité à prendre en compte les effets inertiels dans les WIFF mésoscopiques quand les deux mécanismes d'WIFF prédominent à des fréquences similaires. Néanmoins, nos résultats montrent aussi que l'erreur associée est relativement faible, même à des fréquences à laquelle sont attendus les deux effets d'inertie et de diffusion, indiquant que le modèle proposé fournit une solution qui est remarquablement proche de sa référence numérique.

Muuttuvat tilat lasten- ja nuortenkirjallisuuden tutkijoiden pohdittavana : : metamorphic spaces in children's literature and culture 17.-19.1.2008 Jyväskylä

Ajankohtaista

Evoluutio Rocks!

Kirjallisuusarvostelu

Reactive dissolution of sedimentary rocks in flue gas desulfurization : modeling and experimental investigation

Svavel förekommer i kol och olja och oxideras vid förbränning till svaveldioxid (SO2). Årligen utsläpps stora mängder svaveldioxid som åstadkommer sura regn, minskning av stratosfäriskt ozon och sjukdomar. Av dessa orsaker är det nödvändigt att minska utsläppen av svaveldioxid. Den teknologi som rör planering av reaktorer och processer för rökgasavsvavling (FGD) har utvecklats kraftigt och idag använder man olika typer av lösningar. De vanligaste alternativen för FGD är våtskrubber- och semitorra skrubbersystem (Spray Dry Scrubbers, SDS) och injektionsprocesser för absorbenter. SDS-processer har studerats bland annat av Ahlbeck [4] och Klingspor [5] och injektionsprocesser av Alvfors [6]. Kalksten, som i huvusak består av kalciumkarbonat, används i rökgasavsvavling på grund av sin förmåga att binda svavel i form av sulfatsalter. Den vanligaste rökgasavsvavlingsmetoden är våtskrubbning, där det sedimentära stenmaterialets upplösningshastighet är en av de faktorer som påverkar resultatet mest. Utvärdering av kalkstensreaktivitet är därför speciellt viktig vid planering och drift av anläggningar för rökgasavsvavling. Målsättningen med detta arbete var att modellera upplösningen av olika typers kalksten för att få en kvantitativ utvärdering av kvaliteten på de analyserade proverna. Därtill testades även karbonatbiprodukter från stålindustrin för att utvärdera möjligheter att använda andra råmaterial. Det transienta förloppet har analyserats, varvid upplösningshastigheten modellerades bl.a. i avseende på tid och pH. Under arbetets gång har antalet empiriska korrelationer minskats till fördel för fysikaliska modeller av diffusiva och konvektiva masstransportfenomen. En målsättning var att skapa en effektiv och snabb metod för att testa olika absorbenter för rökgasavsvavling under transienta förlopp. I arbetet användes PSD-analys, gjordes pH-mätningar och andra utvärderingar av de fysikaliska parametrar som ingår i beräkningarna. On-line mätningar för de icke-stationära variablerna tid och pH ger möjlighet att eliminera osäkerheter. Vissa modeller kan vara komplicerade. En modell för upplösningshastigheten med mer detaljerad utvärdering av parametrar och färre approximationer är därför nödvändig då man vill utvärdera reaktionshastigheten för fasta partiklar i sur miljö. Arbetet utfördes under fyra år och fem peer review-artiklar ingår i avhandlingen.

The Sottunga-Jurmo shear zone : structure and deformation history of a crustal-scale ductile shear zone in SW Finland

The aim of this study is to gain a better understanding of the structure and the deformation history of a NW-SE trending regional, crustal-scale shear structure in the Åland archipelago, SW Finland, called the Sottunga-Jurmo shear zone (SJSZ). Approaches involving e.g. structural geology, geochronology, geochemistry and metamorphic petrology were utilised in order to reconstruct the overall deformation history of the study area. The study therefore describes several features of the shear zone including structures, kinematics and lithologies within the study area, the ages of the different deformation phases (ductile to brittle) within the shear zone, as well as some geothermobarometric results. The results indicate that the SJSZ outlines a major crustal discontinuity between the extensively migmatized rocks NE of the shear zone and the unmigmatised, amphibolite facies rocks SW of the zone. The main SJSZ shows overall dextral lateral kinematics with a SW-side up vertical component and deformation partitioning into pure shear and simple shear dominated deformation styles that was intensified toward later stages of the deformation history. The deformation partitioning resulted in complex folding and refolding against the SW margin of the SJSZ, including conical and sheath folds, and in a formation of several minor strike-slip shear zones both parallel and conjugate to the main SJSZ in order to accommodate the regional transpressive stresses. Different deformation phases within the study area were dated by SIMS (zircon U-Pb), ID-TIMS (titanite U-Pb) and 40Ar/39Ar (pseudotachylyte wholerock) methods. The first deformation phase within the ca. 1.88 Ga rocks of the study area is dated at ca. 1.85 Ga, and the shear zone was reactivated twice within the ductile regime (at ca. 1.83 Ga and 1.79 Ga), during which the strain was successively increasingly partitioned into the main SJSZ and the minor shear zones. The age determinations suggest that the orogenic processes within the study area did not occur in a temporal continuum; instead, the metamorphic zircon rims and titanites show distinct, 10-20 Ma long breaks in deformation between phases of active deformation. The results of this study further imply slow cooling of the rocks through 600-700ºC so that at 1.79 Ga, 2 the temperature was still at least 600ºC. The highest recorded metamorphic pressures are 6.4-7.1 kbar. At the late stages or soon after the last ductile phase (ca. 1.79 Ga), relatively high-T mylonites and ultramylonites were formed, witnessing extreme deformation partitioning and high strain rates. After the rocks reached lower amphibolite facies to amphibolite-greenschist facies transitional conditions (ca. 500-550ºC), they cooled rapidly, probably due to crustal uplift and exhumation. The shear zone was reactivated at least once within the semi-brittle to brittle regime between ca. 1.79 Ga and 1.58 Ga, as evidenced by cataclasites and pseudotachylytes. In summary, the results of this study suggest that the Sottunga-Jurmo shear zone (and the South Finland shear zone) defines a major crustal discontinuity, and played a central role in accommodating the regional stresses during and after the Svecofennian orogeny.