ALBERT

Description: The Arabian Sea (AS) hosts one of the most intense oxygen minimum zones (OMZs) in the world. Observations suggest a decline in O2 in the northern AS over the recent decades accompanied by an intensification of the suboxic conditions there. Over the same period, the local sea surface temperature has risen significantly, particularly over the Arabian Gulf (also known as Persian Gulf, hereafter the Gulf), while summer monsoon winds may have intensified. Here, we simulate the evolution of dissolved oxygen in the AS from 1982 through 2010 and explore its controlling factors, with a focus on changing atmospheric conditions. To this end, we use a set of eddy-resolving hindcast simulations forced with winds and heat and freshwater fluxes from an atmospheric reanalysis. We find a significant deoxygenation in the northern AS, with O2 inventories north of 20∘ N dropping by over 6 % per decade between 100 and 1000 m. These changes cause an expansion of the OMZ volume north of 20∘ N at a rate of 0.6 % per decade as well as an increase in the volume of suboxia and the rate of denitrification by 14 and 15 % per decade, respectively. We also show that strong interannual and decadal variability modulate dissolved oxygen in the northern AS, with most of the O2 decline taking place in the 1980s and 1990s. Using a set of sensitivity simulations we demonstrate that deoxygenation in the northern AS is essentially caused by reduced ventilation induced by the recent fast warming of the sea surface, including in the Gulf, with a contribution from concomitant summer monsoon wind intensification. This is because, on the one hand, surface warming enhances vertical stratification and increases Gulf water buoyancy, thus inhibiting vertical mixing and ventilation of the thermocline. On the other hand, summer monsoon wind intensification causes a rise in the thermocline depth in the northern AS that lowers O2 levels in the upper ocean. Our findings confirm that the AS OMZ is strongly sensitive to upper-ocean warming and concurrent changes in the Indian monsoon winds. Finally, our results also demonstrate that changes in the local climatic forcing play a key role in regional dissolved oxygen changes and hence need to be properly represented in global models to reduce uncertainties in future projections of deoxygenation.

Description: McMurdo Sound sea ice can generally be partitioned into two regimes: (1) a stable fast-ice cover, forming south of approximately 77.6∘ S around March–April and then breaking out the following January–February, and (2) a more dynamic region north of 77.6∘ S that the McMurdo Sound and Ross Sea polynyas regularly impact. In 2019, a stable fast-ice cover formed unusually late due to repeated break-out events. We analyse the 2019 sea-ice conditions and relate them to a modified storm index (MSI), a proxy for southerly wind events. We find there is a strong correlation between the timing of break-out events and several unusually large MSI events.

Description: The global distribution of cropping intensity (CI) is essential to our understanding of agricultural land use management on Earth. Optical remote sensing has revolutionized our ability to map CI over large areas in a repeated and cost-efficient manner. Previous studies have mainly focused on investigating the spatiotemporal patterns of CI ranging from regions to the entire globe with the use of coarse-resolution data, which are inadequate for characterizing farming practices within heterogeneous landscapes. To fill this knowledge gap, in this study, we utilized multiple satellite data to develop a global, spatially continuous CI map dataset at 30 m resolution (GCI30). Accuracy assessments indicated that GCI30 exhibited high agreement with visually interpreted validation samples and in situ observations from the PhenoCam network. We carried out both statistical and spatial comparisons of GCI30 with six existing global CI estimates. Based on GCI30, we estimated that the global average annual CI during 2016–2018 was 1.05, which is close to the mean (1.09) and median (1.07) CI values of the existing six global CI estimates, although the spatial resolution and temporal coverage vary significantly among products. A spatial comparison with two satellite-based land surface phenology products further suggested that GCI30 was not only capable of capturing the overall pattern of global CI but also provided many spatial details. GCI30 indicated that single cropping was the primary agricultural system on Earth, accounting for 81.57 % (12.28×106 km2) of the world's cropland extent. Multiple-cropping systems, on the other hand, were commonly observed in South America and Asia. We found large variations across countries and agroecological zones, reflecting the joint control of natural and anthropogenic drivers on regulating cropping practices. As the first global-coverage, fine-resolution CI product, GCI30 is expected to fill the data gap for promoting sustainable agriculture by depicting worldwide diversity of agricultural land use intensity. The GCI30 dataset is available on Harvard Dataverse: https://doi.org/10.7910/DVN/86M4PO (Zhang et al., 2020).

Description: This article summarises the results of an analysis of solar radio bursts (SRBs) detected by the Compound Astronomical Low-cost Low-frequency Instrument for Spectroscopy and Transportable Observatory (CALLISTO) spectrometer hosted by the University of Rwanda. The data analysed were detected during the first year (2014–2015) of the instrument operation. Using quick plots provided by the e-CALLISTO website, a total of 201 intense and well-separated solar radio bursts detected by the CALLISTO station located in Rwanda, are found consisting of 4 type II, 175 type III and 22 type IV radio bursts. It is found that all analysed type II and ∼ 37 % of type III bursts are associated with impulsive solar flares, while the minority (∼ 13 %) of type IV radio bursts are associated with solar flares. Furthermore, all type II radio bursts are associated with coronal mass ejections (CMEs), ∼ 44 % of type III bursts are associated with CMEs, and the majority (∼ 82 %) of type IV bursts were accompanied by CMEs. With aid of the atmospheric imaging assembly (AIA) images on board the Solar Dynamics Observatory (SDO), the location of open magnetic field lines of non-flare-associated type III radio bursts are shown. The same images are used to show the magnetic loops in the solar corona for type IV radio bursts observed in the absence of solar flares and/or CMEs. Findings from this study indicate that analysis of SRBs that are observed from the ground can provide a significant contribution to the early diagnosis of solar transients phenomena, such as solar flares and CMEs, which are major drivers of potential space weather hazards.

Description: Two feature-based verification methods, thus far only used for the diagnostic evaluation of atmospheric models, have been applied to compare ∼7 km resolution pre-operational analyses of chlorophyll-a (Chl-a) concentrations to a 1 km gridded satellite-derived Chl-a concentration product. The aim of this study was to assess the value of applying such methods to ocean models. Chl-a bloom objects were identified in both data sets for the 2019 bloom season (1 March to 31 July). These bloom objects were analysed as discrete (2-D) spatial features, but also as space–time (3-D) features, providing the means of defining the onset, duration and demise of distinct bloom episodes and the season as a whole. The new feature-based verification methods help reveal that the model analyses are not able to represent small coastal bloom objects, given the coarser definition of the coastline, also wrongly producing more bloom objects in deeper Atlantic waters. Model analyses' concentrations are somewhat higher overall. The bias manifests itself in the size of the model analysis bloom objects, which tend to be larger than the satellite-derived bloom objects. The onset of the bloom season is delayed by 26 d in the model analyses, but the season also persists for another month beyond the diagnosed end. The season was diagnosed to be 119 d long in the model analyses, compared to 117 d from the satellite product. Geographically, the model analyses and satellite-derived bloom objects do not necessarily exist in a specific location at the same time and only overlap occasionally.

Description: The source apportionment of aerosol iron (Fe), including natural and combustion Fe, is an important issue because aerosol Fe can enhance oceanic primary production in the surface ocean. Based on our previous finding that combustion Fe emitted by evaporation processes has Fe isotope ratios (δ56Fe) that are approximately 4 ‰ lower than those of natural Fe, this study aimed to distinguish aerosol Fe sources over the northwestern Pacific using two size-fractionated marine aerosols. The δ56Fe values of fine and coarse particles from the eastern or northern Pacific were found to be similar to each other, ranging from 0.0 ‰ to 0.4 ‰. Most of them were close to the crustal average, suggesting the dominance of natural Fe. On the other hand, particles from the direction of East Asia demonstrated lower δ56Fe values in fine particles (−0.5 ‰ to −2.2 ‰) than in coarse particles (on average −0.02 ± 0.12 ‰). The correlations between the δ56Fe values and the enrichment factors of lead and vanadium suggested that the low δ56Fe values obtained were due to the presence of combustion Fe. The δ56Fe values of the soluble component of fine particles in this region were lower than the total, indicating the preferential dissolution of combustion Fe. In addition, we found a negative correlation between the δ56Fe value and the fractional Fe solubility in air masses from the direction of East Asia. These results suggest that the presence of combustion Fe is an important factor in controlling the fractional Fe solubility in air masses from the direction of East Asia, whereas other factors are more important in the other areas. By assuming typical δ56Fe values for combustion and natural Fe, the contribution of combustion Fe to the total (acid-digested) Fe in aerosols was estimated to reach up to 50 % of fine and 21 % of bulk (coarse + fine) particles in air masses from the direction of East Asia, whereas its contribution was small in the other areas. The contribution of combustion Fe to the soluble Fe component estimated for one sample was approximately twice as large as the total, indicating the importance of combustion Fe as a soluble Fe source despite lower emissions than the natural. These isotope-based estimates were compared with those estimated using an atmospheric chemical transport model (IMPACT), in which the fractions of combustion Fe in fine particles, especially in air masses from the direction of East Asia, were consistent with each other. In contrast, the model estimated a relatively large contribution from combustion Fe in coarse particles, probably because of the different characteristics of combustion Fe that are included in the model calculation and the isotope-based estimation. This highlights the importance of observational data on δ56Fe for size-fractionated aerosols to scale the combustion Fe emission by the model. The average deposition fluxes of soluble Fe to the surface ocean were 1.4 and 2.9 nmol m−2 d−1 from combustion and natural aerosols, respectively, in air masses from the direction of East Asia, which suggests that combustion Fe could be an important Fe source to the surface seawater among other Fe sources. Distinguishing Fe sources using the δ56Fe values of marine aerosols and seawater is anticipated to lead to a more quantitative understanding of the Fe cycle in the atmosphere and surface ocean.

Description: The importance of Antarctic sea ice and Southern Ocean warming has come into the focus of polar research during the last couple of decades. Especially around West Antarctica, where warm water masses approach the continent and where sea ice has declined, the distribution and evolution of sea ice play a critical role in the stability of nearby ice shelves. Organic geochemical analyses of marine seafloor surface sediments from the Antarctic continental margin allow an evaluation of the applicability of biomarker-based sea-ice and ocean temperature reconstructions in these climate-sensitive areas. We analysed highly branched isoprenoids (HBIs), such as the sea-ice proxy IPSO25 and phytoplankton-derived HBI-trienes, as well as phytosterols and isoprenoidal glycerol dialkyl glycerol tetraethers (GDGTs), which are established tools for the assessment of primary productivity and ocean temperatures respectively. The combination of IPSO25 with a phytoplankton marker (i.e. the PIPSO25 index) permits semi-quantitative sea-ice reconstructions and avoids misleading over- or underestimations of sea-ice cover. Comparisons of the PIPSO25-based sea-ice distribution patterns and TEX86L- and RI-OH′-derived ocean temperatures with (1) sea-ice concentrations obtained from satellite observations and (2) instrument measurements of sea surface and subsurface temperatures corroborate the general capability of these proxies to determine oceanic key variables properly. This is further supported by model data. We also highlight specific aspects and limitations that need to be taken into account for the interpretation of such biomarker data and discuss the potential of IPSO25 as an indicator for the former occurrence of platelet ice and/or the export of ice-shelf water.

Description: Forest fires modify soil organic carbon and suppress soil respiration for many decades after the initial disturbance. The associated changes in soil autotrophic and heterotrophic respiration from the time of the forest fire, however, are less well characterized. The FireResp model predicts soil autotrophic and heterotrophic respiration parameterized with a novel dataset across a fire chronosequence in the Yukon and Northwest Territories of Canada. The dataset consisted of soil incubation experiments and field measurements of soil respiration and soil carbon stocks. The FireResp model contains submodels that consider a Q10 (exponential) model of respiration compared to models of heterotrophic respiration using Michaelis–Menten kinetics parameterized with soil microbial carbon. For model evaluation we applied the Akaike information criterion and compared predicted patterns in components of soil respiration across the chronosequence. Parameters estimated with data from the 5 cm soil depth had better model–data comparisons than parameters estimated with data from the 10 cm soil depth. The model–data fit was improved by including parameters estimated from soil incubation experiments. Models that incorporated microbial carbon with Michaelis–Menten kinetics reproduced patterns in autotrophic and heterotrophic soil respiration components across the chronosequence. Autotrophic respiration was associated with aboveground tree biomass at more recently burned sites, but this association was less robust at older sites in the chronosequence. Our results provide support for more structured soil respiration models than standard Q10 exponential models.

Description: A new luminescence erosion meter has huge potential for inferring erosion rates on sub-millennial scales for both steady and transient states of erosion, which is not currently possible with any existing techniques capable of measuring erosion. This study applies new rock luminescence techniques to a well-constrained scenario provided by the Beinn Alligin rock avalanche, NW Scotland. Boulders in this deposit are lithologically consistent and have known cosmogenic nuclide ages and independently derived Holocene erosion rates. We find that luminescence-derived exposure ages for the Beinn Alligin rock avalanche were an order of magnitude younger than existing cosmogenic nuclide exposure ages, suggestive of high erosion rates (as supported by field evidence of quartz grain protrusions on the rock surfaces). Erosion rates determined by luminescence were consistent with independently derived rates measured from boulder edge roundness. Inversion modelling indicates a transient state of erosion reflecting the stochastic nature of erosional processes over the last ∼4.5 kyr in the wet, temperate climate of NW Scotland. Erosion was likely modulated by known fluctuations in moisture availability and to a lesser extent temperature, which controlled the extent of chemical weathering of these highly lithified rocks prior to erosion. The use of a multi-elevated temperature, post-infra-red, infra-red stimulated luminescence (MET-pIRIR) protocol (50, 150 and 225 ∘C) was advantageous as it identified samples with complexities that would not have been observed using only the standard infra-red stimulated luminescence (IRSL) signal measured at 50 ∘C, such as that introduced by within-sample variability (e.g. surficial coatings). This study demonstrates that the luminescence erosion meter can infer accurate erosion rates on sub-millennial scales and identify transient states of erosion (i.e. stochastic processes) in agreement with independently derived erosion rates for the same deposit.

Description: The photo-oxidation of myrcene, a monoterpene species emitted by plants, was investigated at atmospheric conditions in the outdoor simulation chamber SAPHIR (Simulation of Atmospheric PHotochemistry In a Large Reaction Chamber). The chemical structure of myrcene consists of one moiety that is a conjugated π system (similar to isoprene) and another moiety that is a triple-substituted olefinic unit (similar to 2-methyl-2-butene). Hydrogen shift reactions of organic peroxy radicals (RO2) formed in the reaction of isoprene with atmospheric OH radicals are known to be of importance for the regeneration of OH. Structure–activity relationships (SARs) suggest that similar hydrogen shift reactions like in isoprene may apply to the isoprenyl part of RO2 radicals formed during the OH oxidation of myrcene. In addition, SAR predicts further isomerization reactions that would be competitive with bimolecular RO2 reactions for chemical conditions that are typical for forested environments with low concentrations of nitric oxide. Assuming that OH peroxy radicals can rapidly interconvert by addition and elimination of O2 like in isoprene, bulk isomerization rate constants of 0.21 and 0.097 s−1 (T=298 K) for the three isomers resulting from the 3′-OH and 1-OH addition, respectively, can be derived from SAR. Measurements of radicals and trace gases in the experiments allowed us to calculate radical production and destruction rates, which are expected to be balanced. The largest discrepancies between production and destruction rates were found for RO2. Additional loss of organic peroxy radicals due to isomerization reactions could explain the observed discrepancies. The uncertainty of the total radical (ROx=OH+HO2+RO2) production rates was high due to the uncertainty in the yield of radicals from myrcene ozonolysis. However, results indicate that radical production can only be balanced if the reaction rate constant of the reaction between hydroperoxy (HO2) and RO2 radicals derived from myrcene is lower (0.9 to 1.6×10-11 cm3 s−1) than predicted by SAR. Another explanation of the discrepancies would be that a significant fraction of products (yield: 0.3 to 0.6) from these reactions include OH and HO2 radicals instead of radical-terminating organic peroxides. Experiments also allowed us to determine the yields of organic oxidation products acetone (yield: 0.45±0.08) and formaldehyde (yield: 0.35±0.08). Acetone and formaldehyde are produced from different oxidation pathways, so that yields of these compounds reflect the branching ratios of the initial OH addition to myrcene. Yields determined in the experiments are consistent with branching ratios expected from SAR. The yield of organic nitrate was determined from the gas-phase budget analysis of reactive oxidized nitrogen in the chamber, giving a value of 0.13±0.03. In addition, the reaction rate constant for myrcene + OH was determined from the measured myrcene concentration, yielding a value of (2.3±0.3)×10-10 cm3 s−1.