Search NASA⌕ Search

Engineering topics

Peter R Colarco

Publications and source records attributed to Peter R Colarco.

At least 19 records

Benchmarking GOCART-2G in the Goddard Earth Observing System (GEOS)

The Goddard Chemistry Aerosol Radiation and Transport (GOCART) model, which controls the sources sinks and chemistry within the Goddard Earth Observing System, recently underwent a major refactoring and update to the representation of physical processes. This paper serves to document code changes that were included in GOCART 2nd Generation (GOCART-2G) and establishes a benchmark simulation that is to be used for future development of the system. The code refactoring increases flexibility such multiple instances of an aerosol species can be run and interact with radiation and cloud microphysics, in addition to the output of multiple wavelength aerosol optical properties in support of data assimilation. From a science perspective, a new radiatively active tracer, brown carbon, was added to distinguish smoke from other sources of organic aerosol thereby improving optical properties entering the radiative calculations. A four-year benchmark simulation was evaluated using in situ and space borne measurements to develop a baseline and prioritize future development. A comparison of simulated aerosol optical depth between GOCART-2G and MODIS retrievals indicates the model captures the overall spatial pattern and seasonal cycle of aerosol optical depth but overestimates aerosol extinction over dusty regions and underestimates aerosol extinction over northern hemisphere boreal forests, requiring further tuning of emissions. This MODIS-based analysis is corroborated by comparisons to MISR and selected AERONET stations. Despite the underestimate of aerosol optical depth in biomass burning regions in GEOS, there is an overestimate in the surface mass of organic carbon in the United States, especially during the summer months.

Allison B Collow↗

Change in Tropospheric Ozone in the Recent Decades and its Contribution to Global Total Ozone

Tropospheric ozone is a key chemically active trace gas and radiative forcer. Understanding its long-term changes is important to properly interpret observed changes in total column ozone and stratospheric ozone recovery. We investigate global and regional tropospheric ozone changes and their impact on total column ozone during 2005-2018 using satellite measurements and the NASA Goddard Earth Observing System Chemistry Climate Model (GEOSCCM). Global total ozone increased ~4 DU during 2005-2018 (+0.28±0.06 DU yr(exp -1)) as inferred from Ozone Monitoring Instrument (OMI). Consistent with previous studies, the OMI/MLS (Microwave Limb Sounder) derived global tropospheric ozone increased 2.2 DU during this period, 60% of the global total column ozone increase. While GEOSCCM reproduces reasonably well the total column increase, it overestimates the stratospheric ozone increase and underestimates the tropospheric ozone increase. We find that the tropospheric ozone increases are likely attributed to a growth of regional emissions of key ozone precursors, especially volatile organic compounds (VOCs) as reflected by the positive trends in formaldehyde (CH2O). Although carbon monoxide (CO) has been decreasing everywhere around the globe, it has relatively small impact on the tropospheric ozone trend. Trends in nitrogen dioxide (NO2) vary with regions, and these changes counteract or reinforce the positive effects of CH2O on the tropospheric ozone increases. The model underestimates the observed tropospheric ozone increase, especially over the US and Europe, because of underestimated NO2 emissions changes used in the model. The stratospheric ozone contribution increases during this period in the Northern Hemisphere and contributes to the tropospheric ozone increase.

Atmospheric chemistry↗

Dynamical Perturbation of the Stratosphere by a Pyrocumulonimbus Injection of Carbonaceous Aerosols

The Pacific Northwest pyrocumulonimbus Event (PNE) took place in British Columbia during the evening and night-time hours between the 12th and 13th of August 2017. Several pyroconvective clouds erupted in this occasion, and released in the upper troposphere/lower stratosphere unprecedented amounts of carbonaceous aerosols (300 ktn). Only a few years later, an even larger pyroCb injection took place over Australia. This event, named the Australian New Year (ANY) event, injected up to 1100 ktn of aerosol between December 29th 2019 and January 4th 2020. Such large injections of carbonaceous aerosol modify the stratospheric radiative budgets, locally perturbing stratospheric temperatures and winds. In this study, we use the Goddard Earth Observing System Chemistry Climate Model (GEOS CCM) to study the pertubations on the stratospheric meteorology induced by an aerosol injection of the magnitude of the PNE. Our simulations include the radiative interactions of aerosols, so that their impact on temperatures and winds are explicitly simulated. We show how the presence of the carbonaceous aerosols from the pyroCb causes the formation and maintenance of a synoptic scale stratospheric anticyclone. We follow this disturbance considering the potential vorticity anomaly and the brown carbon aerosol loading and we describe its dynamical and thermodynamical structure and its evolution in time. The analysis presented here shows that the simulated anticyclone undergoes daily expansion-compression cycles governed by the radiative heating, which are directly related to the vertical motion of the plume, and that the aerosol radiative heating is essential in maintaining the anticyclone itself.

Pyrocumulonimbus↗

A Coupled Evaluation of Operational MODIS and Model Aerosol Products for Maritime Environments Using Sun Photometry: Evaluation of the Fine and Coarse Mode

Although satellite retrievals and data assimilation have progressed to where there is a good skill for monitoring maritime Aerosol Optical Depth (AOD), there remains uncertainty in achieving further degrees of freedom, such as distinguishing fine and coarse mode dominated species in maritime environments (e.g., coarse mode sea salt and dust versus fine mode terrestrial anthropogenic emissions, biomass burning, and maritime secondary production). For the years 2016 through 2019, we performed an analysis of 550 nm total AOD550, fine mode AOD (FAOD550; also known as FM AOD in the literature), coarse mode AOD (CAOD550), and fine mode fraction (η550) between Moderate Resolution Spectral Imaging Radiometer (MODIS) V6.1 MOD/MYD04 dark target aerosol retrievals and the International Cooperative for Aerosol Prediction (ICAP) core four multi-model consensus (C4C) of analyses/short term forecasts that assimilate total MODIS AOD550. Differences were adjudicated by the global shipboard Maritime Aerosol Network (MAN) and selected island AERONET sun photometer observations with the application of the spectral deconvolution algorithm (SDA). Through a series of conditional and regional analyses, we found divergence included regions of terrestrial influence and latitudinal dependencies in the remote oceans. Notably, MODIS and the C4C and its members, while having good correlations overall, have a persistent +0.04 to +0.02 biases relative to MAN and AERONET for typical AOD550 values (84th% < 0.28), with the C4C underestimating significant events thereafter. Second, high biases in AOD550 are largely associated with the attribution of the fine mode in satellites and models alike. Thus, both MODIS and C4C members are systematically overestimating AOD550 and FAOD550 but perform better in characterizing the CAOD550. Third, for MODIS, findings are consistent with previous reports of a high bias in the retrieved Ångström Exponent, and we diagnosed both the optical model and cloud masking as likely causal factors for the AOD550 and FAOD550 high bias, whereas for the C4C, it is likely from secondary overproduction and perhaps numerical diffusion. Fourth, while there is no wind-speed-dependent bias for surface winds <12 m s−1, the C4C and MODIS AOD550s also overestimate CAOD550 and FAOD550, respectively, for wind speeds above 12 m/s. Finally, sampling bias inherent in MAN, as well as other circumstantial evidence, suggests biases in MODIS are likely MODIS and the C4C products have their own strengths and challenges for a given climate application and discuss needed research. even larger than what was diagnosed here. We conclude with a discussion on how

Aerosols↗

The Response of the Amazon Ecosystem to the Photosynthetically Active Radiation Fields: Integrating Impacts of Biomass Burning Aerosol and Clouds in the NASA GEOS Earth System Model

The Amazon experiences fires every year, and the resulting biomass burning aerosols, together with cloud particles, influence the penetration of sunlight through the atmosphere, increasing the ratio of diffuse to direct photosynthetically active radiation (PAR) reaching the vegetation canopy and thereby potentially increasing ecosystem productivity. In this study, we use the NASA Goddard Earth Observing System (GEOS) model with coupled aerosol, cloud, radiation, and ecosystem modules to investigate the impact of Amazon biomass burning aerosols on ecosystem productivity, as well as the role of the Amazon’s clouds in tempering this impact. The study focuses on a seven-year period (2010-2016) during which the Amazon experienced a variety of dynamic environments (e.g., La Niña, normal years, and El Niño). The direct radiative impact of biomass burning aerosols on ecosystem productivity—called here the aerosol diffuse radiation fertilization effect —is found to increase Amazonian Gross Primary Production (GPP) by 2.6% via a 3.8% increase in diffuse PAR (DFPAR) despite a 5.4% decrease in direct PAR (DRPAR) on multiyear average during burning seasons. On a monthly basis, this increase in GPP can be as large as 9.9% (occurring in August 2010). Consequently, the net primary production (NPP) in Amazon is increased by 1.5%, or ~92 Tg C a-1– equivalent to ~37% of the average carbon lost due to Amazon fires over the seven years considered. Clouds, however, strongly regulate the effectiveness of the aerosol diffuse radiation fertilization effect. The efficiency of this fertilization effect is the highest in cloud-free conditions and linearly decreases with increasing cloud amount until the cloud fraction reaches ~0.8, at which point the aerosol-influenced light changes from being a stimulator to an inhibitor of plant growth. Nevertheless, interannual changes in the overall strength of the aerosol diffuse radiation fertilization effect are primarily controlled by the large interannual changes in biomass burning aerosols rather than by changes in cloudiness during the studied period.

biomass burning aerosols on ecosystem productivity↗

The Long-term Transport and Radiative Impacts of the 2017 British Columbia Pyrocumulonimbus Smoke Aerosols in the Stratosphere

Interactions of meteorology with wildfires in British Columbia, Canada during August 2017 led to three major pyrocumulonimbus (pyroCb) events that resulted in the injection of large amounts of smoke aerosols and other combustion products at the local upper troposphere and lower stratosphere (UTLS). These plumes of UTLS smoke with elevated values of aerosol extinction and backscatter compared to the background state were readily tracked by multiple satellite-based instruments as they spread across the Northern Hemisphere (NH). The plumes were observed in the lower stratosphere for about 8-10 months following the fire injections, with a stratospheric aerosol e-folding time of about 5 months. To investigate the radiative impacts of these events on the Earth system, we performed a number of simulations with the Goddard Earth Observing System (GEOS) atmospheric general circulation model (AGCM). Observations from multiple remote-sensing instruments were used to calibrate the injection parameters (location, amount, composition and heights) and optical properties of the smoke aerosols in the model. The resulting simulations of three-dimensional smoke transport were evaluated for a year from the day of injections using daily observations from OMPS-LP (Ozone Mapping Profiler Suite Limb Profiler). The model simulated rate of ascent, hemispheric spread and residence time (or e-folding time) of the smoke aerosols in the stratosphere are in close agreement with OMPS-LP observations. We found that both aerosol self-lofting and the large-scale atmospheric motion play important roles in lifting the smoke plumes from near the tropopause altitudes (~12 km) to about 22-23 km into the atmosphere. Further, our estimations of the radiative impacts of the pyroCb-emitted smoke aerosols showed that the smoke caused an additional warming of the atmosphere by about 0.6-1 W/m2 (zonal mean) that persisted for about 2-3 months after the injections in regions north of 40oN. The surface experienced a comparable magnitude of cooling. The atmospheric warming is mainly located in the stratosphere, coincident with the location of the smoke plumes, leading to an increase in zonal mean shortwave (SW) heating rates of 0.02-0.04 K/day during September 2017.

pyrocumulonimbus (pyroCb) events↗

Improved Representation of the Global Dust Cycle Using Observational Constraints on Dust Properties and Abundance

Even though desert dust is the most abundant aerosol by mass in Earth’s atmosphere, atmospheric models 25 struggle to accurately represent its spatial and temporal distribution. These model errors are partially caused byfundamental difficulties in simulating dust emission in coarse-resolution models and in accurately representing dust microphysical properties. Here we mitigate these problems by developing a new methodology that yields an improved representation of the global dust cycle. We present an analytical framework that uses inverse modeling to integrate an ensemble of global model simulations with observational constraints on the dust size distribution, 30 extinction efficiency, and regional dust aerosol optical depth. We then compare the inverse model results against independent measurements of dust surface concentration and deposition flux and find that errors are reduced by approximately a factor of two relative to current model simulations of the Northern Hemisphere dust cycle. The inverse model results show smaller improvements in the less dusty Southern Hemisphere, most likely because both the model simulations and the observational constraints used in the inverse model are less accurate. On a global 35 basis, we find that the emission flux of dust with geometric diameter up to 20 μm (PM20) is approximately 5,000 Tg/year, which is greater than most models account for. This larger PM20 dust flux is needed to match observational constraints showing a large atmospheric loading of coarse dust. We obtain gridded data sets of dust emission, vertically integrated loading, dust aerosol optical depth, (surface) concentration, and wet and dry deposition fluxes that are resolved by season and particle size. As our results indicate that this data set is more accurate than current 40 model simulations and the MERRA-2 dust reanalysis product, it can be used to improve quantifications of dust impacts on the Earth system.

aerosols, dust↗

Contribution of the World's Main Dust Source Regions to the Global Cycle of Desert Dust

Even though desert dust is the most abundant aerosol by mass in Earth's atmosphere, the relative contributions of the world's major source regions to the global dust cycle remain poorly constrained. This problem hinders accounting for the potentially large impact of regional differences in dust properties on clouds, the Earth's energy balance, and terrestrial and marine biogeochemical cycles. Here, we constrain the contribution of each of the world's main dust source regions to the global dust cycle. We use an analytical framework that integrates an ensemble of global aerosol model simulations with observationally informed constraints on the dust size distribution, extinction efficiency, and regional dust aerosol optical depth (DAOD). We obtain a dataset that constrains the relative contribution of nine major source regions to size-resolved dust emission, atmospheric loading, DAOD, concentration, and deposition flux. We find that the 22–29 Tg (1 standard error range) global loading of dust with a geometric diameter up to 20 µm is partitioned as follows: North African source regions contribute ∼ 50 % (11–15 Tg), Asian source regions contribute ∼ 40 % (8–13 Tg), and North American and Southern Hemisphere regions contribute ∼ 10 % (1.8–3.2 Tg). These results suggest that current models on average overestimate the contribution of North African sources to atmospheric dust loading at ∼ 65 %, while underestimating the contribution of Asian dust at ∼ 30 %. Our results further show that each source region's dust loading peaks in local spring and summer, which is partially driven by increased dust lifetime in those seasons. We also quantify the dust deposition flux to the Amazon rainforest to be ∼ 10 Tg yr−1, which is a factor of 2–3 less than inferred from satellite data by previous work that likely overestimated dust deposition by underestimating the dust mass extinction efficiency. The data obtained in this paper can be used to obtain improved constraints on dust impacts on clouds, climate, biogeochemical cycles, and other parts of the Earth system.

Dust↗

Improved representation of the global dust cycle using observational constraints on dust properties and abundance

Even though desert dust is the most abundant aerosol by mass in Earth's atmosphere, atmospheric models struggle to accurately represent its spatial and temporal distribution. These model errors are partially caused by fundamental difficulties in simulating dust emission in coarse-resolution models and in accurately representing dust microphysical properties. Here we mitigate these problems by developing a new methodology that yields an improved representation of the global dust cycle. We present an analytical framework that uses inverse modeling to integrate an ensemble of global model simulations with observational constraints on the dust size distribution, extinction efficiency, and regional dust aerosol optical depth. We then compare the inverse model results against independent measurements of dust surface concentration and deposition flux and find that errors are reduced by approximately a factor of 2 relative to current model simulations of the Northern Hemisphere dust cycle. The inverse model results show smaller improvements in the less dusty Southern Hemisphere, most likely because both the model simulations and the observational constraints used in the inverse model are less accurate. On a global basis, we find that the emission flux of dust with a geometric diameter up to 20 µm (PM20) is approximately 5000 Tg yr−1, which is greater than most models account for. This larger PM20 dust flux is needed to match observational constraints showing a large atmospheric loading of coarse dust. We obtain gridded datasets of dust emission, vertically integrated loading, dust aerosol optical depth, (surface) concentration, and wet and dry deposition fluxes that are resolved by season and particle size. As our results indicate that this dataset is more accurate than current model simulations and the MERRA-2 dust reanalysis product, it can be used to improve quantifications of dust impacts on the Earth system.

Dust cycle↗

Chemical Transport Models Often Underestimate Inorganic Atmospheric Aerosol Acidity in Remote Regions of the Atmosphere

The inorganic fraction of fine particles affects numerous physicochemical processes in the atmosphere. However, there is large uncertainty in its burden and composition due to limited global measurements. Here, we present observations from eleven different aircraft campaigns from around the globe and investigate how aerosol pH and ammonium balance change from polluted to remote regions, such as over the oceans. Both parameters show increasing acidity with remoteness, at all altitudes, with pH decreasing from about 3 to about -1 and ammonium balance decreasing from almost 1 to nearly 0. We compare these observations against nine widely used chemical transport models and find that the simulations show more scatter (generally R(exp 2) < 0.50) and typically predict less acidic aerosol in the most remote regions. These differences in observations and predictions are likely to result in underestimating the model-predicted direct radiative cooling effect for sulfate, nitrate, and ammonium aerosol by 15-39%.

Chemical transport model↗

Revisiting the Relationship between Atlantic Dust and Tropical Cyclone Activity using Aerosol Optical Depth Reanalyses: 2003-2018

Previous studies have noted a relationship between African dust and Atlantic tropical cyclone (TC) activity. However, due to the limitations of past dust analyses, the strength of this relationship remains uncertain. The emergence of aerosol reanalyses, including the Navy Aerosol Analysis and Prediction System (NAAPS) aerosol optical depth (AOD) reanalysis, NASA Modern-Era Retrospective analysis for Research and Applications, Version 2 (MERRA-2), and ECMWF Copernicus Atmosphere Monitoring Service reanalysis (CAMSRA), enables an investigation of the relationship between African dust and TC activity over the tropical Atlantic and Caribbean in a consistent temporal and spatial manner for 2003–2018. Although June–July–August (JJA) 550 nm dust AOD (DAOD) from all three reanalysis products correlates significantly over the tropical Atlantic and Caribbean, the difference in DAOD magnitude between products can be as large as 60 % over the Caribbean and 20 % over the tropical North Atlantic. Based on the three individual reanalyses, we have created an aerosol multi-reanalysis consensus (MRC). The MRC presents overall better root mean square error over the tropical Atlantic and Caribbean compared to individual reanalyses when verified with ground-based AErosol RObotic NETwork (AERONET) AOD measurements. Each of the three individual reanalyses and the MRC have significant negative correlations between JJA Caribbean DAOD and seasonal Atlantic accumulated cyclone energy (ACE), while the correlation between JJA tropical North Atlantic DAOD and seasonal ACE is weaker. Possible reasons for this regional difference are provided. A composite analysis of 3 high-JJA-Caribbean-DAOD years versus 3 low-JJA-Caribbean-DAOD years reveals large differences in overall Atlantic TC activity. We also show that JJA Caribbean DAOD is significantly correlated with large-scale fields associated with variability in interannual Atlantic TC activity including zonal wind shear, mid-level moisture, and sea surface temperature (SST), as well as the El Niño–Southern Oscillation (ENSO) and the Atlantic Meridional Mode (AMM), implying confounding effects of these factors on the dust–TC relationship. We find that seasonal Atlantic DAOD and the AMM, the leading mode of coupled Atlantic variability, are inversely related and intertwined in the dust–TC relationship. Overall, DAOD in both the tropical Atlantic and Caribbean is negatively correlated with Atlantic hurricane frequency and intensity, with stronger correlations in the Caribbean than farther east in the tropical North Atlantic.

Dust particles↗

Apparent dust size discrepancy in aerosol reanalysis in north African dust after long-range transport

North African dust reaches the southeastern United States every summer. Size-resolved dust mass measurements taken in Miami, Florida, indicate that more than one-half of the surface dust mass concentrations reside in particles with geometric diameters less than 2.1 µm, while vertical profiles of micropulse lidar depolarization ratios show dust reaching above 4 km during pronounced events. These observations are compared to the representation of dust in the Modern-Era Retrospective analysis for Research and Applications, version 2 (MERRA-2) aerosol reanalysis and closely related Goddard Earth Observing System model version 5 (GEOS-5) Forward Processing (FP) aerosol product, both of which assimilate satellite-derived aerosol optical depths using a similar protocol and inputs. These capture the day-to-day variability in aerosol optical depth well, in a comparison to an independent sun-photometer-derived aerosol optical depth dataset. Most of the modeled dust mass resides in diameters between 2 and 6 µm, in contrast to the measurements. Model-specified mass extinction efficiencies equate light extinction with approximately 3 times as much aerosol mass, in this size range, compared to the measured dust sizes. GEOS-5 FP surface-layer sea salt mass concentrations greatly exceed observed values, despite realistic winds and relative humidities. In combination, these observations help explain why, despite realistic total aerosol optical depths, (1) free-tropospheric model volume extinction coefficients are lower than those retrieved from the micro-pulse lidar, suggesting too-low model dust loadings in the free troposphere, and (2) model dust mass concentrations near the surface can be higher than those measured. The modeled vertical distribution of dust, when captured, is reasonable. Large, aspherical particles exceeding the modeled dust sizes are also occasionally present, but dust particles with diameters exceeding 10 µm contribute little to the measured total dust mass concentrations after such long-range transport. Remaining uncertainties warrant a further integrated assessment to confirm this study's interpretations.

MEERA-2 aerosols reanalysis↗

Optimized Profile Retrievals of Aerosol Microphysical Properties from Simulated Spaceborne Multiwavelength Lidar

This work is an expanded study of one previously published onretrievals of aerosol microphysical properties from space-borne multiwavelengthlidarmeasurements. The earlier studiesand this one weredone in the framework of the NASA Aerosol-Clouds-Ecosystems (now the Aerosol Clouds Convection and Precipitation) NASA mission. The focus here is on the capabilities of a simulated spacebornemultiwavelengthlidar system for retrieving aerosol complex refractive index (m = mr+ imi) and spectral single scattering albedo (SSA(λ)), although other bulk parameters such as effective (reff) radius and particle volume (V) and surface (S) concentrations are also studied. The novelty presented here is the use of recently published, case dependent optimized-constraints on the microphysical retrievals using three backscattering coefficients (β) at 355, 532 and 1064 nm and two extinction coefficients (α) at 355 and 532 nm, typically known as the stand-alone 3β+2α lidar inversion. Case-dependent optimized-constraints (CDOC) limit the ranges of refractive index, both real (mr) and imaginary (mi) parts, and of radii that are permitted in the retrievals. Such constraints are selected directly from the 3β+2α41measurements through an analysis of the relationship between spectral dependence of aerosol extinction-to-backscatter ratios (LR) and the Ångström exponent of extinction. The analyses presented here for different sets ofsize distributions and refractive indices reveal that the direct determination of CDOCareonly feasible for cases where the uncertaintiesin the input optical data areless than 15 %.Forthe same simulated spacebornesystem and yield than in Whiteman et al., (2018), we demonstrated that the use of CDOC as essential for the retrievals of refractive index and also largely improved retrieval of bulk parameters. A discussion of the global representativeness of CDOC is presented using simulated lidar data from a 24-hour satellite track using GEOS model output to initialize the lidar simulator.We found that CDOCare representative of many aerosol mixtures in spite of some outliers (e.g. highly hydrated particles) associatedwith the assumptions of bimodal size distributions and of the same refractive index for fine and coarse modes. Moreover, sensitivity tests performed using synthetic data reveal that retrievals of imaginary refractive index (mi) and SSA are extremely sensitive to β(355).

NASA Aerosol-Clouds-Ecosystems↗

Total Dust Deposition Flux During Precipitation in Toyama, Japan, in the Spring of 2009: A Sensitivity Analysis with the NASA GEOS-5 Model

We compared the observed total dust deposition fluxes during precipitation (TDP) mainly at Toyama in Japan during the period January - April 2009 with results available from four NASA GEOS-5 global model experiments. The modeled results were obtained from three previous experiments and carried out in one experiment, which were all driven by assimilated meteorology and simulating aerosol distributions for the time period. We focus mainly on the observations of two distinct TDP events, which were reported in Osada et al. (2011), at Toyama, Japan, in February (Event B) and March 2009 (Event C). Although all of our GEOS-5 simulations captured aspects of the observed TDP, we found that our low horizontal spatial resolution control experiment performed generally the worst. The other three experiments were run at a higher spatial resolution, with the first differing only in that respect from the control, the second adding imposed a prescribed corrected precipitation product, and the final experiment adding as well assimilation of aerosol optical depth based on MODIS observations. During Event C, the increased horizontal resolution could increase TDP with precipitation increase. There was no significant improvement, however, due to the imposition of the corrected precipitation product. The simulation that incorporated aerosol data assimilation performed was by far the best for this event, but even so could only reproduce less than half of the observed TDP despite the significantly increased atmospheric dust mass concentrations. All three of the high spatial resolution experiments had higher simulated precipitation at Toyama than was observed and that in the lower resolution control run. During Event B, the aerosol data assimilation run did not perform appreciably better than the other higher resolution simulations, suggesting that upstream conditions (i.e., upstream cloudiness), or vertical or horizontal misplacement of the dust plume did not allow for significant improvement in the simulated aerosol distributions. Furthermore, a detailed comparison of observed hourly precipitation and surface particulate mass concentration data suggests that the observed TDP during Event B was highly dependent on short periods of weak precipitation correlated with elevated dust surface concentrations, important details possibly not captured well in a current global model.

Dust; Wet deposition; Numerical modeling; Japan; P↗

The GOddard SnoW Impurity Module (GOSWIM) for the NASA GEOS-5 Earth System Model: Preliminary Comparisons with Observations in Sapporo, Japan

The snow darkening module evaluating dust, black carbon, and organic carbon depositions on mass and albedo has been developed for the NASA Goddard Earth Observing System, Version 5 (GEOS-5) Earth System Model, as the GOddard SnoW Impurity Module (GOSWIM). GOSWIM consists of the snow albedo scheme from a previous study (Yasunari et al. 2011) with updates and a newly developed mass concentration scheme, using aerosol depositions from the chemical transport model (GOCART) in GEOS-5. Compared to observations at Sapporo, the numerical experiments, forced by observation-based meteorology and aerosol depositions from GOES-5, better simulated the seasonal migration of snow depth, albedos, and impurities of dust, BC, and OC in the snow surface. However, the magnitude of the impurities is underestimated, compared to the sporadic snow impurity measurements. Increasing the deposition rates of dust and BC could explain the differences on the snow darkening effect between observation and simulation. Ignoring BC deposition can possibly lead to an extension of snow cover duration in Sapporo for four days. Comparing the off-line GOSWIM and the GEOS-5 global simulations, we found that determining better local precipitation and deposition rates of the aerosols are key factors in generating better GOSWIM snow darkening simulation in NASA GEOS-5.

land surface processes↗

Impacts of aerosol size-dependent below-cloud scavenging on tropospheric aerosol in the NASA GEOS model

Wet deposition is the major sink for a suite of tropospheric aerosols, and thus largely affects the global distribution and lifetime of aerosols as well as the formation of clouds in the atmosphere. In this study, we implement the aerosol size-dependent below-cloud scavenging (BCS; washout) parameterization of Croft et al. (Atmos. Chem. Phys., 2009; hereafter referred to as C09) in the GOCART bulk aerosol scheme coupled with the NASA Goddard Earth Observing System (GEOS) model, and evaluate the simulated tropospheric aerosols against surface, airborne, and satellite observations. The current below-cloud scavenging in GOCART employs a first-order removal function with a constant scavenging coefficient that is independent of aerosol sizes and precipitation fluxes. The C09 BCS parameterization instead calculates aerosol below-cloud scavenging coefficients based on a size-dependent collision efficiency between aerosol particles and raindrops (or snow crystals). A look-up table of compiled collision efficiencies provides scavenging coefficients for sixty aerosol sizes and nine rain flux rates. We determine instantaneous coefficients according to the C09 BCS look-up table, precipitation fluxes, and a log-normal size distribution of aerosols. The implemented BCS parameterization tends to result in a large diversity of scavenging coefficient values with higher magnitudes, leading to more efficient scavenging. The model-simulated global mean tropospheric Pb-210 aerosol lifetime decreases from 7.4 days to 5.8 days. The separate consideration of rain and snow conditions in the scavenging coefficient lookup table leads to a bi-mode distribution of scavenging coefficients. These changes also cause increased mass fractions of smaller aerosols compared to the current simulation. We examine the impacts of the C09 BCS parameterization on the simulated aerosol vertical profiles, tropospheric lifetimes, and deposition fluxes through comparisons with NASA aircraft measurements (ATom), surface observations of aerosol concentrations (IMPROVE) and wet deposition fluxes (NADP and APQMP), and MODIS aerosol optical depth (AOD) retrievals.

Aerosols↗