The western and northern margin off svalbard     

Eirik Sundvor  Olav Eldholm   《Tectonophysics》1979,59(1-4)

Recent geophysical measurements, including multi-channel seismic reflection, on the Svalbard passive margin have revealed that it has undergone a complex geological history which largely reflects the plate tectonic evolution of the Greenland Sea and the Arctic Ocean. The western margin (75–80°N) is of a sheared-rifted type, along which the rifted margin developed subsequent to a change in the pole of plate rotation about 36 m.y. B.P. The north-trending Hornsund Fault on the central shelf and the eastern escarpment of the Knipovich Ridge naturally divide the margin into three structural units. These main marginal structures strike north, paralleling the regional onshore fault trends. This trend also parallels the direction of Early Tertiary plate motion between Svalbard and Greenland. Thus, the western Svalbard margin was initially a zone of shear, and the shear movements have affected the adjacent continental crust. Although, the nature and location of the continent—ocean crustal transition is somewhat uncertain, it is unlikely to lie east of the Hornsund Fault. The northern margin, including the Yermak marginal plateau, is terminated to the west by the Spitsbergen Fracture Zone system. This margin is of a rifted type and the preliminary analysis indicates that the main part of the investigated area is underlain by continental crust.  相似文献   

Monotonicity for MPFA methods on triangular grids     

Eirik Keilegavlen  Ivar Aavatsmark 《Computational Geosciences》2011,15(1):3-16

Flow in a porous medium can be described by a set of non-linear partial differential equations. The pressure variable satisfies a maximum principle, which guarantees that the solution will have no oscillations. A discretisation of the pressure equation should preserve this monotonicity property. Whether a numerical method is monotone will depend both on the medium and on the grid. We study monotonicity of Multi-point Flux Approximation methods on triangular grids. We derive necessary conditions for monotonicity on uniform grids. Further, we study the robustness of the methods on rough grids, and quantify the violations of the maximum principle. These investigations are done for single phase flow, however, they are supported by two-phase simulations.  相似文献   

Heterogeneity preserving upscaling for heat transport in fractured geothermal reservoirs     

Anna Nissen  Eirik Keilegavlen  Tor Harald Sandve  Inga Berre  Jan Martin Nordbotten 《Computational Geosciences》2018,22(2):451-467

In simulation of fluid injection in fractured geothermal reservoirs, the characteristics of the physical processes are severely affected by the local occurence of connected fractures. To resolve these structurally dominated processes, there is a need to develop discretization strategies that also limit computational effort. In this paper, we present an upscaling methodology for geothermal heat transport with fractures represented explicitly in the computational grid. The heat transport is modeled by an advection-conduction equation for the temperature, and solved on a highly irregular coarse grid that preserves the fracture heterogeneity. The upscaling is based on different strategies for the advective term and the conductive term. The coarse scale advective term is constructed from sums of fine scale fluxes, whereas the coarse scale conductive term is constructed based on numerically computed basis functions. The method naturally incorporates the coupling between solution variables in the matrix and in the fractures, respectively, via the discretization. In this way, explicit transfer terms that couple fracture and matrix solution variables are avoided. Numerical results show that the upscaling methodology performs well, in particular for large upscaling ratios, and that it is applicable also to highly complex fracture networks.  相似文献   

Non-hydrostatic pressure in σ-coordinate ocean models     

Eirik Keilegavlen  Jarle Berntsen 《Ocean Modelling》2009,28(4):240-249

Numerical ocean modelling is computationally very demanding. Traditionally, the hydrostatic approximation has been applied to reduce the computational burden. This approximation is valid in large scale studies with coarse grid resolution. With faster computers and gradually smaller grid sizes, we may expect that more studies will be performed with non-hydrostatic ocean models. In recent papers several methods for including non-hydrostatic pressure in σ-coordinate models have been suggested. In this paper the sensitivity of the non-hydrostatic pressure field, the velocity fields, and the density fields to changes in the method for computing non-hydrostatic pressure in σ-coordinate ocean models is addressed.The first test case used involves the propagation and breaking of an internal wave at an incline in a tank. The other test case concerns tidally driven flow over a sill in a stratified fjord. The results from our numerical exercises suggest that the velocity and density fields are very robust to the model choices investigated here. The differences between the model results are of the same order as the uncertainty due to the internal pressure gradient error, and they are smaller than an estimate of the uncertainty due to subgrid scale closure.  相似文献   

Coronal heating and solar wind acceleration; gyrotropic electron-proton solar wind     

Endeve  Eirik  Leer  Egil 《Solar physics》2001,200(1-2):235-250

In coronal holes the electron (proton) density is low, and heating of the proton gas produces a rapidly increasing proton temperature in the inner corona. In models with a reasonable electron density in the upper transition region the proton gas becomes collisionless some 0.2 to 0.3 solar radii into the corona. In the collisionless region the proton heat flux is outwards, along the temperature gradient. The thermal coupling to electrons is weak in coronal holes, so the heat flux into the transition region is too small to supply the energy needed to heat the solar wind plasma to coronal temperatures. Our model studies indicate that in models with proton heating the inward heat conduction may be so inefficient that some of the energy flux must be deposited in the transition region to produce the proton fluxes that are observed in the solar wind. If we allow for coronal electron heating, the energy that is needed in the transition region to heat the solar wind to coronal temperatures, may be supplied by heat conduction from the corona.  相似文献   

Rotational dynamics of subsolar sublimating triaxial comets     

Eirik Mysen 《Planetary and Space Science》2004,52(10):897-907

A model with subsolar water sublimation on a triaxial, ellipsoidal comet nucleus is presented for the calculation of reactive torques. The resulting differential equations describing the comet's rotation are then Hamiltonian, and gravity-gradients are trivial to include. While effects derived from a weak perturbing function are neither able to change the rotational excitation nor the spin magnitude of the nucleus, it is shown how the spin orientation of comets can change significantly over an orbital run. However, of the four comets studied, 1P, 19P, 46P and the Rosetta target 67P, 19P and 46P were the only objects clearly exhibiting this feature, thereby confirming a technique used to derive the consequences of a more elaborate model of sublimation induced torques. In particular, the rotational parameters of 67P were seen to be very stable, indicating that a highly kinematical model of its rotation for the mapping of the comet's gravitational field during the Rosetta mission can be used. The model's hierarchy with 1P/Halley as the object with highest excitation probability, is consistent with observations.  相似文献   

Inversion of inline and broadside marine controlled‐source electromagnetic data with constraints derived from seismic data          下载免费PDF全文

Torgeir Wiik  Janniche Irén Nordskag  Eirik Øverland Dischler  Anh Kiet Nguyen 《Geophysical Prospecting》2015,63(6):1371-1382

We present a structural smoothing regularization scheme in the context of inversion of marine controlled‐source electromagnetic data. The regularizing hypothesis is that the electrical parameters have a structure similar to that of the elastic parameters observed from seismic data. The regularization is split into three steps. First, we ensure that our inversion grid conforms with the geometry derived from seismic. Second, we use a seismic stratigraphic attribute to define a spatially varying regularization strength. Third, we use an indexing strategy on the inversion grid to define smoothing along the seismic geometry. Enforcing such regularization in the inversion will encourage an inversion result that is more intuitive for the interpreter to deal with. However, the interpreter should also be aware of the bias introduced by using seismic data for regularization. We illustrate the method using one synthetic example and one field data example. The results show how the regularization works and that it clearly enforces the structure derived from seismic data. From the field data example we find that the inversion result improves when the structural smoothing regularization is employed. Including the broadside data improves the inversion results even more, due to a better balancing between the sensitivities for the horizontal and vertical resistivities.  相似文献   

U–Th–Pb–REE systematics of organic-rich shales from the ca. 2.15 Ga Sengoma Argillite Formation, Botswana: Evidence for oxidative continental weathering during the Great Oxidation Event     

Gerald D. Pollack  Eirik J. Krogstad  Andrey Bekker   《Chemical Geology》2009,260(3-4):172-185

The redox-sensitive geochemical behavior of uranium permits the use of Th/U ratios as a geochemical proxy for the oxidation state of the atmosphere during deposition. Due to the effects of post-depositional uranium mobility on Th/U ratios during events involving oxygenated fluids, direct measurements of Th/U ratios are often misleading even for drill core samples. Because both of these elements radioactively decay and produce lead isotopes, the Pb isotope composition may reflect the depositional Th/U ratio, although the Th/U ratios induced by changes shortly after deposition may not be distinguished from the true depositional Th/U ratios. In order to effectively evaluate the time-integrated Th/U ratio (κ_a), values for the initial depositional Pb isotope composition must be determined or accepted from the models for the whole Earth.While the timing for the rise of atmospheric oxygen is reasonably well constrained now, its effect on continental weathering and ocean redox state remains poorly constrained and debated. The ca. 2.15 Ga Sengoma Argillite Formation of Botswana contains organic-rich shales deposited during the Great Oxidation Event. The slope of the ²⁰⁷Pb/²⁰⁴Pb–²⁰⁶Pb/²⁰⁴Pb array of shales from the Sengoma Argillite Formation corresponds to a Pb–Pb age that is within analytical error of the depositional age and is, therefore, inferred to be the time by which the time-integrated thorogenic and uranogenic lead growth started. The time-integrated lead growth corresponds to an average κ_a of 2.63 (± 0.62, 1σ) for the organic-rich shales of the Sengoma Argillite Formation. This is lower than Th/U ratios measured in Archean shale suites or estimated for the Archean–Proterozoic average upper continental crust [Taylor, S.R. and McLennan, S.M., 1985. The Continental Crust: Its Composition and Evolution. Blackwell, Oxford, 312 pp.], which indicates that these samples were enriched in uranium with respect to thorium (and perhaps lead) at the time of deposition. In the modern ocean, sediments are enriched in uranium under reducing conditions by reduction of the water-soluble uranyl ion, which is delivered to the ocean by oxidative weathering of continental crust. To evaluate the potential mobility of Th, U, and Pb during post-depositional processes, the concentrations of the rare earth elements (REE) were also determined. Interelement ratios of the largely immobile REE (in this study, La/Nd and Gd/Er) can be used as a proxy for the Th/U ratio, as the geochemical behavior of the lanthanide and actinide elements is similar under a variety of conditions. Furthermore, close similarity in the chondrite-normalized REE patterns and small range in La/Nd and Gd/Er ratios in studied samples indicate that variations in κ_a values are not likely to have been controlled by mixing of one or more REE-, Th-, and U-rich heavy minerals from the multiple detrital sources. Our study of shales from the ca. 2.15 Ga Sengoma Argillite Formation indicates that decoupling of U from Th, most likely related to the oxidative continental weathering, began by 2.15 Ga, at the latest.  相似文献   

Past and future climate variations in the Norwegian Arctic: overview and novel analyses   总被引：1，自引：1，他引：1  

Eirik J. Førland  Inger Hanssen-Bauer 《Polar research》2003,22(2):113-124

Sparse stations and serious measuring problems hamper analyses of climatic conditions in the Arctic. This paper presents a discussion of measuring problems in the Arctic and gives an overview of observed past and projected future climate variations in Svalbard and Jan Mayen. Novel analyses of temperature conditions during precipitation and trends in fractions of solid/liquid precipitation at the Arctic weather stations are also outlined. Analyses based on combined and homogenized series from the regular weather stations in the region indicate that the measured annual precipitation has increased by more than 2.5% per decade since the measurements started in the beginning of the 20th century. The annual temperature has increased in Svalbard and Jan Mayen during the latest decades, but the present level is still lower than in the 1930s. Downscaled scenarios for Svalbard Airport indicate a further increase in temperature and precipitation. Analyses based on observations of precipitation types at the regular weather stations demonstrate that the annual fraction of solid precipitation has decreased at all stations during the latest decades. The reduced fraction of solid precipitation implies that the undercatch of the precipitation gauges is reduced. Consequently, part of the observed increase in the annual precipitation is fictitious and is due to a larger part of the "true" precipitation being caught by the gauges. With continued warming in the region, this virtual increase will be measured in addition to an eventual real increase.  相似文献   

Geological events during the early formation of a passive margin     

Olav Eldholm  Eirik Sundvor   《Tectonophysics》1979,59(1-4)

The concept of plate tectonics implies that the normal sea floor spreading stage is preceded by a sequence of events associated with the break-up of continental crust. Thus, evidence of the early development of “non-failed” rifts is to be found at passive continental margins. Of special interest is the question of the extent of the continental crust and the structural and compositional changes associated with the change in crustal type. In addressing these topics, we have focused attention on the Norwegian margin between the Jan Mayen and Senja fracture zones (66°–70°N) in an attempt to understand its history of rifting and early sea floor spreading. p ]The southern part of this rifted margin is characterized by a wide shelf and the marginal Vøring Plateau interrupts a gentle slope at a level of about 1500 m. However, the margin becomes progressively narrower towards the north and a typical narrow shelf and steep slope emerge off the Lofo—tenVesterålen Islands (Fig. 1). In a reconstructed pre-opening configuration (Talwani and Eldholm, 1977) the narrowest part of the juxtaposed EastGreenland margin is found in the south and a wide shelf and slope corresponds to the Lofoten-Vesterålen margin.The most prominent structural element is a buried basement high underneath the Vøring Plateau. The high is bounded landward by the Vøring Plateau Escarpment, a major structural boundary which defines typical changes in the geophysical parameters. These are: (1) a sudden increase of depth to acoustic basement; (2) changes in the velocity-depth function; (3) a gravity gradient; and (4) a magnetic edge anomaly separating sea-floor spreading type anomalies from a quiet zone on the landward side (Talwani and Eldholm, 1972). These observations were interpreted in terms of a sharp ocea—ncontinent crustal transition along the escarpment with sea-floor spreading commencing between anomaly 24 and 25 time (56–58 m.y. B.P.). Alternatively, the concept of ancient oceanic crust landward of this escarpment and the possible existence of continental crust under the outer basement high have been argued and we refer to Eldholm et al. (1979) for a detailed discussion.  相似文献