Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Statistical forecasting”

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 541 records · Page 30
All sourcesNASA NTRSDOE OSTI
NASA NTRS2015

Evaluation of the 7-km GEOS-5 Nature Run

This report documents an evaluation by the Global Modeling and Assimilation Office (GMAO) of a two-year 7-km-resolution non-hydrostatic global mesoscale simulation produced with the Goddard Earth Observing System (GEOS-5) atmospheric general circulation model. The simulation was produced as a Nature Run for conducting observing system simulation experiments (OSSEs). Generation of the GEOS-5 Nature Run (G5NR) was motivated in part by the desire of the OSSE community for an improved high-resolution sequel to an existing Nature Run produced by the European Centre for Medium-Range Weather Forecasts (ECMWF), which has served the community for several years. The intended use of the G5NR in this context is for generating simulated observations to test proposed observing system designs regarding new instruments and their deployments. Because NASA's interest in OSSEs extends beyond traditional weather forecasting applications, the G5NR includes, in addition to standard meteorological components, a suite of aerosol types and several trace gas concentrations, with emissions downscaled to 10 km using ancillary information such as power plant location, population density and night-light information. The evaluation exercise described here involved more than twenty-five GMAO scientists investigating various aspects of the G5NR performance, including time mean temperature and wind fields, energy spectra, precipitation and the hydrological cycle, the representation of waves, tropical cyclones and midlatitude storms, land and ocean surface characteristics, the representation and forcing effects of clouds and radiation, dynamics of the stratosphere and mesosphere, and the representation of aerosols and trace gases. Comparisons are made with observational data sets when possible, as well as with reanalyses and other long model simulations. The evaluation is broad in scope, as it is meant to assess the overall realism of basic aspects of the G5NR deemed relevant to the conduct of OSSEs. However, because of the relatively short record and other practical considerations, these comparisons cannot provide a definitive, statistically sound assessment of all model deficiencies, or guarantee the G5NR's suitability for all OSSE applications. Differences between the observed and simulated behavior also must be judged in the context of basic internal atmospheric variability which can introduce variations that are not necessarily controlled by the prescribed sea surface temperatures used in generating the G5NR. The results show that the G5NR performs well as measured by the majority of metrics applied in this evaluation. Particular benefits derived from the 7-km resolution of G5NR include realistic representations of extreme weather events in both the tropics and extratropics including tropical cyclones, Nor'easters and mesoscale convective complexes; improved representation of the diurnal cycle of precipitation over land; well-resolved surface-atmosphere interactions such as katabatic wind flows over Antarctica and Greenland; and resolution of orographically generated gravity waves that propagate into the upper atmosphere and influence the large scale circulation. Obvious deficiencies in the G5NR include a "splitting" of the inter-tropical convergence zone, which leads to a weaker-than-observed Hadley circulation and related deficiencies in the depiction of stationary wave patterns. Also, while the G5NR captures global cloud features and radiative effects well in general, close comparison with observations reveals higher-than-observed cloud brightness, likely due to an overabundance of cloud condensate; less distinct cloud minima in subtropical subsidence zones, consistent with a weak Hadley circualtion; and too few near-coastal marine stratocumulus clouds.

GEOS-5↗
NASA NTRS2007

Electrification in Hurricanes over the Tropical Americas: Implication for Stratospheric Water Vapor

This study explores the relation between lightning activity and water vapor in the Tropical Tropopause Layer (TTL) over hurricane systems in the Tropical Americas. The hypothesis herein is that hurricanes that exhibit enhanced lightning activity are associated with stronger updrafts that can transport more moisture directly into the TTL (and subsequently into the tropical stratosphere) or even directly into the tropical stratosphere over this region. The TTL over the Tropical Americas, which includes the Caribbean and Gulf of Mexico, is of particular interest, because summertime cold point tropopause is the lowest in height and thus the warmest in temperature over the tropics. The latter condition implies higher saturation values and thus potential for more water vapor to enter the stratosphere. Climate forecast is very sensitive to stratospheric water vapor abundance, because of the key role that water vapor plays in regulating the chemical and radiative properties of the stratosphere. Given the potential for increases in hurricane intensity and frequency under predicted warmer conditions, it becomes essential to understand the effect of hurricanes on stratospheric water vapor. In this study, we use a combination of ground and space-borne observations as well as trajectory calculations. The observations include: cloud-to-ground (CG) lightning data from the U.S. National Lightning Detection Network (NLDN), geostationary infrared observations from the National Climatic Data Center Hurricane Satellite (HURSAT) data set, cloud properties from Aqua-MODIS, and water vapor from Aura-MLS. We analyze hurricanes from the 2005 season when Aura-MLS data are available, namely: Dennis, Emily, Katrina, Rita, and Wilma. Our analysis consists of examining CG lightning, cloud-top properties, and TTL water vapor (i.e., 100 and 147 mb) over the hurricane while it remains over water in the Tropical Americas region. We investigate daily as well as diurnal statistical properties. The hurricanes analyzed in this study showed that lightning activity is negatively correlated with minimum infrared brightness temperature and positively correlated with 100-mb water vapor. An examination of the maxima in water vapor observed over the hurricane not only shows larger magnitudes, but also larger differences between water vapor averages and water vapor maxima over the hurricane as lightning activity increases. Trajectory calculations are performed using the Flextra model in order to investigate the fate of the moister air masses found in the TTL.

Pittman, Jasna V.↗
NASA NTRS2012

Influences of the El Nino Southern Oscillation and the Pacific Decadal Oscillation on the Timing of the North American Spring

Detrended, modelled first leaf dates for 856 sites across North America for the period 1900-2008 are used to examine how the El Nino Southern Oscillation (ENSO) and the Pacific Decadal Oscillation (PDO) separately and together might influence the timing of spring. Although spring (mean March through April) ENSO and PDO signals are apparent in first leaf dates, the signals are not statistically significant (at a 95% confidence level (p <0.05)) for most sites. The most significant ENSO/PDO signal in first leaf dates occurs for El Nino and positive PDO conditions. An analysis of the spatial distributions of first leaf dates for separate and combined ENSO/PDO conditions features a northwest-southeast dipole that is significantly (at p <0.05) different than the distributions for neutral conditions. The nature of the teleconnection between Pacific SST's and first leaf dates is evident in comparable composites for detrended sea level pressure (SLP) in the spring months. During positive ENSO/PDO, there is an anomalous flow of warm air from the southwestern US into the northwestern US and an anomalous northeasterly flow of cold air from polar regions into the eastern and southeastern US. These flow patterns are reversed during negative ENSO/PDO. Although the magnitudes of first leaf date departures are not necessarily significantly related to ENSO and PDO, the spatial patterns of departures are significantly related to ENSO and PDO. These significant relations and the long-lived persistence of SSTs provide a potential tool for forecasting the tendencies for first leaf dates to be early or late.

McCabe, Gregory J.↗
NASA NTRS1987

Statistics and dynamics of persistent anomalies

Persistent anomalies with recurrent spatial patterns play an important role in the atmosphere's low-frequency variability. A connection between statistical and dynamical methods of description and prediction of persistent anomalies is established by computing and analyzing the empirical orthogonal functions (EOFs) in a simple deterministic model, on the one hand, and in Southern Hemisphere geopotential heights, on the other. The dynamical model is governed by the fully nonlinear, equivalent-barotropic vorticity equation on the sphere, with simplified forcing, dissipation and topography. The Southern Hemisphere data consist in gridded daily maps of 500 mb heights from June 1972 to July 1983. Two types of persistent anomalies appear in this time series, both having a strong wavenumber-three component; they differ by the value of the constant phase of this wave and by the strength of the wavenumber-one component. The first two EOFs bear a striking resemblance to these two patterns. It is concluded that the dynamical interpretation of EOFs is their pointing from the time mean to the most populated regions of the system's phase space. Pursuing this interpretation, a Markov-chain formulation of transitions from one persistent anomaly regime to another is introduced, and the implications for long-range forecasting are discussed.

Mo, Kingtse C.↗
NASA NTRS2015

Rao-Blackwellization for Adaptive Gaussian Sum Nonlinear Model Propagation

When dealing with imperfect data and general models of dynamic systems, the best estimate is always sought in the presence of uncertainty or unknown parameters. In many cases, as the first attempt, the Extended Kalman filter (EKF) provides sufficient solutions to handling issues arising from nonlinear and non-Gaussian estimation problems. But these issues may lead unacceptable performance and even divergence. In order to accurately capture the nonlinearities of most real-world dynamic systems, advanced filtering methods have been created to reduce filter divergence while enhancing performance. Approaches, such as Gaussian sum filtering, grid based Bayesian methods and particle filters are well-known examples of advanced methods used to represent and recursively reproduce an approximation to the state probability density function (pdf). Some of these filtering methods were conceptually developed years before their widespread uses were realized. Advanced nonlinear filtering methods currently benefit from the computing advancements in computational speeds, memory, and parallel processing. Grid based methods, multiple-model approaches and Gaussian sum filtering are numerical solutions that take advantage of different state coordinates or multiple-model methods that reduced the amount of approximations used. Choosing an efficient grid is very difficult for multi-dimensional state spaces, and oftentimes expensive computations must be done at each point. For the original Gaussian sum filter, a weighted sum of Gaussian density functions approximates the pdf but suffers at the update step for the individual component weight selections. In order to improve upon the original Gaussian sum filter, Ref. [2] introduces a weight update approach at the filter propagation stage instead of the measurement update stage. This weight update is performed by minimizing the integral square difference between the true forecast pdf and its Gaussian sum approximation. By adaptively updating each component weight during the nonlinear propagation stage an approximation of the true pdf can be successfully reconstructed. Particle filtering (PF) methods have gained popularity recently for solving nonlinear estimation problems due to their straightforward approach and the processing capabilities mentioned above. The basic concept behind PF is to represent any pdf as a set of random samples. As the number of samples increases, they will theoretically converge to the exact, equivalent representation of the desired pdf. When the estimated qth moment is needed, the samples are used for its construction allowing further analysis of the pdf characteristics. However, filter performance deteriorates as the dimension of the state vector increases. To overcome this problem Ref. [5] applies a marginalization technique for PF methods, decreasing complexity of the system to one linear and another nonlinear state estimation problem. The marginalization theory was originally developed by Rao and Blackwell independently. According to Ref. [6] it improves any given estimator under every convex loss function. The improvement comes from calculating a conditional expected value, often involving integrating out a supportive statistic. In other words, Rao-Blackwellization allows for smaller but separate computations to be carried out while reaching the main objective of the estimator. In the case of improving an estimator's variance, any supporting statistic can be removed and its variance determined. Next, any other information that dependents on the supporting statistic is found along with its respective variance. A new approach is developed here by utilizing the strengths of the adaptive Gaussian sum propagation in Ref. [2] and a marginalization approach used for PF methods found in Ref. [7]. In the following sections a modified filtering approach is presented based on a special state-space model within nonlinear systems to reduce the dimensionality of the optimization problem in Ref. [2]. First, the adaptive Gaussian sum propagation is explained and then the new marginalized adaptive Gaussian sum propagation is derived. Finally, an example simulation is presented.

state estimation↗
NASA NTRS1995

Sea-air of CO2 in the North Pacific using shipboard and satellite data

A method has been developed to produce high-resolution maps of pCO2 in surface water for the North Pacific using satellite sea surface temperature (SST) data and statistical relationships between measured pCO2 and temperature. In the subtropical North Pacific the pCO in seawater is controlled primarily by temperature. Accordingly, pCO2 values that are calculated from the satellite SST data have good agreement with the measured values (rms deviation of +/- microatm). In the northwestern subpolar region the pCO2 is controlled not only by temperature, but also by significant seasonal changes in the total CO2 concentration, which are caused by seasonal changes in primary production, mixing with subsurface waters and sea-air exchange. Consequently, the parameterization of oceanic p CO2 based on SST data alone is not totally successful in the northwestern region (rms deviation of +/- 40 microatm). The use of additional satellite products, such as wind and ocean color data, as planned for a future study, is considered necessary to account for the pCO2 variability caused by seasonal changes in the total CO2 concentration. The net CO2 flux for the area of the North Pacific included in this study (north of 10 deg N) has been calculated using the monthly pCO2 distributions computed, and monthly wind speeds from the European Centre for Medium-Range Weather Forecasts. The region is found to be a net source to the atmosphere of 1.9 x 10(exp 12) to 5.8 x 10(exp 12) moles of CO2 per year (or 0.02-0.07 Gt C/yr), most of the outflux occurring in the subtropics.

Stephens, Mark P.↗
NASA NTRS2012

Incorporation of the Mass Concentration and the New Snow Albedo Schemes into the Global Forecasting Model, GEOS-5 and the Impact of the New Schemes over Himalayan Glaciers

Recently the issue on glacier retreats comes up and many factors should be relevant to the issue. The absorbing aerosols such as dust and black carbon (BC) are considered to be one of the factors. After they deposited onto the snow surface, it will reduce snow albedo (called snow darkening effect) and probably contribute to further melting of glacier. The Goddard Earth Observing System version 5 (GEOS-5) has developed at NASA/GSFC. However, the original snowpack model used in the land surface model in the GEOS-5 did not consider the snow darkening effect. Here we developed the new snow albedo scheme which can consider the snow darkening effect. In addition, another scheme on calculating mass concentrations on the absorbing aerosols in snowpack was also developed, in which the direct aerosol depositions from the chemical transport model in the GEOS-5 were used. The scheme has been validated with the observed data obtained at backyard of the Institute of Low Temperature Science, Hokkaido University, by Dr. Teruo Aoki (Meteorological Research Institute) et aL including me. The observed data was obtained when I was Ph.D. candidate. The original GEOS-5during 2007-2009 over the Himalayas and Tibetan Plateau region showed more reductions of snow than that of the new GEOS-5 because the original one used lower albedo settings. On snow cover fraction, the new GEOS-5 simulated more realistic snow-covered area comparing to the MODIS snow cover fraction. The reductions on snow albedo, snow cover fraction, and snow water equivalent were seen with statistically significance if we consider the snow darkening effect comparing to the results without the snow darkening effect. In the real world, debris cover, inside refreezing process, surface flow of glacier, etc. affect glacier mass balance and the simulated results immediately do not affect whole glacier retreating. However, our results indicate that some surface melting over non debris covered parts of the glacier would be explained by the snow darkening effect. Further discussion and observations are necessary to assess the glacier issue.

Yasunari, Teppei↗
NASA NTRS2015

Combined Effect of El Nino Southern Oscillation and Atlantic Multidecadal Oscillation on Lake Chad Level Variability Region

In this study, the combined effect of the Atlantic Multidecadal Oscillation (AMO) and El Niño Southern Oscillation (ENSO) on the Lake Chad (LC) level variability is explored. Our results show that the lake level at the Bol monitoring station has a statistically significant correlation with precipitation (R2 = 0.6, at the 99.5% confidence level). The period between the late 1960s and early 1970s marked a turning point in the response of the regional rainfall to climatic drivers, thereby severely affecting the LC level. Our results also suggest that the negative impact of the cold phase of AMO on Sahel precipitation masks and supersedes the positive effect of La Niña in the early the 1970s. The drop in the size of LC level from 282.5 m in the early 1960s to about 278.1 m in 1983/1984 was the largest to occur within the period of study (1900-2010) and coincides with the combined cold phase of AMO and strong El Niño phase of ENSO. Further analyses show that the current warm phase of AMO and increasing La Niña episodes appear to be playing a major role in the increased precipitation in the Sahel region. The LC level is responding to this increase in precipitation by a gradual recovery, though it is still below the levels of the 1960s. This understanding of the AMO-ENSO-rainfall-LC level association will help in forecasting the impacts of similar combined episodes in the future. These findings also have implications for long-term water resources management in the LC region.

ENSO↗
NASA NTRS2019

Validation of Regional CO2 Concentrations in the ECMWF Real-Time Analysis and Carbon-Tracker Reanalysis with Airborne Observations from ACT-A Field Campaign

Through verifying against hundreds of hours of airborne in-situ measurements from the NASA-sponsored Atmospheric Carbon and Transport – America (ACT-A) field campaign, this study systematically examines the regional uncertainties and biases of the carbon dioxide (CO2) concentrations from two of the state-of-the-art global analysis products, namely the real-time analysis from the European Center (EC) for Medium Range Forecasting and NOAA’s near real-time Carbon Tracker (CT) reanalysis. It is found that both the EC and CT-NRT analyses agree reasonably well with the independent ACT-A flight-level CO2 measurements in the free troposphere but the uncertainties are considerably larger in the boundary layer during both the summer months of 2016 and the winter months of 2017. There are also strong variabilities in accuracy and bias between seasons, and across three different subregions in the United States (Mid-Atlantic, Midwest and South). Overall, the analysis uncertainties of the EC and CT-NRT analyses in terms of root-mean square deviations against airborne data are comparable to each other, both of which are between 1-2 ppm in the free troposphere but can be as large as 10 ppm near the surface, which are grossly consistent with the difference between the two analyses. The current study not only provides systematic uncertainty estimates for both analysis products over North America but also demonstrated that these two independent estimates can be used to approximate the overall regional CO2 analysis uncertainties. Both statistics are important in future studies in quantifying the uncertainties of regional carbon concentration and flux estimates, as well as in assessing the impact of regional transport through more refined regional modeling and analysis systems.

Chen, Hans W.↗
NASA NTRS1998

T wave alternans as a predictor of recurrent ventricular tachyarrhythmias in ICD recipients: prospective comparison with conventional risk markers

INTRODUCTION: The current standard for arrhythmic risk stratification is electrophysiologic (EP) testing, which, due to its invasive nature, is limited to patients already known to be at high risk. A number of noninvasive tests, such as determination of left ventricular ejection fraction (LVEF) or heart rate variability, have been evaluated as additional risk stratifiers. Microvolt T wave alternans (TWA) is a promising new risk marker. Prospective evaluation of noninvasive risk markers in low- or moderate-risk populations requires studies involving very large numbers of patients, and in such studies, documentation of the occurrence of ventricular tachyarrhythmias is difficult. In the present study, we identified a high-risk population, recipients of an implantable cardioverter defibrillator (ICD), and prospectively compared microvolt TWA with invasive EP testing and other risk markers with respect to their ability to predict recurrence of ventricular tachyarrhythmias as documented by ICD electrograms. METHODS AND RESULTS: Ninety-five patients with a history of ventricular tachyarrhythmias undergoing implantation of an ICD underwent EP testing, assessment of TWA, as well as determination of LVEF, baroreflex sensitivity, signal-averaged ECG, analysis of 24-hour Holter monitoring, and QT dispersion from the 12-lead surface ECG. The endpoint of the study was first appropriate ICD therapy for electrogram-documented ventricular fibrillation or tachycardia during follow-up. Kaplan-Meier survival analysis revealed that TWA (P < 0.006) and LVEF (P < 0.04) were the only significant univariate risk stratifiers. EP testing was not statistically significant (P < 0.2). Multivariate Cox regression analysis revealed that TWA was the only statistically significant independent risk factor. CONCLUSIONS: Measurement of microvolt TWA compared favorably with both invasive EP testing and other currently used noninvasive risk assessment methods in predicting recurrence of ventricular tachyarrhythmias in ICD recipients. This study suggests that TWA might also be a powerful tool for risk stratification in low- or moderate-risk patients, and needs to be prospectively evaluated in such populations.

Non-NASA Center↗
NASA NTRS2008

A Goddard Multi-Scale Modeling System with Unified Physics

Numerical cloud resolving models (CRMs), which are based the non-hydrostatic equations of motion, have been extensively applied to cloud-scale and mesoscale processes during the past four decades. Recent GEWEX Cloud System Study (GCSS) model comparison projects have indicated that CRMs agree with observations in simulating various types of clouds and cloud systems from different geographic locations. Cloud resolving models now provide statistical information useful for developing more realistic physically based parameterizations for climate models and numerical weather prediction models. It is also expected that Numerical Weather Prediction (NWP) and regional scale model can be run in grid size similar to cloud resolving model through nesting technique. Current and future NASA satellite programs can provide cloud, precipitation, aerosol and other data at very fine spatial and temporal scales. It requires a coupled global circulation model (GCM) and cloud-scale model (termed a szrper-parameterization or multi-scale modeling -framework, MMF) to use these satellite data to improve the understanding of the physical processes that are responsible for the variation in global and regional climate and hydrological systems. The use of a GCM will enable global coverage, and the use of a CRM will allow for better and more sophisticated physical parameterization. NASA satellite and field campaign can provide initial conditions as well as validation through utilizing the Earth Satellite simulators. At Goddard, we have developed a multi-scale modeling system with unified physics. The modeling system consists a coupled GCM-CRM (or MMF); a state-of-the-art weather research forecast model (WRF) and a cloud-resolving model (Goddard Cumulus Ensemble model). In these models, the same microphysical schemes (2ICE, several 3ICE), radiation (including explicitly calculated cloud optical properties), and surface models are applied. In addition, a comprehensive unified Earth Satellite simulator has been developed at GSFC, which is designed to fully utilize the multi-scale modeling system. A brief review of the multi-scale modeling system with unified physics/simulator and examples is presented in this article.

Tao, W.K.↗
NASA NTRS2008

Completion of the Edward Air Force Base Statistical Guidance Wind Tool

The goal of this task was to develop a GUI using EAFB wind tower data similar to the KSC SLF peak wind tool that is already in operations at SMG. In 2004, MSFC personnel began work to replicate the KSC SLF tool using several wind towers at EAFB. They completed the analysis and QC of the data, but due to higher priority work did not start development of the GUI. MSFC personnel calculated wind climatologies and probabilities of 10-minute peak wind occurrence based on the 2-minute average wind speed for several EAFB wind towers. Once the data were QC'ed and analyzed the climatologies were calculated following the methodology outlined in Lambert (2003). The climatologies were calculated for each tower and month, and then were stratified by hour, direction (10" sectors), and direction (45" sectors)/hour. For all climatologies, MSFC calculated the mean, standard deviation and observation counts of the Zminute average and 10-minute peak wind speeds. MSFC personnel also calculated empirical and modeled probabilities of meeting or exceeding specific 10- minute peak wind speeds using PDFs. The empirical PDFs were asymmetrical and bounded on the left by the 2- minute average wind speed. They calculated the parametric PDFs by fitting the GEV distribution to the empirical distributions. Parametric PDFs were calculated in order to smooth and interpolate over variations in the observed values due to possible under-sampling of certain peak winds and to estimate probabilities associated with average winds outside the observed range. MSFC calculated the individual probabilities of meeting or exceeding specific 10- minute peak wind speeds by integrating the area under each curve. The probabilities assist SMG forecasters in assessing the shuttle FR for various Zminute average wind speeds. The A M ' obtained the processed EAFB data from Dr. Lee Bums of MSFC and reformatted them for input to Excel PivotTables, which allow users to display different values with point-click-drag techniques. The GUI was created from the PivotTables using VBA code. It is run through a macro within Excel and allows forecasters to quickly display and interpret peak wind climatology and probabilities in a fast-paced operational environment. The GUI was designed to look and operate exactly the same as the KSC SLF tool since SMG forecasters were already familiar with that product. SMG feedback was continually incorporated into the GUI ensuring the end product met their needs. The final version of the GUI along with all climatologies, PDFs, and probabilities has been delivered to SMG and will be put into operational use.

Dreher, Joseph G.↗
NASA NTRS1995

Estimating the urban bias of surface shelter temperatures using upper-air and satellite data. Part 1: Development of models predicting surface shelter temperatures

Multiple regression techniques were used to predict surface shelter temperatures based on the time period 1986-89 using upper-air data from the European Centre for Medium-Range Weather Forecasts (ECMWF) to represent the background climate and site-specific data to represent the local landscape. Global monthly mean temperature models were developed using data from over 5000 stations available in the Global Historical Climate Network (GHCN). Monthly maximum, mean, and minimum temperature models for the United States were also developed using data from over 1000 stations available in the U.S. Cooperative (COOP) Network and comparative monthly mean temperature models were developed using over 1150 U.S. stations in the GHCN. Three-, six-, and full-variable models were developed for comparative purposes. Inferences about the variables selected for the various models were easier for the GHCN models, which displayed month-to-month consistency in which variables were selected, than for the COOP models, which were assigned a different list of variables for nearly every month. These and other results suggest that global calibration is preferred because data from the global spectrum of physical processes that control surface temperatures are incorporated in a global model. All of the models that were developed in this study validated relatively well, especially the global models. Recalibration of the models with validation data resulted in only slightly poorer regression statistics, indicating that the calibration list of variables was valid. Predictions using data from the validation dataset in the calibrated equation were better for the GHCN models, and the globally calibrated GHCN models generally provided better U.S. predictions than the U.S.-calibrated COOP models. Overall, the GHCN and COOP models explained approximately 64%-95% of the total variance of surface shelter temperatures, depending on the month and the number of model variables. In addition, root-mean-square errors (rmse's) were over 3 C for GHCN models and over 2 C for COOP models for winter months, and near 2 C for GHCN models and near 1.5 C for COOP models for summer months.

Epperson, David L.↗
NASA NTRS2008

The Global Modeling and Assimilation Office (GMAO) 4d-Var and its Adjoint-based Tools

The fifth generation of the Goddard Earth Observing System (GEOS-5) Data Assimilation System (DAS) is a 3d-var system that uses the Grid-point Statistical Interpolation (GSI) system developed in collaboration with NCEP, and a general circulation model developed at Goddard, that includes the finite-volume hydrodynamics of GEOS-4 wrapped in the Earth System Modeling Framework and physical packages tuned to provide a reliable hydrological cycle for the integration of the Modern Era Retrospective-analysis for Research and Applications (MERRA). This MERRA system is essentially complete and the next generation GEOS is under intense development. A prototype next generation system is now complete and has been producing preliminary results. This prototype system replaces the GSI-based Incremental Analysis Update procedure with a GSI-based 4d-var which uses the adjoint of the finite-volume hydrodynamics of GEOS-4 together with a vertical diffusing scheme for simplified physics. As part of this development we have kept the GEOS-5 IAU procedure as an option and have added the capability to experiment with a First Guess at the Appropriate Time (FGAT) procedure, thus allowing for at least three modes of running the data assimilation experiments. The prototype system is a large extension of GEOS-5 as it also includes various adjoint-based tools, namely, a forecast sensitivity tool, a singular vector tool, and an observation impact tool, that combines the model sensitivity tool with a GSI-based adjoint tool. These features bring the global data assimilation effort at Goddard up to date with technologies used in data assimilation systems at major meteorological centers elsewhere. Various aspects of the next generation GEOS will be discussed during the presentation at the Workshop, and preliminary results will illustrate the discussion.

Todling, Ricardo↗
NASA NTRS2011

Magnetic Clouds at/near the 2007 - 2009 Solar Minimum: Frequency of Occurrence and Some Unusual Properties

Magnetic clouds (MCs) have been identified for the period 2007 2009 (at/near the recent solar minimum) from Wind data, then confirmed through MC parameter fitting using a force-free model. A dramatic increase in the frequency of occurrence of these events took place from the two early years of 2007 (with five MCs) and 2008 (one MC) compared to 2009 (12 MCs). This pattern approximately mirrors the occurrence-frequency profile that was observed over a three-year interval 12 years earlier, with eight events in 1995, four in 1996, and 17 in 1997, but decreased overall by a factor of 0.62 in number. However, the average estimated axial field strength taken over all of the 18 events of 2007 - 2009 (called the "recent period" here) was only 11.0 nT, whereas |BO| for the 29 events of 1995 - 1997 (called the "earlier period" ) was 16.5 nT. This 33% average drop in |BO| is more or less consistent with the decreased three-year average interplanetary magnetic field intensity between these two periods, which shows a 23% drop. In the earlier period, the MCs were clearly of mixed types but predominantly of the South-to-North type, whereas those in the recent period are almost exclusively the North-to-South type; this change is consistent with global solar field changes predicted by Bothmer and Rust (Geophys. Monogr. Ser. 99, 139, 1997). As we have argued in earlier work (Lepping and Wu, J. Geophys. Res. 112, A10103, 2007), this change should make it possible to carry out (accurate short-term) magnetic storm forecasting by predicting the latter part of an MC from the earlier part, using a good MC parameter-fitting model with real-time data from a spacecraft at L1, for example. The recent set s average duration is 15.2 hours, which is a 27% decrease compared to that of the earlier set, which had an average duration of 20.9 hours. In fact, all physical aspects of the recent MC set are shown to drop with respect to the earlier set; e.g., as well as the average internal magnetic field drop, the recent set had a somewhat low average speed of 379 km/s (5% drop), and the average diameter had a 24% drop. Hence, compared to the earlier set, the recent set consists of events that are smaller, slightly slower, and weaker in every respect (and fewer in number), but in a relative sense the two three-year sets have similar frequency-of-occurrence profiles. It is also interesting that the two sets have almost the same average axial inclinations, i.e., axial latitude approx. = 31deg (in GSE). These MC characteristics are compared to relevant solar features and their changes. A preliminary assessment of the statistics on possible shocks and pressure pulses upstream of these recent MCs yields the following: About 28% of the MCs, at most, had shocks, and 33% had shocks and/or pressure pulses. These are low values, since typically the percentage of cases with shocks is about 50%, and the percentage with shocks and/or pressure pulses is usually about 75%.

Lepping. R. P.↗
NASA NTRS2011

Analysis of Ultra High Resolution Sea Surface Temperature Level 4 Datasets

Sea surface temperature (SST) studies are often focused on improving accuracy, or understanding and quantifying uncertainties in the measurement, as SST is a leading indicator of climate change and represents the longest time series of any ocean variable observed from space. Over the past several decades SST has been studied with the use of satellite data. This allows a larger area to be studied with much more frequent measurements being taken than direct measurements collected aboard ship or buoys. The Group for High Resolution Sea Surface Temperature (GHRSST) is an international project that distributes satellite derived sea surface temperatures (SST) data from multiple platforms and sensors. The goal of the project is to distribute these SSTs for operational uses such as ocean model assimilation and decision support applications, as well as support fundamental SST research and climate studies. Examples of near real time applications include hurricane and fisheries studies and numerical weather forecasting. The JPL group has produced a new 1 km daily global Level 4 SST product, the Multiscale Ultrahigh Resolution (MUR), that blends SST data from 3 distinct NASA radiometers: the Moderate Resolution Imaging Spectroradiometer (MODIS), the Advanced Very High Resolution Radiometer (AVHRR), and the Advanced Microwave Scanning Radiometer ? Earth Observing System(AMSRE). This new product requires further validation and accuracy assessment, especially in coastal regions.We examined the accuracy of the new MUR SST product by comparing the high resolution version and a lower resolution version that has been smoothed to 19 km (but still gridded to 1 km). Both versions were compared to the same data set of in situ buoy temperature measurements with a focus on study regions of the oceans surrounding North and Central America as well as two smaller regions around the Gulf Stream and California coast. Ocean fronts exhibit high temperature gradients (Roden, 1976), and thus satellite data of SST can be used in the detection of these fronts. In this case, accuracy is less of a concern because the primary focus is on the spatial derivative of SST. We calculated the gradients for both versions of the MUR data set and did statistical comparisons focusing on the same regions.

sea surface temperature SST)↗
NASA NTRS2019

Impacts of the Mount Pinatubo eruption on ENSO in the GEOS seasonal-to-subseasonal forecasting system

The eruption of Mount Pinatubo in June 1991 introduced a perturbation of the Earth's global energy budget by increasing the stratospheric aerosol loading by an order of magnitude, with effects on the global climate. In this presentation we analyze the effects of the Mt. Pinatubo eruption on the seasonal forecast performed with Goddard Earth Observing System Seasonal-to Subseasonal (GEOS-S2S) system, an Earth System Model that includes an interactive ocean and a bulk aerosol model coupled to radiation. We performed 10-member ensembles for the year after the eruption (June 1991-May 1992) at ~0.5 horizontal resolution, with and without the inclusion of the Mt. Pinatubo eruption. In GEOS-S2S, the eruption leads to ta strengthening of El Nino peaking in January 1992. The strengthening is mainly due to the weakening of the trade winds, which is caused by a attening of the temperature gradient across the Pacic due to a differential response to the volcanic forcing between the central and eastern Pacic (ocean-dynamical thermostat). This response largely depends on the assumed size for the volcanic aerosols. Indeed, we performed simulations assuming a volcanic aerosol effective radius of 0.35 m (similar to tropospheric aerosol, and the default in GEOS) and 0.6 m (closer to observations of volcanic aerosol from Pinatubo-sized eruptions). We nd that in the latter case the tropical radiative forcing is lower, since smaller aerosols scatter shortwave radiation more eciently than larger ones. Accordingly, the impact on ENSO is not statistically signicant when a larger and more realistic particle radius is assumed.

Aquila, Valentina↗
NASA NTRS1999

Characterization of Climate Change and Variability with GPS

We compared zonal mean specific humidity derived from the 21 June-4 July 1995 Global Positioning System (GPS)/MET occultation observations with that derived from the European Center for Medium-Range Weather Forecasts (ECMWF) global analyses. The GPS/MET results indicate a drier troposphere, especially near the subtropical tradewind inversion. A small, moist bias in the GPS/MET upper northern-hemisphere troposphere compared to ECMWF may be due to a small radiosonde temperature bias. A diagram shows the difference (g/kg) between the GPS/MET zonal mean specific humidity and that for June-August derived from 1963-1973 radiosondes. Although the observing period is short, GPS and ECMWF results both indicate a significantly wetter boundary layer at most latitudes consistent with decadal trends observed in radiosonde data. GPS/MET results exhibit higher tropical convective available potential energy (CAPE), suggesting a more vigorous tropical Hadley circulation. Drier, free troposphere air in the descending branches of the Hadley circulation is due in part to a moist radiosonde bias but may also reflect some negative moisture feedback. Using 1992-1997 ground GPS observations and recent advancements in GPS technology, we removed an apparent altimetric drift (-1.2 +/- 0.4 mm/yr) due to columnar water vapor from the Topography (Ocean) Experiment (TOPEX) microwave radiometer, which brought the TOPEX mean sea level change estimates into better agreement with historical tide gauge records, suggesting global mean sea level is rising at a rate of 1.5-2.0 mm/yr. We can also discern a statistically significant increase of 0.2 +/- 0.1 kg/square m/yr in mean columnar water vapor over the ocean from 1992-1997. Optimal fingerprinting can be used for the detection and attribution of tropospheric warming due to an anthropogenic greenhouse. Optimal fingerprinting distinguishes between different types of signals according to their spatial and temporal patterns, while minimizing the influence of natural climate variability. S. Leroy concludes that the signal-to-noise ratio of global warming detection increases by unity approximately every 10 years if a single oceanic region is chosen. Less time for detection is likely when many global regions are considered simultaneously. GPS occultation constellations allow the possibility of detecting small changes in upper air temperature with inconsequential calibration errors, making occultation an ideal data type for global warming detection studies. Our initial study of a 22-GHz satellite-satellite occultation system predicts upper troposphere moisture sensitivities of 3-5 ppmv and 1-2 percent in the middle and lower troposphere. Additional information contained in original.

Kursinski, R.↗