Refine
Document Type
- Doctoral Thesis (9)
- Article (3)
- Master's Thesis (1)
Is part of the Bibliography
- yes (13)
Keywords
- Diffusion (13) (remove)
During lower sea levels in glacial periods, deep permafrost formed on large continental shelf areas of the Arctic Ocean. Subsequent sea level rise and coastal erosion created subsea permafrost, which generally degrades after inundation under the influence of a complex suite of marine, near-shore processes. Global warming is especially pronounced in the Arctic, and will increase the transition to and the degradation of subsea permafrost, with implications for atmospheric climate forcing, offshore infrastructure, and aquatic ecosystems.
This thesis combines new geophysical, borehole observational and modelling approaches to enhance our understanding of subsea permafrost dynamics. Three specific areas for advancement were identified: (I) sparsity of observational data, (II) lacking implementation of salt infiltration mechanisms in models, and (III) poor understanding of the regional differences in key driving parameters. This study tested the combination of spectral ratios of the ambient vibration seismic wavefield, together with estimated shear wave velocity from seismic interferometry analysis, for estimating the thickness of the unfrozen sediment overlying the ice-bonded permafrost offshore. Mesoscale numerical calculations (10^1 to 10^2 m, thousands of years) were employed to develop and solve the coupled heat diffusion and salt transport equations including phase change effects. Model soil parameters were constrained by borehole data, and the impact of a variety of influences during the transgression was tested in modelling studies. In addition, two inversion schemes (particle swarm optimization and a least-square method) were used to reconstruct temperature histories for the past 200-300 years in the Laptev Sea region in Siberia from two permafrost borehole temperature records. These data were evaluated against larger scale reconstructions from the region.
It was found (I) that peaks in spectral ratios modelled for three-layer, one-dimensional systems corresponded with thaw depths. Around Muostakh Island in the central Laptev Sea seismic receivers were deployed on the seabed. Derived depths of the ice-bonded permafrost table were between 3.7-20.7 m ± 15 %, increasing with distance from the coast. (II) Temperatures modelled during the transition to subsea permafrost resembled isothermal conditions after about 2000 years of inundation at Cape Mamontov Klyk, consistent with observations from offshore boreholes. Stratigraphic scenarios showed that salt distribution and infiltration had a large impact on the ice saturation in the sediments. Three key factors were identified that, when changed, shifted the modelled permafrost thaw depth most strongly: bottom water temperatures, shoreline retreat rate and initial temperature before inundation. Salt transport based on diffusion and contribution from arbitrary density-driven mechanisms only accounted for about 50 % of observed thaw depths at offshore sites hundreds to thousands of years after inundation. This bias was found consistently at all three sites in the Laptev Sea region. (III) In the temperature reconstructions, distinct differences in the local temperature histories between the western Laptev Sea and the Lena Delta sites were recognized, such as a transition to warmer temperatures a century later in the western Laptev Sea as well as a peak in warming three decades later. The local permafrost surface temperature history at Sardakh Island in the Lena Delta was reminiscent of the circum-Arctic regional average trends. However, Mamontov Klyk in the western Laptev Sea was consistent to Arctic trends only in the most recent decade and was more similar to northern hemispheric mean trends. Both sites were consistent with a rapid synoptic recent warming.
In conclusion, the consistency between modelled response, expected permafrost distribution, and observational data suggests that the passive seismic method is promising for the determination of the thickness of unfrozen sediment on the continental Arctic shelf. The quantified gap between currently modelled and observed thaw depths means that the impact of degradation on climate forcing, ecosystems, and infrastructure is larger than current models predict. This discrepancy suggests the importance of further mechanisms of salt penetration and thaw that have not been considered – either pre-inundation or post-inundation, or both. In addition, any meaningful modelling of subsea permafrost would have to constrain the identified key factors and their regional differences well. The shallow permafrost boreholes provide missing well-resolved short-scale temperature information in the coastal permafrost tundra of the Arctic. As local differences from circum-Arctic reconstructions, such as later warming and higher warming magnitude, were shown to exist in this region, these results provide a basis for local surface temperature record parameterization of climate and, in particular, permafrost models. The results of this work bring us one step further to understanding the full picture of the transition from terrestrial to subsea permafrost.
The two hallmark features of Brownian motion are the linear growth < x2(t)> = 2Ddt of the mean squared displacement (MSD) with diffusion coefficient D in d spatial dimensions, and the Gaussian distribution of displacements. With the increasing complexity of the studied systems deviations from these two central properties have been unveiled over the years. Recently, a large variety of systems have been reported in which the MSD exhibits the linear growth in time of Brownian (Fickian) transport, however, the distribution of displacements is pronouncedly non-Gaussian (Brownian yet non-Gaussian, BNG). A similar behaviour is also observed for viscoelastic-type motion where an anomalous trend of the MSD, i.e., <x2(t)> ~ ta, is combined with a priori unexpected non-Gaussian distributions (anomalous yet non-Gaussian, ANG). This kind of behaviour observed in BNG and ANG diffusions has been related to the presence of heterogeneities in the systems and a common approach has been established to address it, that is, the random diffusivity approach.
This dissertation explores extensively the field of random diffusivity models. Starting from a chronological description of all the main approaches used as an attempt of describing BNG and ANG diffusion, different mathematical methodologies are defined for the resolution and study of these models. The processes that are reported in this work can be classified in three subcategories, i) randomly-scaled Gaussian processes, ii) superstatistical models and iii) diffusing diffusivity models, all belonging to the more general class of random diffusivity models. Eventually, the study focuses more on BNG diffusion, which is by now well-established and relatively well-understood. Nevertheless, many examples are discussed for the description of ANG diffusion, in order to highlight the possible scenarios which are known so far for the study of this class of processes.
The second part of the dissertation deals with the statistical analysis of random diffusivity processes. A general description based on the concept of moment-generating function is initially provided to obtain standard statistical properties of the models. Then, the discussion moves to the study of the power spectral analysis and the first passage statistics for some particular random diffusivity models. A comparison between the results coming from the random diffusivity approach and the ones for standard Brownian motion is discussed. In this way, a deeper physical understanding of the systems described by random diffusivity models is also outlined.
To conclude, a discussion based on the possible origins of the heterogeneity is sketched, with the main goal of inferring which kind of systems can actually be described by the random diffusivity approach.
Humankind and their environment need to be protected from the harmful effects of spent nuclear fuel, and therefore disposal in deep geological formations is favoured worldwide. Suitability of potential host rocks is evaluated, among others, by the retention capacity with respect to radionuclides. Safety assessments are based on the quantification of radionuclide migration lengths with numerical simulations as experiments cannot cover the required temporal (1 Ma) and spatial scales (>100 m).
Aim of the present thesis is to assess the migration of uranium, a geochemically complex radionuclide, in the potential host rock Opalinus Clay. Radionuclide migration in clay formations is governed by diffusion due to their low permeability and retarded by sorption. Both processes highly depend on pore water geochemistry and mineralogy that vary between different facies. Diffusion is quantified with the single-component (SC) approach using one diffusion coefficient for all species and the process-based multi-component (MC) option. With this, each species is assigned its own diffusion coefficient and the interaction with the diffuse double layer is taken into account. Sorption is integrated via a bottom-up approach using mechanistic surface complexation models and cation exchange. Therefore, reactive transport simulations are conducted with the geochemical code PHREEQC to quantify uranium migration, i.e. diffusion and sorption, as a function of mineralogical and geochemical heterogeneities on the host rock scale.
Sorption processes are facies dependent. Migration lengths vary between the Opalinus Clay facies by up to 10 m. Thereby, the geochemistry of the pore water, in particular the partial pressure of carbon dioxide (pCO2), is more decisive for the sorption capacity than the amount of clay minerals. Nevertheless, higher clay mineral quantities compensate geochemical variations. Consequently, sorption processes must be quantified as a function of pore water geochemistry in contact with the mineral assemblage.
Uranium diffusion in the Opalinus Clay is facies independent. Speciation is dominated by aqueous ternary complexes of U(VI) with calcium and carbonate. Differences in the migration lengths between SC and MC diffusion are with +/-5 m negligible. Further, the application of the MC approach highly depends on the quality and availability of the underlying data. Therefore, diffusion processes can be adequately quantified with the SC approach using experimentally determined diffusion coefficients.
The hydrogeological system governs pore water geochemistry within the formation rather than the mineralogy. Diffusive exchange with the adjacent aquifers established geochemical gradients over geological time scales that can enhance migration by up to 25 m. Consequently, uranium sorption processes must be quantified following the identified priority: pCO2 > hydrogeology > mineralogy.
The presented research provides a workflow and orientation for other potential disposal sites with similar pore water geochemistry due to the identified mechanisms and dependencies. With a maximum migration length of 70 m, the retention capacity of the Opalinus Clay with respect to uranium is sufficient to fulfill the German legal minimum requirement of a thickness of at least 100 m.