ALBERT

Description: Deposition to the sea surface is a major atmospheric loss pathway for many important trace gases, such as sulfur dioxide (SO2). The air–sea transfer of SO2 is controlled entirely on the atmospheric side of the air–sea interface due to high effective solubility and other physical–chemical properties. There have been few direct field measurements of such fluxes due to the challenges associated with making fast-response measurements of highly soluble trace gases at very low ambient levels. In this study, we report direct eddy covariance air–sea flux measurements of SO2, sensible heat, water vapor, and momentum. The measurements were made over shallow coastal waters from the Scripps Pier, La Jolla, CA, using negative ion chemical ionization mass spectrometry as the SO2 sensor. The observed transfer velocities for SO2, sensible heat, water vapor, and momentum and their wind speed dependences indicate that SO2 fluxes can be reliably measured using this approach. As expected, the transfer velocities for SO2, sensible heat, and water vapor are lower than that for momentum, demonstrating the contribution of molecular diffusion to the overall air-side resistance to gas transfer. Furthermore, transfer velocities of SO2 were lower than those of sensible heat and water vapor when observed simultaneously. This result is attributed to diffusive resistance in the interfacial layer of the air–sea interface.

Description: Low concentrations of ice-nucleating particles (INPs) are thought to be important for the properties of mixed-phase clouds, but their detection is challenging. Hence, there is a need for instruments where INP concentrations of less than 0.01 L−1 can be routinely and efficiently determined. The use of larger volumes of suspension in drop assays increases the sensitivity of an experiment to rarer INPs or rarer active sites due to the increase in aerosol or surface area of particulates per droplet. Here we describe and characterise the InfraRed-Nucleation by Immersed Particles Instrument (IR-NIPI), a new immersion freezing assay that makes use of IR emissions to determine the freezing temperature of individual 50 µL droplets each contained in a well of a 96-well plate. Using an IR camera allows the temperature of individual aliquots to be monitored. Freezing temperatures are determined by detecting the sharp rise in well temperature associated with the release of heat caused by freezing. In this paper we first present the calibration of the IR temperature measurement, which makes use of the fact that following ice nucleation aliquots of water warm to the ice–liquid equilibrium temperature (i.e. 0 ∘C when water activity is ∼1), which provides a point of calibration for each individual well in each experiment. We then tested the temperature calibration using ∼100 µm chips of K-feldspar, by immersing these chips in 1 µL droplets on an established cold stage (µL-NIPI) as well as in 50 µL droplets on IR-NIPI; the results were consistent with one another, indicating no bias in the reported freezing temperature. In addition we present measurements of the efficiency of the mineral dust NX-illite and a sample of atmospheric aerosol collected on a filter in the city of Leeds. NX-illite results are consistent with literature data, and the atmospheric INP concentrations were in good agreement with the results from the µL-NIPI instrument. This demonstrates the utility of this approach, which offers a relatively high throughput of sample analysis and access to low INP concentrations.

Description: The European Slope Current provides a shelf-edge conduit for Atlantic Water, a substantial fraction of which is destined for the northern North Sea, with implications for regional hydrography and ecosystems. Drifters drogued at 50 m in the European Slope Current at the Hebridean shelf break follow a wide range of pathways, indicating highly variable Atlantic inflow to the North Sea. Slope Current pathways, timescales and transports over 1988–2007 are further quantified in an eddy-resolving ocean model hindcast. Particle trajectories calculated with model currents indicate that Slope Current water is largely recruited from the eastern subpolar North Atlantic. Observations of absolute dynamic topography and climatological density support theoretical expectations that Slope Current transport is to first order associated with meridional density gradients in the eastern subpolar gyre, which support a geostrophic inflow towards the slope. In the model hindcast, Slope Current transport variability is dominated by abrupt 25–50 % reductions of these density gradients over 1996–1998. Concurrent changes in wind forcing, expressed in terms of density gradients, act in the same sense to reduce Slope Current transport. This indicates that coordinated regional changes of buoyancy and wind forcing acted together to reduce Slope Current transport during the 1990s. Particle trajectories further show that 10–40 % of Slope Current water is destined for the northern North Sea within 6 months of passing to the west of Scotland, with a general decline in this percentage over 1988–2007. Salinities in the Slope Current correspondingly decreased, evidenced in ocean analysis data. Further to the north, in the Atlantic Water conveyed by the Slope Current through the Faroe–Shetland Channel (FSC), salinity is observed to increase over this period while declining in the hindcast. The observed trend may have broadly compensated for a decline in the Atlantic inflow, limiting salinity changes in the northern North Sea during this period. Proxies for both Slope Current transport and Atlantic inflow to the North Sea are sought in sea level height differences across the FSC and between Shetland and the Scottish mainland (Wick). Variability of Slope Current transport on a wide range of timescales, from seasonal to multi-decadal, is implicit in sea level differences between Lerwick (Shetland) and Tórshavn (Faroes), in both tide gauge records from 1957 and a longer model hindcast spanning 1958–2012. Wick–Lerwick sea level differences in tide gauge records from 1965 indicate considerable decadal variability in the Fair Isle Current transport that dominates Atlantic inflow to the northwest North Sea, while sea level differences in the hindcast are dominated by strong seasonal variability. Uncertainties in the Wick tide gauge record limit confidence in this proxy.

Description: Surface ozone (O3) pollution levels are strongly correlated with daytime surface temperatures, especially in highly polluted regions. This correlation is nonlinear and occurs through a variety of temperature-dependent mechanisms related to O3 precursor emissions, lifetimes, and reaction rates, making the reproduction of temperature sensitivities – and the projection of associated human health risks – a complex problem. Here we explore the summertime O3–temperature relationship in the United States and Europe using the chemical transport model GEOS-Chem. We remove the temperature dependence of several mechanisms most frequently cited as causes of the O3–temperature “climate penalty”, including PAN decomposition, soil NOx emissions, biogenic volatile organic compound (VOC) emissions, and dry deposition. We quantify the contribution of each mechanism to the overall correlation between O3 and temperature both individually and collectively. Through this analysis we find that the thermal decomposition of PAN can explain, on average, 20 % of the overall O3–temperature correlation in the United States. The effect is weaker in Europe, explaining 9 % of the overall O3–temperature relationship. The temperature dependence of biogenic emissions contributes 3 % and 9 % of the total O3–temperature correlation in the United States and Europe on average, while temperature-dependent deposition (6 % and 1 %) and soil NOx emissions (10 % and 7 %) also contribute. Even considered collectively these mechanisms explain less than 46 % of the modeled O3–temperature correlation in the United States and 36 % in Europe. We use commonality analysis to demonstrate that covariance with other meteorological phenomena such as stagnancy and humidity can explain the bulk of the remainder of the O3–temperature correlation. Thus, we demonstrate that the statistical correlation between O3 and temperature alone may greatly overestimate the direct impacts of temperature on O3, with implications for the interpretation of policy-relevant metrics such as climate penalty.

Description: We report on the appearance of a 2 km wide, 70 m deep circular depression located 50 km inland of the southwestern margin of the Greenland Ice Sheet that provides the first direct evidence for concentrated, long-term storage, and sudden release, of meltwater at the bed. Drainage of the lake may have been triggered by the recent increase in meltwater runoff. The abundance of such lakes and their potential importance to the ice sheet's hydrologic system and flow regime remain unknown.

Description: The Antarctic outer coastal margin (i.e., the coastline itself, or the terminus/front of ice shelves, whichever is adjacent to the ocean) is a key interface between the ice-sheet and terrestrial environments and the Southern Ocean. Its physical configuration (including both length scale of variation and orientation/aspect) has direct bearing on several closely associated cryospheric, biological, oceanographical and ecological processes, yet no study has quantified the coastal complexity or orientation of Antarctica’s coastal margin. This first-of-a-kind characterisation of Antarctic coastal complexity aims to address this knowledge gap. We quantify and investigate the physical configuration and complexity of Antarctica’s circumpolar outer coastal margin using a novel, technique based on ∼ 40,000 random points selected along a vector coastline derived from the MODIS Mosaic of Antarctica dataset. At each point, a complexity metric is calculated at length scales from 1 to 256 km, giving a multiscale estimate of the magnitude and direction of undulation or complexity at each point location along the entire coastline. Using a cluster analysis to determine characteristic complexity ‘signatures’ for random nodes, the coastline is found to comprise three basic groups or classes: (i) low complexity at all scales; (ii) most complexity at shorter scales; and (iii) most complexity at longer scales. These classes are somewhat heterogeneously distributed throughout the continent. We also consider bays and peninsulas separately and characterise their multi-scale orientation. This unique dataset and its summary analysis have numerous applications for both geophysical and biological studies and will contribute to Antarctic research requiring quantitative information on, and related to, coastal complexity and configuration. All these data are referenced by https://doi.org/10.26179/5d1af0ba45c03 and are available free of charge at https://data.antarctica.gov.au.

Description: A unique stratigraphic sequence of fossil leaves of Eotrigonobalanus furcinervis (extinct trees of the beech family, Fagaceae) from central Germany has been used to derive an atmospheric pCO2 record with multiple data points spanning the late middle to late Eocene, two sampling levels which may be earliest Oligocene, and two samples from later in the Oligocene. Using the inverse relationship between the density of stomata and pCO2, we show that pCO2 decreased continuously from the late middle to late Eocene, reaching a relatively stable low value before the end of the Eocene. Based on the subsequent records, pCO2 in parts of the Oligocene was similar to latest Eocene values. These results suggest that a decrease in pCO2 preceded the large shift in marine oxygen isotope records that characterizes the Eocene–Oligocene transition and that when a certain threshold of pCO2 change was crossed, the cumulative effects of this and other factors resulted in rapid temperature decline, ice build up on Antarctica and hence a change of climate mode.

Description: A unique stratigraphic sequence of fossil leaves of Eotrigonobalanus furcinervis (extinct trees of the beech family, Fagaceae) from central Germany has been used to derive an atmospheric pCO2 record with multiple data points spanning the late middle to late Eocene, two sampling levels which may be earliest Oligocene, and two samples from later in the Oligocene. Using the inverse relationship between the density of stomata and pCO2, we show that pCO2 decreased continuously from the late middle to late Eocene, reaching a relatively stable low value before the end of the Eocene. Based on the subsequent records, pCO2 in parts of the Oligocene was similar to latest Eocene values. These results show that a decrease in pCO2 preceded the large shift in marine oxygen isotope records that characterizes the Eocene–Oliogocene transition. This may be related to the "hysteresis effect" previously proposed – where a certain threshold of pCO2 change was crossed before the cumulative effects of this and other factors resulted in rapid temperature decline, ice build up on Antarctica and hence a change of climate mode.

Description: We present an approach which we call PSyKAl that is designed to achieve portable performance for parallel, finite-difference Ocean models. In PSyKAl the code related to the underlying science is formally separated from code related to parallelisation and single-core optimisations. This separation of concerns allows scientists to code their science independently of the underlying hardware architecture and for optimisation specialists to be able to tailor the code for a particular machine independently of the science code. We have taken the free-surface part of the NEMO ocean model and created a new, shallow-water model named NEMOLite2D. In doing this we have a code which is of a manageable size and yet which incorporates elements of full ocean models (input/output, boundary conditions, etc.). We have then manually constructed a PSyKAl version of this code and investigated the transformations that must be applied to the middle/PSy layer in order to achieve good performance, both serial and parallel. We have produced versions of the PSy layer parallelised with both OpenMP and OpenACC; in both cases we were able to leave the natural-science parts of the code unchanged while achieving good performance on both multi-core CPUs and GPUs. In quantifying whether or not the obtained performance is `good' we also consider the limitations of the basic roofline model and improve on it by generating kernel-specific CPU ceilings.

Description: Weather and climate models are complex pieces of software which include many individual components, each of which is evolving under the pressure to exploit advances in computing to enhance some combination of a range of possible improvements (higher spatio/temporal resolution, increased fidelity in terms of resolved processes, more quantification of uncertainty etc). However, after many years of a relatively stable computing environment with little choice in processing architecture or programming paradigm (basically X86 processors using MPI for parallelism), the existing menu of processor choices includes significant diversity, and more is on the horizon. This computational diversity, coupled with ever increasing software complexity, leads to the very real possibility that weather and climate modelling will arrive at a chasm which will separate scientific aspiration from our ability to develop and/or rapidly adapt codes to the available hardware. In this paper we review the hardware and software trends which are leading us towards this chasm, before describing current progress in addressing some of the tools which we may be able to use to bridge the chasm. This brief introduction to current tools and plans is followed by a discussion outlining the scientific requirements for quality model codes which have satisfactory performance and portability, while simultaneously supporting productive scientific evolution. We assert that the existing method of incremental model improvements employing small steps which adjust to the changing hardware environment is likely to be inadequate for crossing the chasm between aspiration and hardware at a satisfactory pace, in part because institutions cannot have all the relevant expertise in house. Instead, we outline a methodology based on large community efforts in engineering and standardisation, one which will depend on identifying a taxonomy of key activities – perhaps based on existing efforts to develop domain specific languages, identify common patterns in weather and climate codes, and develop community approaches to commonly needed tools, libraries etc – and then collaboratively building up those key components. Such a collaborative approach will depend on institutions, projects and individuals adopting new interdependencies and ways of working.