Seasonal Prediction Potential for Springtime Dustiness in the United States

Most dust forecast models focus on short, subseasonal lead times, that is, 3 to 6 days, and the skill of seasonal prediction is not clear. In this study we examine the potential of seasonal dust prediction in the United States using an observation‐constrained regression model and key variables predicted by a seasonal prediction model developed at National Oceanic and Atmospheric Administration Geophysical Fluid Dynamics Laboratory, the Forecast‐Oriented Low Ocean Resolution (FLOR) model. Our method shows skillful predictions of spring dustiness 3 to 6 months in advance. It is found that the regression model explains about 71% of the variances of dust event frequency over the Great Plains and 63% over the southwestern United States in March‐May from 2004 to 2016 using predictors from FLOR initialized on 1 December. Variations in springtime dustiness are dominated by springtime climatic factors rather than wintertime factors. Findings here will help development of a seasonal dust prediction system and hazard prevention.NASA (NNH14ZDA001N-ACMAP, NNH16ZDA001N-MAP)Princeton University's Cooperative Institute for Climate ScienceNOA

Aerosols, Chemistry, and Radiative Forcing: A 3-D Model Analysis of Satellite and ACE-Asia data (ACMAP)

We propose a research project to incorporate a global 3-D model and satellite data into the multi-national Aerosol Characterization Experiment-Asia (ACE-Asia) mission. Our objectives are (1) to understand the physical, chemical, and optical properties of aerosols and the processes that control those properties over the Asian-Pacific region, (2) to investigate the interaction between aerosols and tropospheric chemistry, and (3) to determine the aerosol radiative forcing over the Asia-Pacific region. We will use the Georgia TecWGoddard Global Ozone Chemistry Aerosol Radiation and Transport (GOCART) model to link satellite observations and the ACE-Asia measurements. First, we will use the GOCART model to simulate aerosols and related species, and evaluate the model with satellite and in-situ observations. Second, the model generated aerosol vertical profiles and compositions will be used to validate the satellite products; and the satellite data will be used for during- and post- mission analysis. Third, we will use the model to analyze and interpret both satellite and ACE- Asia field campaign data and investigate the aerosol-chemistry interactions. Finally, we will calculate aerosol radiative forcing over the Asian-Pacific region, and assess the influence of Asian pollution in the global atmosphere. We propose a research project to incorporate a global 3-D model and satellite data int

Do MODIS-defined dust sources have a geomorphological signature?

The preferential dust source (PDS) scheme enables large-scale mapping of geomorphology in terms of importance for dust emissions but has not been independently tested other than at local scales. We examine the PDS qualitative conceptual model of surface emissivity alongside a quantitative measurement of dust loading from Moderate Resolution Imaging Spectroradiometer (MODIS) Deep Blue Collection 6 for the Chihuahuan Desert. The predicted ranked importance of each geomorphic type for dust emissions is compared with the actual ranked importance as determined from the satellite-derived dust loading. For this region, the predicted variability and magnitude of dust emissions from most surface types present coincides with the observed characteristics demonstrating the significance of geomorphological controls on emission. The exception is for areas of low magnitude but persistent emissions such as alluvial surfaces where PDS overpredicts dustiness. As PDS is a good predictor of emissions and incorporates surface dynamics it could improve models of future dust emissions

Impact of preindustrial to present-day changes in short-lived pollutant emissions on atmospheric composition and climate forcing

We describe and evaluate atmospheric chemistry in the newly developed Geophysical Fluid Dynamics Laboratory chemistry-climate model (GFDL AM3) and apply it to investigate the net impact of preindustrial (PI) to present (PD) changes in short-lived pollutant emissions (ozone precursors, sulfur dioxide, and carbonaceous aerosols) and methane concentration on atmospheric composition and climate forcing. The inclusion of online troposphere-stratosphere interactions, gas-aerosol chemistry, and aerosol-cloud interactions (including direct and indirect aerosol radiative effects) in AM3 enables a more complete representation of interactions among short-lived species, and thus their net climate impact, than was considered in previous climate assessments. The base AM3 simulation, driven with observed sea surface temperature (SST) and sea ice cover (SIC) over the period 1981–2007, generally reproduces the observed mean magnitude, spatial distribution, and seasonal cycle of tropospheric ozone and carbon monoxide. The global mean aerosol optical depth in our base simulation is within 5% of satellite measurements over the 1982–2006 time period. We conduct a pair of simulations in which only the short-lived pollutant emissions and methane concentrations are changed from PI (1860) to PD (2000) levels (i.e., SST, SIC, greenhouse gases, and ozone-depleting substances are held at PD levels). From the PI to PD, we find that changes in short-lived pollutant emissions and methane have caused the tropospheric ozone burden to increase by 39% and the global burdens of sulfate, black carbon, and organic carbon to increase by factors of 3, 2.4, and 1.4, respectively. Tropospheric hydroxyl concentration decreases by 7%, showing that increases in OH sinks (methane, carbon monoxide, nonmethane volatile organic compounds, and sulfur dioxide) dominate over sources (ozone and nitrogen oxides) in the model. Combined changes in tropospheric ozone and aerosols cause a net negative top-of-the-atmosphere radiative forcing perturbation (−1.05 W m−2) indicating that the negative forcing (direct plus indirect) from aerosol changes dominates over the positive forcing due to ozone increases, thus masking nearly half of the PI to PD positive forcing from long-lived greenhouse gases globally, consistent with other current generation chemistry-climate models

Why Is Improvement of Earth System Models so Elusive? Challenges and Strategies from Dust Aerosol Modeling

Past decades have seen an accelerating increase in computing efficiency, while climate models are representing a rapidly widening set of physical processes. Yet simulations of some fundamental aspects of climate like precipitation or aerosol forcing remain highly uncertain and resistant to progress. Dust aerosol modeling of soil particles lofted by wind erosion has seen a similar conflict between increasing model sophistication and remaining uncertainty. Dust aerosols perturb the energy and water cycles by scattering radiation and acting as ice nuclei, while mediating atmospheric chemistry and marine photosynthesis (and thus the carbon cycle). These effects take place across scales from the dimensions of an ice crystal to the planetary-scale circulation that disperses dust far downwind of its parent soil. Representing this range leads to several modeling challenges. Should we limit complexity in our model, which consumes computer resources and inhibits interpretation? How do we decide if a process involving dust is worthy of inclusion within our model? Can we identify a minimal representation of a complex process that is efficient yet retains the physics relevant to climate? Answering these questions about the appropriate degree of representation is guided by model evaluation, which presents several more challenges. How do we proceed if the available observations do not directly constrain our process of interest? (This could result from competing processes that influence the observed variable and obscure the signature of our process of interest.) Examples will be presented from dust modeling, with lessons that might be more broadly applicable. The end result will either be clinical depression or there assuring promise of continued gainful employment as the community confronts these challenges

How Well Does NASA GEOS Model Perform in Simulating Dust Deposition into the Tropical Atlantic Ocean?

Massive dust emitted from North Africa can transport long distances across the tropical Atlantic Ocean, reaching the Americas. Dust deposition along the transit adds microorganisms and essential nutrients to marine ecosystem, which has important implications for biogeochemical cycle and climate. However, assessing the dust-ecosystemclimate interactions has been hindered in part by the paucity of dust deposition measurements and large uncertainties associated with oversimplified representations of dust processes in current models. We have recently produced a unique dataset of seasonal dust deposition flux and dust loss frequency into the tropical Atlantic Ocean at a nominal resolution of 200 km x 500 km by using the decade-long (2007-2016) record of aerosol three-dimensional distribution from four satellite sensors, namely CALIOP, MODIS, MISR, and IASI. On the basis of the ten-year average, the yearly dust deposition into the tropical Atlantic Ocean is estimated at 98-153 Tg. The dust deposition shows large spatial and temporal (on seasonal and interannual scale) variability. The satellite observations also yield an estimate of annual mean dust loss frequency of 0.052 ~ 0.078 d-1, a useful diagnostic that makes it possible to disentangle the dust transport and removal processes from the dust emissions when identifying the major factors contributing to the uncertainties and biases in the model simulated dust deposition. In this study, we use the dataset along with in situ and remote sensing observations to assess how well NASA GEOS model performs in simulating trans-Atlantic dust transport and deposition. We found that the GEOS modeling of dust deposition falls within the range of satellite-based estimates. However, this reasonable agreement in dust deposition is a compensation of the model's underestimate of dust emissions and overestimate of dust removal efficiency. Further, the overestimate of dust removal efficiency results largely from the model's overestimate of rainfall rate. Our results provide insights into the model's deficiencies at process level, which could better guide model improvements

Interpretation of TOMS observations of tropical tropospheric ozone with a global model and in-situ observations

We interpret the distribution of tropical tropospheric ozone columns (TTOCs) from the Total Ozone Mapping Spectrometer (TOMS) by using a global three-dimensional model of tropospheric chemistry (GEOS-CHEM) and additional information from in situ observations. The GEOS-CHEM TTOCs capture 44% of the variance of monthly mean TOMS TTOCs from the convective cloud differential method (CCD) with no global bias. Major discrepancies are found over northern Africa and south Asia where the TOMS TTOCs do not capture the seasonal enhancements from biomass burning found in the model and in aircraft observations. A characteristic feature of these northern tropical enhancements, in contrast to southern tropical enhancements, is that they are driven by the lower troposphere where the sensitivity of TOMS is poor due to Rayleigh scattering. We develop an efficiency correction to the TOMS retrieval algorithm that accounts for the variability of ozone in the lower troposphere. This efficiency correction increases TTOCs over biomass burning regions by 3–5 Dobson units (DU) and decreases them by 2–5 DU over oceanic regions, improving the agreement between CCD TTOCs and in situ observations. Applying the correction to CCD TTOCs reduces by ∼5 DU the magnitude of the “tropical Atlantic paradox” [Thompson et al., 2000], i.e. the presence of a TTOC enhancement over the southern tropical Atlantic during the northern African biomass burning season in December–February. We reproduce the remainder of the paradox in the model and explain it by the combination of upper tropospheric ozone production from lightning NOx, persistent subsidence over the southern tropical Atlantic as part of the Walker circulation, and cross-equatorial transport of upper tropospheric ozone from northern midlatitudes in the African “westerly duct.” These processes in the model can also account for the observed 13–17 DU persistent wave-1 pattern in TTOCs with a maximum over the tropical Atlantic and a minimum over the tropical Pacific during all seasons. The photochemical effects of mineral dust have only a minor role on the modeled distribution of TTOCs, including over northern Africa, due to multiple competing effects. The photochemical effects of mineral dust globally decrease annual mean OH concentrations by 9%. A global lightning NOx source of 6 Tg N yr−1 in the model produces a simulation that is most consistent with TOMS and in situ observations

Retrieving the global distribution of the threshold of wind erosion from satellite data and implementing it into the Geophysical Fluid Dynamics Laboratory land–atmosphere model (GFDL AM4.0/LM4.0)

Dust emission is initiated when surface wind velocities exceed the threshold of wind erosion. Many dust models used constant threshold values globally. Here we use satellite products to characterize the frequency of dust events and land surface properties. By matching this frequency derived from Moderate Resolution Imaging Spectroradiometer (MODIS) Deep Blue aerosol products with surface winds, we are able to retrieve a climatological monthly global distribution of the wind erosion threshold (Vthreshold) over dry and sparsely vegetated surfaces. This monthly two-dimensional threshold velocity is then implemented into the Geophysical Fluid Dynamics Laboratory coupled land–atmosphere model (AM4.0/LM4.0). It is found that the climatology of dust optical depth (DOD) and total aerosol optical depth, surface PM10 dust concentrations, and the seasonal cycle of DOD are better captured over the “dust belt” (i.e., northern Africa and the Middle East) by simulations with the new wind erosion threshold than those using the default globally constant threshold. The most significant improvement is the frequency distribution of dust events, which is generally ignored in model evaluation. By using monthly rather than annual mean Vthreshold, all comparisons with observations are further improved. The monthly global threshold of wind erosion can be retrieved under different spatial resolutions to match the resolution of dust models and thus can help improve the simulations of dust climatology and seasonal cycles as well as dust forecasting

The MONARCH high-resolution reanalysis of desert dust aerosol over Northern Africa, the Middle East and Europe (2007-2016)

One of the challenges in studying desert dust aerosol along with its numerous interactions and impacts is the paucity of direct in-situ measurements, particularly in the areas most affected by dust storms. Satellites typically provide columnintegrated aerosol measurements, but observationally-constrained continuous 3D dust fields are needed to assess dust variability, climate effects and impacts upon a variety of socio-economic sectors. Here, we present a high resolution regional reanalysis data set of desert dust aerosols that covers Northern Africa, the Middle East and Europe along with the Mediterranean sea and parts of Central Asia, and the Atlantic and Indian Oceans between 2007 and 2016. The horizontal resolution is 0.1° latitude × 0.1° longitude, and the temporal resolution is 3 hours. The reanalysis was produced using Local Ensemble Transform Kalman Filter (LETKF) data assimilation in the Multiscale Online Non-hydrostatic AtmospheRe CHemistry model (MONARCH) developed at the Barcelona Supercomputing Center (BSC). The assimilated data are coarse-mode dust optical depth retrieved from the Moderate Resolution Imaging Spectroradiometer (MODIS) Deep Blue Level 2 products. The reanalysis data set consists of upper air (dust mass concentrations and extinction coefficient), surface (dust deposition and solar irradiance fields, among them) and total column (e.g., dust optical depth and load) variables. Some dust variables, such as concentrations and wet and dry deposition, are expressed for a binned size distribution that ranges from 0.2 to 20 μm in particle diameter. Both analysis and first-guess (analysis-initialized simulation) fields are available for the variables that are diagnosed from the state vector.We acknowledge the DustClim project which is part of ERA4CS, an ERA-NET initiated by JPI Climate, and funded by FORMAS (SE), DLR (DE), BMWFW (AT), IFD (DK), MINECO (ES), ANR (FR) with co-funding by the European Union’s Horizon 2020 research and innovation programme (Grant n. 690462)