Search

We investigate 3-D mountain effects on solar flux distributions and their impact on surface hydrology over the western United States, specifically the Rocky Mountains and the Sierra Nevada, using the global CCSM4 (Community Climate System Model version 4; Community Atmosphere Model/Community Land Model – CAM4/CLM4) with a 0.23° × 0.31° resolution for simulations over 6 years. In a 3-D radiative transfer parameterization, we have updated surface topography data from a resolution of 1 km to 90 m to improve parameterization accuracy. In addition, we have also modified the upward-flux deviation (3-D–PP (plane-parallel)) adjustment to ensure that the energy balance at the surface is conserved in global climate simulations based on 3-D radiation parameterization. We show that deviations in the net surface fluxes are not only affected by 3-D mountains but also influenced by feedbacks of cloud and snow in association with the long-term simulations. Deviations in sensible heat and surface temperature generally follow the patterns of net surface solar flux. The monthly snow water equivalent (SWE) deviations show an increase in lower elevations due to reduced snowmelt, leading to a reduction in cumulative runoff. Over higher-elevation areas, negative SWE deviations are found because of increased solar radiation available at the surface. Simulated precipitation increases for lower elevations, while it decreases for higher elevations, with a minimum in April. Liquid runoff significantly decreases at higher elevations after April due to reduced SWE and precipitation.

We systematically evaluate the black carbon (BC) simulations for 2006 over the Tibetan Plateau by a global 3-D chemical transport model (CTM) (GEOS-Chem) driven by GEOS-5 assimilated meteorological fields, using in situ measurements of BC in surface air, BC in snow, and BC absorption aerosol optical depth (AAOD). Using improved anthropogenic BC emission inventories for Asia that account for rapid technology renewal and energy consumption growth (Zhang et al., 2009; Lu et al., 2011) and improved global biomass burning emission inventories that account for small fires (van der Werf et al., 2010; Randerson et al., 2012), we find that model results of both BC in surface air and in snow are statistically in good agreement with observations (biases < 15%) away from urban centers. Model results capture the seasonal variations of the surface BC concentrations at rural sites in the Indo-Gangetic Plain, but the observed elevated values in winter are absent. Modeled surface-BC concentrations are within a factor of 2 of the observations at remote sites. Part of the discrepancy is explained by the deficiencies of the meteorological fields over the complex Tibetan terrain. We find that BC concentrations in snow computed from modeled BC deposition and GEOS-5 precipitation are spatiotemporally consistent with observations (r = 0.85). The computed BC concentrations in snow are a factor of 2–4 higher than the observations at several Himalayan sites because of excessive BC deposition. The BC concentrations in snow are biased low by a factor of 2 in the central plateau, which we attribute to the absence of snow aging in the CTM and strong local emissions unaccounted for in the emission inventories. Modeled BC AAOD is more than a factor of 2 lower than observations at most sites, particularly to the northwest of the plateau and along the southern slopes of the Himalayas in winter and spring, which is attributable in large part to underestimated emissions and the assumption of external mixing of BC aerosols in the model. We find that assuming a 50% increase of BC absorption associated with internal mixing reduces the bias in modeled BC AAOD by 57% in the Indo-Gangetic Plain and the northeastern plateau and to the northeast of the plateau, and by 16% along the southern slopes of the Himalayas and to the northwest of the plateau. Both surface BC concentration and AAOD are strongly sensitive to anthropogenic emissions (from China and India), while BC concentration in snow is especially responsive to the treatment of BC aerosol aging. We find that a finer model resolution (0.5° × 0.667° nested over Asia) reduces the bias in modeled surface-BC concentration from 15 to 2%. The large range and non-homogeneity of discrepancies between model results and observations of BC across the Tibetan Plateau undoubtedly undermine current assessments of the climatic and hydrological impact of BC in the region and thus warrant imperative needs for more extensive measurements of BC, including its concentration in surface air and snow, AAOD, vertical profile and deposition.

The version 6 cloud products of the Atmospheric Infrared Sounder (AIRS) and Advanced Microwave Sounding Unit (AMSU) instrument suite are described. The cloud top temperature, pressure, and height and effective cloud fraction are now reported at the AIRS field-of-view (FOV) resolution. Significant improvements in cloud height assignment over version 5 are shown with FOV-scale comparisons to cloud vertical structure observed by the CloudSat 94 GHz radar and the Cloud-Aerosol LIdar with Orthogonal Polarization (CALIOP). Cloud thermodynamic phase (ice, liquid, and unknown phase), ice cloud effective diameter (De), and ice cloud optical thickness (τ) are derived using an optimal estimation methodology for AIRS FOVs, and global distributions for 2007 are presented. The largest values of τ are found in the storm tracks and near convection in the tropics, while De is largest on the equatorial side of the midlatitude storm tracks in both hemispheres, and lowest in tropical thin cirrus and the winter polar atmosphere. Over the Maritime Continent the diurnal variability of τ is significantly larger than for the total cloud fraction, ice cloud frequency, and De, and is anchored to the island archipelago morphology. Important differences are described between northern and southern hemispheric midlatitude cyclones using storm center composites. The infrared-based cloud retrievals of AIRS provide unique, decadal-scale and global observations of clouds over portions of the diurnal and annual cycles, and capture variability within the mesoscale and synoptic scales at all latitudes.

We investigate 3-D mountains/snow effects on solar flux distributions and their impact on surface hydrology over the western United States, specifically the Rocky Mountains and Sierra Nevada. The Weather Research and Forecasting (WRF) model, applied at a 30 km grid resolution, is used in conjunction with a 3-D radiative transfer parameterization covering a time period from 1 November 2007 to 31 May 2008, during which abundant snowfall occurred. A comparison of the 3-D WRF simulation with the observed snow water equivalent (SWE) and precipitation from Snowpack Telemetry (SNOTEL) sites shows reasonable agreement in terms of spatial patterns and daily and seasonal variability, although the simulation generally has a positive precipitation bias. We show that 3-D mountain features have a profound impact on the diurnal and monthly variation of surface radiative and heat fluxes, and on the consequent elevation-dependence of snowmelt and precipitation distributions. In particular, during the winter months, large deviations (3-D-PP, in which PP denotes the plane-parallel approach) of the monthly mean surface solar flux are found in the morning and afternoon hours due to shading effects for elevations below 2.5 km. During spring, positive deviations shift to the earlier morning. Over mountaintops higher than 3 km, positive deviations are found throughout the day, with the largest values of 40–60 W m−2 occurring at noon during the snowmelt season of April to May. The monthly SWE deviations averaged over the entire domain show an increase in lower elevations due to reduced snowmelt, which leads to a reduction in cumulative runoff. Over higher elevation areas, positive SWE deviations are found because of increased solar radiation available at the surface. Overall, this study shows that deviations of SWE due to 3-D radiation effects range from an increase of 18% at the lowest elevation range (1.5–2 km) to a decrease of 8% at the highest elevation range (above 3 km). Since lower elevation areas occupy larger fractions of the land surface, the net effect of 3-D radiative transfer is to extend snowmelt and snowmelt-driven runoff into the warm season. Because 60–90% of water resources originate from mountains worldwide, the aforementioned differences in simulated hydrology due solely to 3-D interactions between solar radiation and mountains/snow merit further investigation in order to understand the implications of modeling mountain water resources, and these resources' vulnerability to climate change and air pollution.

Atmospheric mineral dust particles exert significant direct radiative forcings and are important drivers of climate and climate change. We used the GEOS-Chem global three-dimensional chemical transport model (CTM) coupled with the Fu-Liou-Gu (FLG) radiative transfer model (RTM) to investigate the dust radiative forcing and heating rate based on different vertical profiles for April 2006. We attempt to actually quantify the sensitivities of radiative forcing to dust vertical profiles, especially the discrepancies between using realistic and climatological vertical profiles. In these calculations, dust emissions were constrained by observations of aerosol optical depth (AOD). The coupled calculations utilizing a more realistic dust vertical profile simulated by GEOS-Chem minimize the physical inconsistencies between 3-D CTM aerosol fields and the RTM. The use of GEOS-Chem simulated vertical profile of dust extinction, as opposed to the FLG prescribed vertical profile, leads to greater and more spatially heterogeneous changes in the estimated radiative forcing and heating rate produced by dust. Both changes can be attributed to a different vertical structure between dust and non-dust source regions. Values of the dust vertically resolved AOD per grid level (VRAOD) are much larger in the middle troposphere, though smaller at the surface when the GEOS-Chem simulated vertical profile is used, which leads to a much stronger heating rate in the middle troposphere. Compared to the FLG vertical profile, the use of GEOS-Chem vertical profile reduces the solar radiative forcing at the top of atmosphere (TOA) by approximately 0.2–0.25 W m−2 over the African and Asian dust source regions. While the Infrared (IR) radiative forcing decreases 0.2 W m−2 over African dust belt, it increases 0.06 W m−2 over the Asian dust belt when the GEOS-Chem vertical profile is used. Differences in the solar radiative forcing at the surface between the use of the GEOS-Chem and FLG vertical profiles are most significant over the Gobi desert with a value of about 1.1 W m−2. The radiative forcing effect of dust particles is more pronounced at the surface over the Sahara and Gobi deserts by using FLG vertical profile, while it is less significant over the downwind area of Eastern Asia.

A three-dimensional variational data assimilation (3-DVAR) algorithm for aerosols in a WRF/Chem model is presented. The WRF/Chem model uses the MOSAIC (Model for Simulating Aerosol Interactions and Chemistry) scheme, which explicitly treats eight major species (elemental/black carbon, organic carbon, nitrate, sulfate, chloride, ammonium, sodium and the sum of other inorganic, inert mineral and metal species) and represents size distributions using a sectional method with four size bins. The 3-DVAR scheme is formulated to take advantage of the MOSAIC scheme in providing comprehensive analyses of species concentrations and size distributions. To treat the large number of state variables associated with the MOSAIC scheme, this 3-DVAR algorithm first determines the analysis increments of the total mass concentrations of the eight species, defined as the sum of the mass concentrations across all size bins, and then distributes the analysis increments over four size bins according to the background error variances. The number concentrations for each size bin are adjusted based on the ratios between the mass and number concentrations of the background state. Additional flexibility is incorporated to further lump the eight mass concentrations, and five lumped species are used in the application presented. The system is evaluated using the analysis and prediction of PM2.5 in the Los Angeles basin during the CalNex 2010 field experiment, with assimilation of surface PM2.5 and speciated concentration observations. The results demonstrate that the data assimilation significantly reduces the errors in comparison with a simulation without data assimilation and improved forecasts of the concentrations of PM2.5 as well as individual species for up to 24 h. Some implementation difficulties and limitations of the system are discussed.

A surface solar radiation parameterization based on deviations between 3-D and conventional plane-parallel radiative transfer models has been incorporated into the Weather Research and Forecasting (WRF) model to understand the solar insolation over mountain/snow areas and to investigate the impact of the spatial and temporal distribution and variation of surface solar fluxes on land-surface processes. Using the Sierra-Nevada in the western United States as a testbed, we show that mountain effect could produce up to −50 to + 50 W m−2 deviations in the surface solar fluxes over the mountain areas, resulting in a temperature increase of up to 1 °C on the sunny side. Upward surface sensible and latent heat fluxes are modulated accordingly to compensate for the change in surface solar fluxes. Snow water equivalent and surface albedo both show decreases on the sunny side of the mountains, indicating more snowmelt and hence reduced snow albedo associated with more solar insolation due to mountain effect. Soil moisture increases on the sunny side of the mountains due to enhanced snowmelt, while decreases on the shaded side. Substantial differences are found in the morning hours from 8–10 a.m. and in the afternoon around 3–5 p.m., while differences around noon and in the early morning and late afternoon are comparatively smaller. Variation in the surface energy balance can also affect atmospheric processes, such as cloud fields, through the modulation of vertical thermal structure. Negative changes of up to −40 g m−2 are found in the cloud water path, associated with reductions in the surface insolation over the cloud region. The day-averaged deviations in the surface solar flux are positive over the mountain areas and negative in the valleys, with a range between −12~12 W m−2. Changes in sensible and latent heat fluxes and surface skin temperature follow the solar insolation pattern. Differences in the domain-averaged diurnal variation over the Sierras show that the mountain area receives more solar insolation during early morning and late afternoon, resulting in enhanced upward sensible heat and latent heat fluxes from the surface and a corresponding increase in surface skin temperature. During the middle of the day, however, the surface insolation and heat fluxes show negative changes, indicating a cooling effect. Hence overall, the diurnal variations of surface temperature and surface fluxes in the Sierra-Nevada are reduced through the interactions of radiative transfer and mountains. The hourly differences of the surface solar insolation in higher elevated regions, however, show smaller magnitude in negative changes during the middle of the day and possibly more solar fluxes received during the whole day.

The climatic effects of dust aerosols in North Africa have been investigated using the atmospheric general circulation model (AGCM) developed at the University of California, Los Angeles (UCLA). The model includes an efficient and physically based radiation parameterization scheme developed specifically for application to clouds and aerosols. Parameterization of the effective ice particle size in association with the aerosol first indirect effect based on ice cloud and aerosol data retrieved from A-Train satellite observations have been employed in climate model simulations. Offline simulations reveal that the direct solar, IR, and net forcings by dust aerosols at the top of the atmosphere (TOA) generally increase with increasing aerosol optical depth. When the dust semi-direct effect is included with the presence of ice clouds, positive IR radiative forcing is enhanced since ice clouds trap substantial IR radiation, while the positive solar forcing with dust aerosols alone has been changed to negative values due to the strong reflection of solar radiation by clouds, indicating that cloud forcing associated with aerosol semi-direct effect could exceed direct aerosol forcing. With the aerosol first indirect effect, the net cloud forcing is generally reduced in the case for an ice water path (IWP) larger than 20 g m−2. The magnitude of the reduction increases with IWP. AGCM simulations show that the reduced ice crystal mean effective size due to the aerosol first indirect effect results in less OLR and net solar flux at TOA over the cloudy area of the North Africa region because ice clouds with smaller size trap more IR radiation and reflect more solar radiation. The precipitation in the same area, however, increases due to the aerosol indirect effect on ice clouds, corresponding to the enhanced convection as indicated by reduced OLR. Adding the aerosol direct effect into the model simulation reduces the precipitation in the normal rainfall band over North Africa, where precipitation is shifted to the south and the northeast produced by the absorption of sunlight and the subsequent heating of the air column by dust particles. As a result, rainfall is drawn further inland to the northeast. This study represents the first attempt to quantify the climate impact of the aerosol indirect effect using a GCM in connection with A-Train satellite data. The parameterization for the aerosol first indirect effect developed in this study can be readily employed for application to other GCMs.

Forest fires are an important source to carbonaceous aerosols in the Western United States (WUS). We quantify the relative contribution of biomass burning to black carbon (BC) in the WUS mountain ranges by analyzing surface BC observations for 2006 from the Interagency Monitoring of PROtected Visual Environment (IMPROVE) network using the GEOS-Chem global chemical transport model. Observed surface BC concentrations show broad maxima during late June to early November. Enhanced potassium concentrations and potassium/sulfur ratios observed during the high-BC events indicate a dominant biomass burning influence during the peak fire season. Model surface BC reproduces the observed day-to day and synoptic variabilities in regions downwind of but near urban centers. Major discrepancies are found at elevated mountainous sites during the July-October fire season when simulated BC concentrations are biased low by a factor of two. We attribute these low biases largely to the underestimated (by more than a factor of two) and temporally misplaced biomass burning emissions of BC in the model. Additionally, we find that the biomass burning contribution to surface BC concentrations in the USA likely was underestimated in a previous study using GEOS-Chem (Park et al., 2003), because of the unusually low planetary boundary layer (PBL) heights in the GEOS-3 meteorological reanalysis data used to drive the model. PBL heights from GEOS-4 and GEOS-5 reanalysis data are comparable to those from the North American Regional Reanalysis (NARR). Model simulations show slightly improved agreements with the observations when driven by GEOS-5 reanalysis data, but model results are still biased low. The use of biomass burning emissions with diurnal cycle, synoptic variability, and plume injection has relatively small impact on the simulated surface BC concentrations in the WUS.

Global observations of cloud and humidity distributions in the upper troposphere within all geophysical conditions are critically important in order to monitor the present climate and to provide necessary data for validation of climate models to project future climate change. Towards this end, tropical oceanic distributions of thin cirrus optical depth (τ), effective diameter (De), and relative humidity with respect to ice (RHi) within cirrus (RHic) are simultaneously derived from the Atmospheric Infrared Sounder (AIRS). Corresponding increases in De and cloud temperature are shown for cirrus with τ>0.25 that demonstrate quantitative consistency to other surface-based, in situ and satellite retrievals. However, inferred cirrus properties are shown to be less certain for increasingly tenuous cirrus. In-cloud supersaturation is observed for 8–12% of thin cirrus and is several factors higher than all-sky conditions; even higher frequencies are shown for the coldest and thinnest cirrus. Spatial and temporal variations in RHic correspond to cloud frequency while regional variability in RHic is observed to be most prominent over the N. Indian Ocean basin. The largest cloud/clear sky RHi anomalies tend to occur in dry regions associated with vertical descent in the sub-tropics, while the smallest occur in moist ascending regions in the tropics. The characteristics of RHic frequency distributions depend on τ and a peak frequency is located between 60–80% that illustrates RHic is on average biased dry. The geometrical thickness of cirrus is typically less than the vertical resolution of AIRS temperature and specific humidity profiles and thus leads to the observed dry bias, shown with coincident cloud vertical structure obtained from the Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observation (CALIPSO). The joint distributions of thin cirrus microphysics and humidity derived from AIRS provide unique and important regional and global-scale insights on upper tropospheric processes not available from surface, in situ, and other contemporary satellite observing platforms.

Topography exerts significant influences on the incoming solar radiation at the land surface. A few stand-alone regional and global atmospheric models have included parameterizations for sub-grid topographic effects on solar radiation. However, nearly all Earth system models (ESMs) that participated in the Coupled Model Intercomparison Project (CMIP6) use a plane-parallel (PP) radiative transfer scheme that assumes that the terrain is flat. In this study, we incorporated a well-validated sub-grid topographic (TOP) parameterization in the Energy Exascale Earth System Model (E3SM) Land Model (ELM) version 1.0 to quantify the effects of sub-grid topography on solar radiation flux, including the shadow effects and multi-scattering between adjacent terrain. We studied the role of sub-grid topography by performing ELM simulations with the PP and TOP schemes over the Tibetan Plateau (TP). Additional ELM simulations were performed at multiple spatial resolutions to investigate the role of spatial scale on sub-grid topographic effects on solar radiation. The Moderate Resolution Imaging Spectroradiometer (MODIS) data was used to compare with the ELM simulations. The results show that topography has non-negligible effects on surface energy budget, snow cover, snow depth, and surface temperature over the TP. The absolute differences in surface energy fluxes for net solar radiation, latent heat flux, and sensible heat flux between TOP and PP exceed 20, 10, and 5 W m−2, respectively. The differences in land surface albedo, snow cover fraction, snow depth, and surface temperature between TOP and PP exceed 0.1, 0.1, 10 cm, and 1 K, respectively. The magnitude of the sub-grid topographic effects is dependent on seasons and elevations and is also sensitive to the spatial scales. Although the sub-grid topographic effects on solar radiation are larger with more spatial details at finer spatial scales, they cannot be simply neglected at coarse spatial scales. When compared to MODIS data, incorporating the sub-grid topographic effects overall reduces the biases of ELM in simulating surface energy balance, snow cover, and surface temperature, especially in the high-elevation and snow-covered regions over the TP. The inclusion of sub-grid topographic effects on solar radiation parameterization in ELM will contribute to advancing our understanding of the role of the surface topography on terrestrial processes over complex terrain.

In response to the coronavirus disease of 2019 (COVID-19), California issued statewide stay-at-home orders, bringing about abrupt and dramatic reductions in air pollutant emissions. This crisis offers us an unprecedented opportunity to evaluate the effectiveness of emission reductions in terms of air quality. Here we use the Weather Research and Forecasting model with Chemistry (WRF-Chem) in combination with surface observations to study the impact of the COVID-19 lockdown measures on air quality in southern California. Based on activity level statistics and satellite observations, we estimate the sectoral emission changes during the lockdown. Due to the reduced emissions, the population-weighted concentrations of fine particulate matter (PM2.5) decrease by 15 % in southern California. The emission reductions contribute 68 % of the PM2.5 concentration decrease before and after the lockdown, while meteorology variations contribute the remaining 32 %. Among all chemical compositions, the PM2.5 concentration decrease due to emission reductions is dominated by nitrate and primary components. For O3 concentrations, the emission reductions cause a decrease in rural areas but an increase in urban areas; the increase can be offset by a 70 % emission reduction in anthropogenic volatile organic compounds (VOCs). These findings suggest that a strengthened control on primary PM2.5 emissions and a well-balanced control on nitrogen oxides and VOC emissions are needed to effectively and sustainably alleviate PM2.5 and O3 pollution in southern California.

The climatic and health effects of aerosols are strongly dependent on the intra-annual variations in their loading and properties. While the seasonal variations of regional aerosol optical depth (AOD) have been extensively studied, understanding the temporal variations in aerosol vertical distribution and particle types is also important for an accurate estimate of aerosol climatic effects. In this paper, we combine the observations from four satellite-borne sensors and several ground-based networks to investigate the seasonal variations of aerosol column loading, vertical distribution, and particle types over three populous regions: the Eastern United States (EUS), Western Europe (WEU), and Eastern and Central China (ECC). In all three regions, column AOD, as well as AOD at heights above 800 m, peaks in summer/spring, probably due to accelerated formation of secondary aerosols and hygroscopic growth. In contrast, AOD below 800 m peaks in winter over WEU and ECC regions because more aerosols are confined to lower heights due to the weaker vertical mixing. In the EUS region, AOD below 800 m shows two maximums, one in summer and the other in winter. The temporal trends in low-level AOD are consistent with those in surface fine particle (PM2.5) concentrations. AOD due to fine particles ( 1.4 µm diameter), generally shows small variability, except that a peak occurs in spring in the ECC region due to the prevalence of airborne dust during this season. When aerosols are classified according to sources, the dominant type is associated with anthropogenic air pollution, which has a similar seasonal pattern as total AOD. Dust and sea-spray aerosols in the WEU region peak in summer and winter, respectively, but do not show an obvious seasonal pattern in the EUS region. Smoke aerosols, as well as absorbing aerosols, present an obvious unimodal distribution with a maximum occurring in summer over the EUS and WEU regions, whereas they follow a bimodal distribution with peaks in August and March (due to crop residue burning) over the ECC region.

A version of the WRF-Chem model with fully coupled aerosol–meteorology–snowpack is employed to investigate the impacts of various aerosol sources on precipitation and snowpack in California. In particular, the impacts of locally emitted anthropogenic and dust aerosols, and aerosols transported from outside California are studied. We differentiate three pathways of aerosol effects: aerosol–radiation interaction (ARI), aerosol–snow interaction (ASI), and aerosol–cloud interaction (ACI). The convection-permitting model simulations show that precipitation, snow water equivalent (SWE), and surface air temperature averaged over the whole domain (34–42° N, 117–124° W, not including ocean points) are reduced when aerosols are included, therefore reducing large biases in these variables due to the absence of aerosol effects in the model. Aerosols affect California water resources through the warming of mountaintops and the reduction of precipitation; however, different aerosol sources play different roles in changing surface temperature, precipitation, and snowpack in California by means of various weights of the three pathways. ARI by all aerosols mainly cools the surface, leading to slightly increased SWE over the mountains. Locally emitted dust aerosols warm the surface of mountaintops through ASI, in which the reduced snow albedo associated with dusty snow leads to more surface absorption of solar radiation and reduced SWE. Transported aerosols and local anthropogenic aerosols play a dominant role in increasing nonprecipitating clouds but reducing precipitation through ACI, leading to reduced SWE and runoff on the Sierra Nevada, as well as the warming of mountaintops associated with decreased SWE and hence lower surface albedo. The average changes in surface temperature from October 2012 to June 2013 are about −0.19 and 0.22 K for the whole domain and over mountaintops, respectively. Overall, the averaged reduction during October to June is about 7 % for precipitation, 3 % for SWE, and 7 % for surface runoff for the whole domain, while the corresponding numbers are 12, 10, and 10 % for the mountaintops. The reduction in SWE is more significant in a dry year, with 9 % for the whole domain and 16 % for the mountaintops. The maximum reduction of ∼ 20 % in precipitation occurs in May and is associated with the maximum aerosol loading, leading to the largest decrease in SWE and surface runoff over that period. It is also found that dust aerosols can cause early snowmelt on the mountaintops and reduced surface runoff after April.

The interactions between aerosols and ice clouds represent one of the largest uncertainties in global radiative forcing from pre-industrial time to the present. In particular, the impact of aerosols on ice crystal effective radius (Rei), which is a key parameter determining ice clouds' net radiative effect, is highly uncertain due to limited and conflicting observational evidence. Here we investigate the effects of aerosols on Rei under different meteorological conditions using 9-year satellite observations. We find that the responses of Rei to aerosol loadings are modulated by water vapor amount in conjunction with several other meteorological parameters. While there is a significant negative correlation between Rei and aerosol loading in moist conditions, consistent with the "Twomey effect" for liquid clouds, a strong positive correlation between the two occurs in dry conditions. Simulations based on a cloud parcel model suggest that water vapor modulates the relative importance of different ice nucleation modes, leading to the opposite aerosol impacts between moist and dry conditions. When ice clouds are decomposed into those generated from deep convection and formed in situ, the water vapor modulation remains in effect for both ice cloud types, although the sensitivities of Rei to aerosols differ noticeably between them due to distinct formation mechanisms. The water vapor modulation can largely explain the difference in the responses of Rei to aerosol loadings in various seasons. A proper representation of the water vapor modulation is essential for an accurate estimate of aerosol–cloud radiative forcing produced by ice clouds.

We develop and examine a microphysics-based black carbon (BC) aerosol aging scheme that accounts for condensation, coagulation, and heterogeneous chemical oxidation processes in a global 3-D chemical transport model (GEOS-Chem) by interpreting the BC measurements from the HIAPER Pole-to-Pole Observations (HIPPO, 2009–2011) using the model. We convert aerosol mass in the model to number concentration by assuming lognormal aerosol size distributions and compute the microphysical BC aging rate (excluding chemical oxidation aging) explicitly from the condensation of soluble materials onto hydrophobic BC and the coagulation between hydrophobic BC and preexisting soluble particles. The chemical oxidation aging is tested in the sensitivity simulation. The microphysical aging rate is ∼ 4 times higher in the lower troposphere over source regions than that from a fixed aging scheme with an e-folding time of 1.2 days. The higher aging rate reflects the large emissions of sulfate–nitrate and secondary organic aerosol precursors hence faster BC aging through condensation and coagulation. In contrast, the microphysical aging is more than 5-fold slower than the fixed aging in remote regions, where condensation and coagulation are weak. Globally, BC microphysical aging is dominated by condensation, while coagulation contribution is largest over eastern China, India, and central Africa. The fixed aging scheme results in an overestimate of HIPPO BC throughout the troposphere by a factor of 6 on average. The microphysical scheme reduces this discrepancy by a factor of ∼ 3, particularly in the middle and upper troposphere. It also leads to a 3-fold reduction in model bias in the latitudinal BC column burden averaged along the HIPPO flight tracks, with largest improvements in the tropics. The resulting global annual mean BC lifetime is 4.2 days and BC burden is 0.25 mg m−2, with 7.3 % of the burden at high altitudes (above 5 km). Wet scavenging accounts for 80.3 % of global BC deposition. We find that, in source regions, the microphysical aging rate is insensitive to aerosol size distribution, condensation threshold, and chemical oxidation aging, while it is the opposite in remote regions, where the aging rate is orders of magnitude smaller. As a result, global BC burden and lifetime show little sensitivity (< 5 % change) to these three factors.

A theoretical black carbon (BC) aging model is developed to account for three typical evolution stages, namely, freshly emitted aggregates, BC coated by soluble material, and BC particles undergoing further hygroscopic growth. The geometric-optics surface-wave (GOS) approach is employed to compute the BC single-scattering properties at each aging stage, which are subsequently compared with laboratory measurements. Theoretical calculations are consistent with measurements in extinction and absorption cross sections for fresh BC aggregates with different BC sizes (i.e., mobility diameters of 155, 245, and 320 nm), with differences of ≤ 25 %. The measured optical cross sections for BC coated by sulfuric acid and for that undergoing further hygroscopic growth are generally captured (differences < 30 %) by theoretical calculations using a concentric core-shell structure, with an overestimate in extinction and absorption of the smallest BC size and an underestimate in scattering of the largest BC size. We find that the absorption and scattering cross sections of fresh BC aggregates vary by 20–40 and 50–65 %, respectively, due to the use of upper (1.95–0.79i) and lower (1.75–0.63i) bounds of BC refractive index, while the variations are < 20 % in absorption and < 50 % in scattering in the case of coated BC particles. Sensitivity analyses of the BC morphology show that the optical properties of fresh BC aggregates are more sensitive to fractal dimension than primary spherule size. The absorption and scattering cross sections of coated BC particles vary by more than a factor of 2 due to different coating structures. We find an increase of 20–250 % in absorption and a factor of 3–15 in scattering during aging, significantly depending on coating morphology and aging stages. This study suggests that an accurate estimate of BC radiative effects requires the incorporation of a dynamic BC aging process that accounts for realistic coating structures in climate models.

We estimate black carbon (BC) emissions in the western United States for July–September 2006 by inverting surface BC concentrations from the Interagency Monitoring of Protected Visual Environments (IMPROVE) network using a global chemical transport model (GEOS-Chem) and its adjoint. Our best estimate of the BC emissions is 49.9 Gg at 2° × 2.5° (a factor of 2.1 increase) and 47.3 Gg at 0.5° × 0.667° (1.9 times increase). Model results now capture the observed major fire episodes with substantial bias reductions (~ 35 % at 2° × 2.5° and ~ 15 % at 0.5° × 0.667°). The emissions are ~ 20–50 % larger than those from our earlier analytical inversions (Mao et al., 2014). The discrepancy is especially drastic in the partitioning of anthropogenic versus biomass burning emissions. The August biomass burning BC emissions are 4.6–6.5 Gg and anthropogenic BC emissions 8.6–12.8 Gg, varying with the model resolution, error specifications, and subsets of observations used. On average both anthropogenic and biomass burning emissions in the adjoint inversions increase 2-fold relative to the respective {a priori} emissions, in distinct contrast to the halving of the anthropogenic and tripling of the biomass burning emissions in the analytical inversions. We attribute these discrepancies to the inability of the adjoint inversion system, with limited spatiotemporal coverage of the IMPROVE observations, to effectively distinguish collocated anthropogenic and biomass burning emissions on model grid scales. This calls for concurrent measurements of other tracers of biomass burning and fossil fuel combustion (e.g., carbon monoxide and carbon isotopes). We find that the adjoint inversion system as is has sufficient information content to constrain the total emissions of BC on the model grid scales.

We estimate biomass burning and anthropogenic emissions of black carbon (BC) in the western US for May–October 2006 by inverting surface BC concentrations from the Interagency Monitoring of PROtected Visual Environment (IMPROVE) network using a global chemical transport model. We first use active fire counts from the Moderate Resolution Imaging Spectroradiometer (MODIS) to improve the spatiotemporal distributions of the biomass burning BC emissions from the Global Fire Emissions Database (GFEDv2). The adjustment primarily shifts emissions from late to middle and early summer (a 33% decrease in September–October and a 56% increase in June–August) and leads to appreciable increases in modeled surface BC concentrations in early and middle summer, especially at the 1–2 and 2–3 km altitude ranges. We then conduct analytical inversions at both 2° × 2.5° and 0.5° × 0.667° (nested over North America) horizontal resolutions. The a posteriori biomass burning BC emissions for July–September are 31.7 Gg at 2° × 2.5° (an increase by a factor of 4.7) and 19.2 Gg at 0.5° × 0.667° (an increase by a factor of 2.8). The inversion results are rather sensitive to model resolution. The a posteriori biomass burning emissions at the two model resolutions differ by a factor of ~6 in California and the Southwest and by a factor of 2 in the Pacific Northwest. The corresponding a posteriori anthropogenic BC emissions are 9.1 Gg at 2° × 2.5° (a decrease of 48%) and 11.2 Gg at 0.5° × 0.667° (a decrease of 36%). Simulated surface BC concentrations with the a posteriori emissions capture the observed major fire episodes at most sites and the substantial enhancements at the 1–2 and 2–3 km altitude ranges. The a posteriori emissions also lead to large bias reductions (by ~30% on average at both model resolutions) in modeled surface BC concentrations and significantly better agreement with observations (increases in Taylor skill scores of 95% at 2° × 2.5° and 42 % at 0.5° × 0.667°).

Satellite observations show that ice cloud effective radius (re) increases with ice water content (IWC) but decreases with aerosol optical thickness (AOT). Using least-squares fitting to the observed data, we obtain an analytical formula to describe the variations of re with IWC and AOT for several regions with distinct characteristics of re-IWC-AOT relationships. As IWC directly relates to convective strength and AOT represents aerosol loading, our empirical formula provides a means to quantify the relative roles of dynamics and aerosols in controlling re in different geographical regions, and to establish a framework for parameterization of aerosol effects on re in climate models.