Modeling Whitecaps on Global Scale

Whitecaps play an important role in the surface-atmosphere interactions across the ocean. They are directly linked to the energy dissipation rate during wave breaking and transfer of heat, momentum, and gas/aerosol exchange from the wind to the sea surface. Although the first models of W were dependent only on wind speeds, a large number of diverse models based on wind and sea state which include wave height, wave age, friction velocity, and stability effects have been proposed since then. However, it is recognized that most of the proposed W models have strong systematic (e.g., zonal bias) and random errors when compared against observations. This is partly due to the differences in environmental conditions, measurement techniques, and geographical locations among these studies. But, some of these biases are linked to the inability of the proposed models to capture the variability in W in certain wind/wave regimes. Despite the knowledge of existing biases, W residual relationships from the models with wind and wave fields remain highly uncertain, with residual trends varying between the published studies. Here, we take advantage of the availability of relatively dense observations of W from WindSat microwave satellite retrievals in combination with the University of Miami wave model which was recently incorporated within the NASA GMAO/GEOS system (GEOS-UMWM). We use Windsat W retrievals to assess and constrain the previously published W models and understand the relationships of residuals from models in different wind/wave regimes. We link these unexplained residual variations to additional factors such as swell index, drag coefficient etc and add information to the existing whitecap models. Since Windsat retrievals cover wide range of environmental conditions, it helps to reduce the uncertainties associated with differences in measurement techniques. Regression of wind-wave fields against all Windsat data points (CTL) results in larger residuals for lower wave age and W is overestimated upto ~4% for wave age < 10 and underestimated by upto ~2% as wave age increases. We attest to this bias by considering two approaches. One is to perform regression separately for different stages of wave development such as developing sea, fully developed, and wind sea regimes thereby understanding the sensitivity of regression coefficients to sea state (EXP1). Another is to derive coefficients of W models in EXP1 as a function of additional wind/wave factors such as swell index, drag coefficient, and mean squared slope, deriving more nonlinear W models (EXP2). EXP2 provides reduction in Root Mean Squared Error (RMSE) by 0.1-0.3%. Sea surface drag has a stronger relationship with regression coefficients compared to swell index.These additional factors provide improved parameterizations in different wind and wave age regimes, with smaller unexplained/residual variations in W that has been a major concern in the W community

The Quick Fire Emissions Dataset (QFED): Documentation of Versions 2.1, 2.2 and 2.4

Biomass burning is an important source of particulates and trace gases and a major element of the terrestrial carbon cycle. Well constrained emissions from vegetation fires are needed to model direct and indirect effects of biomass burning aerosols, to model homogeneous and heterogeneous chemistry in the atmosphere, and to perform credible Earth system analysis, and climate and air pollution studies. To improve the performance of NASA Goddard Earth Observing System Model (GEOS) in the areas of atmospheric constituent modeling with a focus on biomass burning we developed the Quick Fire Emissions Dataset (QFED). The QFED emissions are based on the fire radiative power (top-down) approach and draw on the cloud correction method developed in the Global Fire Assimilation System (GFAS). Location and fire radiative power of fires are obtained from the Moderate Resolution Imaging Spectroradiometer (MODIS) Level 2 fire products (MOD14 and MYD14) and the MODIS Geolocation products (MOD03 and MYD03). QFED strengths are high spatial and temporal resolutions and near-real time availability. Daily mean emissions are available at 0.3125 times 0.25 degrees and in recent versions also at 0.1 times 0.1 degrees. QFED provides emissions of black carbon, organic carbon, sulfur dioxide, carbon monoxide, carbon dioxide, PM2.5, ammonia, nitrogen oxides, methyl ethyl ketone, propylene, ethane, propane, n- and i-butane, acetaldehyde, formaldehyde, acetone and methane. Two QFED product systems are maintained by the NASA Global Modeling and Assimilation Office (GMAO): one that produces near real-time daily emissions used operationally in the GEOS-5 Data Assimilation System, and one that produces an extended historical dataset with daily emissions from March 2000 to the present. The historical dataset also provides monthly mean and monthly climatological emissions

Challenges in Developing Better Observational Constraints and Models for Aerosols : Emerging Ideas for Design and Use of Future Observing Systems

No abstract availabl

Sea State Based Estimation of White Cap Fraction: Implications for Primary Marine Aerosol Fluxes

Oceanic whitecaps (hereafter, W) or the characteristic whiteness of the sea foam is an important feature for predicting exchange of gases, sea spray aerosols (SSAs), heat and momentum transfer between the ocean and the atmosphere at the air-sea interface. Due to its increased surface emission and brightness temperature, whitecaps are critical for satellite retrievals of ocean albedo, ocean color, ocean surface wind vectors from satellite borne radiometer and microwave instruments. Most of the existing models predict W using wind speed and sea surface temperature (SST). However, numerous publications have pointed out that there are large uncertainties in the predicted W and using parameterizations based on wind-wave state can improve the precision of the predicted W. Here, we integrate the University of Miami Wave Model - 2.0 (UMWM) in Goddard Earth Observing System (GEOS) and use wave diagnostics to predict W. We choose the year 2006 for our global UMWM/GEOS runs because of the availability of W dataset from satellite observations. We run UMWM/GEOS at 0.5o x 0.5o by replaying to MERRA2 meteorology and evaluate the wave diagnostics using measurements from fixed buoys and satellite altimeters. We use three different parameterizations for W based on: 1) Reynolds number, 2) wave dissipation energy, and 3) volume of air entrained by breaking waves. We compare our results of W with previous studies and also with the satellite based observational dataset. Predicting W is important for understanding the processes at the air-sea interface. Therefore, this work is a step further in improving the uncertainties in the aerosol and atmospheric chemistry modules of the global models

A New Global Anthropogenic SO2 Emission Inventory for the Last Decade: A Mosaic of Satellite-Derived and Bottom-Up Emissions

Sulfur dioxide (SO2) measurements from the Ozone Monitoring Instrument (OMI) satellite sensor have been used to detect emissions from large point sources. Emissions from over 400 sources have been quantified individually based on OMI observations, accounting for about a half of total reported anthropogenic SO2 emissions. Here we report a newly developed emission inventory, OMI-HTAP, by combining these OMI-based emission estimates and the conventional bottom-up inventory, HTAP, for smaller sources that OMI is not able to detect. OMI-HTAP includes emissions from OMI-detected sources that are not captured in previous leading bottom-up inventories, enabling more accurate emission estimates for regions with such missing sources. In addition, our approach offers the possibility of rapid updates to emissions from large point sources that can be detected by satellites. Our methodology applied to OMI-HTAP can also be used to merge improved satellite-derived estimates with other multi-year bottom-up inventories, which may further improve the accuracy of the emission trends. OMI-HTAP SO2 emissions estimates for Persian Gulf, Mexico, and Russia are 59%, 65%, and 56% larger than HTAP estimates, respectively, in year 2010. We have evaluated the OMI-HTAP inventory by performing simulations with the Goddard Earth Observing System version 5 (GEOS-5) model. The GEOS-5 simulated SO2 concentrations driven by both HTAP and OMI-HTAP were compared against in situ measurements. We focus for the validation on year 2010 for which HTAP is most valid and for which a relatively large number of in situ measurements are available. Results show that the OMI-HTAP inventory improves the agreement between the model and observations, in particular over the US, with the normalized mean bias decreasing from 0.41 (HTAP) to -0.03 (OMI-HTAP) for year 2010. Simulations with the OMI-HTAP inventory capture the worldwide major trends of large anthropogenic SO2 emissions that are observed with OMI. Correlation coefficients of the observed and modelled surface SO2 in 2014 increase from 0.16 (HTAP) to 0.59 (OMI-HTAP) and the normalized mean bias dropped from 0.29 (HTAP) to 0.05 (OMI-HTAP), when we updated 2010 HTAP emissions with 2014 OMI-HTAP emissions in the model

GEOS-5 Aerosol Modeling Data Assimilation: Update on Recent and Future Development

No abstract availabl

Current and Future Applications of the GEOS-5 Aerosol Modeling System

The presentation summarizes current and proposed activities for the GEOS-5 aerosol modeling system. Activities discussed include (i) forecasting and event simulation, (ii) observation simulation, (iii) aerosol-chemistry-climate applications, and (iv) future activities. The document was presented at the 2013 AEROCENTER Annual Meeting held at the GSFC Visitors Center May 31, 2013. The Organizers of the meeting are posting the talks to the public Aerocenter website, after the meeting

Mass concentration estimates of long-range-transported Canadian biomass burning aerosols from a multi-wavelength Raman polarization lidar and a ceilometer in Finland

A quantitative comparison study for Raman lidar and ceilometer observations, and for model simulations of mass concentration estimates of smoke particles is presented. Layers of biomass burning aerosol particles were observed in the lower troposphere, at 2 to 5km height on 4 to 6 June 2019, over Kuopio, Finland. These long-range-transported smoke particles originated from a Canadian wildfire event. The most pronounced smoke plume detected on 5 June was intensively investigated. Optical properties were retrieved from the multi-wavelength Raman polarization lidar PollyXT. Particle linear depolarization ratios (PDRs) of this plume were measured to be 0.08 +/- 0.02 at 355nm and 0.05 +/- 0.01 at 532nm, suggesting the presence of partly coated soot particles or particles that have mixed with a small amount of dust or other non-spherical aerosol type. The layer-mean PDR at 355nm (532nm) decreased during the day from similar to 0.11 (0.06) in the morning to similar to 0.05 (0.04) in the evening; this decrease with time could be linked to the particle aging and related changes in the smoke particle shape properties. Lidar ratios were derived as 47 +/- 5sr at 355nm and 71 +/- 5sr at 532nm. A complete ceilometer data processing for a Vaisala CL51 ceilometer is presented from a sensor-provided attenuated backscatter coefficient to particle mass concentration (including the water vapor correction for high latitude for the first time). Aerosol backscatter coefficients (BSCs) were measured at four wavelengths (355, 532, 1064nm from PollyXT and 910nm from CL51). Two methods, based on a combined lidar and sun-photometer approach, are applied for mass concentration estimations from both PollyXT and the ceilometer CL51 observations. In the first method, no. 1, we used converted BSCs at 532nm (from measured BSCs) by corresponding measured backscatter-related angstrom ngstrom exponents, whereas in the second method, no. 2, we used measured BSCs at each wavelength independently. A difference of similar to 12% or similar to 36% was found between PollyXT and CL51 estimated mass concentrations using method no. 1 or no. 2, showing the potential of mass concentration estimates from a ceilometer. Ceilometer estimations have an uncertainty of similar to 50% in the mass retrieval, but the potential of the data lies in the great spatial coverage of these instruments. The mass retrievals were compared with the Modern-Era Retrospective analysis for Research and Applications, version 2 (MERRA-2) meteorological and aerosol reanalysis. The inclusion of dust (as indicated by MERRA-2 data) in the retrieved mass concentration is negligible considering the uncertainties, which also shows that ceilometer observations for mass retrievals can be used even without exact knowledge of the composition of the smoke-dominated aerosol plume in the troposphere

A new global anthropogenic SO2 emission inventory for the last decade: A mosaic of satellite-derived and bottom-up emissions

Sulfur dioxide (SO2) measurements from the Ozone Monitoring Instrument (OMI) satellite sensor have been used to detect emissions from large point sources. Emissions from over 400 sources have been quantified individually based on OMI observations, accounting for about a half of total reported anthropogenic SO2 emissions. Here we report a newly developed emission inventory, OMI-HTAP, by combining these OMI-based emission estimates and the conventional bottom-up inventory, HTAP, for smaller sources that OMI is not able to detect. OMI-HTAP includes emissions from OMI-detected sources that are not captured in previous leading bottom-up inventories, enabling more accurate emission estimates for regions with such missing sources. In addition, our approach offers the possibility of rapid updates to emissions from large point sources that can be detected by satellites. Our methodology applied to OMI-HTAP can also be used to merge improved satellite-derived estimates with other multi-year bottom-up inventories, which may further improve the accuracy of the emission trends. OMI-HTAP SO2 emissions estimates for Persian Gulf, Mexico, and Russia are 59 %, 65 %, and 56% larger than HTAP estimates in 2010, respectively. We have evaluated the OMI-HTAP inventory by performing simulations with the Goddard Earth Observing System version 5 (GEOS-5) model. The GEOS-5 simulated SO2 concentrations driven by both HTAP and OMI-HTAP were compared against in situ measurements. We focus for the validation on 2010 for which HTAP is most valid and for which a relatively large number of in situ measurements are available. Results show that the OMI-HTAP inventory improves the agreement between the model and observations, in particular over the US, with the normalized mean bias decreasing from 0.41 (HTAP) to -0:03 (OMIHTAP) for 2010. Simulations with the OMI-HTAP inventory capture the worldwide major trends of large anthropogenic SO2 emissions that are observed with OMI. Correlation coefficients of the observed and modeled surface SO2 in 2014 increase from 0.16 (HTAP) to 0.59 (OMI-HTAP) and the normalized mean bias dropped from 0.29 (HTAP) to 0.05 (OMIHTAP), when we updated 2010 HTAP emissions with 2014 OMI-HTAP emissions in the model.JRC.D.6-Knowledge for Sustainable Development and Food Securit