High-resolution IP<Subscript>25</Subscript>-based reconstruction of sea-ice variability in the western North Pacific and Bering Sea during the past 18,000 years

Due to its strong influence on heat and moisture exchange between the ocean and the atmosphere, sea ice is an essential component of the global climate system. In the context of its alarming decrease in terms of concentration, thickness and duration, understanding the processes controlling sea-ice variability and reconstructing paleo-sea-ice extent in polar regions have become of great interest for the scientific community. In this study, for the first time, IP₂₅, a recently developed biomarker sea-ice proxy, was used for a high-resolution reconstruction of the sea-ice extent and its variability in the western North Pacific and western Bering Sea during the past 18,000 years. To identify mechanisms controlling the sea-ice variability, IP₂₅ data were associated with published sea-surface temperature as well as diatom and biogenic opal data. The results indicate that a seasonal sea-ice cover existed during cold periods (Heinrich Stadial 1 and Younger Dryas), whereas during warmer intervals (Bølling-Allerød and Holocene) reduced sea ice or ice-free conditions prevailed in the study area. The variability in sea-ice extent seems to be linked to climate anomalies and sea-level changes controlling the oceanographic circulation between the subarctic Pacific and the Bering Sea, especially the Alaskan Stream injection though the Aleutian passes.     

Postglacial uplift and relative sea level changes in Finnmark,northern Norway

The outer coast of Finnmark in northern Norway is where the former Fennoscandian and Barents Sea ice sheets coalesced. This key area for isostatic modelling and deglaciation history of the ice sheets has abundant raised shorelines, but only a few existing radiocarbon dates constrain their chronology. Here we present three Holocene sea level curves based on radiocarbon dated deposits from isolation basins at the outermost coast of Finnmark; located at the islands Sørøya and Rolvsøya and at the Nordkinn peninsula. We analysed animal and plant remains in the basin deposits to identify the transitions between marine and lacustrine sediments. Terrestrial plant fragments from these transitions were then radiocarbon dated. Radiocarbon dated mollusk shells and marine macroalgae from the lowermost deposits in several basins suggest that the first land at the outer coast became ice free around 14,600 cal yr BP. We find that the gradients of the shorelines are much lower than elsewhere along the Norwegian coast because of substantial uplift of the Barents Sea. Also, the anomalously high elevation of the marine limit in the region can be attributed to uplift of the adjacent seafloor. After the Younger Dryas the coast emerged 1.6–1.0 cm per year until about 9500–9000 cal yr BP. Between 9000 and 7000 cal yr BP relative sea level rose 2–4 m and several of the studied lakes became submerged. At the outermost locality Rolvsøya, relative sea level was stable at the transgression highstand for more than 3000 years, between ca 8000 and 5000 cal yr BP. Deposits in five of the studied lakes were disturbed by the Storegga tsunami ca 8200–8100 cal yr BP.     

Three‐dimensional modeling of problems in poro‐elasticity via a mixed least‐squares method using linear tetrahedral elements

In a previous publication we developed a new mixed least‐squares method for poro‐elasticity. The approximate solution was obtained via a minimization of a least‐squares functional, based upon the equations of equilibrium, the equations of continuity and weak forms of the constitutive relationships for elasticity and Darcy flow. The formulation involved four independent types of variables: displacements, stresses, pore pressures and velocities. All of them were approximated by linear continuous triangles. Encouraged by the computational results, obtained from the two‐dimensional implementation of the method, we extended our formulation to three dimensions. In this paper we present numerical examples for the performance of continuous linear tetrahedra within the context of the mixed least‐squares method. The initial results suggest that the method works well in the nearly and entirely incompressible limits for elasticity. For poro‐elasticity, the obtained pore pressures are stable without exhibiting the oscillations, which are observed when the standard Galerkin formulation is used. Copyright © 2010 John Wiley & Sons, Ltd.     

Nonlinear output constraints handling for production optimization of oil reservoirs

Adjoint-based gradient computations for oil reservoirs have been increasingly used in closed-loop reservoir management optimizations. Most constraints in the optimizations are for the control input, which may either be bound constraints or equality constraints. This paper addresses output constraints for both state and control variables. We propose to use a (interior) barrier function approach, where the output constraints are added as a barrier term to the objective function. As we assume there always exist feasible initial control inputs, the method maintains the feasibility of the constraints. Three case examples are presented. The results show that the proposed method is able to preserve the computational efficiency of the adjoint methods.     

Paired uraninite and molybdenite dating of the K?nigshain granite: implications for the onset of late-Variscan magmatism in the Lausitz Block

We present geochronological data for late-Variscan magmatism in the Lausitz Block of the Saxo-Thuringian Zone, Germany. The Th–U–total Pb age of uraninite and the Re–Os age of molybdenite from the composite biotite–monzogranite pluton of K?nigshain overlap at the 2σ confidence limit: 328.6 ± 1.9 Ma (uraninite), and 327.0 ± 1.3 Ma and 327.6 ± 1.3 Ma (molybdenite), indicating that crystallization of magmatic uraninite and deposition of molybdenite were nearly contemporaneous. These data imply that magmatic processes in this part of the Variscan orogen already started in latest Visean time, about 10 Ma earlier than previously assumed (315–320 Ma). The new ages correspond to ages for plutonic rocks in the Elbe Zone immediately west of the Lausitz (around 335–325 Ma) and the bulk of late-Variscan igneous rocks in the Saxo-Thuringian Zone (335–320 Ma).     

Grid adaptation for the Dirichlet–Neumann representation method and the multiscale mixed finite-element method

A Dirichlet–Neumann representation method was recently proposed for upscaling and simulating flow in reservoirs. The DNR method expresses coarse fluxes as linear functions of multiple pressure values along the boundary and at the center of each coarse block. The number of flux and pressure values at the boundary can be adjusted to improve the accuracy of simulation results and, in particular, to resolve important fine-scale details. Improvement over existing approaches is substantial especially for reservoirs that contain high-permeability streaks or channels. As an alternative, the multiscale mixed finite-element (MsMFE) method was designed to obtain fine-scale fluxes at the cost of solving a coarsened problem, but can also be used as upscaling methods that are flexible with respect to geometry and topology of the coarsened grid. Both methods can be expressed in mixed-hybrid form, with local stiffness matrices obtained as “inner products” of numerically computed basis functions with fine-scale sub-resolution. These basis functions are determined by solving local flow problems with piecewise linear Dirichlet boundary conditions for the DNR method and piecewise constant Neumann conditions for MsMFE. Adding discrete pressure points in the DNR method corresponds to subdividing faces in the coarse grid and hence increasing the number of basis functions in the MsMFE method. The methods show similar accuracy for 2D Cartesian cases, but the MsMFE method is more straightforward to formulate in 3D and implement for general grids.     

Excitation of solar-like oscillations: From PMS to MS stellar models

The amplitude of solar-like oscillations results from a balance between excitation and damping. As in the sun, the excitation is attributed to turbulent motions that stochastically excite thep modes in the uppermost part of the convective zone. We present here a model for the excitation mechanism. Comparisons between modeled amplitudes and helio and stellar seismic constraints are presented and the discrepancies discussed. Finally the possibility and the interest of detecting such stochastically excited modes in pre-main sequence stars are also discussed.     

Accumulation of particulate organic carbon at the Eurasian continental margin during late Quaternary times: controlling mechanisms and paleoenvironmental significance

Data on the amount and composition of organic carbon were determined in sediment cores from the Kara and Laptev Sea continental margin, representing oxygen isotope stages 1–6. The characterization of organic matter is based on hydrogen index (HI) values, n-alkanes and maceral composition, indicating the predominance of terrigenous organic matter through space and time. The variations in the amount and composition of organic carbon are mainly influenced by changes in fluvial sediment supply, Atlantic water inflow, and continental ice sheets. During oxygen isotope stage (OIS) 6, high organic carbon contents in sediments from the Laptev Sea and western East Siberian Sea continental margin were probably caused by the increased glacial erosion and further transport in the eastward-flowing boundary current along the continental margin. During OIS 5 and early OIS 3, some increased amounts of marine organic matter were preserved in sediments east of the Lomonosov Ridge, suggesting an influence of nutrient-rich Pacific waters. During OIS 2, terrigenous organic carbon supply was increased along the Barents and western Kara Sea continental margin caused by extended continental ice sheets in the Barents Sea (Svalbard to Franz Josef Land) area and increased glacial erosion. Along the Laptev Sea continental margin, on the other hand, the supply of terrigenous (organic) matter was significantly reduced due to the lack of major ice sheets and reduced river discharge. Towards the Holocene, the amount of total organic carbon (TOC) increased along the Kara and Laptev Sea continental margin, reaching average values of up to 0.5 g C cm⁻² ky⁻¹. Between about 8 and 10 ka (9 and 11 Cal ka), i.e., during times when the inner shallow Kara and Laptev seas became largely flooded for the first time after the Last Glacial Maximum, maximum supply of terrigenous organic carbon occurred, which is related to an increase in coastal erosion and Siberian river discharge. During the last 8000 years, the increased amount of marine organic carbon preserved in the sediments from the Kara and Laptev Sea continental margin is interpreted as a result of the intensification of Atlantic water inflow along the Eurasian continental margin.     

Extinct isotope heterogeneities in the mantles of Earth and Mars: Implications for mantle stirring rates

Heterogeneities in terrestrial samples for ¹⁸²W/¹⁸³W and ¹⁴²Nd/¹⁴⁴Nd are only preserved in Hadean and Archean rocks while heterogeneities in ¹²⁹Xe/¹³⁰Xe and ¹³⁶Xe/¹³⁰Xe persist to very young mantle‐derived rocks. In contrast, meteorites from Mars show that the Martian mantle preserves heterogeneities in ¹⁸²W/¹⁸³W and ¹⁴²Nd/¹⁴⁴Nd up to the present. As a consequence of the probable “deep magma ocean” core formation process, we assume that the Earth and Mars both had a very early two‐mantle‐reservoir structure with different initial extinct nuclide isotopic compositions (different ¹⁸²W/¹⁸³W, ¹⁴²Nd/¹⁴⁴Nd, ¹²⁹Xe/¹³⁰Xe, ¹³⁶Xe/¹³⁰Xe ratios). Based on this assumption, we developed a simple stochastic model to trace the evolution of a mantle with two initially distinct layers for the extinct isotopes and its development into a heterogeneous mantle by convective mixing and stretching of these two layers. Using the extinct isotope system ¹⁸²Hf‐¹⁸²W, we find that the mantles of Earth and Mars exhibit substantially different mixing or stirring rates. This is consistent with Mars having cooled faster than the Earth due to its smaller size, resulting in less efficient mantle mixing for Mars. Moreover, the mantle stirring rate obtained for Earth using ¹⁸²Hf‐¹⁸²W is consistent with the mantle stirring rate of ~500 Myr constrained by the long‐lived isotope system, ⁸⁷Rb‐⁸⁷Sr and ¹⁴⁷Sm‐¹⁴³Nd. The apparent absence of ¹⁸²W/¹⁸³W isotopic heterogeneity in modern terrestrial rocks is attributed to very active mantle stirring which reduced the ¹⁸²W/¹⁸³W isotopic heterogeneity to a relatively small scale (~83 m for a mantle stirring rate of 500 Myr) compared to the common sampling scale of terrestrial basalts (~30 or 100 km). Our results also support the “deep magma ocean” core formation model as being applicable to both Mars and Earth.     

Isotopic climate record in a Holocene stalagmite from Ursilor Cave (Romania)

The PU‐2 stalagmite from Ursilor Cave provides the first dated Romanian isotope record for the Holocene. The overall growth rate of the speleothem was 3.5 cm kyr^?1, corresponding to a temporal resolution of 142 y between each isotope analysis. The ‘Hendy’ tests indicate that isotopic equilibrium conditions occurred during the formation of PU‐2, and hence that it is suitable for palaeoclimatic studies. The relationship between δ¹⁸O and temperature was found to be positive. This can be interpreted either as rain‐out with distance from the west‐northwest ocean source of evaporation or shifts in air mass source with changing North Atlantic Oscillation indices. Applying five U–Th thermal ionisation mass spectrometric (TIMS) dates to a 17.5 cm isotope profile (δ¹⁸O and δ¹³C) along the stalagmite growth axis enabled a tentative interpretation of the palaeoclimate signal over the past 7.1 kyr. Spikes of depleted isotopic δ¹⁸O values are centred near ca. 7, ca. 5.2 and ca. 4 ka, reflecting cool conditions. The record shows two warm intervals between ca. 3.8 and ca. 3.2 ka (the maximum warmth) and from ca. 2 to ca. 1.4 ka, when the δ¹⁸O values were less negative than present. The ‘Holocene Climate Optimum’ spanning the time interval from ca. 6.8 to ca. 4.4 ka is not well expressed in the PU‐2 stalagmite. Individual spikes of lighter δ¹³C are interpreted as indicative of periods of heavy rainfall, at ca. 7, ca. 5.5, and ca. 3.5 ka. The overall trend to lighter δ¹³C in the PU‐2 stalagmite may reflect a gradual decrease in water–rock interaction. The results demonstrate that the effect of North Atlantic oceanic changes extended to the investigated area. Nevertheless, some differences in temporal correlation and intensity of stable isotopic response to these climatic events have been found, but the exact nature of these differences and the underlying mechanism is yet to be determined. Copyright © 2002 John Wiley & Sons, Ltd.