ALBERT

Description: We present a plate kinematic evolution of the South Atlantic which is based largely on the determination of the equatorial fracture zone trends between the African and South American continental margins. Four main opening phases are dated by oceanic magnetic anomalies, notably MO, A34, and A13, and are correlated with volcanism and tectonic events on land around the South Atlantic Ocean. The Ceara and Sierra Leone rises are probably of oceanic origin and were created 80 m.y. ago or later in their present-day positions with respect to South America and Africa.

Description: This paper concerns the linear response of the ocean to forcing at a specified frequency and wave number in the absence of mean currents. It discusses the details of the forcing function, the general properties of the equations of motion, and possible simplifications of these equations. Two representations for the oceanic response to forcing are described in detail. One solution is in terms of the normal modes of the ocean. The vertical structure of these modes corresponds to that of the barotropic and baroclinic modes; their latitudinal structure corresponds to that of inertia‐gravity and Rossby waves. These waves are eigenfunctions of Laplace's tidal equations (LTE) with the frequency as eigenvalue. The description in terms of vertically standing modes is particularly useful if the forcing is nonlocal, because only these modes can propagate into undisturbed regions. The principal result is that it is extremely difficult for baroclinic (but not barotropic) disturbances to propagate horizontally away from a forced region. Instabilities of the Gulf Stream excite disturbances that are confined to the immediate neighborhood of the current; disturbances due to instabilities of equatorial currents do not propagate far latitudinally. A second representation of the oceanic response to forcing is in terms of vertically propagating, or vertically trapped, latitudinal modes. These modes are eigenfunctions of LTE with the equivalent depth h (not the frequency) as eigenvalue. Both positive and negative eigenvalues h are necessary for completeness. The modes with h 〉 0 consist of an infinite set of inertia‐gravity waves and a finite set of Rossby waves which either propagate vertically or form vertically standing modes. The latitudinally gravest modes are equatorially trapped and have been observed in the Atlantic and Pacific oceans. The modes with h 〈 0 are necessary to describe the oceanic response to nonresonant forcing. In the vertical this response attenuates with increasing distance from the forcing region. Because of the shallowness of the ocean the large eastward traveling atmospheric cyclones in mid‐latitudes and high latitudes force a response down to the ocean floor. Interaction with the bottom topography will result in smaller‐scale disturbances and will affect the frequency spectrum of the response when bottom‐trapped waves are excited.

Description: A tsunami earthquake is defined as a shock which generates extensive tsunamis but relatively weak seismic waves. A comparative study is made for the two recent tsunami earthquakes, and a subduction mechanism near a deep-sea trench is discussed. These two earthquakes occurred at extremely shallow depths far off the coasts of the Kurile Islands and of eastern Hokkaido on October 20, 1963, and on June 10, 1975, respectively. Both can be regarded as an aftershock of the preceding larger events. Their tsunami heights and seismic wave amplitudes are compared with those of the preceding events. The results show that the time constants involved in the tsunami earthquakes are relatively long but not long enough to explain the observed disproportionality between the tsunamis and the seismic waves. The process times are estimated to be less than 100 s. The spatio-temporal characteristics of the two events suggest that they represent a seaward and upward extension of the rupture associated with a great earthquake which did not break the free surface at the coseismic stage. The amplitude and phase spectra of long-period surface waves and the long-period P waveforms indicate that this extension of the rupture did not take place entirely along the lithospheric interface emerging as a trench axis. It rather branched upward from the interface in a complex way through the wedge portion at the leading edge of the continental lithosphere. This wedge portion consists in large part of thick deformable sediments. A large vertical deformation and hence extensive tsunamis result from such a branching process. A shallowest source depth, steepening of rupture surfaces, and a deformable nature of the source region all enhance generation of tsunamis. The wedge portion ruptured by a tsunami earthquake is usually characterized by a very low seismic activity which is presumably due to ductility of the sediments. We suggest that this portion fractures in a brittle way to generate a tsunami earthquake when it is loaded suddenly by the occurrence of a great earthquake and that otherwise it yields slowly. Upward branching of the rupture from the lithospheric interface produces permanent deformation of the free surface which is relative uplift landward and relative subsidence trenchward of the zone of surface break. This surface break zone geomorphologically corresponds to the lower continental slope between the deep-sea terrace and the trench. Such a mode of permanent deformation seems to be consistent with a rising feature of the outer ridge of the deep-sea terrace and a depressional feature of the trench. This consistency implies a causal relationship between great earthquake activities and geomorphological features near the trench.

Description: A tsunami earthquake is defined as a shock which generates extensive tsunamis but relatively weak seismic waves. A comparative study is made for the two recent tsunami earthquakes, and a subduction mechanism near a deep-sea trench is discussed. These two earthquakes occurred at extremely shallow depths far off the coasts of the Kurile Islands and of eastern Hokkaido on October 20, 1963, and on June 10, 1975, respectively. Both can be regarded as an aftershock of the preceding larger events. Their tsunami heights and seismic wave amplitudes are compared with those of the preceding events. The results show that the time constants involved in the tsunami earthquakes are relatively long but not long enough to explain the observed disproportionality between the tsunamis and the seismic waves. The process times are estimated to be less than 100 s. The spatio-temporal characteristics of the two events suggest that they represent a seaward and upward extension of the rupture associated with a great earthquake which did not break the free surface at the coseismic stage. The amplitude and phase spectra of long-period surface waves and the long-period P waveforms indicate that this extension of the rupture did not take place entirely along the lithospheric interface emerging as a trench axis. It rather branched upward from the interface in a complex way through the wedge portion at the leading edge of the continental lithosphere. This wedge portion consists in large part of thick deformable sediments. A large vertical deformation and hence extensive tsunamis result from such a branching process. A shallowest source depth, steepening of rupture surfaces, and a deformable nature of the source region all enhance generation of tsunamis. The wedge portion ruptured by a tsunami earthquake is usually characterized by a very low seismic activity which is presumably due to ductility of the sediments. We suggest that this portion fractures in a brittle way to generate a tsunami earthquake when it is loaded suddenly by the occurrence of a great earthquake and that otherwise it yields slowly. Upward branching of the rupture from the lithospheric interface produces permanent deformation of the free surface which is relative uplift landward and relative subsidence trenchward of the zone of surface break. This surface break zone geomorphologically corresponds to the lower continental slope between the deep-sea terrace and the trench. Such a mode of permanent deformation seems to be consistent with a rising feature of the outer ridge of the deep-sea terrace and a depressional feature of the trench. This consistency implies a causal relationship between great earthquake activities and geomorphological features near the trench.

Description: Published geological and geophysical data are reviewed. The Walvis Ridge is a complex linear feature made up of three parts of unequal lengths and differing basement morphologies: an eastern sector composed of rugged, subparallel basement ridges; a low-lying central sector with subdued basement morphology; and a western sector consisting of seamounts and guyots (including Tristan da Cunha and Gough islands). Rock samples and geophysical data suggest that the Ridge is composed of alkali basalt which becomes progressively older eastwards. Gravity data indicate that at least parts of the ridge are in local isostatic equilibrium. A mantle plume mechanism of formation is rejected in favour of a centre of abnormally high volcanic activity on the spreading ridge axis. The location of this centre, whose relative movement has been southwards, is determined by fracture zones crossing the spreading ridge axis.

Description: Observations of the temporal and spatial distribution of N2O in solution are not yet sufficient to permit quantitative assessment of the role of the ocean in the budget of atmospheric N2O. Consideration of the global nitrogen cycle suggests that the land should be the primary source of N2O. The gas is removed in the atmosphere by photolysis and by reaction with O(1D), and there may be additional sinks in the ocean.

Description: One hundred and five new heat flow measurements in the Gulf of California support the premise that conductive heat loss is not the only mode by which heat is lost from a sea floor spreading center, even in an area with thick sediment cover. Theoretical estimates suggest that the average heat flow in the Guaymas and Farallon basins should be at least 11 μcal/cm2 s (HFU) (325 mW/m2). Outside a 30-km-wide zone centered on the central troughs, the heat flow values measured are reasonably uniform but average only 4.3±0.2 HFU (180±10 mW/m2). Although the high sedimentation rate may depress the measured heat flow, the effect probably does not exceed 15%. Some heat, particularly in the smaller basins, may be lost to the adjacent cooler continental blocks. The discrepancy between the measured and predicted heat losses, which is at least 30%, may be due to the discharge of thermal waters, through the thinner sediment cover in the central troughs or along active faults.

Description: The Cape Verde Islands are emerged portions of a Mesozoic-Cenozoic volcanic accretion in the form of a westward-opening horseshoe along fracture zones converging from the mid-Atlantic ridge toward Africa. An interior abyssal plain slopes westward, increasing in depth from 2.7 to 4.5 km. The plain is underlain by low relief on acoustic basement that is associated with a 300-gamma negative magnetic anomaly. The flanks of the Sal-Maio ridge appear bounded by large-displacement normal faults; superficial slumping is common. The trends of magnetic anomalies are linear N-S north of the islands and less linear within the islands and may change coincident with E-W bathymetric trends south of the islands. A triangular pattern of reversed refraction lines 200–250 km long along the north and east ridges and NW-SE across the interior abyssal plain indicated 2–3 km of semiconsolidated sediments underlain by 3–6 km of basalt and 6–8 km of plutonic rocks. The depth of the Moho is between 16 and 17 km. A deep NW-SE trending fault intersects the Sal-Maio ridge near Boa Vista. The consistent depth to Moho and the regional Bouguer anomaly indicate lack of local relief at the base of the crust. The crustal load of the entire archipelago is regionally adjusted.

Description: An intensive three-dimensional survey of the Antarctic Polar Front was made in the Drake Passage in March 1976. The front, which was imbedded within one of the high-velocity cores of the circumpolar current, is viewed as a water mass boundary demarking the northern extent of near-surface antarctic waters. Within the front, water masses are observed to intrude, one above the other, with characteristic vertical scales of 50–100 m. The intrusions are horizontally anisotropic, being elongated in the alongstream direction and constrained primarily to the upper 800 m of the front. The spatial and temporal persistence of the variability is examined through the analysis of continuous vertical profiles of horizontal velocity, temperature, salinity, and oxygen with discrete sampling of nutrients. Analysis of the velocity data showed the mean current flowing to the NNE with speeds of the order of 30–40 cm s−1 in the upper 600 m, with temporal variability over a 28-hour ‘yo-yo’ due primarily to internal gravity waves. The thermohaline variability was not internal wave induced but rather was associated with nearly isentropic advection of different water masses across the front. Cold fresh and warm salty intrusions did not conserve potential density, however, and double-diffusive transfers are strongly suggested as being crucial to an understanding of the dynamics of the intrusions. Applying a model (Joyce, 1977) for lateral mixing we estimate poleward temperature and salinity fluxes due to interleaving of 0.086°C cm s−1 and 0.069‰ cm s−1, respectively. If these values are typical, interleaving could play a significant role in large-scale balance of salt and, to a lesser extent, heat for the Southern Ocean.

Description: The structural evolution of the northwestern Iberian margin has been reconstructed from the results of IPOD drill site 398, as well as from numerous dredgings and a dense network of seismic profiles. During the Mesozoíc the margin first underwent two consecutive extensional phases interpreted as the result of two episodes of rifting in the Atlantic. Then during Eocene, subsidence was interrupted by compression and related deformation caused by subduction of oceanic sea floor of the Bay of Biscay beneath the Iberian Peninsula. Present day marginal banks are interpreted as blocks of the older passive margin uplifted during early Tertiary as a result of that subduction. Fault escarpments provide opportunities to sample older sediments and basement by dredging.

Description: A model for interstitial silica concentrations is derived, incorporating biological mixing of sediments. This model predicts concentrations and gradients and can account for the observed geographical variations in interstitial silica on the basis of a dynamic balance between solution of silica particles and diffusion from the sediments. The flux of particulate biogenous silica into the sediments is confirmed as an important parameter controlling interstitial silica concentrations. Biological mixing of sea floor sediments also has an important influence on interstitial composition by modifyirig the depth at which dissolving particles react. Faster mixing raises the interstitial concentration. The rate at which siliceous particles dissolve also plays a role; the slower they dissolve, the greater the interstitial silica concentration. Measurements on near‐bottom waters of the Atlantic show no consistent gradients in dissolved silica, but antarctic bottom water seems significantly more variable in the benthic boundary layer than in the water mass above or in the benthic zone of North Atlantic deep water.

Description: Biological mixing in deep‐sea sediments is described in terms of a time‐dependent eddy diffusion model where mixing takes place to a depth L at constant eddy diffusivity D. The differential equation that describes this model has been solved for an impulse source of tracer delivered to the plane surface that forms the top of the mixed layer. The solution then serves as a Green's function, which can be used to determine the distribution of tracer in depth and in time for a surface input of tracer specified as any arbitrary function of time. The characteristic properties of the solution are dependent on the dimensionless parameter D/Lυ, where υ is the sedimentation rate. If D/Lυ is greater than 10, the surface layer becomes homogeneous, and the model is identical to the homogeneous layer model proposed by Berger and Heath (1968). If D/Lυ is less than 0.1, little mixing can take place before the sediments are buried, and so the surface concentration propagates downward into the sediments with little dispersion. For all values of D/Lυ the weighted mean depth of the concentration distribution is the depth at which an impulse source would be found in the sediment if no mixing had taken place. The microtektite data of Glass (1969, 1972) and Glass et al. (1973) indicate that abyssal sediments are mixed from the surface to a maximum mixing depth that ranges between 17 and 40 cm below the surface. Mixing occurs at rates between 1 and 100 cm2 kyr−1. Higher mixing rates may occur nearer the surface, but microtektite distributions cannot be used to estimate these rates in the presence of the deeper, slower mixing. Estimates for D based on dimensional analysis of sediment reworking rates for nearshore organisms (103–106 cm2 kyr−1) are used to predict abyssal mixing rates between 1 and 103 cm2 kyr−1 by invoking the assumption that mixing is proportional to biomass. Plutonium distributions in deep‐sea sediments (Noshkin and Bowen, 1973) indicate abyssal mixing rates ranging from 100 to 400 cm2 kyr−1.

Description: Deep‐sea drilling in the Antarctic region (Deep‐Sea Drilling Project legs 28, 29, 35, and 36) has provided many new data about the development of circum‐Antarctic circulation and the closely related glacial evolution of Antarctica. The Antarctic continent has been in a high‐latitude position since the middle to late Mesozoic. Glaciation commenced much later, in the middle Tertiary, demonstrating that near‐polar position is not sufficient for glacial development. Instead, continental glaciation developed as the present‐day Southern Ocean circulation system became established when obstructing land masses moved aside. During the Paleocene (t = ∼65 to 55 m.y. ago), Australia and Antarctica were joined. In the early Eocene (t = ∼55 m.y. ago), Australia began to drift northward from Antarctica, forming an ocean, although circum‐Antarctic flow was blocked by the continental South Tasman Rise and Tasmania. During the Eocene (t = 55 to 38 m.y. ago) the Southern Ocean was relatively warm and the continent largely nonglaciated. Cool temperate vegetation existed in some regions. By the late Eocene (t = ∼39 m.y. ago) a shallow water connection had developed between the southern Indian and Pacific oceans over the South Tasman Rise. The first major climatic‐glacial threshold was crossed 38 m.y. ago near the Eocene‐Oligocene boundary, when substantial Antarctic sea ice began to form. This resulted in a rapid temperature drop in bottom waters of about 5°C and a major crisis in deep‐sea faunas. Thermohaline oceanic circulation was initiated at this time much like that of the present day. The resulting change in climatic regime increased bottom water activity over wide areas of the deep ocean basins, creating much sediment erosion, especially in western parts of oceans. A major (∼2000 m) and apparently rapid deepening also occurred in the calcium carbonate compensation depth (CCD). This climatic threshold was crossed as a result of the gradual isolation of Antarctica from Australia and perhaps the opening of the Drake Passage. During the Oligocene (t = 38 to 22 m.y. ago), widespread glaciation probably occurred throughout Antarctica, although no ice cap existed. By the middle to late Oligocene (t = ∼30 to 25 m.y. ago), deep‐seated circum‐Antarctic flow had developed south of the South Tasman Rise, as this had separated sufficiently from Victoria Land, Antarctica. Major reorganization resulted in southern hemisphere deep‐sea sediment distribution patterns. The next principal climatic threshold was crossed during the middle Miocene (t = 14 to 11 m.y. ago) when the Antarctic ice cap formed. This occurred at about the time of closure of the Australian‐Indonesian deep‐sea passage. During the early Miocene, calcareous biogenic sediments began to be displaced northward by siliceous biogenic sediments with higher rates of sedimentation reflecting the beginning of circulation related to the development of the Antarctic Convergence. Since the middle Miocene the East Antarctic ice cap has remained a semipermanent feature exhibiting some changes in volume. The most important of these occurred during the latest Miocene (t = ∼5 m.y. ago) when ice volumes increased beyond those of the present day. This event was related to global climatic cooling, a rapid northward movement of about 300 km of the Antarctic Convergence, and a eustatic sea level drop that may have been partly responsible for the isolation of the Mediterranean basin. Northern hemisphere ice sheet development began about 2.5–3 m.y. ago, representing the next major global climatic threshold, and was followed by the well‐known major oscillations in northern ice sheets. In the Southern Ocean the Quaternary marks a peak in activity of oceanic circulation as reflected by widespread deep‐sea erosion, very high biogenic productivity at the Antarctic Convergence and resulting high rates of biogenic sedimentation, and maximum northward distribution of ice‐rafted debris.

Description: The olivine melilitite diatemes of the Swabian Alb, frequently compared with kimberlite diatremes, are discussed in terms of hydrogeological setting, internal structure and juvenile fraction. The hydrogeological conditions of the Swabian Alb at the time of diatreme emplacement were characterized by copious amounts of groundwater within the sedimentary cover of the basement. Subsequently to the eruptions groundwater accumulated within the maars of the larger diatremes forming fresh‐water lakes as also happened nearby in the Steinheim and Ries impact craters. The diatremes reveal subsidence structures composed of large wall‐rock blocks, subaerially deposited pyroclastic beds, and well‐bedded reworked pyroclastic debris which accumulated on the floor of the fresh‐water crater lakes. The latter fact implies availability of groundwater at the time the diatremes formed. The juvenile fraction is developed in the shape of spherical to ovoid nucleated autoliths of ash to lapilli size that are macroscopically nearly devoid of vesicles. The autoliths are interpreted as the product of water vapor explosions which took place when rising olivine melilitite magma contacted groundwater and was fragmented into magma droplets. The droplets were rapidly chilled and thus preserved their shape. Because of the hydrogeological data, the diatreme structure, and the chilled nature of the autoliths a phreatomagmatic origin of the Swabian diatremes is suggested.

Description: The olivine melilitite diatemes of the Swabian Alb, frequently compared with kimberlite diatremes, are discussed in terms of hydrogeological setting, internal structure and juvenile fraction. The hydrogeological conditions of the Swabian Alb at the time of diatreme emplacement were characterized by copious amounts of groundwater within the sedimentary cover of the basement. Subsequently to the eruptions groundwater accumulated within the maars of the larger diatremes forming fresh‐water lakes as also happened nearby in the Steinheim and Ries impact craters. The diatremes reveal subsidence structures composed of large wall‐rock blocks, subaerially deposited pyroclastic beds, and well‐bedded reworked pyroclastic debris which accumulated on the floor of the fresh‐water crater lakes. The latter fact implies availability of groundwater at the time the diatremes formed. The juvenile fraction is developed in the shape of spherical to ovoid nucleated autoliths of ash to lapilli size that are macroscopically nearly devoid of vesicles. The autoliths are interpreted as the product of water vapor explosions which took place when rising olivine melilitite magma contacted groundwater and was fragmented into magma droplets. The droplets were rapidly chilled and thus preserved their shape. Because of the hydrogeological data, the diatreme structure, and the chilled nature of the autoliths a phreatomagmatic origin of the Swabian diatremes is suggested.

Description: The current geochronological state of the art for applying the radiocarbon (14C) method to deep-sea sediment archives lacks key information on sediment bioturbation. Here, we apply a sediment accumulation model that simulates the sedimentation and bioturbation of millions of foraminifera, whereby realistic 14C activities (i.e. from a 14C calibration curve) are assigned to each single foraminifera based on its simulation time step. We find that the normal distribution of 14C age typically used to represent discrete-depth sediment intervals (based on the reported laboratory 14C age and measurement error) is unlikely to be a faithful reflection of the actual 14C age distribution for a specific depth interval. We also find that this deviation from the actual 14C age distribution is greatly amplified during the calibration process. Specifically, we find a systematic underestimation of total geochronological error in many cases (by up to thousands of years), as well as the generation of age–depth artefacts in downcore calibrated median age. Even in the case of “perfect” simulated sediment archive scenarios, whereby sediment accumulation rate (SAR), bioturbation depth, reservoir age and species abundance are all kept constant, the 14C measurement and calibration processes generate temporally dynamic median age–depth artefacts on the order of hundreds of years – whereby even high SAR scenarios (40 and 60 cm kyr−1) are susceptible. Such age–depth artefacts can be especially pronounced during periods corresponding to dynamic changes in the Earth's Δ14C history, when single foraminifera of varying 14C activity can be incorporated into single discrete-depth sediment intervals. For certain lower-SAR scenarios, we find that downcore discrete-depth true median age can systematically fall outside the calibrated age range predicted by the 14C measurement and calibration processes, thus leading to systematically inaccurate age estimations. In short, our findings suggest the possibility of 14C-derived age–depth artefacts in the literature. Furthermore, since such age–depth artefacts are likely to coincide with large-scale changes in global Δ14C, which themselves can coincide with large-scale changes in global climate (such as the last deglaciation), 14C-derived age–depth artefacts may have been previously incorrectly attributed to changes in SAR coinciding with global climate. Our study highlights the need for the development of improved deep-sea sediment 14C calibration techniques that include an a priori representation of bioturbation for multi-specimen samples.

Description: The present study investigates the response of the high-latitude carbon cycle to changes in atmospheric greenhouse gas (GHG) concentrations in idealized climate change scenarios. To this end we use an adapted version of JSBACH – the land surface component of the Max Planck Institute for Meteorology Earth System Model (MPI-ESM) – that accounts for the organic matter stored in the permafrost-affected soils of the high northern latitudes. The model is run under different climate scenarios that assume an increase in GHG concentrations, based on the Shared Socioeconomic Pathway 5 and the Representative Concentration Pathway 8.5, which peaks in the years 2025, 2050, 2075 or 2100, respectively. The peaks are followed by a decrease in atmospheric GHGs that returns the concentrations to the levels at the beginning of the 21st century, reversing the imposed climate change. We show that the soil CO2 emissions exhibit an almost linear dependence on the global mean surface temperatures that are simulated for the different climate scenarios. Here, each degree of warming increases the fluxes by, very roughly, 50 % of their initial value, while each degree of cooling decreases them correspondingly. However, the linear dependence does not mean that the processes governing the soil CO2 emissions are fully reversible on short timescales but rather that two strongly hysteretic factors offset each other – namely the net primary productivity and the availability of formerly frozen soil organic matter. In contrast, the soil methane emissions show a less pronounced increase with rising temperatures, and they are consistently lower after the peak in the GHG concentrations than prior to it. Here, the net fluxes could even become negative, and we find that methane emissions will play only a minor role in the northern high-latitude contribution to global warming, even when considering the high global warming potential of the gas. Finally, we find that at a global mean temperature of roughly 1.75 K (±0.5 K) above pre-industrial levels the high-latitude ecosystem turns from a CO2 sink into a source of atmospheric carbon, with the net fluxes into the atmosphere increasing substantially with rising atmospheric GHG concentrations. This is very different from scenario simulations with the standard version of the MPI-ESM, in which the region continues to take up atmospheric CO2 throughout the entire 21st century, confirming that the omission of permafrost-related processes and the organic matter stored in the frozen soils leads to a fundamental misrepresentation of the carbon dynamics in the Arctic.

Description: A realistic simulation of the surface mass balance (SMB) is essential for simulating past and future ice-sheet changes. As most state-of-the-art Earth system models (ESMs) are not capable of realistically representing processes determining the SMB, most studies of the SMB are limited to observations and regional climate models and cover the last century and near future only. Using transient simulations with the Max Planck Institute ESM in combination with an energy balance model (EBM), we extend previous research and study changes in the SMB and equilibrium line altitude (ELA) for the Northern Hemisphere ice sheets throughout the last deglaciation. The EBM is used to calculate and downscale the SMB onto a higher spatial resolution than the native ESM grid and allows for the resolution of SMB variations due to topographic gradients not resolved by the ESM. An evaluation for historical climate conditions (1980–2010) shows that derived SMBs compare well with SMBs from regional modeling. Throughout the deglaciation, changes in insolation dominate the Greenland SMB. The increase in insolation and associated warming early in the deglaciation result in an ELA and SMB increase. The SMB increase is caused by compensating effects of melt and accumulation: the warming of the atmosphere leads to an increase in melt at low elevations along the ice-sheet margins, while it results in an increase in accumulation at higher levels as a warmer atmosphere precipitates more. After 13 ka, the increase in melt begins to dominate, and the SMB decreases. The decline in Northern Hemisphere summer insolation after 9 ka leads to an increasing SMB and decreasing ELA. Superimposed on these long-term changes are centennial-scale episodes of abrupt SMB and ELA decreases related to slowdowns of the Atlantic meridional overturning circulation (AMOC) that lead to a cooling over most of the Northern Hemisphere.

Description: Variability of surface water masses of the Laptev and the East Siberian seas in August–September 2018 is studied using in situ and satellite data. In situ data were collected during the ARKTIKA-2018 expedition and then complemented with satellite-derived sea surface temperature (SST), salinity (SSS), sea surface height, wind speed, and sea ice concentration. The estimation of SSS fields is challenging in high-latitude regions, and the precision of soil moisture and ocean salinity (SMOS) SSS retrieval is improved by applying a threshold on SSS weekly error. For the first time in this region, the validity of DMI (Danish Meteorological Institute) SST and SMOS SSS products is thoroughly studied using ARKTIKA-2018 expedition continuous thermosalinograph measurements and conductivity–temperature–depth (CTD) casts. They are found to be adequate to describe large surface gradients in this region. Surface gradients and mixing of the river and the sea water in the ice-free and ice-covered areas are described with a special attention to the marginal ice zone at a synoptic scale. We suggest that the freshwater is pushed northward, close to the marginal ice zone (MIZ) and under the sea ice, which is confirmed by the oxygen isotope analysis. The SST-SSS diagram based on satellite estimates shows the possibility of investigating the surface water mass transformation at a synoptic scale and reveals the presence of river water on the shelf of the East Siberian Sea. The Ekman transport is calculated to better understand the pathway of surface water displacement on the shelf and beyond.

Description: Although numerous pollen records are available worldwide in various databases, their use for synthesis works is limited as the chronologies are, as yet, not harmonized globally, and temporal uncertainties are unknown. We present a chronology framework named LegacyAge 1.0 that includes harmonized chronologies of 2831 palynological records (out of 3471 available records), downloaded from the Neotoma Paleoecology Database (last access: April 2021) and 324 additional Asian records. All chronologies use the Bayesian framework implemented in Bacon version 2.5.3. Optimal parameter settings of priors (accumulation.shape, memory.strength, memory.mean, accumulation.rate, thickness) were identified based on previous experiences or iteratively after preliminary model inspection. The most common control points for the chronologies are radiocarbon dates (86.1 %), calibrated by the latest calibration curves (IntCal20 and SHcal20 for the terrestrial radiocarbon dates in the northern and southern hemispheres; Marine20 for marine materials). The original literature was consulted when dealing with obvious outliers and inconsistencies. Several major challenges when setting up the chronologies included the waterline issue (18.8 % of records), reservoir effect (4.9 %), and sediment deposition discontinuity (4.4 %). Finally, we numerically compare the LegacyAge 1.0 chronologies to the original ones and show that the chronologies of 95.4 % of records could be improved according to our assessment. Our chronology framework and revised chronologies provide the opportunity to make use of the ages and age uncertainties in synthesis studies of, for example, pollen-based vegetation and climate change. The LegacyAge 1.0 dataset and R code used are open-access and available at PANGAEA (https://doi.pangaea.de/10.1594/PANGAEA.933132) and Github (https://github.com/LongtermEcology/LegacyAge-1.0), respectively.

Description: Geological records show that vast proglacial lakes existed along the land terminating margins of palaeo ice sheets in Europe and North America. Proglacial lakes impact ice sheet dynamics by imposing marine-like boundary conditions at the ice margin. These lacustrine boundary conditions include changes in the ice sheet’s geometry, stress balance and frontal ablation and therefore affect the entire ice sheet’s mass balance. This interaction, however, has not been rigorously implemented in ice sheet models. In this study, the implementation of an adaptive lake boundary into the Parallel Ice Sheet Model (PISM) is described and applied to the glacial retreat of the Laurentide Ice Sheet (LIS). The results show that the presence of proglacial lakes locally enhances the ice flow. Along the continental ice margin, ice streams and ice lobes can be observed. Lacustrine terminating ice streams cause immense thinning of the ice sheet’s interior and thus play a significant role in the demise of the LIS. Due to the presence of lakes, a process similar to the marine ice sheet instability causes the collapse of the ice saddle over Hudson Bay, which blocked drainage via the Hudson Strait. In control experiments without a lake model, Hudson Bay is still glaciated at the end of the simulation. Future studies should target the development of parametrizations that better describe the glacial-lacustrine interactions.

Description: Nitrogen fixers, or diazotrophs, play a key role in the carbon and nitrogen cycle of the world oceans, but the controlling mechanisms are not comprehensively understood yet. The present study compares two paradigms on the ecological niche of diazotrophs in an Earth System Model (ESM). In our standard model configuration, which is representative for most of the state-of-the-art pelagic ecosystem models, diazotrophs take advantage of zooplankton featuring a lower food preference for diazotrophs than for ordinary phytoplankton. We compare this paradigm with the idea that diazotrophs are more competitive under oligotrophic conditions, characterized by low (dissolved, particulate, organic and inorganic) phosphorous availability. Both paradigms are supported by observational evidence and lead to a similar good agreement to the most recent and advanced observation-based nitrogen fixation estimate in our ESM framework. Further, we illustrate that the similarity between the two paradigms breaks in a RCP 8.5 anthropogenic emission scenario. We conclude that a more advanced understanding of the ecological niche of diazotrophs is mandatory for assessing the cycling of essential nutrients, especially under changing environmental conditions. Our results call for more in-situ measurements of cyanobacteria biomass if major controls of nitrogen fixation in the oceans are to be dissected.

Description: The availability of dissolved iron (dFe) exerts an important control on primary production. Recent ocean observation programs have provided information on dFe in many parts of the ocean, but knowledge is still limited concerning the rates of processes that control the concentrations and cycling of dFe in the ocean and hence the role of dFe as a determinant of global primary production. We constructed a three-dimensional gridded dataset of oceanic dFe concentrations by using both observations and a simple model of the iron cycle, and estimated the difference of processes among the ocean basins in controlling the dFe distributions. A Green's function approach was used to integrate the observations and the model. The reproduced three-dimensional dFe distribution indicated that iron influx from aeolian dust and from shelf sediment were 7.6 Gmol yr and 4.4 Gmol yr in the Atlantic Ocean and 0.4 Gmol yr and 4.1 Gmol yr in the Pacific Ocean. The residence times were estimated to be 12.2 years in the Atlantic and 80.4 years in the Pacific. These estimates imply large differences in the cycling of dFe between the two ocean basins that would need to be taken into consideration when projecting future iron biogeochemical cycling under different climate change scenarios. Although there is some uncertainty in our estimates, global estimates of iron cycle characteristics based on this approach can be expected to enhance our understanding of the material cycle and hence of the current and future rates of marine primary production.