Numerical modelling of permafrost in bedrock in northern Fennoscandia during the Holocene

The occurrence of permafrost in bedrock in northern Fennoscandia and its dependence on past and presently ongoing climatic variations was investigated with one- (1D) and two-dimensional (2D) numerical models by solving the transient heat conduction equation with latent heat effects included. The study area is characterized by discontinuous permafrost occurrences such as palsa mires and local mountain permafrost. The ground temperature changes during the Holocene were constructed using climatic proxy data. This variation was used as a forcing function at the ground surface in the calculations. Several versions of the present ground temperature were applied, resulting in different subsurface freezing–thawing conditions in the past depending on the assumed porosity and geothermal conditions.Our results suggest that in high altitude areas with a cold climate (present mean annual ground temperature between 0°C and −3°C), there may have been considerable variations in permafrost thickness (ranging from 0 to 150 m), as well as periods of no permafrost at all. The higher is the porosity of bedrock filled with ice, the stronger is the retarding effect of permafrost against climatic variations.Two-dimensional models including topographic effects with altitude-dependent ground temperatures and slope orientation and inclination dependent solar radiation were applied to a case of mountain permafrost in Ylläs, western Finnish Lapland, where bedrock permafrost is known to occur in boreholes to a depth of about 60 m. Modelling suggests complicated changes in permafrost thickness with time as well as contrasting situations on southern and northern slopes of the mountain.Extrapolating the climatic warming of the last 200 years to the end of the next century when the anticipated increase in the annual average air temperature is expected to be about 2 K indicates that the permafrost occurrences in bedrock in northern Fennoscandia would be thawing rapidly in low-porosity formations. However, already a porosity of 5% filled with ice would retard the thawing considerably.     

An Example of Island-Arc Petrogenesis: Geochemistry and Petrology of the Southern Luzon Arc, Philippines

Volcanism throughout the Luzon arc is associated with eastwardsubduction of the South China Sea floor along the Manila Trench.The southern section of the arc, the focus of this study, extendsfrom the Lingayen-Dingalan fault to the small islands just southof Luzon. Two segments appear to exist along this section ofthe arc the northern Bataan and southern Mindoro segments whichare separated by the Macolod Corridor. The volcanic rocks have typical arc phenocryst mineralogies:olivine, clinopyroxene, plagioclase, and titanomagnetite inthe most mafic rocks and clinopyroxene, plagioclase, orthopyroxene,titanomagnetite, ? amphibole in the more felsic samples. Complexzoning, sieve textures, and decoupling of incompatible traceelements suggest that processes such as assimilation have takenplace. The rocks from the study area range from basalts to rhyolitesand show typical calc-alkaline features. The rocks of the MacolodCorridor and Mindoro segment are particularly enriched in largeion lithophile elements (LILE), light rare earth elements (LREE),and radiogenic Sr compared with the Bataan segment. The datafall within the mantle array on Sr-Nd isotopic diagrams andgrade toward higher Sr and lower Nd isotopic values from northto south. A likely source for the volcanics of this study is either amid-ocean ridge basalt (MORB)-type mantle that undergoes higherdegrees of partial melting than regions involved in MORB generationor a previously depleted source. We suggest that the high fieldstrength element (HFSE) anomalies have been derived throughdifferential element partitioning during fluid transport fromthe subducted lithosphere to the mantle wedge. Continental crustal material seems to play a significant roleparticularly in the Macolod Corridor and the Mindoro segment,based on the high LILE, La/Sm ratios, radiogenic Sr isotopes,and ¹⁸O values. The Macolod Corridor and the Mindoro segmenthave undergone source contamination by crustal material fromthe North Palawan-Mindoro crustal block either during the collisionof this block with the Manila Trench or by subduction of sedimentsrich in this crustal material. A similar component has alsobeen detected in the Bataan segment but in minor amounts. Thetrace element and isotopic differences between the northernand southern sections of the arc are interpreted in terms ofvariable composition (i. e., variable amounts of a crustal componentintroduced from the Palawan-Mindoro crustal terrain) of themetasomatic fluids released into the source.     

The relation between molecular clouds and stellar kinematics

The relation between molecular clouds, star clusters, and the stellar component of the galactic disk is investigated. According to Elmegreen (1985) bound stellar systems, e.g., open star clusters, can be formed from molecular cloud of mass 10⁴ M . A close encounter with a giant molecular cloud or massive black hole disrupts such stellar systems and forms superclusters. This explains why some open star clusters are so mass-deficient. Unbound stellar systems, e.g., expanding OB associations, are formed from molecular clouds of mass 10⁵ M . When disruptive O-type stars appear the star formation is halted and the cloud is destroyed. An example of the relict of GMC disruption in the solar vicinity is Gould's belt. The velocity dispersion-versus-age relation is also investigated and explained as a consequence of gravitational scattering of stars on GMC, or massive black holes, or as due to recurrent transient spirals.Paper presented at a Workshop on The Role of Dust in Dense Regions of Interstellar Matter, held at Georgenthal, G.D.R., in March 1986.     

Opening of the Fram Strait gateway: A review of plate tectonic constraints

We have revised the regional crustal structure, oceanic age distribution, and conjugate margin segmentation in and around the Lena Trough, the oceanic part of the Fram Strait between the Norwegian–Greenland Sea and the Eurasia Basin (Arctic Ocean). The Lena Trough started to open after Eurasia–Greenland relative plate motions changed from right-lateral shear to oblique divergence at Chron 13 times (33.3 Ma; earliest Oligocene). A new Bouguer gravity map, supported by existing seismic data and aeromagnetic profiles, has been applied to interpret the continent–ocean transition and the influence of Eocene shear structures on the timing of breakup and initial seafloor spreading. Assuming that the onset of deep-water exchange depended on the formation of a narrow, oceanic corridor, the gateway formed during early Miocene times (20–15 Ma). However, if the initial Lena Trough was blocked by terrigenous sediments or was insufficiently subsided to allow for deep-water circulation, the gateway probably formed with the first well developed magnetic seafloor spreading anomaly around Chron 5 times (9.8 Ma; Late Miocene). Paleoceanographic changes at ODP Site 909 (northern Hovgård Ridge) are consistent with both hypotheses of gateway formation. We cannot rule out that a minor gateway formed across stretched continental crust prior to the onset of seafloor spreading in the Lena Trough. The gravity, seismic and magnetic observations question the prevailing hypotheses on the Yermak Plateau and the Morris Jesup Rise as Eocene oceanic plateaus and the Hovgård Ridge as a microcontinent.     

Determining effective wellbore permeability from a field pressure test: a numerical analysis of detection limits

We propose a simple pressure test that can be used in the field to determine the effective permeability of existing wellbores. Such tests are motivated by the need to understand and quantify leakage risks associated with geological storage of CO₂ in mature sedimentary basins. If CO₂ is injected into a deep geological formation, and the resulting CO₂ plume encounters a wellbore, leakage may occur through various pathways in the “disturbed zone” surrounding the well casing. The effective permeability of this composite zone, on the outside of the well casing, is an important parameter for models of leakage. However, the data that exist on this key parameter do not exist in the open literature, and therefore specific field tests need to be done in order to reduce the uncertainty inherent in the leakage estimates. The test designed and analyzed herein is designed to measure effective wellbore permeability within a low-permeability caprock, bounded above and below by permeable reservoirs, by pressurizing the reservoir below and measuring the response in the reservoir above. Alternatively, a modified test can be performed within the caprock without directly contacting the reservoirs above and below. We use numerical simulation to relate pressure response to effective well permeability and then evaluate the range of detection of the effective permeability based on instrument measurement error and limits on fracture pressure. These results can guide field experiments associated with site characterization and leakage analysis.     

Modeling of water–rock interactions in an aquitard of sandy clay in the coastal area near Beihai,China

Unconsolidated sand, gravel and clay deposits near Beihai and in the Leizhou Peninsula in southern China form an unconfined aquifer, aquitard and a confined aquifer. Water and soil samples were collected from the two aquifers in the coastal Beihai area for the determination of chemical compositions, minerals and soluble ions. Hydrogeochemical modeling of three flow paths through the aquitard are carried out using PHREEQC to determine water–rock interactions along the flow paths. The results indicate that the dissolution of anorthite, fluorite, halite, rhodochrosite and CO₂, and precipitation of potash feldspar and kaolinite may be occurring when groundwater leaks through the aquitard from the unconfined aquifer to the confined aquifer. Cation exchanges between Na and Ca can also happen along the flow paths.     

Reformulated O correction of mass spectrometric stable isotope measurements in carbon dioxide and a critical appraisal of historic ‘absolute’ carbon and oxygen isotope ratios

Mass-spectrometric stable isotope measurements of CO₂ use molecular ion currents at mass-to-charge ratios m/z 44, 45 and 46 to derive the elemental isotope ratios n(¹³C)/n(¹²C) and n(¹⁸O)/n(¹⁶O), abbreviated ¹³C/¹²C and ¹⁸O/¹⁶O, relative to a reference. The ion currents have to be corrected for the contribution of ¹⁷O-bearing isotopologues, the so-called ‘¹⁷O correction’. The magnitude of this correction depends on the calibrated isotope ratios of the reference. Isotope ratio calibrations are difficult and are therefore a matter of debate. Here, I provide a comprehensive evaluation of the existing ¹³C/¹²C (¹³R), ¹⁷O/¹⁶O (¹⁷R) and ¹⁸O/¹⁶O (¹⁸R) calibrations of the reference material Vienna Standard Mean Ocean Water (VSMOW) and CO₂ generated from the reference material Vienna Pee Dee Belemnite (VPDB) by reaction with 100% H₃PO₄ at 25 °C (VPDB-CO₂). I find , ¹⁸R_VSMOW/10⁻⁶ = 2005.20 ± 0.45, ¹³R_VPDB-CO2/10^-6= 11124 ± 45, and ¹⁸R_VPDB-CO2/10^-6=2088.37±0.90. I also rephrase the calculation scheme for the ¹⁷O correction completely in terms of relative isotope ratio differences (δ values). This reveals that only ratios of isotope ratios (namely, ¹⁷R/¹³R and ¹³R¹⁷R/¹⁸R) are required for the ¹⁷O correction. These can be, and have been, measured on conventional stable isotope mass spectrometers. I then show that the remaining error for these ratios of isotope ratios can lead to significant uncertainty in the derived relative ¹³C/¹²C difference, but not for¹⁸O/¹⁶O. Even though inter-laboratory differences can be corrected for by a common ‘ratio assumption set’ and/or normalisation, the ultimate accuracy of the ¹⁷O correction is hereby limited. Errors of similar magnitude can be introduced by the assumed mass-dependent relationship between ¹⁷O/¹⁶O and ¹⁸O/¹⁶O isotope ratios. For highest accuracy in the ¹³C/¹²C ratio, independent triple oxygen isotope measurements are required. Finally, I propose an experiment that allows direct measurement of ¹³R¹⁷R/¹⁸R.     

Deep-seated gravitational slope deformations and their influence on consequent mass movements (case studies from the highest part of the Czech Carpathians)

Gravitational spreading of mountain ridges displays primary disequilibrium of flysch mountain areas of the Czech Carpathians. The progression of various types of mass movements is a product of long-term ridge disintegration and is predisposed by the geological structure of the area and the upper Tertiary-Quaternary morphogenesis of the mountain area. Deep-seated slope deformations are spatially interconnected by the occurrence of some other types of slope deformations (e.g. debris flows, debris slides, slumps, rock avalanches, etc.), which pose a considerable risk for the existence of human society. An important causative factor in these dynamically developing hazardous processes is, among other factors, the way in which land has been used in the last three centuries. Therefore, the occurrence of various types of slope deformations is studied in terms of their relation to deep-seated gravitational deformations and in terms of other limiting factors (structural geological, morphological and climatic factors, manmade impacts, etc.). The paper presents several case studies of slope deformations (Velká Čantoryje Mt, Lysá hora Mt, Ropice Mt and Smrk Mt) in the area of the Outer Carpathians within the territory of the Czech Republic and also adverts to some consequences in terms of the socioeconomic structure of the landscape.