Search NASA⌕ Search

Engineering topics

da Silva, Arlindo M.

Publications and source records attributed to da Silva, Arlindo M..

A Geostationary Instrument Simulator for Aerosol Observing System Simulation Experiments

In the near future, there will be several new instruments measuring atmospheric composition from geostationary orbit over North America, East Asia, and Europe. This constellation of satellites will provide high resolution, time resolved measurements of trace gases and aerosols for monitoring air quality and tracking pollution sources. This paper describes a detailed, fast, and accurate (less than 1.0% uncertainty) method for calculating synthetic top of the atmosphere (TOA) radiances from a global simulation with a mesoscale free running model, the GEOS-5 Nature Run, for remote sensing instruments in geostationary orbit that measure in the ultraviolet-visible spectral range (UV-Vis). Generating these synthetic observations is the first step of an Observing System Simulation Experiment (OSSE), a framework for evaluating the impact of a new observation or algorithm. This paper provides details of the model sampling, aerosol and cloud optical properties, surface reflectance modeling, Rayleigh scattering calculations, and a discussion of the uncertainties of the simulated TOA radiance. An application for the simulated TOA radiance observations is demonstrated in the manuscript. Simulated TEMPO (Tropospheric Emissions: Monitoring of Pollution) and GOES-R (Geostationary Operational Environmental Satellites) observations were used to show how observations from the two instruments could be combined to facilitate aerosol type discrimination. The results demonstrate the viability of a detailed instrument simulator for radiance measurements in the UV-Vis that is capable of accurately simulating high resolution, time-resolved measurements with reasonable computational efficiency.

Castellanos, Patricia↗

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.

Sulfur dioxide (SO2)↗

Preliminary Evaluation of Influence of Aerosols on the Simulation of Brightness Temperature in the NASA's Goddard Earth Observing System Atmospheric Data Assimilation System

This document reports on preliminary results obtained when studying the impact of aerosols on the calculation of brightness temperature (BT) for satellite infrared (IR) instruments that are currently assimilated in a 3DVAR configuration of Goddard Earth Observing System (GEOS)-atmospheric data assimilation system (ADAS). A set of fifteen aerosol species simulated by the Goddard Chemistry Aerosol Radiation and Transport (GOCART) model is used to evaluate the influence of the aerosol fields on the Community Radiative Transfer Model (CRTM) calculations taking place in the observation operators of the Gridpoint Statistical Interpolation (GSI) analysis system of GEOSADAS. Results indicate that taking aerosols into account in the BT calculation improves the fit to observations over regions with significant amounts of dust. The cooling effect obtained with the aerosol-affected BT leads to a slight warming of the analyzed surface temperature (by about 0:5oK) in the tropical Atlantic ocean (off northwest Africa), whereas the effect on the air temperature aloft is negligible. In addition, this study identifies a few technical issues to be addressed in future work if aerosol-affected BT are to be implemented in reanalysis and operational settings. The computational cost of applying CRTM aerosol absorption and scattering options is too high to justify their use, given the size of the benefits obtained. Furthermore, the differentiation between clouds and aerosols in GSI cloud detection procedures needs satisfactory revision.

3DVAR↗

Aerosol Data Assimilation at GMAO

This presentation presents an overview of the aerosol data assimilation work performed at GMAO. The GMAO Forward Processing system and the biomass burning emissions from QFED are first presented. Then, the current assimilation of Aerosol Optical Depth (AOD), performed by means of the analysis splitting method is briefly described, followed by some results on the quality control of observations using a Neural Network trained using AERONET AOD. Some applications are shown such as the Mount Pinatubo eruption in 1991 using the MERRA-2 aerosol dataset. Finally preliminary results on the EnKF implementation for aerosol assimilation are presented.

Splitting↗

Can Aerosol Direct Radiative Effects Account for Analysis Increments of Temperature in the Tropical Atlantic?

In the late 1990's, prior to the launch of the Terra satellite, atmospheric general circulation models (GCMs) did not include aerosol processes because aerosols were not properly monitored on a global scale and their spatial distributions were not known well enough for their incorporation in operational GCMs. At the time of the first GEOS Reanalysis (Schubert et al. 1993), long time series of analysis increments (the corrections to the atmospheric state by all available meteorological observations) became readily available, enabling detailed analysis of the GEOS-1 errors on a global scale. Such analysis revealed that temperature biases were particularly pronounced in the Tropical Atlantic region, with patterns depicting a remarkable similarity to dust plumes emanating from the African continent as evidenced by TOMS aerosol index maps. Yoram Kaufman was instrumental encouraging us to pursue this issue further, resulting in the study reported in Alpert et al. (1998) where we attempted to assess aerosol forcing by studying the errors of a the GEOS-1 GCM without aerosol physics within a data assimilation system. Based on this analysis, Alpert et al. (1998) put forward that dust aerosols are an important source of inaccuracies in numerical weather-prediction models in the Tropical Atlantic region, although a direct verification of this hypothesis was not possible back then. Nearly 20 years later, numerical prediction models have increased in resolution and complexity of physical parameterizations, including the representation of aerosols and their interactions with the circulation. Moreover, with the advent of NASA's EOS program and subsequent satellites, atmospheric aerosols are now monitored globally on a routine basis, and their assimilation in global models are becoming well established. In this talk we will reexamine the Alpert et al. (1998) hypothesis using the most recent version of the GEOS-5 Data Assimilation System with assimilation of aerosols. We will explicitly calculate the impact of aerosols on the temperature analysis increments in the tropical Atlantic and assess the extent to which inclusion of atmospheric aerosols have reduced these increments.

General Circulation Models↗

The Quick Fire Emissions Dataset (QFED): Documentation of Versions 2.1, 2.2 and 2.4: Technical Report Series on Global Modeling and Data Assimilation - Volume 38

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.

Koster, Randal D.↗

File Specification for the 7-km GEOS-5 Nature Run, Ganymed Release Non-Hydrostatic 7-km Global Mesoscale Simulation

This document describes the gridded output files produced by a two-year global, non-hydrostatic mesoscale simulation for the period 2005-2006 produced with the non-hydrostatic version of GEOS-5 Atmospheric Global Climate Model (AGCM). In addition to standard meteorological parameters (wind, temperature, moisture, surface pressure), this simulation includes 15 aerosol tracers (dust, sea-salt, sulfate, black and organic carbon), O3, CO and CO2. This model simulation is driven by prescribed sea-surface temperature and sea-ice, daily volcanic and biomass burning emissions, as well as high-resolution inventories of anthropogenic sources. A description of the GEOS-5 model configuration used for this simulation can be found in Putman et al. (2014). The simulation is performed at a horizontal resolution of 7 km using a cubed-sphere horizontal grid with 72 vertical levels, extending up to to 0.01 hPa (approximately 80 km). For user convenience, all data products are generated on two logically rectangular longitude-latitude grids: a full-resolution 0.0625 deg grid that approximately matches the native cubed-sphere resolution, and another 0.5 deg reduced-resolution grid. The majority of the full-resolution data products are instantaneous with some fields being time-averaged. The reduced-resolution datasets are mostly time-averaged, with some fields being instantaneous. Hourly data intervals are used for the reduced-resolution datasets, while 30-minute intervals are used for the full-resolution products. All full-resolution output is on the model's native 72-layer hybrid sigma-pressure vertical grid, while the reduced-resolution output is given on native vertical levels and on 48 pressure surfaces extending up to 0.02 hPa. Section 4 presents additional details on horizontal and vertical grids. Information of the model surface representation can be found in Appendix B. The GEOS-5 product is organized into file collections that are described in detail in Appendix C. Additional details about variables listed in this file specification can be found in a separate document, the GEOS-5 File Specification Variable Definition Glossary. Documentation about the current access methods for products described in this document can be found on the GEOS-5 Nature Run portal: http://gmao.gsfc.nasa.gov/projects/G5NR. Information on the scientific quality of this simulation will appear in a forthcoming NASA Technical Report Series on Global Modeling and Data Assimilation to be available from http://gmao.gsfc.nasa.gov/pubs/tm/.

GEOS-5↗

Evaluating the Impact of Aerosols on Numerical Weather Predictions

This presentation focused on proposed Aerosol related case studies, and details were given about participating centers and modeling systems, as well as webpage and analysis tools. Preliminary results were given with regards to a case study on the dust storm over Egypt and another case study on the extreme pollution event in China.

Aerosol↗

Monte Carlo Bayesian Inference on a Statistical Model of Sub-gridcolumn Moisture Variability Using High-resolution Cloud Observations : Sensitivity Tests and Results - Part II

Part I presented a Monte Carlo Bayesian method for constraining a complex statistical model of GCM sub-gridcolumn moisture variability using high-resolution MODIS cloud data, thereby permitting large-scale model parameter estimation and cloud data assimilation. This part performs some basic testing of this new approach, verifying that it does indeed significantly reduce mean and standard deviation biases with respect to the assimilated MODIS cloud optical depth, brightness temperature and cloud top pressure, and that it also improves the simulated rotational-Ramman scattering cloud optical centroid pressure (OCP) against independent (non-assimilated) retrievals from the OMI instrument. Of particular interest, the Monte Carlo method does show skill in the especially difficult case where the background state is clear but cloudy observations exist. In traditional linearized data assimilation methods, a subsaturated background cannot produce clouds via any infinitesimal equilibrium perturbation, but the Monte Carlo approach allows finite jumps into regions of non-zero cloud probability. In the example provided, the method is able to restore marine stratocumulus near the Californian coast where the background state has a clear swath. This paper also examines a number of algorithmic and physical sensitivities of the new method and provides guidance for its cost-effective implementation. One obvious difficulty for the method, and other cloud data assimilation methods as well, is the lack of information content in the cloud observables on cloud vertical structure, beyond cloud top pressure and optical thickness, thus necessitating strong dependence on the background vertical moisture structure. It is found that a simple flow-dependent correlation modification due to Riishojgaard (1998) provides some help in this respect, by better honoring inversion structures in the background state.

Monte Carlo Bayesian inference↗