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: This report summarizes the first WP6 consultation on ocean liming, which focused on the life-cycle assessment (LCA) of lime produced within existing industrial processes for carbon dioxide removal through ocean alkalinity enhancement (see Deliverable 6.2 for further details on the LCA).

Description: This document is the first of three periodic reports on available data.

Description: This report provides an overview on the cooperation activities with relevant EU projects and other initiatives on negative emission technologies in the first phase of the OceanNETs project (01.07.2020 – 31.12.2021).

Description: This report provides an overview of OceanNETs engagement with relevant stakeholders and the OceanNETs stakeholder reference group in the first 18 month of the project (01.07.2020 – 31.12.2021).

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: Research undertaken in Task 2.2 identified a range of governance challenges to ocean-based NETs related to the global ocean governance framework, e.g., linked to the transboundary nature of the ocean, potential effects of ocean-based NETs on the ocean’s condition and marine ecosystem services, as well as the many unknowns and uncertainties linked to NET-deployment. The fragmented approaches and frameworks in place to govern the global ocean further complicate comprehensive governance of these emerging technologies. This deliverable presents results from a workshop that explored how oceanbased NETs should be governed to best confront these challenges and integrate international climate targets as well as global goals for ocean and biodiversity conservation, in addition to global ambitions towards sustainable development. The workshop is part of research undertaken by Task 2.2 to assess how ocean-based NETs are addressed by the current global ocean governance framework and develop governance scenarios and recommendations to policy makers for a “good governance” of NETs in the ocean.

Description: This working paper presents first insights on lay public perceptions of marine carbon dioxide removal (CDR) approaches. In seven focus groups, three in Germany and four in Norway (including one pilot) the researchers asked members of the lay public to share their views of the ocean and the effects of climate change, four CDR approaches, as well as their reflections on responsible research and innovation (RRI) of marine CDR. The four CDR methods were ocean iron fertilization, ocean alkalinity enhancement, artificial upwelling, and blue carbon management through restoration of coastal and marine ecosystems. In addition, respondents were asked to compare the four approaches. Our findings indicate that the public will be very supportive of blue carbon management irrespective of its actual carbon sequestration potential, due in part to the perceived bad state of marine ecosystems worldwide. Participants were skeptical whether any of the CDR approaches could have relevant effect on carbon sequestration and long-term storage; they reasoned about issues such as the ability to scale up treatments in time and space, unforeseen or unforeseeable effects on ecosystems in time and space, and the role of industry in the implementation process. They argued that despite the potential availability of marine CDR, industry and the general public should stop polluting behaviors and practices. Nevertheless, the participants universally agreed that further research on all four CDR methods should be pursued to better understand effects on climate, ecosystems, local communities, and the economy.

Description: This is the stylized description of our ocean liming case study, which we are using the introduce our life-cycle assessment to stakeholders in our consultation process

Description: Any integration of extra carbon dioxide removal (CDR) via terrestrial or marine sink enhancement into climate policies requires accounting for their effectiveness in reducing atmospheric carbon concentration and translating this information into the amount of carbon credits (to be used in official and voluntary emission trading schemes). Here, we assess accounting schemes in their appropriateness of assigning carbon credits. We discuss the role of temporary carbon storage and present the various ccounting methods for carbon credit assignment. We explain how we have implemented the methods numerically and analyse carbon assignments across the different accounting schemes, using stylized, model-based ocean sink enhancement experiments.

Description: This report summarizes OceanNETs stakeholder engagement activities in the Canary Islands, prior and during the mesocosm study in ocean alkalinity enhancement carried out on the island of Gran Canaria in September-October 2021. It also presents ideas for future stakeholder engagement in the islands.

Description: Realistic alkalinization scenarios, under the global context, are proposed and examined, which can be extent to include spatial considerations and specific technical and regulatory constraints. Results provide a set of stylistic projections of total mineral (carbonate and silicate) addition, with its temporal timeframe spanning from as early as 2025 up to 2100). Among others, these estimates can be used to constrain model simulations that will be carried out in Work Package 4.

Description: We applied a 1-D plankton ecosystem-biogeochemical model to assess the impacts of ocean alkalinity enhancement (OAE) on seasonal changes in biogeochemistry and plankton dynamics. Depending on deployment scenarios, OAE should theoretically have variable effects on pH and seawater pCO2, which might in turn affect (a) plankton growth conditions and (b) the efficiency of carbon dioxide removal (CDR) via OAE. Thus, a major focus of our work is how different magnitudes and temporal frequencies of OAE might affect seasonal response patterns of net primary productivity (NPP), ecosystem functioning and biogeochemical cycling. With our study we aimed at identifying a parameterization of how magnitude and frequency of OAE affect net growth rates, so that these effects could be employed for Earth System Modell applications. So far we learned that a meaningful response parameterization has to resolve positive and negative anomalies that covary with temporal shifts. As to the intricacy of the response patterns, the derivation of such parameterization is work in progress. However, our study readily provides valuable insights to how OAE can alter plankton dynamics and biogeochemistry. Our modelling study first focuses at a local site where time series data are available (European Station for Time series in the Ocean Canary Islands ESTOC), including measurements of pH, concentrations of total alkalinity, dissolved inorganic carbon (DIC), chlorophyll-a and dissolved inorganic nitrogen (DIN). These observational data were made available by Andres Cianca (personal communication, PLOCAN, Spain), Melchor Gonzalez and Magdalena Santana Casiano (personal communication, Universidad de Las Palmas de Gran Canaria). The choice of this location was underpinned by the fact that the first OAE mesocosm experiment was conducted on the Canary Island Gran Canaria, which will facilitate synthesizing our modelling approach with experimental findings. For our simulations at the ESTOC site in the Subtropical North Atlantic we found distinct, non-linear responses of NPP to different temporal modes of alkalinity deployment. In particular, phytoplankton bloom patterns displayed pronounced temporal phase shifts and changes in their amplitude. Notably, our simulations suggest that OAE can have a slightly stimulating effect on NPP, which is however variable, depending on the magnitude of OAE and the temporal mode of alkalinity addition. Furthermore, we find that increasing alkalinity perturbations can lead to a shift in phytoplankton community composition (towards coccolithophores), which even persists after OAE has stopped. In terms of CDR, we found that a decrease in efficiency with increasing magnitude of alkalinity addition, as well as substantial differences related to the timing of addition. Altogether, our results suggest that annual OAE during the right season (i.e. physical and biological conditions), could be a reasonable compromise in terms of logistical feasibility, efficiency of CDR and side-effects on marine biota. With respect to transferability to global models, the complex, non-linear responses of biological processes to OAE identified in our simulations do not allow for simple parameterizations that can easily adapted. Dedicated future work is required to transfer the observed responses at small spatiotemporal scales to the coarser resolution of global models.

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.

Description: Net-zero climate policies foresee deployment of atmospheric carbon dioxide removal wit geological, terrestrial, or marine carbon storage. While terrestrial and geological storage would be governed under the framework of national property rights, marine storage implies that carbon is transferred from one global common, the atmosphere, to another global common, the ocean, in particular if storage exceeds beyond coastal applications. This paper investigates the option of carbon dioxide removal (CDR) and storage in different (marine) reservoir types in an analytic climate-economy model, and derives implications for optimal mitigation efforts and CDR deployment. We show that the introduction of CDR lowers net energy input and net emissions over the entire time path. Furthermore, CDR affects the Social Cost of Carbon (SCC) via changes in total economic output but leaves the analytic structure of the SCC unchanged. In the first years after CDR becomes available the SCC is lower and in later years it is higher compared to a standard climate-economy model. Carbon dioxide emissions are first higher and then lower relative to a world without CDR. The paper provides the basis for the analysis of decentralized and potentially non-cooperative CDR policies.

Description: This study uses an existing perturbed parameter ensemble (PPE) of simulated ocean CO2 removal (CDR) to better determine sustainable pathways of ocean-based NET deployment and to provide information to constrain the design of subsequent modelling experiments. The results show that ocean alkalinity enhancement (OAE) can only help meet SDG13 (Climate Action) when other ambitious mitigation efforts are taken. This reinforces that OAE is not a substitute for emissions reduction, but could contribute to meeting our climate goals (if other factors suggest OAE is worth doing). For SDG14 (Life Below Water), the results suggest OEA can contribute to limiting or even reversing ocean acidification. Meeting many other SDG14 objectives is closely linked to also meeting SDG13. A key recommendation is therefore, that subsequent simulations in OceanNETs should only use SDG13 compatible baseline scenarios, unless there is some specific need for process understanding at higher levels of climate change. The analysis has also determined that the idealized CDR in the PPE is not suitable for determining many socio-economic constraints and the implications that these have for meeting the SDGs. Another key recommendation is therefore, that subsequent simulations within OceanNETs should use more realistic scenarios of CDR deployment.

Description: In part one of this deliverable, an ensemble of 14 CMIP6 Earth System Models is evaluated regarding their performance in simulating alkalinity and related parameters. The majority of the models and the multi-model-mean underestimate surface alkalinity compared to climatological observations. Alkalinity biases stemming from the parametrization of calcium carbonate formation and dissolution can be as big as biases stemming from model physics. In part two, we test the sensitivity of parametrizations concerning the carbonate chemistry in the FESOM2.1-REcoM3 and give recommendations for addressing alkalinity biases.

Description: This is the data management plan for the research project OceanNETs. It compiles OceanNETs research data output and describes the data handling during and after the projects duration with the aim to make OceanNETs research data FAIR – sustainably available for the scientific community. This data management plan is a living document; it will be continuously developed in close cooperation with the consortium members throughout the project duration. Version 2

Description: The potential biogeochemical and ecological impacts of ocean alkalinity enhancement were tested in a 5-weeks mesocosm experiment conducted in the subtropical, oligotrophic waters off Gran Canaria in September/October 2021. In the nine mesocosms, each with a volume of about 10 m3 inhabiting a natural plankton community, alkalinity enhancement was achieved through addition of a mix of sodium bicarbonate and sodium carbonate, simulating CO2-equilibrated alkalinization in a gradient from control up to twice the natural alkalinity. The response of the enclosed plankton community to the alkalinity addition was monitored in over 50 parameters which were sampled or measured in situ daily or every second day. In addition to the mesocosm experiment, a series of side experiments were conducted, focusing on individual aspects of mineral dissolution, secondary precipitation and biological responses at the primary producer level. This campaign, in which 47 scientists from 6 nations participated, generated the most comprehensive data set collected so far on the ecological and biogeochemical impacts of ocean alkalinity enhancement.

Description: Emissions trading systems (ETS) and markets usually do not allow for the inclusion of carbon dioxide removal (CDR) activities and if they do, removal activities are primarily restricted to afforestation. The New Zealand emission trading system (NZ ETS), for examples, integrates afforestation, and the California Low-Fuel Standard, the Quebec ETS and the Chinese ETS permit the restricted inclusion of afforestation offsets. Furthermore, the California Low-Carbon Fuel Standard System allows for the inclusion of removal via Direct Air Capture. In combination with the 45Q tax credit program, the largest incentives for CDR via Negative Emissions Technologies (NETs) are currently provided in the US. However, both do not yet allow for the inclusion of ocean-based carbon removal. Hence, we provide first a brief overview about the NZ ETS and its inclusion of afforestation, pointing out that the concept will likely not be applicable to ocean-based CDR with the potential exemption of blue carbon projects. Second, we discuss the 45Q tax credit program, the California Low-Fuel Standard System, and the California Compliance Offset Scheme. Third, we provide an overview about the company database related to ocean-based carbon removal. Fourth, we briefly look at the voluntary carbon market, providing some insights for carbon removal accounting.

Description: This deliverable synthesizes the results on public perceptions of marine CDR methods from the first two years of OceanNETs. The purpose is to inform the other work packages in OceanNETs and stakeholders about our results in a timely and brief manner about the ways members of the public view marine CDR specifically but also in the broader context of net-zero targets and climate policy. The deliverable summarises results of two studies: (1) focus groups held in Germany and Norway that covered ocean fertilization, ocean alkalinity enhancement, artificial upwelling and blue carbon management and (2) a deliberative survey in Norway that covered ocean alkalinity enhancement, macroalgae farming with BECCS or biomass sinking and land-based BECCS and enhanced weathering as terrestrial approaches for comparison. Participants in both studies emphasise the importance of reducing emissions and changing consumptions patterns. They hardly discuss the need to remove CO2 from the atmosphere to reach the Paris climate goal and the concept of negative emissions seems difficult for them to engage with. Among the methods, participants prefer ecosystem-based approaches like mangrove or seagrass restoration over other methods like alkalinity enhancement or ocean fertilization. Participants are concerned about the actual feasibility of deployment at a relevant removal scale and for a longer period. Connected to this are concerns about the controllability of the deployment and the methods’ impact, like difficulties to control negative environmental effects from biomass sinking at the seafloor. They also question the buildup of additional infrastructure or additional interventions into nature on top of already existing human interference. The opportunity to deliberate the methods increases participants’ certainty about their assessment but only slightly changes the direction of the assessment.