Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “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 523 records · Page 29

Linking Excessive Heat with Daily Heat-Related Mortality over the Coterminous United States

In the United States, extreme heat is the most deadly weather-related hazard. In the face of a warming climate and urbanization, which contributes to local-scale urban heat islands, it is very likely that extreme heat events (EHEs) will become more common and more severe in the U.S. This research seeks to provide historical and future measures of climate-driven extreme heat events to enable assessments of the impacts of heat on public health over the coterminous U.S. We use atmospheric temperature and humidity information from meteorological reanalysis and from Global Climate Models (GCMs) to provide data on past and future heat events. The focus of research is on providing assessments of the magnitude, frequency and geographic distribution of extreme heat in the U.S. to facilitate public health studies. In our approach, long-term climate change is captured with GCM outputs, and the temporal and spatial characteristics of short-term extremes are represented by the reanalysis data. Two future time horizons for 2040 and 2090 are compared to the recent past period of 1981- 2000. We characterize regional-scale temperature and humidity conditions using GCM outputs for two climate change scenarios (A2 and A1B) defined in the Special Report on Emissions Scenarios (SRES). For each future period, 20 years of multi-model GCM outputs are analyzed to develop a 'heat stress climatology' based on statistics of extreme heat indicators. Differences between the two future and the past period are used to define temperature and humidity changes on a monthly time scale and regional spatial scale. These changes are combined with the historical meteorological data, which is hourly and at a spatial scale (12 km) much finer than that of GCMs, to create future climate realizations. From these realizations, we compute the daily heat stress measures and related spatially-specific climatological fields, such as the mean annual number of days above certain thresholds of maximum and minimum air temperatures, heat indices, and a new heat stress variable developed as part of this research that gives an integrated measure of heat stress (and relief) over the course of a day. Comparisons are made between projected (2040 and 2090) and past (1990) heat stress statistics. Outputs are aggregated to the county level, which is a popular scale of analysis for public health interests. County-level statistics are made available to public health researchers by the Centers for Disease Control and Prevention (CDC) via the Wide-ranging Online Data for Epidemiologic Research (WONDER) system. This addition of heat stress measures to CDC WONDER allows decision and policy makers to assess the impact of alternative approaches to optimize the public health response to EHEs. Through CDC WONDER, users are able to spatially and temporally query public health and heat-related data sets and create county-level maps and statistical charts of such data across the coterminous U.S.

Quattrochi, Dale A.↗

ldentifying Episodes of Earth Science Phenomena Using a Big-Data Technology

A significant portion of Earth Science investigations is phenomenon- (or event-) based, such as the studies of Rossby waves, volcano eruptions, tsunamis, mesoscale convective systems, and tropical cyclones. However, except for a few high-impact phenomena, e.g. tropical cyclones, comprehensive records are absent for the occurrences or events of these phenomena. Phenomenon-based studies therefore often focus on a few prominent cases while the lesser ones are overlooked. Without an automated means to gather the events, comprehensive investigation of a phenomenon is at least time-consuming if not impossible. We have constructed a prototype Automated Event Service (AES) system that is used to methodically mine custom-defined events in the reanalysis data sets of atmospheric general circulation models. Our AES will enable researchers to specify their custom, numeric event criteria using a user-friendly web interface to search the reanalysis data sets. Moreover, we have included a social component to enable dynamic formation of collaboration groups for researchers to cooperate on event definitions of common interest and for the analysis of these events. An Earth Science event (ES event) is defined here as an episode of an Earth Science phenomenon (ES phenomenon). A cumulus cloud, a thunderstorm shower, a rogue wave, a tornado, an earthquake, a tsunami, a hurricane, or an El Nino, is each an episode of a named ES phenomenon, and, from the small and insignificant to the large and potent, all are examples of ES events. An ES event has a duration (often finite) and an associated geo-location as a function of time; it's therefore an entity embedded in four-dimensional (4D) spatiotemporal space. Earth Science phenomena with the potential to cause massive economic disruption or loss of life often rivet the attention of researchers. But, broader scientific curiosity also drives the study of phenomena that pose no immediate danger, such as land/sea breezes. Due to Earth System's intricate dynamics, we are continuously discovering novel ES phenomena. We generally gain understanding of a given phenomenon by observing and studying individual events. This process usually begins by identifying the occurrences of these events. Once representative events are identified or found, we must locate associated observed or simulated data prior to commencing analysis and concerted studies of the phenomenon. Knowledge concerning the phenomenon can accumulate only after analysis has started. However, as mentioned previously, comprehensive records only exist for a very limited set of high-impact phenomena; aside from these, finding events and locating associated data currently may take a prohibitive amount of time and effort on the part of an individual investigator. The reason for the lack of comprehensive records for most of the ES phenomena is mainly due to the perception that they do not pose immediate and/or severe threat to life and property. Thus they are not consistently tracked, monitored, and catalogued. Many phenomena even lack precise and/or commonly accepted criteria for definitions. Moreover, various Earth Science observations and data have accumulated to a previously unfathomable volume; NASA Earth Observing System Data Information System (EOSDIS) alone archives several petabytes (PB) of satellite remote sensing data, which are steadily increasing. All of these factors contribute to the difficulty of methodically identifying events corresponding to a given phenomenon and significantly impede systematic investigations. We have not only envisioned AES as an environment for identifying customdefined events but also aspired for it to be an interactive environment with quick turnaround time for revisions of query criteria and results, as well as a collaborative environment where geographically distributed experts may work together on the same phenomena. A Big Data technology is thus required for the realization of such a system. In the following, we first introduce the technology selected for AES in the next section. We then demonstrate the utility of AES using a use case, Blizzard, before we conclude.

Kuo, Kwo-Sen↗

The Impact of the Atlantic Cold Tongue on West African Monsoon Onset in Regional Model Simulations for 1998-2002

The Atlantic cold tongue (ACT) develops during spring and early summer near the Equator in the Eastern Atlantic Ocean and Gulf of Guinea. The hypothesis that the ACT accelerates the timing of West African monsoon (WAM) onset is tested by comparing two regional climate model (RM3) simulation ensembles. Observed sea surface temperatures (SST) that include the ACT are used to force a control ensemble. An idealized, warm SST perturbation is designed to represent lower boundary forcing without the ACT for the experiment ensemble. Summer simulations forced by observed SST and reanalysis boundary conditions for each of five consecutive years are compared to five parallel runs forced by SST with the warm perturbation. The article summarizes the sequence of events leading to the onset of the WAM in the Sahel region. The representation of WAM onset in RM3 simulations is examined and compared to Tropical Rainfall Measuring Mission (TRMM), Global Precipitation Climatology Project (GPCP) and reanalysis data. The study evaluates the sensitivity of WAM onset indicators to the presence of the ACT by analysing the differences between the two simulation ensembles. Results show that the timing of major rainfall events and therefore theWAM onset in the Sahel are not sensitive to the presence of the ACT. However, the warm SST perturbation does increase downstream rainfall rates over West Africa as a consequence of enhanced specific humidity and enhanced northward moisture flux in the lower troposphere.

Regional model simulation↗

Estimating Planetary Boundary Layer Heights from NOAA Profiler Network Wind Profiler Data

An algorithm was developed to estimate planetary boundary layer (PBL) heights from hourly archived wind profiler data from the NOAA Profiler Network (NPN) sites located throughout the central United States. Unlike previous studies, the present algorithm has been applied to a long record of publicly available wind profiler signal backscatter data. Under clear conditions, summertime averaged hourly time series of PBL heights compare well with Richardson-number based estimates at the few NPN stations with hourly temperature measurements. Comparisons with clear sky reanalysis based estimates show that the wind profiler PBL heights are lower by approximately 250-500 m. The geographical distribution of daily maximum PBL heights corresponds well with the expected distribution based on patterns of surface temperature and soil moisture. Wind profiler PBL heights were also estimated under mostly cloudy conditions, and are generally higher than both the Richardson number based and reanalysis PBL heights, resulting in a smaller clear-cloudy condition difference. The algorithm presented here was shown to provide a reliable summertime climatology of daytime hourly PBL heights throughout the central United States.

Wind Profiler↗

Development of the GEOS-5 Atmospheric General Circulation Model: Evolution from MERRA to MERRA2.

The Modern-Era Retrospective Analysis for Research and Applications-2 (MERRA2) version of the GEOS-5 (Goddard Earth Observing System Model - 5) Atmospheric General Circulation Model (AGCM) is currently in use in the NASA Global Modeling and Assimilation Office (GMAO) at a wide range of resolutions for a variety of applications. Details of the changes in parameterizations subsequent to the version in the original MERRA reanalysis are presented here. Results of a series of atmosphere-only sensitivity studies are shown to demonstrate changes in simulated climate associated with specific changes in physical parameterizations, and the impact of the newly implemented resolution-aware behavior on simulations at different resolutions is demonstrated. The GEOS-5 AGCM presented here is the model used as part of the GMAO's MERRA2 reanalysis, the global mesoscale "nature run", the real-time numerical weather prediction system, and for atmosphere-only, coupled ocean-atmosphere and coupled atmosphere-chemistry simulations. The seasonal mean climate of the MERRA2 version of the GEOS-5 AGCM represents a substantial improvement over the simulated climate of the MERRA version at all resolutions and for all applications. Fundamental improvements in simulated climate are associated with the increased re-evaporation of frozen precipitation and cloud condensate, resulting in a wetter atmosphere. Improvements in simulated climate are also shown to be attributable to changes in the background gravity wave drag, and to upgrades in the relationship between the ocean surface stress and the ocean roughness. The series of "resolution aware" parameters related to the moist physics were shown to result in improvements at higher resolutions, and result in AGCM simulations that exhibit seamless behavior across different resolutions and applications.

AGCM↗

Decadal Prediction Skill in the GEOS-5 Forecast System

A suite of decadal predictions has been conducted with the NASA Global Modeling and Assimilation Office's (GMAO's) GEOS-5 Atmosphere-Ocean general circulation model. The hind casts are initialized every December 1st from 1959 to 2010, following the CMIP5 experimental protocol for decadal predictions. The initial conditions are from a multivariate ensemble optimal interpolation ocean and sea-ice reanalysis, and from GMAO's atmospheric reanalysis, the modern-era retrospective analysis for research and applications. The mean forecast skill of a three-member-ensemble is compared to that of an experiment without initialization but also forced with observed greenhouse gases. The results show that initialization increases the forecast skill of North Atlantic sea surface temperature compared to the uninitialized runs, with the increase in skill maintained for almost a decade over the subtropical and mid-latitude Atlantic. On the other hand, the initialization reduces the skill in predicting the warming trend over some regions outside the Atlantic. The annual-mean Atlantic meridional overturning circulation index, which is defined here as the maximum of the zonally-integrated overturning stream function at mid-latitude, is predictable up to a 4-year lead time, consistent with the predictable signal in upper ocean heat content over the North Atlantic. While the 6- to 9-year forecast skill measured by mean squared skill score shows 50 percent improvement in the upper ocean heat content over the subtropical and mid-latitude Atlantic, prediction skill is relatively low in the sub-polar gyre. This low skill is due in part to features in the spatial pattern of the dominant simulated decadal mode in upper ocean heat content over this region that differ from observations. An analysis of the large-scale temperature budget shows that this is the result of a model bias, implying that realistic simulation of the climatological fields is crucial for skillful decadal forecasts.

Decadal Prediction↗

Air Pollution Over the Ganges Basin and Northwest Bay of Bengal in the Early Postmonsoon Season Based on NASA MERRAero Data

The MERRA Aerosol Reanalysis (MERRAero) has been recently developed at NASA's Global Modeling Assimilation Office. This reanalysis is based on a version of the Goddard Earth Observing System-5 (GEOS-5) model radiatively coupled with Goddard Chemistry, Aerosol, Radiation, and Transport aerosols, and it includes assimilation of bias-corrected aerosol optical thickness (AOT) from the Moderate Resolution Imaging Spectroradiometer (MODIS) sensor on both Terra and Aqua satellites. In October over the period 2002-2009, MERRAero showed that AOT was lower over the east of the Ganges basin than over the northwest of the Ganges basin: this was despite the fact that the east of the Ganges basin should have produced higher anthropogenic aerosol emissions because of higher population density, increased industrial output, and transportation. This is evidence that higher aerosol emissions do not always correspond to higher AOT over the areas where the effects of meteorological factors on AOT dominate those of aerosol emissions. MODIS AOT assimilation was essential for correcting modeled AOT mainly over the northwest of the Ganges basin, where AOT increments were maximal. Over the east of the Ganges basin and northwest Bay of Bengal (BoB), AOT increments were low and MODIS AOT assimilation did not contribute significantly to modeled AOT. Our analysis showed that increasing AOT trends over northwest BoB (exceeding those over the east of the Ganges basin) were reproduced by GEOS-5, not because of MODIS AOT assimilation butmainly because of the model capability of reproducing meteorological factors contributing to AOT trends. Moreover, vertically integrated aerosol mass flux was sensitive to wind convergence causing aerosol accumulation over northwest BoB.

AOT↗

MERRA-2: File Specification

The second Modern-Era Retrospective analysis for Research and Applications (MERRA-2) is a NASA atmospheric reanalysis that begins in 1980. It replaces the original MERRA reanalysis (Rienecker et al., 2011) using an upgraded version of the Goddard Earth Observing System Model, Version 5 (GEOS-5) data assimilation system. The file collections for MERRA-2 are described in detail in this document, including some important changes from those of the MERRA dataset (Lucchesi, 2012).

MERRA-2↗

Solar Signals in CMIP-5 Simulations: The Stratospheric Pathway

The 11 year solar-cycle component of climate variability is assessed in historical simulations of models taken from the Coupled Model Intercomparison Project, phase 5 (CMIP-5). Multiple linear regression is applied to estimate the zonal temperature, wind and annular mode responses to a typical solar cycle, with a focus on both the stratosphere and the stratospheric influence on the surface over the period approximately 1850-2005. The analysis is performed on all CMIP-5 models but focuses on the 13 CMIP-5 models that resolve the stratosphere (high-top models) and compares the simulated solar cycle signature with reanalysis data. The 11 year solar cycle component of climate variability is found to be weaker in terms of magnitude and latitudinal gradient around the stratopause in the models than in the reanalysis. The peak in temperature in the lower equatorial stratosphere (approximately 70 hPa) reported in some studies is found in the models to depend on the length of the analysis period, with the last 30 years yielding the strongest response. A modification of the Polar Jet Oscillation (PJO) in response to the 11 year solar cycle is not robust across all models, but is more apparent in models with high spectral resolution in the short-wave region. The PJO evolution is slower in these models, leading to a stronger response during February, whereas observations indicate it to be weaker. In early winter, the magnitude of the modeled response is more consistent with observations when only data from 1979-2005 are considered. The observed North Pacific high-pressure surface response during the solar maximum is only simulated in some models, for which there are no distinguishing model characteristics. The lagged North Atlantic surface response is reproduced in both high- and low-top models, but is more prevalent in the former. In both cases, the magnitude of the response is generally lower than in observations.

stratosphere↗

Assessing the Ability of Instantaneous Aircraft and Sonde Measurements to Characterize Climatological Means and Long-Term Trends in Tropospheric Composition

Over four decades of measurements exist that sample the 3-D composition of reactive trace gases in the troposphere from approximately weekly ozone sondes, instrumentation on civil aircraft, and individual comprehensive aircraft field campaigns. An obstacle to using these data to evaluate coupled chemistry-climate models (CCMs)the models used to project future changes in atmospheric composition and climateis that exact space-time matching between model fields and observations cannot be done, as CCMs generate their own meteorology. Evaluation typically involves averaging over large spatiotemporal regions, which may not reflect a true average due to limited or biased sampling. This averaging approach generally loses information regarding specific processes. Here we aim to identify where discrete sampling may be indicative of long-term mean conditions, using the GEOS-Chem global chemical-transport model (CTM) driven by the MERRA reanalysis to reflect historical meteorology from 2003 to 2012 at 2o by 2.5o resolution. The model has been sampled at the time and location of every ozone sonde profile available from the Would Ozone and Ultraviolet Radiation Data Centre (WOUDC), along the flight tracks of the IAGOSMOZAICCARABIC civil aircraft campaigns, as well as those from over 20 individual field campaigns performed by NASA, NOAA, DOE, NSF, NERC (UK), and DLR (Germany) during the simulation period. Focusing on ozone, carbon monoxide and reactive nitrogen species, we assess where aggregates of the in situ data are representative of the decadal mean vertical, spatial and temporal distributions that would be appropriate for evaluating CCMs. Next, we identically sample a series of parallel sensitivity simulations in which individual emission sources (e.g., lightning, biogenic VOCs, wildfires, US anthropogenic) have been removed one by one, to assess where and when the aggregated observations may offer constraints on these processes within CCMs. Lastly, we show results of an additional 31-year simulation from 1980-2010 of GEOS-Chem driven by the MACCity emissions inventory and MERRA reanalysis at 4o by 5o. We sample the model at every WOUDC sonde and flight track from MOZAIC and NASA field campaigns to evaluate which aggregate observations are statistically reflective of long-term trends over the period.

Atmospheric composition↗

Covariability in the Monthly Mean Convective and Radiative Diurnal Cycles in the Amazon

The diurnal cycle of convective clouds greatly influences the radiative energy balance in convectively active regions of Earth, through both direct presence, and the production of anvil and stratiform clouds. Previous studies show that the frequency and properties of convective clouds can vary on monthly timescales as a result of variability in the monthly mean atmospheric state. Furthermore, the radiative budget in convectively active regions also varies by up to 7 Wm-2 in convectively active regions. These facts suggest that convective clouds connect atmospheric state variability and radiation variability beyond clear sky effects alone. Previous research has identified monthly covariability between the diurnal cycle of CERES-observed top-of-atmosphere radiative fluxes and multiple atmospheric state variables from reanalysis over the Amazon region. ASVs that enhance (reduce) deep convection, such as CAPE (LTS), tend to shift the daily OLR and cloud albedo maxima earlier (later) in the day by 2-3 hr. We first test the analysis method using multiple reanalysis products for both the dry and wet seasons to further investigate the robustness of the preliminary results. We then use CloudSat data as an independent cloud observing system to further evaluate the relationships of cloud properties to variability in radiation and atmospheric states. While CERES can decompose OLR variability into clear sky and cloud effects, it cannot determine what variability in cloud properties lead to variability in the radiative cloud effects. Cloud frequency, cloud top height, and cloud microphysics all contribute to the cloud radiative effect, all of which are observable by CloudSat. In addition, CloudSat can also observe the presence and variability of deep convective cores responsible for the production of anvil clouds. We use these capabilities to determine the covariability of convective cloud properties and the radiative diurnal cycle.

Dodson, Jason B.↗

Investigating Ozone Sources in California Using AJAX Airborne Measurements and Models: Implications for Stratospheric Intrusion and Long Range Transport

High ozone concentrations at low altitudes near the surface were detected from airborne Alpha Jet Atmospheric eXperiment (AJAX) measurements on May 30, 2012. We investigate the causes of the elevated ozone concentrations using the airborne measurements and various models. GEOSchem and WRF-STILT model simulations show that the contribution from local sources is small. From MERRA reanalysis, it is found that high potential vorticity (PV) is observed at low altitudes. This high PV appears to be only partially coming through the stratospheric intrusions because the air inside the high PV region is moist, which shows that mixing appears to be enhanced in the low altitudes. Considering that diabatic heating can also produce high PV in the lower troposphere, high ozone is partially coming through stratospheric intrusion, but this cannot explain the whole ozone concentration in the target areas of the western U.S. A back-trajectory model is utilized to see where the air masses originated. The air masses of the target areas came from the lower stratosphere (LS), upper (UT), mid- (MT), and lower troposphere (LT). The relative number of trajectories coming from LS and UT is low (7.7% and 7.6%, respectively) compared to that from LT (64.1%), but the relative ozone concentration coming from LS and UT is high (38.4% and 20.95%, respectively) compared to that from LT (17.7%). The air mass coming from LT appears to be mostly coming from Asia. Q diagnostics show that there is sufficient mixing along the trajectory to indicate that ozone from the different origins is mixed and transported to the western U.S. This study shows that high ozone concentrations can be detected by airborne measurements, which can be analyzed by integrated platforms such as models, reanalysis, and satellite data.

gas composition↗

Investigation of Ozone Sources in California Using AJAX Airborne Measurements and Models: Implications for Stratospheric Intrusion and Long Range Transport

High ozone concentrations at low altitudes near the surface were detected from airborne Alpha Jet Atmospheric eXperiment (AJAX) measurements on May 30, 2012. We investigate the causes of the elevated ozone concentrations using the airborne measurements and various models. GEOS-chem and WRF-STILT model simulations show that the contribution from local sources is small. From MERRA reanalysis, it is found that high potential vorticity (PV) is observed at low altitudes. This high PV appears to be only partially coming through the stratospheric intrusions because the air inside the high PV region is moist, which shows that mixing appears to be enhanced in the low altitudes. Considering that diabatic heating can also produce high PV in the lower troposphere, high ozone is partially coming through stratospheric intrusion, but this cannot explain the whole ozone concentration in the target areas of the western U.S. A back-trajectory model is utilized to see where the air masses originated. The air masses of the target areas came from the lower stratosphere (LS), upper (UT), mid- (MT), and lower troposphere (LT). The relative number of trajectories coming from LS and UT is low (7.7 and 7.6, respectively) compared to that from LT (64.1), but the relative ozone concentration coming from LS and UT is high (38.4 and 20.95, respectively) compared to that from LT (17.7). The air mass coming from LT appears to be mostly coming from Asia. Q diagnostics show that there is sufficient mixing along the trajectory to indicate that ozone from the different origins is mixed and transported to the western U.S. This study shows that high ozone concentrations can be detected by airborne measurements, which can be analyzed by integrated platforms such as models, reanalysis, and satellite data.

low altitude↗

Enabling Research Tools for Sustained Climate Assessment

The U.S. Global Change Research Program Sustained Assessment process benefits from long-term investments in Earth science research that enable the scientific community to conduct assessment-relevant science. To this end, NASA initiated several research programs over the past five years to support the Earth observation community in developing indicators, datasets, research products, and tools to support ongoing and future National Climate Assessments. These activities complement NASA's ongoing Earth science research programs. One aspect of the assessment portfolio funds four "enabling tools" projects at NASA research centers. Each tool leverages existing capacity within the center, but has developed tailored applications and products for National Climate Assessments. The four projects build on the capabilities of a global atmospheric reanalysis (MERRA-2), a continental U.S. land surface reanalysis (NCA-LDAS), the NASA Earth Exchange (NEX), and a Regional Climate Model Evaluation System (RCMES). Here, we provide a brief overview of each enabling tool, highlighting the ways in which it has advanced assessment science to date. We also discuss how the assessment community can access and utilize these tools for National Climate Assessments and other sustained assessment activities.

Research↗

GEOS Atmospheric Model: Challenges at Exascale

The Goddard Earth Observing System (GEOS) model at NASA's Global Modeling and Assimilation Office (GMAO) is used to simulate the multi-scale variability of the Earth's weather and climate, and is used primarily to assimilate conventional and satellite-based observations for weather forecasting and reanalysis. In addition, assimilations coupled to an ocean model are used for longer-term forecasting (e.g., El Nino) on seasonal to interannual times-scales. The GMAO's research activities, including system development, focus on numerous time and space scales, as detailed on the GMAO website, where they are tabbed under five major themes: Weather Analysis and Prediction; Seasonal-Decadal Analysis and Prediction; Reanalysis; Global Mesoscale Modeling, and Observing System Science. A brief description of the GEOS systems can also be found at the GMAO website. GEOS executes as a collection of earth system components connected through the Earth System Modeling Framework (ESMF). The ESMF layer is supplemented with the MAPL (Modeling, Analysis, and Prediction Layer) software toolkit developed at the GMAO, which facilitates the organization of the computational components into a hierarchical architecture. GEOS systems run in parallel using a horizontal decomposition of the Earth's sphere into processing elements (PEs). Communication between PEs is primarily through a message passing framework, using the message passing interface (MPI), and through explicit use of node-level shared memory access via the SHMEM (Symmetric Hierarchical Memory access) protocol. Production GEOS weather prediction systems currently run at 12.5-kilometer horizontal resolution with 72 vertical levels decomposed into PEs associated with 5,400 MPI processes. Research GEOS systems run at resolutions as fine as 1.5 kilometers globally using as many as 30,000 MPI processes. Looking forward, these systems can be expected to see a 2 times increase in horizontal resolution every two to three years, as well as less frequent increases in vertical resolution. Coupling these resolution changes with increases in complexity, the computational demands on the GEOS production and research systems should easily increase 100-fold over the next five years. Currently, our 12.5 kilometer weather prediction system narrowly meets the time-to-solution demands of a near-real-time production system. Work is now in progress to take advantage of a hybrid MPI-OpenMP parallelism strategy, in an attempt to achieve a modest two-fold speed-up to accommodate an immediate demand due to increased scientific complexity and an increase in vertical resolution. Pursuing demands that require a 10- to 100-fold increases or more, however, would require a detailed exploration of the computational profile of GEOS, as well as targeted solutions using more advanced high-performance computing technologies. Increased computing demands of 100-fold will be required within five years based on anticipated changes in the GEOS production systems, increases of 1000-fold can be anticipated over the next ten years.

ESMF↗

Martian Polar Vortices: Comparison of Reanalyses

The structure and evolution of the Martian polar vortices is examined using two recently available reanalysis systems: version 1.0 of the Mars Analysis Correction Data Assimilation (MACDA) and a preliminary version of the Ensemble Mars Atmosphere Reanalysis System (EMARS). There is quantitative agreement between the reanalyses in the lower atmosphere, where Mars Global Surveyor (MGS) Thermal Emission Spectrometer (TES) data are assimilated, but there are differences at higher altitudes reflecting differences in the free-running general circulation model simulations used in the two reanalyses. The reanalyses show similar potential vorticity (PV) structure of the vortices: There is near-uniform small PV equatorward of the core of the westerly jet, steep meridional PV gradients on the polar side of the jet core, and a maximum of PV located off of the pole. In maps of 30 sol mean PV, there is a near-continuous elliptical ring of high PV with roughly constant shape and longitudinal orientation from fall to spring. However, the shape and orientation of the vortex varies on daily time scales, and there is not a continuous ring of PV but rather a series of smaller scale coherent regions of high PV. The PV structure of the Martian polar vortices is, as has been reported before, very different from that of Earth's stratospheric polar vortices, but there are similarities with Earth's tropospheric vortices which also occur at the edge of the Hadley Cell, and have near-uniform small PV equatorward of the jet, and a large increase of PV poleward of the jet due to increased stratification.

Waugh, D. W.↗

Evaluation of the Sensitivity of the Amazonian Diurnal Cycle to Convective Intensity in Reanalyses

Model parameterizations of tropical deep convection are unable to reproduce the observed diurnal and spatial variability of convection in the Amazon, which contributes to climatological biases in the water cycle and energy budget. Convective intensity regimes are defined using percentiles of daily minimum 3-hourly averaged outgoing longwave radiation (OLR) from Clouds and the Earth's Radiant Energy System (CERES). This study compares the observed spatial variability of convective diurnal cycle statistics for each regime to MERRA-2 and ERA-Interim (ERA) reanalysis data sets. Composite diurnal cycle statistics are computed for daytime hours (06:00-21:00 local time) in the wet season (December-January-February). MERRA-2 matches observations more closely than ERA for domain averaged composite diurnal statistics-specifically precipitation. However, ERA reproduces mesoscale features of OLR and precipitation phase associated with topography and the propagation of the coastal squall line. Both reanalysis models are shown to underestimate extreme convection.

Itterly, Kyle F.↗

Saharan Dust as a Causal Factor of Significant Cloud Cover Along the Saharan Air Layer in the Atlantic Ocean

The tropical Atlantic is frequently affected by Saharan dust intrusions. Based on MODIS cloud fraction (CF) data during the ten-year study period, we found that these dust intrusions contribute to significant cloud cover along the Saharan Air Layer (SAL). Below the temperature inversion at the SAL's base, the presence of large amounts of settling dust particles, together with marine aerosols, produces meteorological conditions suitable for the formation of shallow stratocumulus clouds. The significant cloud fraction along the SAL together with clouds over the Atlantic Inter-tropical Convergence Zone contributes to the 20% hemispheric CF asymmetry between the tropical North and South Atlantic. This leads to the imbalance in strong solar radiation, which reaches the sea surface between the tropical North and South Atlantic, and, consequently, affects climate formation in the tropical Atlantic. Therefore, despite the fact that, over the global ocean, there is no noticeable hemispheric asymmetry in cloud fraction, over the significant area such as the tropical Atlantic the hemispheric asymmetry in CF takes place. Saharan dust is also the major contributor to hemispheric aerosol asymmetry over the tropical Atlantic. The NASA GEOS-5 model with aerosol data assimilation was used to extend the MERRA reanalysis with five atmospheric aerosol species (desert dust, sulfates, organic carbon, black carbon, and sea-salt). The obtained ten-year (2002 - 2012) MERRA-driven aerosol reanalysis dataset (aka MERRAero) showed that, over the tropical Atlantic, dust and carbonaceous aerosols were distributed asymmetrically relative to the equator, while other aerosol species were distributed more symmetrically.

CF↗