Search NASA⌕ Search

SEARCH · Search NASA

Results for “model reanalysis”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 217 records · Page 12

Near Real Time Air Quality Forecasts Using the NASA GEOS Model

This presentation provides an overview of NASA's Global Modeling and Assimilation Office (GMAO) high-resolution global forecast and reanalysis products for weather, aerosols, and air quality. The NASA Global Earth Observing System (GEOS) model radiatively coupled to GOCART aerosol module assimilates 2-dimensional column-integrated aerosol optical depth (AOD) at one wavelength (550 nm) in order to constrain the model's background AOD in order to have the best historical estimate and forecasts of AOD and particulate matter. Furthermore, the GEOS model has been expanded to provide global near-real-time 5-day forecasts of atmospheric composition at unprecedented horizontal resolution of 0.25 degrees (~25 km). This composition forecast system (GEOS-CF) combines the operational GEOS weather forecasting model with the state-of-the-science GEOS-Chem chemistry module (version 12) to provide detailed analysis of a wide range of air pollutants such as ozone, carbon monoxide, nitrogen oxides, and fine particulate matter (PM2.5).

K Emma Knowland↗

Soil Moisture Active Passive (SMAP) Project Assessment Report for Version 4 of the L4_SM Data Product

This report provides an assessment of Version 4 of the SMAP Level 4 Surface and Root Zone Soil Moisture (L4_SM) product, released on 14 June 2018. The assessment includes comparisons of L4_SM soil moisture and temperature estimates with in situ measurements from core validation sites and sparse networks. The assessment further includes a global evaluation of the internal diagnostics from the ensemble-based data assimilation system that is used to generate the L4_SM product, including observation-minus-forecast (O-F) brightness temperature residuals and soil moisture analysis increments.Together, the core validation site comparisons and the statistics of the assimilation diagnostics areconsidered primary validation methodologies for the L4_SM product. Comparisons against in situ measurements from regional-scale sparse networks are considered a secondary validation methodology because such in situ measurements are subject to upscaling errors from the point-scale to the grid-cell scale of the data product.The Version 4 L4_SM product benefits from an improved land surface modeling system and from retrospective surface meteorological forcing data that are as consistent as possible with the present-day datain terms of their climatology. Specifically, the model changes include revised parameters and parameterizations for (i) the surface energy balance, (ii) recharge from below of the model's surface excess reservoir, and (iii) the snow depletion curve. Updated ancillary inputs include improved datasets for landcover, topography, and vegetation height. The Version 4 algorithm further includes a revised approach to precipitation corrections that improves the precipitation climatology in Africa and the high-latitudes. Moreover, for system calibration the model is forced retrospectively with MERRA-2 reanalysis data, which are more consistent with the near-real time GEOS forward processing (FP) data used during the SMAP period than the retrospective GEOS data that were available for previous L4_SM versions. An analysis of the time-average surface and root zone soil moisture shows that the global pattern ofarid and humid regions is captured by the Version 4 L4_SM estimates. Owing to the changes in the landsurface modeling system, surface soil moisture is typically drier by several volumetric percent in Version 4 compared to Version 3, whereas root zone soil moisture is wetter in Version 4 in some regions and drierin others. Because of these climatological differences, the Version 3 and Version 4 products should not be combined into a single dataset for use in applications.Results from the core validation site comparisons indicate that Version 4 of the L4_SM data product meets the self-imposed L4_SM accuracy requirement, which is formulated in terms of the RMSE after removal of the long-term mean difference (ubRMSE). The overall ubRMSE of the 3-hourly L4_SM dataat the 9 km scale is 0.039 m3 m-3 for surface soil moisture and 0.029 m3 m-3 for root zone soil moisture,below the 0.04 m3 m-3 requirement. The L4_SM estimates are an improvement over estimates from a model-only Nature Run version 7.2 (NRv7.2), which demonstrates the beneficial impact of the SMAP brightness temperature data. Overall, L4_SM surface and root zone soil moisture estimates are more skillful than NRv7.2 estimates, with statistically significant improvements at the 5% level for surface soil moisture R and anomaly R values. Results from comparisons of the L4_SM product to i

Reichle, Rolf H.↗

Overview of the SPARC Reanalysis Intercomparison Project (S-RIP) during 2013-2020

The Stratosphere-troposphere Processes And their Role in Climate (SPARC) project is one of the four core projects of the World Climate Research Programme (WCRP). Researchers interested in SPARC use global atmospheric reanalysis products to understand a wide range of processes and variability in the atmosphere, to validate chemistry climate models, and to investigate and identify climate change. The SPARC Reanalysis Intercomparison Project (S-RIP) was initiated in 2011 and officially started in 2013 to conduct a coordinated intercomparison of all major global atmospheric reanalysis data sets. The S-RIP has been aiming at writing up an assessment report in the SPARC report series (to be published by the end of 2020) (1) on overall quality of temperature, winds, ozone, and water vapor data, (2) on more process- and region-oriented evaluation of the Brewer–Dobson circulation, extratropical stratosphere-troposphere coupling, extratropical upper troposphere and lower stratosphere, the tropical tropopause layer, the quasi-biennial oscillation, polar processes, and the upper stratosphere and lower mesosphere, and (3) with a coordinated description of the reanalysis systems. We also have an inter-journal special issue on "The SPARC Reanalysis Intercomparison Project (S-RIP)" in Atmospheric Chemistry and Physics (ACP) and Earth System Science Data (ESSD). In the presentation, we will discuss key findings and recommendations as well as the evaluation of this rst phase of the S-RIP activity.

Masatomo Fujiwara↗

Long-Term Ozone Variability and Trends from Reanalyses: Can It Be Done?

Stratospheric ozone concentrations have begun to show early signs of recovery following the implementation of the Montreal Protocol and its amendments as well as in response to decreasing upper-stratospheric temperatures. Secular trends in stratospheric ozone are modulated by considerable interannual variability and systematic changes in transport patterns that are expected under increasing concentrations of greenhouse gases, especially in the lower stratosphere. These factors necessitate the continued close monitoring of stratospheric ozone in upcoming decades, with a special focus on the lower stratosphere.As highly resolved data sets combining a plethora of observations with model simulations atmospheric reanalyses are, in principle, well suited for the task. All major reanalyses generate ozone output. However, significant spurious discontinuities that arise from step changes in the observing systems prevent a straightforward analysis of ozone trends and long-term variability. Building on our recent work, in this presentation we will demonstrate that trend detection is nonetheless possible using the ozone record from NASA's MERRA-2 (Modern-Era Retrospective Analysis for Research and Applications, Version 2) reanalysis bias-corrected using a chemistry model simulation as a transfer function. Next, we will outline several strategies to reduce artificial discontinuities in the ozone record in future NASA reanalyses. This discussion will be illustrated by an example of joint assimilation of bias-corrected ozone profiles from the Microwave Limb Sounder (MLS) on the Aura satellite (2004 to present) and the Ozone Mapping Profiler Suite Limb Profiler (OMPS-LP) sensors that are expected to operate on future NOAA platforms.

Wargan, Krzysztof↗

NASA's High-Resolution GEOS Forecasting and Reanalysis Products: A Unified Tool from Local to Global Scales

NASA's GMAO produces high-resolution global forecasts for weather, aerosols, and air quality. The NASA Global Earth Observing System (GEOS) model has been expanded to provide global near-real-time 5-day forecasts of atmospheric chemical composition at unprecedented horizontal resolution of 0.25 degrees (~25 km), freely available to the public. This composition forecast system (GEOS-CF) combines the operational GEOS weather forecasting model with the state-of-the-science GEOS-Chem chemistry module to provide detailed analysis of a wide range of air pollutants such as ozone, carbon monoxide, nitrogen oxides, and fine particulate matter (PM2.5). GEOS-CF also assimilated satellite observations into the system for improved representation of weather and smoke. The assimilation system is currently being expanded to include chemically reactive trace gases. While the main focus of this new product is on tropospheric air quality information, the GEOS-Chem chemistry model used in this system includes the unified tropospheric stratospheric chemistry mechanism for improved forecasts of total column ozone during anomalous dynamical and chemical events. I will discuss current capabilities of the GEOS Constituent Data Assimilation System (CoDAS) to improve atmospheric composition modeling and possible future directions for GEOS-CF and reanalysis products. In addition, I will show how machine learning techniques can be used to correct for sub-grid-scale variability, which further improves model estimates at a given observation site.

Co-DAS↗

Theoretical and software considerations for general dynamic analysis using multilevel substructured models

The dynamic analysis of complex structural systems using the finite element method and multilevel substructured models is presented. The fixed-interface method is selected for substructure reduction because of its efficiency, accuracy, and adaptability to restart and reanalysis. This method is extended to reduction of substructures which are themselves composed of reduced substructures. The implementation and performance of the method in a general purpose software system is emphasized. Solution algorithms consistent with the chosen data structures are presented. It is demonstrated that successful finite element software requires the use of software executives to supplement the algorithmic language. The complexity of the implementation of restart and reanalysis porcedures illustrates the need for executive systems to support the noncomputational aspects of the software. It is shown that significant computational efficiencies can be achieved through proper use of substructuring and reduction technbiques without sacrificing solution accuracy. The restart and reanalysis capabilities and the flexible procedures for multilevel substructured modeling gives economical yet accurate analyses of complex structural systems.

Schmidt, R. J.↗

Land Surface Precipitation and Hydrology in MERRA-2

The Modern-Era Retrospective Analysis for Research and Applications, version 2 (MERRA-2), provides global, 1-hourly estimates of land surface conditions for 1980-present at 50-km resolution. Outside of the high latitudes, MERRA-2 uses observations-based precipitation data products to correct the precipitation falling on the land surface. This paper describes the precipitation correction method and evaluates the MERRA-2 land surface precipitation and hydrology. Compared to monthly GPCPv2.2 observations, the corrected MERRA-2 precipitation (M2CORR) is better than the precipitation generated by the atmospheric models within the cyclingMERRA-2 system and the earlier MERRA reanalysis. Compared to 3-hourlyTRMM observations, the M2CORR diurnal cycle has better amplitude but less realistic phasing than MERRA-2 model-generated precipitation. Because correcting the precipitation within the coupled atmosphere-land modeling system allows the MERRA-2 near-surface air temperature and humidity to respond to the improved precipitation forcing, MERRA-2 provides more self-consistent surface meteorological data than were available from the earlier, offline MERRA-Land reanalysis. Overall, MERRA-2 land hydrology estimates are better than those of MERRA-Land and MERRA. A comparison against GRACE satellite observations of terrestrial water storage demonstrates clear improvements in MERRA-2 over MERRA in South America and Africa but also reflects known errors in the observations used to correct the MERRA-2 precipitation. The MERRA-2 and MERRA-Land surface and root zone soil moisture skill vs. in situ measurements is slightly higher than that of ERA-Interim Land and higher than that of MERRA (significantly for surface soil moisture). Snow amounts from MERRA-2 have lower bias and correlate better against reference data than do those of MERRA-Land and MERRA, with MERRA-2 skill roughly matching that of ERA-Interim Land. Seasonal anomaly R values against naturalized stream flow measurements in the United States are, on balance, highest for MERRA-2 and ERA-Interim Land, somewhat lower for MERRA-Land, and lower still for MERRA.

MERRA-↗

Comparing Flood-Producing Catastrophic Hurricanes in GEOS: A Study of Juan 1985 and Harvey 2017

Hurricanes Juan (1985) and Harvey (2017) bear remarkable similarities because of the exceptional amount of accumulated precipitation over coastal regions in the Gulf, consequent to the peculiar and rare situation in which a hurricane stalls for several days on a coastal area, displaying minimal progress. This is particularly dangerous because: a) terrestrial drag affect only part of the circulation so that weakening is not as fast as when there is rapid motion inland; b) a part of the circulation remains over water and as such low-level moisture flux convergence is maintained; c) precipitation lasts for a very long time over the same area; d) moisture recycling from flooded areas further contributes to maintaining flood conditions. From the forecast perspective, hurricanes with track singularities such as stalling points, loops, or multiple landfalls over the same region are exceptionally difficult to predict because they result from small variations of opposing mechanisms and are therefore intrinsically unstable. Intensity fluctuations and track singularities are strongly coupled. The representation of fluxes undergoing differential drag (ocean vs. land) inside the same circulation is critically dependent upon resolution. Recent modeling development efforts within the NASA Global Modeling and Assimilation Office (GMAO) have brought the horizontal resolution of the global cubed-sphere Goddard Earth Observing System (GEOS) atmospheric model down to 12 km. This has led to innovative findings in the field of tropical cyclone (TC) modeling within a global framework. The GEOS atmospheric model and data assimilation system (DAS), used to produce the MERRA-2 (Modern Era Retrospective analysis for Research and Applications) reanalysis, has a "Replay" capability that adds a forcing term to the model equations to constrain it to follow a specified trajectory. This capability is used to produce global, downscaled versions of the 50 km MERRA-2 reanalysis at a 12.5 km resolution. The full DAS and replay modes both use the Incremental Analysis Update or IAU for applying an analysis increment to the model's background state four times each day at 00Z, 06Z, 12Z and 18Z. In a full DAS mode, a six-hour free run of the model is made (Predictor Step), writing states that will be used as backgrounds or first-guesses for the analysis centered on 00Z. The Analysis system then produces its best estimate of the 00Z state by blending backgrounds with observations. The difference between this best estimate and the model's original first guess at 00Z is treated as an analysis increment and added uniformly as a forcing term to the model during a second, now forced, run of the 12Z to 03Z window, which produces the assimilated results. Replay mode is very similar, except only the 00Z background is needed from the predictor segment, which is thus half as long, and instead of doing an analysis, a previously analyzed state at 00Z is read in. The increment is then computed and applied to the corrector step in the same way as in the original DAS. When downscaling horizontal spatial filters can be applied to the "analysis" increment used to force the Corrector step, so as to reproduce the full meteorological state or selectively choose or discard certain forcings. In this presentation, full global hindcast simulations of Juan and operational forecasts of Harvey are compared.

Flooding↗

NASA's Modern Era Retrospective-Analysis for Research and Applications (MERRA): Early Results and Future Directions

This talk will review the status and progress of the NASA/Global Modeling and Assimilation Office (GMAO) atmospheric global reanalysis project called the Modern Era Retrospective-Analysis for Research and Applications (MERRA). An overview of NASA's emerging capabilities for assimilating a variety of other Earth Science observations of the land, ocean, and atmospheric constituents will also be presented. MERRA supports NASA Earth science by synthesizing the current suite of research satellite observations in a climate data context (covering the period 1979-present), and by providing the science and applications communities with of a broad range of weather and climate data with an emphasis on improved estimates of the hydrological cycle. MERRA is based on a major new version of the Goddard Earth Observing System Data Assimilation System (GEOS-5), that includes the Earth System Modeling Framework (ESMF)-based GEOS-5 atmospheric general circulation model and the new NOAA National Centers for Environmental Prediction (NCEP) unified grid-point statistical interpolation (GST) analysis scheme developed as a collaborative effort between NCEP and the GMAO. In addition to MERRA, the GMAO is developing new capabilities in aerosol and constituent assimilation, ocean, ocean biology, and land surface assimilation. This includes the development of an assimilation capability for tropospheric air quality monitoring and prediction, the development of a carbon-cycle modeling and assimilation system, and an ocean data assimilation system for use in coupled short-term climate forecasting.

Schubert, Siegfried↗

Evaluation of the Antarctic Boundary Layer Thermodynamic Structure In MERRA2 Using Dropsonde Observations from the Concordiasi Campaign

Recent high-resolution dropsonde observations from the 2010 Concordiasi field campaign in austral spring season show that surface-based inversions (SBIs) over Antarctica are frequently eroded, with well-mixed boundary layers occurring 33% and 18% of the time in West and East Antarctica, respectively. In this study, using the dropsonde observations, we evaluate the performance of the Modern-Era Retrospective analysis for Research and Applications, version 2 (MERRA2) in representing the Antarctic boundary layer thermodynamic structure. Results show that MERRA2 has a good overall representation of the Antarctic surface stability and correctly predicts 82% of the SBIs. However, an under prediction of less stable boundary layer occurrence, especially over the elevated East Antarctic plateau, is favored during conditions of increased lower tropospheric stability associated with model dynamics, indicating difficulty in parameterizing turbulence in very stable boundary layers. In addition, a lower tropospheric cool bias (first model level and above) is observed in the MERRA2 reanalysis, especially over West Antarctica, which amplifies in the boundary layer during mixed conditions. The near-surface cold bias is most pronounced when the model fails to predict mixed layers over West Antarctica, and is expected to negatively impact the representation of surface energy budget and melt processes.

Manisha Ganeshan↗

Report of the 4th World Climate Research Programme International Conference on Reanalyses

The 4th WCRP International Conference on Reanalyses provided an opportunity for the international community to review and discuss the observational and modelling research, as well as process studies and uncertainties associated with reanalysis of the Earth System and its components. Characterizing the uncertainty and quality of reanalyses is a task that reaches far beyond the international community of producers, and into the interdisciplinary research community, especially those using reanalysis products in their research and applications. Reanalyses have progressed greatly even in the last 5 years, and newer ideas, projects and data are coming forward. While reanalysis has typically been carried out for the individual domains of atmosphere, ocean and land, it is now moving towards coupling using Earth system models. Observations are being reprocessed and they are providing improved quality for use in reanalysis. New applications are being investigated, and the need for climate reanalyses is as strong as ever. At the heart of it all, new investigators are exploring the possibilities for reanalysis, and developing new ideas in research and applications. Given the many centres creating reanalyses products (e.g. ocean, land and cryosphere research centres as well as NWP and atmospheric centers), and the development of new ideas (e.g. families of reanalyses), the total number of reanalyses is increasing greatly, with new and innovative diagnostics and output data. The need for reanalysis data is growing steadily, and likewise, the need for open discussion and comment on the data. The 4th Conference was convened to provide a forum for constructive discussion on the objectives, strengths and weaknesses of reanalyses, indicating potential development paths for the future.

Bosilovich, Michael G.↗

Fresh Water Content Variability in the Arctic Ocean

Arctic Ocean model simulations have revealed that the Arctic Ocean has a basin wide oscillation with cyclonic and anticyclonic circulation anomalies (Arctic Ocean Oscillation; AOO) which has a prominent decadal variability. This study explores how the simulated AOO affects the Arctic Ocean stratification and its relationship to the sea ice cover variations. The simulation uses the Princeton Ocean Model coupled to sea ice. The surface forcing is based on NCEP-NCAR Reanalysis and its climatology, of which the latter is used to force the model spin-up phase. Our focus is to investigate the competition between ocean dynamics and ice formation/melt on the Arctic basin-wide fresh water balance. We find that changes in the Atlantic water inflow can explain almost all of the simulated fresh water anomalies in the main Arctic basin. The Atlantic water inflow anomalies are an essential part of AOO, which is the wind driven barotropic response to the Arctic Oscillation (AO). The baroclinic response to AO, such as Ekman pumping in the Beaufort Gyre, and ice meldfreeze anomalies in response to AO are less significant considering the whole Arctic fresh water balance.

Hakkinen, Sirpa↗

An Inversion Analysis of Recent Variability in Natural CO2 Fluxes Using GOSAT and In Situ Observations

About one-half of the global CO2 emissions from fossil fuel combustion and deforestation accumulates in the atmosphere, where it contributes to global warming. The rest is taken up by vegetation and the ocean. The precise contribution of the two sinks, and their location and year-to-year variability are, however, not well understood. We use two different approaches, batch Bayesian synthesis inversion and variational data assimilation, to deduce the global spatiotemporal distributions of CO2 fluxes during 2009-2010. One of our objectives is to assess different sources of uncertainties in inferred fluxes, including uncertainties in prior flux estimates and observations, and differences in inversion techniques. For prior constraints, we utilize fluxes and uncertainties from the CASA-GFED model of the terrestrial biosphere and biomass burning driven by satellite observations and interannually varying meteorology. We also use measurement-based ocean flux estimates and two sets of fixed fossil CO2 emissions. Here, our inversions incorporate column CO2 measurements from the GOSAT satellite (ACOS retrieval, filtered and bias-corrected) and in situ observations (individual flask and afternoon-average continuous observations) to estimate fluxes in 108 regions over 8-day intervals for the batch inversion and at 3 x 3.75 weekly for the variational system. Relationships between fluxes and atmospheric concentrations are derived consistently for the two inversion systems using the PCTM atmospheric transport model driven by meteorology from the MERRA reanalysis. We compare the posterior fluxes and uncertainties derived using different data sets and the two inversion approaches, and evaluate the posterior atmospheric concentrations against independent data including aircraft measurements. The optimized fluxes generally resemble those from other studies. For example, the results indicate that the terrestrial biosphere is a net CO2 sink, and a GOSAT-only inversion suggests a shift in the global sink from the tropics south to the north relative to the prior and to an in-situ-only inversion. We also find a smaller terrestrial sink in higher-latitude northern regions in boreal summer of 2010 relative to 2009.

Carbon dioxide↗

Global Variability in Sonic Boom Exposure due to Macroscopic Effects

Supersonic flight over land has been prohibited since 1973 due to the loudness of sonic booms. NASA is building the X-59 aircraft as part of its Quesst mission to demonstrate low-loudness shaped sonic booms, or “sonic thumps.” The Quesst mission will gather human perception data via a series of community noise surveys across the USA. The noise dose and perceptual response data will be provided to the International Civil Aviation Organization (ICAO) and the Federal Aviation Administration for use in determining potential future supersonic aircraft noise certification standards, effectively changing the prohibition from a speed limit to a noise limit. These noise regulations must be globally effective, as long travel distances see the largest benefit to supersonic flight. The state of the atmosphere through which a sonic boom travels affects the size of the region exposed to sound, the “carpet width” (CW), as well as the loudness. The focus of this dissertation is to understand and quantify the expected loudness and CW of sonic booms due to the macroscopic atmospheric effects around the world. A pair of large-scale propagation simulation studies were conducted using the NASA PCBoom code to compare predicted sonic boom loudness and CW statistics first across the USA and then across the world. For the USA study, near-field data of the X-59 in steady cruise was propagated at 4 cardinal headings at 138 locations through 5 years of Climate Forecast System Version 2 (CFSv2) atmospheric profiles. Results of a bootstrap forest predictor screening model indicated the importance of climate zone, latitude, ground elevation, season, and heading. It also noted the unimportance of time of day for predicting loudness and CW. The data is visualized in aggregate, and then broken out geographically, by season and heading, and by climate zone. Multiple linear regression models were fit to the data from the 138 locations so that estimates of the loudness and CW can be produced anywhere in the US. The results can aid in planning when and where to fly the X-59. For the global study, near-field data from three aircraft, the X-59 in a quiet and loud configuration, B-58, and Concorde, were propagated at four cardinal headings through data from three atmospheric models, the CFSv2, the Global Forecast System (GFS), and the ECMWF Reanalysis Version 5 (ERA5), at 100 global locations over 1 year. Results of a bootstrap forest predictor screening model indicated the importance of climate zone, ground elevation, season, and heading. Similar to the US study, the model indicated time of day was not an important predictor. The model also indicated that choice of weather model was not important, so the atmospheric model data are effectively interchangeable. The ERA5 model was chosen for use in an extension of the study to include 18 additional locations to ensure sampling of every climate zone. Loudness and CW results are shown in aggregate, and split geographically and by heading, season, and climate. Multiple linear regression models were fit to the data from the 118 locations so that estimates of loudness and CW can be produced around the world. N-waves and shaped booms did not have the same global variability. Koppen-Geiger climate zones were used as the climate zone definition for the global study. These are available as present-day and future climate projections. Making use of the multiple linear regression models, the future climate zones were input to estimate the effect of the changing climate on sonic boom loudness and CW. Results indicate that a changing climate would have little impact on the effectiveness of noise regulations.

X-59↗

An Overview of the GEOS-5 Aerosol Reanalysis

GEOS-5 is the latest version of the NASA Global Modeling and Assimilation Office (GMAO) earth system model. GEOS-5 contains components for atmospheric circulation and composition (including data assimilation), ocean circulation and biogeochemistry, and land surface processes. In addition to traditional meteorological parameters, GEOS-5 includes modules representing the atmospheric composition, most notably aerosols and tropospheric/stratospheric chemical constituents, taking explicit account of the impact of these constituents on the radiative processes of the atmosphere. MERRA is a NASA meteorological reanalysis for the satellite era (1979-present) using GEOS-5. This project focuses on historical analyses of the hydrological cycle on a broad range of weather and climate time scales. As a first step towards an integrated Earth System Analysis (IESA), the GMAO is extending MERRA with reanalyses for other components of the earth system: land, ocean, bio-geochemistry and atmospheric constituents. In this talk we will present results from the MERRA-driven aerosol reanalysis covering the Aqua period (2003-present). The assimilation of Aerosol Optical Depth (AOD) in GEOS-5 involves very careful cloud screening and homogenization of the observing system by means of a Neural Net scheme that translates MODIS radiances into AERONET calibrated AOD. These measurements are further quality controlled using an adaptive buddy check scheme, and assimilated using the Local Displacement Ensemble (LDE) methodology. For this reanalysis, GEOS-5 runs at a nominal 50km horizontal resolution with 72 vertical layers (top at approx. 8Skm). GEOS-5 is driven by daily biomass burning emissions derived from MODIS fire radiative power retrievals. We will present a summary of our efforts to validate such dataset. The GEOS-5 assimilated aerosol fields are first validated by comparison to independent in-situ measurements (AERONET and PM2.5 surface concentrations). In order to asses aerosol absorption on a global scale, we perform a detailed radiative transfer calculation to simulate the UV aerosol index, comparing our results to OMI measurements. By simulating aerosol attenuated backscatter, we use CALIPSO measurements to evaluate the vertical structure of our aerosol estimates, in particular in regions where we have larger discrepancies with OMI. Finally, the consistency of our AOD estimates with estimates from MISR, MODIS/Deep Blue, OMI and PARASOL will be briefly discussed.

da Silva, Arlindo↗

Approximations to eigenvalues of modified general matrices

The reanalysis of non-self-adjoint dynamic models is computationally very expensive in design optimization applications. This paper describes several approximations that can be applied to eigenvalues of non-hermitian matrices to reduce that computational cost. Approximations based on eigenvalue derivatives, generalized Rayleigh quotient and the trace theorem are presented and their accuracy and computational cost are estimated. The accuracy and cost estimates are verified by applying the approximations to random matrices and matrices arising in flutter analysis of compressor blades. Recommendations are made for selection of the best approximation when the derivatives are available and when they are not. In particular, it is concluded that the quadratic approximation for eigenvalues should never be used as higher order approximations are always more accurate as well as more efficient.

Murthy, Durbha V.↗

Diagnosing Warm Season Precipitation Over the GCIP Region from a GCM and Reanalysis

A 45 year simulation using a global general circulation model (GCM), the National Center for Atmospheric Research (NCAR) Community Climate Model v.3 (CCM3), forced with observed sea surface temperatures (SST), and 39 years of global National Centers for Environmental Prediction (NCEP) reanalyses were analyzed to determine Mississippi River basin warm season (May, June, July or MJJ) wet and dry year composites in the water and energy budgets. Years that have increased MJJ soil moisture over the GEWEX (Global Water and Energy Experiment) Continental Interior Project (GCIP) region also have high precipitation, lower surface temperature, decreased Bowen ratio, and reduced 500 hPa geopotential height (essentially reduced MJJ ridging). The reverse is true for years that have reduced MJJ soil moisture. Wet years are also accompanied by a general increase in moisture transport from the Gulf of Mexico into the central U.S. There is some indication (though weaker) that soil moisture may then affect precipitation and other quantities and be affected in turn by 500 hPa geopotential heights. The correlations are somewhat low, however, demonstrating the difficulty in providing definitive physical links between the remote and local effects. Analysis of two individual years with an extreme wet event (1993) and an extreme dry event (1988) yields the same general relationships as with the wet and dry composites. The composites from this study are currently serving as the basis for a series of experiments aimed at determining the predictability of the land surface and remote SST on the Mississippi River basin and other large-scale river basins.

Oglesby, Robert↗

The Use of MERRA-2 Near Surface Meteorology to Understand the Behavior of Planetary Boundary Layer Heights Derived from Wind Profiler Data over the US Great Plains

The atmospheric general circulation model (GCM) that underlies the MERRA-2 reanalysis includes a suite of physical parameterizations that describe the processes that occur in the planetary boundary layer (PBL). The data assimilation system assures that the atmospheric state variables used as input to these parameterizations are constrained to the best fit to all of the available observations. Many studies, however, have shown that the GCM-based estimates of MERRA-2 PBL heights are biased high, and so are not reliable for application related to constituent transport or the carbon cycle. A new 20-year record of PBL heights was derived from Wind Profiler (WP) backscatter data measured at a wide network of stations throughout the US Great Plains and has been validated against independent estimates. The behavior of these PBL heights shows geographical and temporal variations that are difficult to attribute to particular physical processes without additional information that are not part of the observational record. In the present study, we use information on physical processes from MERRA-2 to understand the behavior of the WP derived PBL heights. The behavior of the annual cycle of both MERRA-2 and WP PBL heights shows three classes of behavior: (i) canonical, where the annual cycle follows the annual cycle of the sun, (ii) delayed, where the PBL height reaches its annual maximum after the annual maximum of the solar insolation, and (iii) double maxima, where the PBL height begins to rise with the solar insolation but falls sometimes during the summer and then rises again. Although the magnitude of these types of variations is described by the WP PBL record, the explanation for these behaviors and the relationship to local precipitation, temperature, hydrology and sensible and latent heat fluxes is articulated using information from MERRA-2.

MERRA↗