ALBERT

Description: The observing system design of multidisciplinary field measurements involves a variety of considerations on logistics, safety, and science objectives. Typically, this is done based on investigator intuition and designs of prior field measurements. However, there is potential for considerable increases in efficiency, safety, and scientific success by integrating numerical simulations in the design process. Here, we present a novel numerical simulation–environmental response function (NS–ERF) approach to observing system simulation experiments that aids surface–atmosphere synthesis at the interface of mesoscale and microscale meteorology. In a case study we demonstrate application of the NS–ERF approach to optimize the Chequamegon Heterogeneous Ecosystem Energy-balance Study Enabled by a High-density Extensive Array of Detectors 2019 (CHEESEHEAD19). During CHEESEHEAD19 pre-field simulation experiments, we considered the placement of 20 eddy covariance flux towers, operations for 72 h of low-altitude flux aircraft measurements, and integration of various remote sensing data products. A 2 h high-resolution large eddy simulation created a cloud-free virtual atmosphere for surface and meteorological conditions characteristic of the field campaign domain and period. To explore two specific design hypotheses we super-sampled this virtual atmosphere as observed by 13 different yet simultaneous observing system designs consisting of virtual ground, airborne, and satellite observations. We then analyzed these virtual observations through ERFs to yield an optimal aircraft flight strategy for augmenting a stratified random flux tower network in combination with satellite retrievals. We demonstrate how the novel NS–ERF approach doubled CHEESEHEAD19's potential to explore energy balance closure and spatial patterning science objectives while substantially simplifying logistics. Owing to its modular extensibility, NS–ERF lends itself to optimizing observing system designs also for natural climate solutions, emission inventory validation, urban air quality, industry leak detection, and multi-species applications, among other use cases.

Description: This paper presents a new technique to derive thermospheric temperature from space-based disk observations of far ultraviolet airglow. The technique, guided by findings from principal component analysis of synthetic daytime Lyman–Birge–Hopfield (LBH) disk emissions, uses a ratio of the emissions in two spectral channels that together span the LBH (2,0) band to determine the change in band shape with respect to a change in the rotational temperature of N2. The two-channel-ratio approach limits representativeness and measurement error by only requiring measurement of the relative magnitudes between two spectral channels and not radiometrically calibrated intensities, simplifying the forward model from a full radiative transfer model to a vibrational–rotational band model. It is shown that the derived temperature should be interpreted as a column-integrated property as opposed to a temperature at a specified altitude without utilization of a priori information of the thermospheric temperature profile. The two-channel-ratio approach is demonstrated using NASA GOLD Level 1C disk emission data for the period of 2–8 November 2018 during which a moderate geomagnetic storm has occurred. Due to the lack of independent thermospheric temperature observations, the efficacy of the approach is validated through comparisons of the column-integrated temperature derived from GOLD Level 1C data with the GOLD Level 2 temperature product as well as temperatures from first principle and empirical models. The storm-time thermospheric response manifested in the column-integrated temperature is also shown to corroborate well with hemispherically integrated Joule heating rates, ESA SWARM mass density at 460 km, and GOLD Level 2 column O/N2 ratio.

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: Non-Gaussian forecast error is a challenge for ensemble-based data assimilation (DA), particularly for more nonlinear convective dynamics. In this study, we investigate the degree of the non-Gaussianity of forecast error distributions at 1 km resolution using a 1000-member ensemble Kalman filter, and how it is affected by the DA update frequency and observation number. Regional numerical weather prediction experiments are performed with the SCALE (Scalable Computing for Advanced Library and Environment) model and the LETKF (local ensemble transform Kalman filter) assimilating phased array radar observations every 30 s. The results show that non-Gaussianity develops rapidly within convective clouds and is sensitive to the DA frequency and the number of assimilated observations. The non-Gaussianity is reduced by up to 40 % when the assimilation window is shortened from 5 min to 30 s, particularly for vertical velocity and radar reflectivity.

Description: China has experienced dramatic changes in emissions since 2010, which accelerated following the implementation of the Clean Air Action program in 2013. These changes have resulted in significant air quality improvements that are reflected in observations from both surface networks and satellite observations. Air pollutants, such as PM2.5, surface ozone, and their precursors, have long enough lifetimes in the troposphere to be easily transported downwind. Emission changes in China will thus not only change the domestic air quality but will also affect the air quality in other regions. In this study, we use a global chemistry transport model (CAM-chem) to simulate the influence of Chinese emission changes from 2010 to 2017 on both domestic and foreign air quality. We then quantify the changes in air-pollution-associated (including both PM2.5 and O3) premature mortality burdens at regional and global scales. Within our simulation period, the population-weighted annual PM2.5 concentration in China peaks in 2011 (94.1 µg m−3) and decreases to 69.8 µg m−3 by 2017. These estimated national PM2.5 concentration changes in China are comparable with previous studies using fine-resolution regional models, though our model tends to overestimate PM2.5 from 2013 to 2017 when evaluated with surface observations. Relative to 2010, emission changes in China increased the global PM2.5-associated premature mortality burdens through 2013, among which a majority of the changes (∼ 93 %) occurred in China. The sharp emission decreases after 2013 generated significant benefits for human health. By 2017, emission changes in China reduced premature deaths associated with PM2.5 by 108 800 (92 800–124 800) deaths per year globally, relative to 2010, among which 92 % were realized in China. In contrast, the population-weighted, annually averaged maximum daily 8 h ozone concentration peaked in 2014 and did not reach 2010 levels by 2017. As such, O3 generated nearly 8500 (6500–9900) more premature deaths per year in 2017 compared to 2010. Downwind regions, such as South Korea, Japan, and the United States, generally experienced O3 improvements following 2013 due to the decreased export of ozone and its precursors. Overall, we conclude that the sharp emission reductions in China over the past decade have generated substantial benefits for air quality that have reduced premature deaths associated with air pollution at a global scale.

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: 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 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: 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: Aeromagnetic exploration is an important method of geophysical exploration. We study the compensation method of the towed bird system and establish the towed bird interference model. Due to the geomagnetic gradient changing greatly, the geomagnetic gradient is considered in the towed bird interference model. In this paper, we model the geomagnetic field gradient and analyze the influence of the towed bird system on the aeromagnetic compensation results. Finally, we apply the ridge regression method to solve the problem. We verify the feasibility of this compensation method through actual flight tests and further improve the data quality of the towed bird interference.