ALBERT

Description: The University of Victoria Earth system climate model of intermediate complexity has been a useful tool in recent assessments of long-term climate changes including paleo-climate modelling. Since the last official release of the UVic ESCM 2.9, and the two official updates during the last decade, a lot of model development has taken place in multiple groups. The new version 2.10 of the University of Victoria Earth System Climate Model (UVic ESCM), to be used in the 6th phase of the coupled model intercomparison project (CMIP6), presented here combines and brings together multiple model developments and new components that have taken place since the last official release of the model. To set the foundation of its use, we here describe the UVic ESCM 2.10 and evaluate results from transient historical simulations against observational data. We find that the UVic ESCM 2.10 is capable of reproducing well changes in historical temperature and carbon fluxes, as well as the spatial distribution of many ocean tracers, including temperature, salinity, phosphate and nitrate. This is connected to a good representation of ocean physical properties. For the moment, there remain biases in ocean alkalinity and dissolved inorganic carbon, which will be addressed in the next updates to the model.

Description: Combining proxy information and climate model simulations reconciles these sources of information about past climates. This, in turn, strengthens our understanding of past climatic changes. The analogue or proxy surrogate reconstruction method is a computationally cheap data assimilation approach, which searches in a pool of simulated climate states the best fit to proxy data. We use the approach to reconstruct European summer mean temperature from the 13th century until present using the Euro 2k set of proxy records and a pool of global climate simulation output fields. Our focus is on quantifying the uncertainty of the reconstruction, because previous applications of the analogue method rarely provided uncertainty ranges. We show several ways of estimating reconstruction uncertainty for the analogue method, which take into account the non-climate part of the variability in each proxy record. In general, our reconstruction agrees well at multi-decadal timescales with the Euro 2k reconstruction, which was conducted with two different statistical methods and no information from model simulations. In both methodological approaches, the decades around the year 1600 CE were the coldest. However, the approaches disagree on the warmest pre-industrial periods. The reconstructions from the analogue method also represent the local variations of the observed proxies. The diverse uncertainty estimates obtained from our analogue approaches can be locally larger or smaller than the estimates from the Euro 2k effort. Local uncertainties of the temperature reconstructions tend to be large in areas that are poorly covered by the proxy records. Uncertainties highlight the ambiguity of field-based reconstructions constrained by a limited set of proxies.

Description: The Parallel Ice Sheet Model (PISM) is applied to the Antarctic Ice Sheet over the last two glacial cycles (≈210 000 years) with a resolution of 16 km. An ensemble of 256 model runs is analyzed in which four relevant model parameters have been systematically varied using full-factorial parameter sampling. Parameters and plausible parameter ranges have been identified in a companion paper (Albrecht et al., 2020) and are associated with ice dynamics, climatic forcing, basal sliding and bed deformation and represent distinct classes of model uncertainties. The model is scored against both modern and geologic data, including reconstructed grounding-line locations, elevation–age data, ice thickness, surface velocities and uplift rates. An aggregated score is computed for each ensemble member that measures the overall model–data misfit, including measurement uncertainty in terms of a Gaussian error model (Briggs and Tarasov, 2013). The statistical method used to analyze the ensemble simulation results follows closely the simple averaging method described in Pollard et al. (2016). This analysis reveals clusters of best-fit parameter combinations, and hence a likely range of relevant model and boundary parameters, rather than individual best-fit parameters. The ensemble of reconstructed histories of Antarctic Ice Sheet volumes provides a score-weighted likely range of sea-level contributions since the Last Glacial Maximum (LGM) of 9.4±4.1 m (or 6.5±2.0×106km3 ), which is at the upper range of most previous studies. The last deglaciation occurs in all ensemble simulations after around 12 000 years before present and hence after the meltwater pulse 1A (MWP1a). Our ensemble analysis also provides an estimate of parametric uncertainty bounds for the present-day state that can be used for PISM projections of future sea-level contributions from the Antarctic Ice Sheet.

Description: Uncertainties in projections from Earth system models (ESMs) are associated to a large degree with the imperfect representation of the marine plankton ecosystem, in particular the physiology of primary and secondary producers. Here we describe the implementation of an optimality-based plankton-ecosystem model (OPEM) with variable C:N:P stoichiometry in the University of Victoria ESM (UVic) and the behaviour of two calibrated reference configurations, which differ in the assumed temperature dependence of diazotrophs. Predicted tracer distributions of oxygen and dissolved inorganic nutrients are similar to those of an earlier fixed-stoichiometry model (Keller et al., 2012). Compared to the classic fixed-stoichiometry model, OPEM is closer to recent satellite-based estimates of net community production (NCP), despite overestimating net primary production (NPP), can better reproduce deep-ocean gradients in the NO3:PO43− ratio, and partially explains observed patterns of particulate C:N:P in the surface ocean. Allowing diazotrophs to grow (but not necessarily fix N2) at similar temperatures as other phytoplankton results in a better representation of surface Chl and NPP in the Arctic and Antarctic Oceans. Deficiencies of our calibrated OPEM configurations may serve as a magnifying glass for shortcomings in global biogeochemical models and hence guide future model development. The overestimation of NPP at low latitudes indicates the need for improved representations of temperature effects on biotic processes, as well as phytoplankton community composition, which may be represented by locally-varying parameters based on suitable trade-offs. Discrepancies between observed and predicted vertical gradients in particulate C:N:P ratios suggest the need to include preferential P remineralisation, which could also benefit the representation of N2 fixation. While OPEM yields a much improved distribution of surface N* (NO3− 16·PO43− + 2.9 mmol m−3), it still fails to reproduce observed N* in the Arctic, possibly related to a mis-representation of the phytoplankton community there and the lack of benthic denitrification in the model. Coexisting ordinary and diazotrophic phytoplankton can exert strong control on N* in our simulations, which questions the interpretation of N* as reflecting the balance of N2 fixation and denitrification.

Description: Simulations of the glacial–interglacial history of the Antarctic Ice Sheet provide insights into dynamic threshold behavior and estimates of the ice sheet's contributions to global sea-level changes for the past, present and future. However, boundary conditions are weakly constrained, in particular at the interface of the ice sheet and the bedrock. Also climatic forcing covering the last glacial cycles is uncertain, as it is based on sparse proxy data. We use the Parallel Ice Sheet Model (PISM) to investigate the dynamic effects of different choices of input data, e.g., for modern basal heat flux or reconstructions of past changes of sea level and surface temperature. As computational resources are limited, glacial-cycle simulations are performed using a comparably coarse model grid of 16 km and various parameterizations, e.g., for basal sliding, iceberg calving, or for past variations in precipitation and ocean temperatures. In this study we evaluate the model's transient sensitivity to corresponding parameter choices and to different boundary conditions over the last two glacial cycles and provide estimates of involved uncertainties. We also discuss isolated and combined effects of climate and sea-level forcing. Hence, this study serves as a “cookbook” for the growing community of PISM users and paleo-ice sheet modelers in general. For each of the different model uncertainties with regard to climatic forcing, ice and Earth dynamics, and basal processes, we select one representative model parameter that captures relevant uncertainties and motivates corresponding parameter ranges that bound the observed ice volume at present. The four selected parameters are systematically varied in a parameter ensemble analysis, which is described in a companion paper.

Description: The project MarParCloud (Marine biological production, organic aerosol Particles and marine Clouds: a process chain) aims to improve our understanding of the genesis, modification and impact of marine organic matter (OM) from its biological production, to its export to marine aerosol particles and, finally, to its ability to act as ice-nucleating particles (INPs) and cloud condensation nuclei (CCN). A field campaign at the Cape Verde Atmospheric Observatory (CVAO) in the tropics in September–October 2017 formed the core of this project that was jointly performed with the project MARSU (MARine atmospheric Science Unravelled). A suite of chemical, physical, biological and meteorological techniques was applied, and comprehensive measurements of bulk water, the sea surface microlayer (SML), cloud water and ambient aerosol particles collected at a ground-based and a mountain station took place. Key variables comprised the chemical characterization of the atmospherically relevant OM components in the ocean and the atmosphere as well as measurements of INPs and CCN. Moreover, bacterial cell counts, mercury species and trace gases were analyzed. To interpret the results, the measurements were accompanied by various auxiliary parameters such as air mass back-trajectory analysis, vertical atmospheric profile analysis, cloud observations and pigment measurements in seawater. Additional modeling studies supported the experimental analysis. During the campaign, the CVAO exhibited marine air masses with low and partly moderate dust influences. The marine boundary layer was well mixed as indicated by an almost uniform particle number size distribution within the boundary layer. Lipid biomarkers were present in the aerosol particles in typical concentrations of marine background conditions. Accumulation- and coarse-mode particles served as CCN and were efficiently transferred to the cloud water. The ascent of ocean-derived compounds, such as sea salt and sugar-like compounds, to the cloud level, as derived from chemical analysis and atmospheric transfer modeling results, denotes an influence of marine emissions on cloud formation. Organic nitrogen compounds (free amino acids) were enriched by several orders of magnitude in submicron aerosol particles and in cloud water compared to seawater. However, INP measurements also indicated a significant contribution of other non-marine sources to the local INP concentration, as (biologically active) INPs were mainly present in supermicron aerosol particles that are not suggested to undergo strong enrichment during ocean–atmosphere transfer. In addition, the number of CCN at the supersaturation of 0.30 % was about 2.5 times higher during dust periods compared to marine periods. Lipids, sugar-like compounds, UV-absorbing (UV: ultraviolet) humic-like substances and low-molecular-weight neutral components were important organic compounds in the seawater, and highly surface-active lipids were enriched within the SML. The selective enrichment of specific organic compounds in the SML needs to be studied in further detail and implemented in an OM source function for emission modeling to better understand transfer patterns, the mechanisms of marine OM transformation in the atmosphere and the role of additional sources. In summary, when looking at particulate mass, we see oceanic compounds transferred to the atmospheric aerosol and to the cloud level, while from a perspective of particle number concentrations, sea spray aerosol (i.e., primary marine aerosol) contributions to both CCN and INPs are rather limited.

Description: The Ligurian Basin is located in the Mediterranean Sea to the north-west of Corsica at the transition from the Western Alpine orogen to the Apennine system and was generated by the south-eastward trench retreat of the Apennines–Calabrian subduction zone. Late-Oligocene-to-Miocene rifting caused continental extension and subsidence, leading to the opening of the basin. Yet it remains unclear if rifting caused continental break-up and seafloor spreading. To reveal its lithospheric architecture, we acquired a 130 km long seismic refraction and wide-angle reflection profile in the Ligurian Basin. The seismic line was recorded in the framework of SPP2017 4D-MB, a Priority Programme of the German Research Foundation (DFG) and the German component of the European AlpArray initiative, and trends in a NE–SW direction at the centre of the Ligurian Basin, roughly parallel to the French coastline. The seismic data were recorded on the newly developed GEOLOG recorder, designed at GEOMAR, and are dominated by sedimentary refractions and show mantle Pn arrivals at offsets of up to 70 km and a very prominent wide-angle Mohorovičić discontinuity (Moho) reflection. The main features share several characteristics (e.g. offset range, continuity) generally associated with continental settings rather than documenting oceanic crust emplaced by seafloor spreading. Seismic tomography results are complemented by gravity data and yield a ∼ 6–8 km thick sedimentary cover and the seismic Moho at 11–13 km depth below the sea surface. Our study reveals that the oceanic domain does not extend as far north as previously assumed. Whether Oligocene–Miocene extension led to extremely thinned continental crust or exhumed subcontinental mantle remains unclear. A low grade of mantle serpentinisation indicates a high rate of syn-rift sedimentation. However, rifting failed before oceanic spreading was initiated, and continental crust thickens towards the NE within the northern Ligurian Basin.

Description: We present a new near-global coupled biogeochemical ocean-circulation model configuration. The configuration features a horizontal discretization with a grid spacing of less than 11 km in the Southern Ocean and gradually coarsens in meridional direction to more than 200 km at 64∘ N, where the model is bounded by a solid wall. The underlying code framework is the Geophysical Fluid Dynamics Laboratory (GFDL)'s Modular Ocean Model coupled to the Biogeochemistry with Light, Iron, Nutrients and Gases (BLING) ecosystem model of Galbraith et al. (2010). The configuration is unique in that it features both a relatively equilibrated oceanic carbon inventory and an eddying ocean circulation based on a realistic model geometry/bathymetry – a combination that has been precluded by prohibitive computational cost in the past. Results from a simulation with climatological forcing and a sensitivity experiment with increasing winds suggest that the configuration is sufficiently equilibrated to explore Southern Ocean carbon uptake dynamics on decadal timescales. The configuration is dubbed MOMSO, a Modular Ocean Model Southern Ocean configuration.

Description: Author Posting. © American Geophysical Union, 2020. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Journal of Geophysical Research: Oceans 125(5), (2020): e2019JC016007, doi:10.1029/2019JC016007.

Description: Benthic inputs of nutrients help support primary production in the Chukchi Sea and produce nutrient‐rich water masses that ventilate the halocline of the western Arctic Ocean. However, the complex biological and redox cycling of nutrients and trace metals make it difficult to directly monitor their benthic fluxes. In this study, we use radium‐228, which is a soluble radionuclide produced in sediments, and a numerical model of an inert, generic sediment‐derived tracer to study variability in sediment inputs to the Chukchi Sea. The 228Ra observations and modeling results are in general agreement and provide evidence of strong benthic inputs to the southern Chukchi Sea during the winter, while the northern shelf receives higher concentrations of sediment‐sourced materials in the spring and summer due to continued sediment‐water exchange as the water mass traverses the shelf. The highest tracer concentrations are observed near the shelfbreak and southeast of Hanna Shoal, a region known for high biological productivity and enhanced benthic biomass.

Description: This study presents data from multiple Arctic expeditions over the past two decades, and we are indebted to the captains, crews, and scientific parties that made this data collection possible. This work was funded by NSF awards OCE‐1458305 to M. Charette, OCE‐1458424 to W. Moore, OCE‐1434085 to D. Kadko, PLR‐1504333 to R. Pickart, and OPP‐1822334 to M. Spall. Funding was also provided by National Oceanic and Atmospheric Administration Grant NA14‐OAR4320158 to R. Pickart. L. Kipp was supported by an Ocean Frontier Institute Postdoctoral Fellowship. Radium data used in this manuscript are available in Table S1.

Description: 2020-10-27

Description: The eastern tropical South Pacific (ETSP) represents one of the most productive areas in the ocean that is characterized by a pronounced oxygen minimum zone (OMZ). Particulate organic matter (POM) that sinks out of the euphotic zone is supplied to the anoxic sediments and utilized by microbial communities. The degradation of POM is associated with dissolved organic matter (DOM) production and reworking. The release of recalcitrant DOM to the overlying waters may represent an important organic matter escape mechanism from remineralization within sediments but received little attention in OMZ regions so far. Here, we combine measurements of dissolved organic carbon (DOC) and dissolved organic nitrogen (DON) with DOM optical properties in the form of chromophoric (CDOM) and fluorescent (FDOM) DOM from pore waters and near-bottom waters of the ETSP off Peru. We evaluate diffusion–driven fluxes and net in situ fluxes of DOC and DON in order to investigate processes affecting DOM cycling at the sediment–water interface along a transect 12° S. To our knowledge, these are the first data for sediment release of DON and pore water CDOM and FDOM for the ETSP off Peru. Pore-water DOC and DON accumulated with increasing sediment depth, suggesting an imbalance between DOM production and remineralization within sediments. High DON accumulation resulted in very low pore water DOC / DON ratios (〉 1) which could be caused by either an "imbalance" in DOC and DON remineralization, or to the presence of an additional nitrogen source. Diffusion driven fluxes of DOC and DON exhibited high spatial variability. They varied from 0.2–0.1 mmol m−2 d−1 to 2.52–1.3 mmol m−2 d−1 and from −0.042–0.02 mmol m−2 d−1 to 3.32–1.7 mmol m−2 d−1, respectively. Generally low net in situ DOC and DON fluxes as well as steepening of spectral slope (S) of CDOM and accumulation of humic-like FDOM at the near-bottom waters over time indicated active microbial DOM utilization at the sediment–water interface, potentially stimulated by nitrate (NO3−) and nitrite (NO2−). The microbial DOC utilization rates, estimated in our study, may be sufficient to support denitrification rates of 0.2–1.4 mmol m−2 d−1, suggesting that sediment release of DOM contributes substantially to nitrogen loss processes in the ETSP off Peru.