SEARCH · Engineering Papers
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.
Factors affecting thermal structure of ocean, error sources in surface temperature, and calculated energy exchange effects at surface to forecast ocean temperature structure
A nine-layer, primitive equation (PE) model of the global atmosphere developed at the Goddard Institute for Space Studies (GISS) has been used to generate six 48-hr forecasts during December 1972 and January 1973. Operational analyses north of 18 N and experimental global analyses made available by the National Meteorological Center (NMC), NOAA, were used as the initial conditions; the operational analyses were used to verify the forecasts at 12-hr intervals over the Northern Hemisphere north of 22 N. The combined analyses were used to verify the forecasts in the global domain. RMS errors of the sea-level pressure, 1000-mb heights, and vector geostrophic winds, and 500-mb heights and vector geostrophic winds indicate that the GISS forecasts have skill comparable to those made by operational PE models.
A new method for performing Monte Carlo forecasts is introduced. The method, called dynamic stratification, selects initial perturbations based on a stratification of the error distribution. A simple implementation is presented in which the error distribution used for the stratification is estimated from a linear model derived from a large ensemble of 12-h forecasts with the full dynamic model. The stratification thus obtained is used to choose a small subsample of initial states with which to perform the dynamical Monte Carlo forecasts. Several test cases are studied using a simple two-level general circulation model with uncertain initial conditions. It is found that the method provides substantial reductions in the sampling error of the forecast mean and variance when compared to the more traditional approach of choosing the initial perturbations at random. The degree of improvement, however, is sensitive to the nature of the initial error distribution and to the base state. In practice the method may be viable only if the computational burden involved in obtaining an adequate estimate of the error distribution is shared with the data-assimilation procedure.
The Space Shuttle is an extremely weather sensitive vehicle with very restrictive constraints for both launches and landings. The most important difference between Shuttle and normal aircraft landings is that the Shuttle has no go-around capability once it begins its decent into the earth's atmosphere. The de-orbit burn decision is generally made approximately 90 minutes before landing requiring a forecast with little room for error. Because of the Shuttle's rapid re-entry to earth, the pilot must be able to see all runway and visual navigation aids from high altitude to land the Shuttle. In addition, the heat resistant tiles which are used to protect the Shuttle during its re-entry into the earth's atmosphere are extremely sensitive to any type of precipitation. Extensive damage to these tiles could occur if the Shuttle passes through any cloud that contains precipitation size particles. To help guard against changing weather conditions or any type of weather problems that might occur prior to landing, flight rules have been developed as guidelines for all landings. Although the rules vary depending on the location of the landing (Kennedy Space Center or Edwards AFB), length of mission, and weight of vehicle, most of the rules can be condensed into 4 major groupings. These are: (1) Cloud ceilings should not be less than 3048 m (10,000 feet), (2) Visibility should not be less than 13 km (7 nm), (3) Cross-wind no greater than 5-8 m/s (10-15 knots); and (4) No showers or thunderstorms at or within 56 km (30 nm) of the Shuttle Landing Facility. This study consisted of developing a climatological database of the Shuttle Landing Facility (SLF) surface observations and performing an analysis of observed conditions one and two hours subsequent to given conditions at the SLF to help analyze the 0.2 cloud cover rule. Particular emphasis was placed on Shuttle landing weather violations and the amounts of cloud cover below 3048 m (10,000 ft.). This analysis has helped to determine the best and worst times to land the Shuttle at KSC. In addition, nomograms have been developed to help forecasters make cloud cover forecasts for End of Mission (EOM) and Return to Launch Site (RTLS) at KSC. Results of categorizing this data by month, season, time of day, and surface and upper-air wind direction are presented.
In a brief evaluation of software-related considerations, it is found that suitable approaches for software development depend to a large degree on the characteristics of the particular project involved. An analysis is conducted of development problems in an environment in which ground support software is produced for spacecraft control. The amount of work involved is in the range from 6 to 10 man-years. Attention is given to a general project summary, a programmer/analyst survey, a component summary, a component status report, a resource summary, a change report, a computer program run analysis, aspects of data collection on a smaller scale, progress forecasting, problems of overhead, and error analysis.
The effects of tropical wind data from the FGGE and tropical latent heating on numerical modeling of ultralong waves are considered in a two-part study. The model studied is the global fourth-order GLAS general circulation model, an energy-conserving format with horizontal differences calculated with fourth-order accuracy. Data assimilation experiments were performed with and without the wind data, with account taken of eastward and northward wind components, the geopotential height, and the relative humidity, all over pressure surfaces. The initial conditions were used to generate six pairs of forecasts, and the tropical wind error decreased after two days of prediction when the initial conditions contained the wind data. The deviations from the measured planetary wave data were attributed to differences in the initial rotational wind field, which varied on a three-day basis. The latent heat initial data had a five-day period and extended its influence beyond the tropical zone. The tropical heat sources sustained the tropical westerlies in the GLAS model, and removal of the tropical heat sources reversed the wind to easterlies.
Space Shuttle landings at Kennedy Space Center (KSC) are of special concern to NASA's landing community because of Florida's rapidly changing weather conditions. Since a large number of Shuttle landing attempts occur in the morning hours (just after sunrise) fog and stratus development are a problem. The deorbit burn decision for a landing at KSC is typically made 90 minutes before Shuttle touchdown. In that 90 minutes weather conditions can change very rapidly. Fog to the west of KSC an advect in and reduce visibility to less than 7 miles. The most important difference between Shuttle and normal aircraft landings is that the Shuttle has no go-around capability requiring a forecast with little room for error. To help guard against rapidly changing weather conditions, flight rules have been developed as guidelines for all landings. This paper concerns fog development that would affect less than 7-statue mile visibility rule which is in effect for End-Of-Mission (EOM) Shuttle landings at KSC (Rule 4-64(A)). Data used for this analysis included hourly surface observations at the X68 Shuttle Landing Facility (SLF) and upper-air observations form the CCAFS (Cape Canaveral Air Force Station-72794) rawinsonde site for the five year period, 1986 to 1990. This investigation focused on rapidly developing fog or stratus that developed between decision time and landing.
The Madden-Julian Oscillation (MJO) is known as the dominant mode of tropical intraseasonal variability and has an important role in the coupled-atmosphere system. This study used twin numerical model experiments to investigate the influence of the MJO activity on weather predictability in the midlatitudes of the Northern Hemisphere during boreal winter. The National Aeronautics and Space Administration (NASA) Goddard laboratory for the Atmospheres (GLA) general circulation model was first used in a 10-yr simulation with fixed climatological SSTs to generate a validation data set as well as to select initial conditions for active MJO periods and Null cases. Two perturbation numerical experiments were performed for the 75 cases selected [(4 MJO phases + Null phase) _ 15 initial conditions in each]. For each alternative initial condition, the model was integrated for 90 days. Mean anomaly correlations in the midlatitudes of the Northern Hemisphere (2O deg N_60 deg.N) and standardized root-mean-square errors were computed to validate forecasts and control run. The analyses of 500-hPa geopotential height, 200-hPa Streamfunction and 850-hPa zonal wind component systematically show larger predictability during periods of active MJO as opposed to quiescent episodes of the oscillation.
It is known that the General Circulation Models (GCMs) have sufficient resolution to accurately simulate hurricane near-eye structure and intensity. To overcome this limitation, the mesoscale-resolving finite-element GCM (fvGCM) has been experimentally deployed on the NASA Columbia supercomputer, and its performance is evaluated choosing hurricane Katrina as an example in this study. On late August 2005 Katrina underwent two stages of rapid intensification and became the sixth most intense hurricane in the Atlantic. Six 5-day simulations of Katrina at both 0.25 deg and 0.125 deg show comparable track forecasts, but the 0,125 deg runs provide much better intensity forecasts, producing center pressure with errors of only +/- 12 hPa. The 0.125 deg simulates better near-eye wind distributions and a more realistic average intensification rate. A convection parameterization (CP) is one of the major limitations in a GCM, the 0.125 deg run with CP disabled produces very encouraging results.
Solar activity predictions using the data assimilation approach have demonstrated great potential to build reliable long-term forecasts of solar activity. In particular, it has been shown that the Ensemble Kalman Filter (EnKF) method applied to a non-linear dynamo model is capable of predicting solar activity up to one sunspot cycle ahead in time, as well as estimating the properties of the next cycle a few years before it begins. These developments assume an empirical relationship between the mean toroidal magnetic field flux and the sunspot number. Estimated from the sunspot number series, variations of the toroidal field have been used to assimilate the data into the Parker-Kleeorin-Ruzmakin (PKR) dynamo model by applying the EnKF method. The dynamo model describes the evolution of the toroidal and poloidal components of the magnetic field and the magnetic helicity. Full-disk magnetograms provide more accurate and complete input data by constraining both the toroidal and poloidal global field components, but these data are available only for the last four solar cycles. In this presentation, using the available magnetogram data, we discuss development of the methodology and forecast quality criteria (including forecast uncertainties and sources of errors). We demonstrate the influence of limited time series observations on the accuracy of solar activity predictions. We present EnKF predictions of the upcoming Solar Cycle 25 based on both the sunspot number series and observed magnetic fields and discuss the uncertainties and potential of the data assimilation approach.
Soil Moisture Active Passive (SMAP) mission L-band brightness temperature (Tb) observations are routinely assimilated into the Catchment land surface model to generate Level-4 Soil Moisture (L4_SM) estimates of global surface and root-zone soil moisture at 9-km, 3-hourly resolution with ~2.5-day latency. The Catchment model in the L4_SM algorithm is driven with ¼-degree, hourly surface meteorological forcing data from the Goddard Earth Observing System (GEOS). Outside of Africa and the high latitudes, GEOS precipitation is corrected using Climate Prediction Center Unified (CPCU) gauge-based, ½-degree, daily precipitation. L4_SM soil moisture was previously shown to improve over land model-only estimates that use CPCU precipitation but no Tb assimilation (CPCU_SIM). Here, we additionally examine the skill of model-only (CTRL) and Tb assimilation-only (SMAP_DA) estimates derived without CPCU precipitation. Soil moisture is assessed versus in situ measurements in well-instrumented regions and globally through the Instrumental Variable (IV) method using independent soil moisture retrievals from the Advanced Scatterometer. At the in situ locations, SMAP_DA and CPCU_SIM have comparable soil moisture skill improvements relative to CTRL for the unbiased root-mean-square error (surface and root-zone) and correlation metrics (root-zone only). In the global average, SMAP Tb assimilation increases the surface soil moisture anomaly correlation by 0.10-0.11 compared to an increase of 0.02-0.03 from the CPCU-based precipitation corrections. The contrast is particularly strong in central Australia, where CPCU is known to have errors and observation-minus-forecast Tb residuals are larger when CPCU precipitation is used. Validation versus streamflow measurements in the contiguous U.S. reveals that CPCU precipitation provides most of the skill gained in L4_SM runoff estimates over CTRL.
Despite long-standing efforts, hydrologists still lack robust tools for calibrating land surface model (LSM) streamflow estimates within ungauged basins. Using surface soil moisture estimates from the Soil Moisture Active Passive Level 4 Soil Moisture (L4_SM) product, precipitation observations, and streamflow gauge measurements for 617 medium-scale (200-10,000 km2) basins in the contiguous United States, we measure the temporal (Spearman) rank correlation between antecedent (i.e., pre-storm) surface soil moisture (ASM) and the storm-scale runoff coefficient (RC; the fraction of storm-scale precipitation accumulation converted into streamflow). In humid and semi-humid basins, this rank correlation is shown to be sufficiently strong to allow for the substitution of storm-scale RC observations (available only in basins that are both lightly regulated and gauged) with high-quality ASM values (available quasi-globally from L4_SM) in streamflow calibration procedures. Using this principle, we define a new, basin-wise LSM streamflow calibration approach based on L4_SM alone and successfully apply it to identify LSM configurations that produce a high rank correlation with observed RC. However, since the approach cannot detect RC bias, it is less successful in identifying LSM configurations with low mean-absolute error. Plain Text Summary Accurately forecasting the fraction of rainfall that runs off into streams, as opposed to infiltrates into the soil, is critical for flash-flood prediction, water-resource monitoring, and tracking the transport of nutrients from agricultural fields into local waterways. Such forecasting is typically performed by hydrologic models that attempt to represent the physical processes responsible for surface runoff generation. However, to provide accurate streamflow forecasts, these models typically need to be calibrated against actual streamflow observations. This is problematic given the relatively poor, and declining, global availability of stream gauges. This paper presents a novel model calibration strategy that uses soil moisture from remote sensing and numerical modeling in place of streamflow observations during calibration. This transition has significant practical advantages because, unlike streamflow observations, the soil moisture data are continuously available across space. Our results demonstrate that this new approach can significantly improve hydrologic models within humid and semi-humid basins lacking sufficient ground-based instrumentation for traditional streamflow calibration.
Marine air temperature observations are collected by volunteer observing ships carrying automatic meteorological sensor packages. Close proximity of the package mount to a dark colored ship's deck is known to introduce a warm bias on sunny afternoons in summer under low wind conditions. This observation set is used as an influential constraint in the MERRA-2 reanalysis, as well as several other reanalyses, in the poorly sampled surface marine boundary layer. Here we explore the consequences of bias in the marine air temperature observation set during the 36 year period 1980-2015. Marine air temperature observations have mainly been collected along major shipping routes, with the densest coverage in the North Atlantic and North Pacific. In the first part of this presentation we review the historical observation coverage and compare air temperature and contemporaneous SST measurements to identify systematic differences. The comparison suggests the presence of diurnal bias in the air temperature observations peaking in subtropics and midlatitudes latitudes in mid-afternoon local time. In the second part we explore the impact of the marine air temperature observations on MERRA-2 through analysis of the terms in the temperature forecast tendency equation, the equation which describes the contributions to changes in temperature during each analysis cycle. Examination of the temperature forecast tendency equation allows us to separate systematic errors due to such processes and radiation, mixing, and advection and allows us to identify the major source of systematic forecast tendencies. The results show that the marine air temperature measurements introduce an error that can reach 0.4K/dy which also peaks in local afternoon at subtropical latitudes. This bias varies with season and has gradually reduced from the beginning of the MERRA-2 reanalysis in 1980. In the final part we conduct numerical experiments with a version of the MERRA-2 analysis system, comparing experiments when the marine air temperature observations are included or withheld from the analysis system.
Department of Mathematical Sciences, University of Alberta, Edmonton, Canada This paper describes the fundamental theory of the ensemble canonical correlation (ECC) algorithm for the seasonal climate forecasting. The algorithm is a statistical regression sch eme based on maximal correlation between the predictor and predictand. The prediction error is estimated by a spectral method using the basis of empirical orthogonal functions. The ECC algorithm treats the predictors and predictands as continuous fields and is an improvement from the traditional canonical correlation prediction. The improvements include the use of area-factor, estimation of prediction error, and the optimal ensemble of multiple forecasts. The ECC is applied to the seasonal forecasting over various parts of the world. The example presented here is for the North America precipitation. The predictor is the sea surface temperature (SST) from different ocean basins. The Climate Prediction Center's reconstructed SST (1951-1999) is used as the predictor's historical data. The optimally interpolated global monthly precipitation is used as the predictand?s historical data. Our forecast experiments show that the ECC algorithm renders very high skill and the optimal ensemble is very important to the high value.
We give an overview of the research at NASA in assimilating tropical rainfall and total precipitable water (TPW) retrievals derived from the TRMM Microwave Imager (TMI) and the Special Sensor Microwave/ Imager (SSM/I) instruments. Global analyses currently contain order-one errors in the primary fields of the hydrological cycle such as precipitation, evaporation, moisture, and the related cloud fields, especially in the tropics. We show that an effective strategy to assimilate tropical rainfall data is to use observations to compensate for errors in moisture tendencies produced by the assimilation model. Results show that assimilating TMI and SSM/I surface rainrates and TPW estimates improves shortrange forecasts and reduces state-dependent systematic errors in the hydrological cycle and related climate parameters such as cloud, radiation, and the large-scale circulation in the tropics. The study provides a demonstration of the potential of using rainfall and moisture observations derived from passive microwave instruments to improve the quality of 4-dimensional global datasets for climate analysis and weather forecasting applications.
The hazards of landfalling hurricanes are well known, but progress on improving the intensity forecasts of these deadly storms at landfall has been slow. Many cite a lack of high-resolution data sets taken inside the core of a hurricane, and the lack of reliable measurements in extreme conditions near the surface of hurricanes, as possible reasons why even the most state-of-the-art forecasting models cannot seem to forecast intensity changes better. The Hurricane Imaging Radiometer (HIRAD) is a new airborne microwave remote sensor for observing hurricanes, and is operated and researched by NASA Marshall Space Flight Center in partnership with the NOAA Atlantic Oceanographic and Meteorological Laboratory/Hurricane Research Division, the University of Central Florida, the University of Michigan, and the University of Alabama in Huntsville. This instrument?s purpose is to study the wind field of a hurricane, specifically observing surface wind speeds and rain rates, in what has traditionally been the most difficult areas for other instruments to study; the high wind and heavy rain regions. Dr. T. N. Krishnamurti has studied various data assimilation techniques for hurricane and monsoon rain rates, and this study builds off of results obtained from utilizing his style of physical initializations of rainfall observations, but obtaining reliable observations in heavy rain regions has always presented trouble to our research of high-resolution rainfall forecasting. Reliable data from these regions at such a high resolution and wide swath as HIRAD provides is potentially very valuable to mesoscale forecasting of hurricane intensity. This study shows how the data assimilation technique of Ensemble Kalman Filtering (EnKF) in the Weather Research and Forecasting (WRF) model can be used to incorporate wind, and later rain rate, data into a mesoscale model forecast of hurricane intensity. The study makes use of an Observing System Simulation Experiment (OSSE) with a simulated HIRAD dataset sampled during a hurricane and uses EnKF to forecast the track and intensity prediction of the hurricane. Comparisons to truth and error metrics are used to assess the model?s forecast performance.
Free, equatorially trapped sinusoidal wave solutions to a linear model on an equatorial beta plane are used to fit the Geosat altimetric sea level observations in the tropical Pacific Ocean. The Kalman filter technique is used to estimate the wave amplitude and phase from the data. The estimation is performed at each time step by combining the model forecast with the observation in an optimal fashion utilizing the respective error covariances. The model error covariance is determined such that the performance of the model forecast is optimized. It is found that the dominant observed features can be described qualitatively by basin-scale Kelvin waves and the first meridional-mode Rossby waves. Quantitatively, however, only 23 percent of the signal variance can be accounted for by this simple model.
We describe a variational continuous assimilation (VCA) algorithm for assimilating tropical rainfall data using moisture and temperature tendency corrections as the control variable to offset model deficiencies. For rainfall assimilation, model errors are of special concern since model-predicted precipitation is based on parameterized moist physics, which can have substantial systematic errors. This study examines whether a VCA scheme using the forecast model as a weak constraint offers an effective pathway to precipitation assimilation. The particular scheme we exarnine employs a '1+1' dimension precipitation observation operator based on a 6-h integration of a column model of moist physics from the Goddard Earth Observing System (GEOS) global data assimilation system DAS). In earlier studies, we tested a simplified version of this scheme and obtained improved monthly-mean analyses and better short-range forecast skills. This paper describes the full implementation ofthe 1+1D VCA scheme using background and observation error statistics, and examines how it may improve GEOS analyses and forecasts of prominent tropical weather systems such as hurricanes. Parallel assimilation experiments with and without rainfall data for Hurricanes Bonnie and Floyd show that assimilating 6-h TMI and SSM/I surfice rain rates leads to more realistic storm features in the analysis, which, in turn, provide better initial conditions for 5-day storm track prediction and precipitation forecast. These results provide evidence that addressing model deficiencies in moisture tendency may be crucial to making effective use of precipitation information in data assimilation.