ALBERT

Description: The aerosol-driven radiative effects on marine low-level cloud represent a large uncertainty in climate simulations, in particular over the Southern Ocean, which is also an important region for sea spray aerosol production. Observations of sea spray aerosol organic enrichment and the resulting impact on water uptake over the remote Southern Hemisphere are scarce, and therefore the region is under-represented in existing parameterisations. The Surface Ocean Aerosol Production (SOAP) voyage was a 23 d voyage which sampled three phytoplankton blooms in the highly productive water of the Chatham Rise, east of New Zealand. In this study we examined the enrichment of organics to nascent sea spray aerosol and the modifications to sea spray aerosol water uptake using in situ chamber measurements of seawater samples taken during the SOAP voyage. Primary marine organics contributed up to 23 % of the sea spray mass for particles with diameter less than approximately 1 µm and up to 79 % of the particle volume for 50 nm diameter sea spray. The composition of the submicron organic fraction was consistent throughout the voyage and was largely composed of a polysaccharide-like component, characterised by very low alkane-to-hydroxyl-concentration ratios of approximately 0.1–0.2. The enrichment of organics was compared to the output from the chlorophyll-a-based sea spray aerosol parameterisation suggested by Gantt et al. (2011) and the OCEANFILMS (Organic Compounds from Ecosystems to Aerosols: Natural Films and Interfaces via Langmuir Molecular Surfactants) models. OCEANFILMS improved on the representation of the organic fraction predicted using chlorophyll a, in particular when the co-adsorption of polysaccharides was included; however, the model still under-predicted the proportion of polysaccharides by an average of 33 %. Nascent 50 nm diameter sea spray aerosol hygroscopic growth factors measured at 90 % relative humidity averaged 1.93±0.08 and did not decrease with increasing sea spray aerosol organic fractions. The observed hygroscopicity was greater than expected from the assumption of full solubility, particularly during the most productive phytoplankton bloom (B1), during which organic fractions were greater than approximately 0.4. The water uptake behaviour observed in this study is consistent with that observed for other measurements of phytoplankton blooms and can be partially attributed to the presence of sea salt hydrates, which lowers the sea spray aerosol hygroscopicity when the organic enrichment is low. The inclusion of surface tension effects only marginally improved the modelled hygroscopicity, and a significant discrepancy between the observed and modelled hygroscopicity at high organic volume fractions remained. The findings from the SOAP voyage highlight the influence of biologically sourced organics on sea spray aerosol composition; these data improve the capacity to parameterise sea spray aerosol organic enrichment and water uptake.

Description: In this study we use ozone and stratospheric ozone tracer simulations from the high-resolution (0.5∘×0.5∘) Goddard Earth Observing System, Version 5 (GEOS-5), in a replay mode to study the impact of stratospheric ozone on tropospheric ozone interannual variability (IAV). We use these simulations in conjunction with ozonesonde measurements from 1990 to 2016 during the winter and spring seasons. The simulations include a stratospheric ozone tracer (StratO3) to aid in the evaluation of the impact of stratospheric ozone IAV on the IAV of tropospheric ozone at different altitudes and locations. The model is in good agreement with the observed interannual variation in tropospheric ozone, except for the post-Pinatubo period (1992–1994) over the region of North America. Ozonesonde data show a negative ozone anomaly in 1992–1994 following the Pinatubo eruption, with recovery thereafter. The simulated anomaly is only half the magnitude of that observed. Our analysis suggests that the simulated stratosphere–troposphere exchange (STE) flux deduced from the analysis might be too strong over the North American (50–70∘ N) region after the Mt. Pinatubo eruption in the early 1990s, masking the impact of lower stratospheric ozone concentration on tropospheric ozone. European ozonesonde measurements show a similar but weaker ozone depletion after the Mt. Pinatubo eruption, which is fully reproduced by the model. Analysis based on the stratospheric ozone tracer identifies differences in strength and vertical extent of stratospheric ozone impact on the tropospheric ozone interannual variation (IAV) between North America and Europe. Over North American stations, the StratO3 IAV has a significant impact on tropospheric ozone from the upper to lower troposphere and explains about 60 % and 66 % of the simulated ozone IAV at 400 hPa and ∼11 % and 34 % at 700 hPa in winter and spring, respectively. Over European stations, the influence is limited to the middle to upper troposphere and becomes much smaller at 700 hPa. The Modern-Era Retrospective analysis for Research and Applications, Version 2 (MERRA-2), assimilated fields exhibit strong longitudinal variations over Northern Hemisphere (NH) mid-high latitudes, with lower tropopause height and lower geopotential height over North America than over Europe. These variations associated with the relevant variations in the location of tropospheric jet flows are responsible for the longitudinal differences in the stratospheric ozone impact, with stronger effects over North America than over Europe.

Description: Anthropogenic emissions of CO2 to the atmosphere have modified the carbon cycle for more than 2 centuries. As the ocean stores most of the carbon on our planet, there is an important task in unraveling the natural and anthropogenic processes that drive the carbon cycle at different spatial and temporal scales. We contribute to this by designing a global monthly climatology of total dissolved inorganic carbon (TCO2), which offers a robust basis in carbon cycle modeling but also for other studies related to this cycle. A feedforward neural network (dubbed NNGv2LDEO) was configured to extract from the Global Ocean Data Analysis Project version 2.2019 (GLODAPv2.2019) and the Lamont–Doherty Earth Observatory (LDEO) datasets the relations between TCO2 and a set of variables related to the former's variability. The global root mean square error (RMSE) of mapping TCO2 is relatively low for the two datasets (GLODAPv2.2019: 7.2 µmol kg−1; LDEO: 11.4 µmol kg−1) and also for independent data, suggesting that the network does not overfit possible errors in data. The ability of NNGv2LDEO to capture the monthly variability of TCO2 was testified through the good reproduction of the seasonal cycle in 10 time series stations spread over different regions of the ocean (RMSE: 3.6 to 13.2 µmol kg−1). The climatology was obtained by passing through NNGv2LDEO the monthly climatological fields of temperature, salinity, and oxygen from the World Ocean Atlas 2013 and phosphate, nitrate, and silicate computed from a neural network fed with the previous fields. The resolution is 1∘×1∘ in the horizontal, 102 depth levels (0–5500 m), and monthly (0–1500 m) to annual (1550–5500 m) temporal resolution, and it is centered around the year 1995. The uncertainty of the climatology is low when compared with climatological values derived from measured TCO2 in the largest time series stations. Furthermore, a computed climatology of partial pressure of CO2 (pCO2) from a previous climatology of total alkalinity and the present one of TCO2 supports the robustness of this product through the good correlation with a widely used pCO2 climatology (Landschützer et al., 2017). Our TCO2 climatology is distributed through the data repository of the Spanish National Research Council (CSIC; https://doi.org/10.20350/digitalCSIC/10551, Broullón et al., 2020).

Description: Uranium–lead (U–Pb) geochronology was conducted by laser ablation – inductively coupled plasma mass spectrometry (LA-ICPMS) on 7175 detrital zircon grains from 29 samples from the Coconino Sandstone, Moenkopi Formation, and Chinle Formation. These samples were recovered from ∼ 520 m of drill core that was acquired during the Colorado Plateau Coring Project (CPCP), located in Petrified Forest National Park (Arizona). A sample from the lower Permian Coconino Sandstone yields a broad distribution of Proterozoic and Paleozoic ages that are consistent with derivation from the Appalachian and Ouachita orogens, with little input from local basement or Ancestral Rocky Mountain sources. Four samples from the Holbrook Member of the Moenkopi Formation yield a different set of Precambrian and Paleozoic age groups, indicating derivation from the Ouachita orogen, the East Mexico arc, and the Permo-Triassic arc built along the Cordilleran margin. A total of 23 samples from the Chinle Formation contain variable proportions of Proterozoic and Paleozoic zircon grains but are dominated by Late Triassic grains. LA-ICPMS ages of these grains belong to five main groups that correspond to the Mesa Redondo Member, Blue Mesa Member and lower part of the Sonsela Member, upper part of the Sonsela Member, middle part of the Petrified Forest Member, and upper part of the Petrified Forest Member. The ages of pre-Triassic grains also correspond to these chronostratigraphic units and are interpreted to reflect varying contributions from the Appalachian orogen to the east, Ouachita orogen to the southeast, Precambrian basement exposed in the ancestral Mogollon Highlands to the south, East Mexico arc, and Permian–Triassic arc built along the southern Cordilleran margin. Triassic grains in each chronostratigraphic unit also have distinct U and thorium (Th) concentrations, which are interpreted to reflect temporal changes in the chemistry of arc magmatism. Comparison of our LA-ICPMS ages with available chemical abrasion thermal ionization mass spectrometry (CA-TIMS) ages and new magnetostratigraphic data provides new insights into the depositional history of the Chinle Formation, as well as methods utilized to determine depositional ages of fluvial strata. For parts of the Chinle Formation that are dominated by fine-grained clastic strata (e.g., mudstone and siltstone), such as the Blue Mesa Member and Petrified Forest Member, all three chronometers agree (to within ∼ 1 Myr), and robust depositional chronologies have been determined. In contrast, for stratigraphic intervals dominated by coarse-grained clastic strata (e.g., sandstone), such as most of the Sonsela Member, the three chronologic records disagree due to recycling of older zircon grains and variable dilution of syn-depositional-age grains. This results in LA-ICPMS ages that significantly predate deposition and CA-TIMS ages that range between the other two chronometers. These complications challenge attempts to establish a well-defined chronostratigraphic age model for the Chinle Formation.

Description: Determining the sensitivity of the Greenland Ice Sheet (GrIS) to Holocene climate changes is a key prerequisite for understanding the future response of the ice sheet to global warming. In this study, we present new information on the Holocene glacial history of the GrIS in Inglefield Land, north Greenland. We use 10Be and in situ 14C exposure dating to constrain the timing of deglaciation in the area and radiocarbon dating of reworked molluscs and wood fragments to constrain when the ice sheet retreated behind its present-day extent. The 10Be ages are scattered ranging from ca. 92.7 to 6.8 ka, whereas the in situ 14C ages range from ca. 14.2 to 6.7 ka. Almost half of the apparent 10Be ages predate the Last Glacial Maximum and up to 89 % are to some degree affected by nuclide inheritance. Based on the few reliable 10Be ages, the in situ 14C ages and existing radiocarbon ages from Inglefield Land, we find that the deglaciation along the coast commenced at ca. 8.6–8.3 ka cal BP in the western part and ca. 7.9 ka in the central part, following the opening of Nares Strait and arrival of warm waters. The ice margin reached its present-day position at ca. 8.2 ka at the Humboldt Glacier and ca. 6.7 ka in the central part of Inglefield Land. Radiocarbon ages of reworked molluscs and wood fragments show that the ice margin was behind its present-day extent from ca. 5.8 to 0.5 ka cal BP. After 0.5 ka cal BP, the ice advanced towards its Little Ice Age position. Our results emphasize that the slowly eroding and possibly cold-based ice in north Greenland makes it difficult to constrain the deglaciation history based on 10Be ages alone unless they are paired with in situ 14C ages. Further, combining our findings with those of recently published studies reveals distinct differences between deglaciation patterns of northwest and north Greenland. Deglaciation of the land areas in northwest Greenland occurred earlier than in north Greenland, and periods of restricted ice extent were longer, spanning the Middle and Late Holocene. Overall, this highlights past ice sheet sensitivity to Holocene climate changes in an area where little information was available just a few years ago.

Description: This is the first of two papers that document the creation of the New European Wind Atlas (NEWA). It describes the sensitivity analysis and evaluation procedures that formed the basis for choosing the final setup of the mesoscale model simulations of the wind atlas. The suitable combination of model setup and parameterizations, bound by practical constraints, was found for simulating the climatology of the wind field at turbine-relevant heights with the Weather Research and Forecasting (WRF) model. Initial WRF model sensitivity experiments compared the wind climate generated by using two commonly used planetary boundary layer schemes and were carried out over several regions in Europe. They confirmed that the most significant differences in annual mean wind speed at 100 m a.g.l. (above ground level) mostly coincide with areas of high surface roughness length and not with the location of the domains or maximum wind speed. Then an ensemble of more than 50 simulations with different setups for a single year was carried out for one domain covering northern Europe for which tall mast observations were available. We varied many different parameters across the simulations, e.g. model version, forcing data, various physical parameterizations, and the size of the model domain. These simulations showed that although virtually every parameter change affects the results in some way, significant changes in the wind climate in the boundary layer are mostly due to using different physical parameterizations, especially the planetary boundary layer scheme, the representation of the land surface, and the prescribed surface roughness length. Also, the setup of the simulations, such as the integration length and the domain size, can considerably influence the results. We assessed the degree of similarity between winds simulated by the WRF ensemble members and the observations using a suite of metrics, including the Earth Mover's Distance (EMD), a statistic that measures the distance between two probability distributions. The EMD was used to diagnose the performance of each ensemble member using the full wind speed and direction distribution, which is essential for wind resource assessment. We identified the most realistic ensemble members to determine the most suitable configuration to be used in the final production run, which is fully described and evaluated in the second part of this study (Dörenkämper et al., 2020).

Description: This is the second of two papers that document the creation of the New European Wind Atlas (NEWA). In Part 1, we described the sensitivity experiments and accompanying evaluation done to arrive at the final mesoscale model setup used to produce the mesoscale wind atlas. In this paper, Part 2, we document how we made the final wind atlas product, covering both the production of the mesoscale climatology generated with the Weather Research and Forecasting (WRF) model and the microscale climatology generated with the Wind Atlas Analysis and Applications Program (WAsP). The paper includes a detailed description of the technical and practical aspects that went into running the mesoscale simulations and the downscaling using WAsP. We show the main results from the final wind atlas and present a comprehensive evaluation of each component of the NEWA model chain using observations from a large set of tall masts located all over Europe. The added value of the WRF and WAsP downscaling of wind climatologies is evaluated relative to the performance of the driving ERA5 reanalysis and shows that the WRF downscaling reduces the mean wind speed bias and spread relative to that of ERA5 from -1.50±1.30 to 0.02±0.78 m s−1. The WAsP downscaling has an added positive impact relative to that of the WRF model in simple terrain. In complex terrain, where the assumptions of the linearized flow model break down, both the mean bias and spread in wind speed are worse than those from the raw mesoscale results.

Description: We examine the relative contribution of processes controlling the interannual variability (IAV) of tropospheric ozone over four sub-regions of the southern hemispheric tropospheric ozone maximum (SHTOM) over a 20-year period. Our study is based on hindcast simulations from the National Aeronautics and Space Administration Global Modeling Initiative chemistry transport model (NASA GMI-CTM) of tropospheric and stratospheric chemistry, driven by assimilated Modern Era Retrospective Analysis for Research and Applications (MERRA) meteorological fields. Our analysis shows that over SHTOM region, the IAV of the stratospheric contribution is the most important factor driving the IAV of upper tropospheric ozone (270 hPa), where ozone has a strong radiative effect. Over the South Atlantic region, the contribution from surface emissions to the IAV of ozone exceeds that from stratospheric input at and below 430 hPa. Over the South Indian Ocean, the IAV of stratospheric ozone makes the largest contribution to the IAV of ozone with little or no influence from surface emissions at 270 and 430 hPa in austral winter. Over the tropical South Atlantic region, the contribution from IAV of stratospheric input dominates in austral winter at 270 hPa and drops to less than half but is still significant at 430 hPa. Emission contributions are not significant at these two levels. The IAV of lightning over this region also contributes to the IAV of ozone in September and December. Over the tropical southeastern Pacific, the contribution of the IAV of stratospheric input is significant at 270 and 430 hPa in austral winter, and emissions have little influence.

Description: Past studies have suggested that ozone in the troposphere has increased globally throughout much of the 20th century due to increases in anthropogenic emissions and transport. We show, by combining satellite measurements with a chemical transport model, that during the last four decades tropospheric ozone does indeed indicate increases that are global in nature, yet still highly regional. Satellite ozone measurements from Nimbus-7 and Earth Probe Total Ozone Mapping Spectrometer (TOMS) are merged with ozone measurements from the Aura Ozone Monitoring Instrument/Microwave Limb Sounder (OMI/MLS) to determine trends in tropospheric ozone for 1979–2016. Both TOMS (1979–2005) and OMI/MLS (2005–2016) depict large increases in tropospheric ozone from the Near East to India and East Asia and further eastward over the Pacific Ocean. The 38-year merged satellite record shows total net change over this region of about +6 to +7 Dobson units (DU) (i.e., ∼15 %–20 % of average background ozone), with the largest increase (∼4 DU) occurring during the 2005–2016 Aura period. The Global Modeling Initiative (GMI) chemical transport model with time-varying emissions is used to aid in the interpretation of tropospheric ozone trends for 1980–2016. The GMI simulation for the combined record also depicts the greatest increases of +6 to +7 DU over India and East Asia, very similar to the satellite measurements. In regions of significant increases in tropospheric column ozone (TCO) the trends are a factor of 2–2.5 larger for the Aura record when compared to the earlier TOMS record; for India and East Asia the trends in TCO for both GMI and satellite measurements are ∼+3 DU decade−1 or greater during 2005–2016 compared to about +1.2 to +1.4 DU decade−1 for 1979–2005. The GMI simulation and satellite data also reveal a tropospheric ozone increases in ∼+4 to +5 DU for the 38-year record over central Africa and the tropical Atlantic Ocean. Both the GMI simulation and satellite-measured tropospheric ozone during the latter Aura time period show increases of ∼+3 DU decade−1 over the N Atlantic and NE Pacific.

Description: The primary instrument on the Greenhouse gases Observing SATellite (GOSAT) is the Thermal And Near infrared Sensor for carbon Observations (TANSO) Fourier transform spectrometer (FTS). TANSO-FTS uses three short-wave infrared (SWIR) bands to retrieve total columns of CO2 and CH4 along its optical line of sight and one thermal infrared (TIR) channel to retrieve vertical profiles of CO2 and CH4 volume mixing ratios (VMRs) in the troposphere. We examine version 1 of the TANSO-FTS TIR CH4 product by comparing co-located CH4 VMR vertical profiles from two other remote-sensing FTS systems: the Canadian Space Agency's Atmospheric Chemistry Experiment FTS (ACE-FTS) on SCISAT (version 3.5) and the European Space Agency's Michelson Interferometer for Passive Atmospheric Sounding (MIPAS) on Envisat (ESA ML2PP version 6 and IMK-IAA reduced-resolution version V5R_CH4_224/225), as well as 16 ground stations with the Network for the Detection of Atmospheric Composition Change (NDACC). This work follows an initial inter-comparison study over the Arctic, which incorporated a ground-based FTS at the Polar Environment Atmospheric Research Laboratory (PEARL) at Eureka, Canada, and focuses on tropospheric and lower-stratospheric measurements made at middle and tropical latitudes between 2009 and 2013 (mid-2012 for MIPAS). For comparison, vertical profiles from all instruments are interpolated onto a common pressure grid, and smoothing is applied to ACE-FTS, MIPAS, and NDACC vertical profiles. Smoothing is needed to account for differences between the vertical resolution of each instrument and differences in the dependence on a priori profiles. The smoothing operators use the TANSO-FTS a priori and averaging kernels in all cases. We present zonally averaged mean CH4 differences between each instrument and TANSO-FTS with and without smoothing, and we examine their information content, their sensitive altitude range, their correlation, their a priori dependence, and the variability within each data set. Partial columns are calculated from the VMR vertical profiles, and their correlations are examined. We find that the TANSO-FTS vertical profiles agree with the ACE-FTS and both MIPAS retrievals' vertical profiles within 4 % (± ∼  40 ppbv) below 15 km when smoothing is applied to the profiles from instruments with finer vertical resolution but that the relative differences can increase to on the order of 25 % when no smoothing is applied. Computed partial columns are tightly correlated for each pair of data sets. We investigate whether the difference between TANSO-FTS and other CH4 VMR data products varies with latitude. Our study reveals a small dependence of around 0.1 % per 10 degrees latitude, with smaller differences over the tropics and greater differences towards the poles.