Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Weather 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 613 records · Page 34

Tropical impacts of SST forcing: A case study for 1987 versus 1988

The response of the NASA/Goddard Institute for Space Studies general circulation model (GCM) to large tropical sea surface temperature (SST) anomalies is investigated by evaluating model simulations of the particularly contrasting summer monsoon seasons 1987 and 1988. These years are representative of the warm and cold phases, respectively, of a recent El Nino-Southern Oscillation (ENSO) event. An ensemble averaging the results of three simulations was considered for each season, using monthly mean observed SST anomalies for June-August 1987 and 1988 as lower boundary forcing. Consistent with the European Center for Medium Weather Forecast (ECMWF)-analyzed winds, the simulations based on 1988 as compared to 1987 SST exhibit stronger upper-tropospheric irrotational circulation between the monsoon regions and the Southern Hemispheric sub-tropical anticyclones, a stronger Pacific Walker cell and a weaker subtropical westerly jet over the South Pacific. In the same vein, the modeled precipitation, indicating a more northerly position of the Pacific Inter-Tropical Convergence Zone (ITCZ) in 1988 compared with 1987, is supported by satellite observations of outgoing longwave radiation and highly reflective clouds.

Druyan, Leonard M.↗

Intraseasonal behavior of clouds, temperature, and motion in the tropics

The spectral character of tropical convection is investigated in an 11-yr record of outgoing longwave radiation from the Advanced Very High Resolution Radiometer (AVHRR) to identify interaction with the tropical circulation. Along the equator in the eastern hemisphere, the space-time spectrum of convection possesses a broad peak at wavenumbers 1-3 and eastward periods of 35-95 days. Significantly broader than the dynamical signal of the Madden-Julian oscillation (MJO), this quasi-discrete convective signal is associatd with a large-scale anomaly that propagates across and modulates time mean or 'climatological convection' over the equatorial Indian Ocean and western Pacific. Outside that region the convective signal is small, even though, under amplified conditions, coherence can be found east of the date line and in the subtropics. Having a zonal scale of approximately wavenumber 2, anomalous convection propagates eastward at some 5 m/s and suppresses as well as reinforces climatological in the eastern hemisphere. The convective signal amplifies to a seasonal maximum near vernal equinox and, to a weaker degree, again near autumnal equinox, when climatological convection and warm sea surface temperature (SST) cross the equator. Contemporaneous records of motion from European Center for medium-Range Weather Forecasts (ECMWF) analyses and tropospheric-mean temperature from Microwave Sounding Unit reveal an anomalous component of the tropical circulation that coexists with the convective signal and embodies many of the established properties of the MJO. In the eastern hemisphere, subtropical Rossby gyres and zonal Kelvin structure along the equator flank the convective anomaly as it tracks eastward, giving the anomalous circulation to form of a 'forced response.' In the western hemisphere, the dynamical signal is composed chiefly of wavenumber-1 Kelvin structure, which as the form of a 'propagating response' that is excited in and radiates away from anomalous convection at some 10 m/s. Kelvin structure comprising the propagating response appears in 850-mb and 200-mb zonal winds even when the convective signal is absent, albeit with much smaller amplitude. In contrast, the signal in 1000-mb convergence appears only when accompanied by anomalous convection, which suggests that convergence in the boundary layer is instrumental in achieving strong interaction with the convective pattern.

Salby, Murry L.↗

Estimating central pressures of oceanic midlatitude cyclones

A method of determining surface pressures in oceanic storm systems using Earth Remote-Sensing Satellite (ERS-1) scatterometer data is employed to determine the lowest pressure in 25 storms. This method uses the surface winds as a lower boundary condition on a planetary boundary layer model to determine gradient winds and, thereby, pressure gradients. An optimization scheme referenced to a pressure outside the storm provides a pressure field and an estimate of the low pressure. The values are compared to European Center for Medium-Range Weather Forecasts (ECMWF) analyses in each case; there is good agreement, with some expected differences.

Brown, R. A.↗

January and July global distributions of atmospheric heating for 1986, 1987, and 1988

Three-dimensional global distributions of atmospheric heating are estimated for January and July of the 3-year period 1986-88 from the European Center for Medium Weather Forecasts (ECMWF) Tropical Ocean Global Atmosphere (TOGA) assimilated datasets. Emphasis is placed on the interseasonal and interannual variability of heating both locally and regionally. Large fluctuations in the magnitude of heating and the disposition of maxima/minima in the Tropics occur over the 3-year period. This variability, which is largely in accord with anomalous precipitation expected during the El Nino-Southern Oscillation (ENSO) cycle, appears realistic. In both January and July, interannual differences of 1.0-1.5 K/day in the vertically averaged heating occur over the tropical Pacific. These interannual regional differences are substantial in comparison with maximum monthly averaged heating rates of 2.0-2.5 K/day. In the extratropics, the most prominent interannual variability occurs along the wintertime North Atlantic cyclone track. Vertical profiles of heating from selected regions also reveal large interannual variability. Clearly evident is the modulation of the heating within tropical regions of deep moist convection associated with the evolution of the ENSO cycle. The heating integrated over continental and oceanic basins emphasizes the impact of land and ocean surfaces on atmospheric energy balance and depicts marked interseasonal and interannual large-scale variability.

Schaack, Todd K.↗

Precipitation recycling in the Amazon basin

Precipitation recycling is the contribution of evaporation within a region to precipitation in that same region. The recycling rate is a diagnostic measure of the potential for interactions between land surface hydrology and regional climate. In this paper we present a model for describing the seasonal and spatial variability of the recycling process. The precipitation recycling ratio, rho, is the basic variable in describing the recycling process. Rho is the fraction of precipitation at a certain location and time which is contributed by evaporation within the region under study. The recycling model is applied in studyiing the hydrologic cycle in the Amazon basin. It is estimated that about 25% of all the rain that falls in the Amazon basin is contributed by evaporation within the basin. This estimate is based on analysis of a data set supplied by the European Centre for Medium-range Weather Forecasts (ECMWF). The same analysis is repeated using a different data set from the Geophysical Fluid Dynamics Laboratory (GFDL). Based on this data set, the recycling ratio is estimated to be 35%. The seasonal variability of the recycling ratio is small compared with the yearly average. The new estimates of the recycling ratio are compared with results of previous studies, and the differences are explained.

Eltahir, E. A. B.↗

Temperature-moisture biases in ECMWF analyses based on clear sky longwave simulations constrained by SSMI and MSU measurements and comparisons to ERBE estimates

Clear sky longwave radiation fluxes for the summer of 1988 and winter of 1989 have been simulated with a radiative transfer model that includes detailed treatment of atmospheric gas absorption. The input data to the model are humidity and temperature profiles from European Center for Medium-Range Weather Forecasts (ECMWF) analyses and surface temperature measurements from the International Satellite Cloud Climatology Project. To reduce the inherent uncertainties in humidity profiles, the ECMWF moisture fields have been adjusted, based on the climatological relationship between the moisture profile and total precipitable water (PW), so that the ECMWF PW is equal to that derived from Special Sensor Microwave Imager (SSMI) data. SSMI and ECMWF PW patterns are generally similar, but significant differences in magnitude are found in the low latitudes, particularly over the subtropical oceans in the Southern Hemisphere. At the same time the tropospheric brightness temperature (T(sub B)) at 53.74 GHz has been computed using a microwave transfer model with ECMWF temperature and humidity fields as inputs. The comparison of computed T(sub B) with the Microwave Sounding Unit (MSU) channel 2 T(sub B) revealed differences larger than 1.5 K over most of the oceans, suggesting that the ECMWF model atmosphere has a cold bias. Maximum biases greater than 3 K are found in low latitudes. After removing the bias from the ECMWF temperature field, a new set of temperature profiles yielding brightness temperature equal to MSU channel 2 T(sub B) has been obtained. Simulation results clearly demonstrate that inclusion of satellite estimates of PW and T(sub B) enhances the accuracy of the clear sky LW flux simulation, substantially reducing the differences from Earth Radiation Budget Experiment (ERBE) estimates. Global averages indicate that the model-derived top-of-atmosphere (TOA) clear sky fluxes from the adjusted moisture and temperature profiles are in very good agreement with satellite-measured ERBE values, differing by only 1.4 W/sq m. By contrast, the calculation using only ECMWF analyses without any adjustments gives a significant underestimation of TOA clear sky fluxes up to 10.7 W/sq m compared to the satellite measurements, suggesting that the constraint method can significantly improve the accuracy in the radiation budget simulation. The results also indicate that the ECMWF-based simulation errors are, at least for the two seasons studied, mainly attributed to the error in the temperature field rather than to errors in the moisture field.

Sohn, Byung-Ju↗

Thermal remote sensing of surface soil water content with partial vegetation cover for incorporation into climate models

This study outlines a method for the estimation of regional patterns of surface moisture availability (M(sub 0)) and fractional vegetation (Fr) in the presence of spatially variable vegetation cover. The method requires relating variations in satellite-derived (NOAA, Advanced Very High Resolution Radiometer (AVHRR)) surface radiant temperature to a vegetation index (computed from satellite visible and near-infrared data) while coupling this association to an inverse modeling scheme. More than merely furnishing surface soil moisture values, the method constitues a new conceptual and practical approach for combining thermal infrared and vegetation index measurements for incorporating the derived values of M(sub 0) into hydrologic and atmospheric prediction models. Application of the technique is demonstrated for a region in and around the city of Newcastle upon Tyne situated in the northeast of England. A regional estimate of M(sub 0) is derived and is probabbly good for fractional vegetation cover up to 80% before errors in the estimated soil water content become unacceptably large. Moreover, a normalization scheme is suggested from which a nomogram, `universal triangle,' is constructed and is seen to fit the observed data well. The universal triangle also simplifies the inclusion of remotely derived M(sub 0) in hydrology and meteorological models and is perhaps a practicable step toward integrating derived data from satellite measurements in weather forecasting.

Gillies, Robert R.↗

Rossby-gravity waves in tropical total ozone

Randel (1992) has recently reported evidence for Rossby-gravity waves in tropical data fields produced by the European Centre for Medium Range Weather Forecasts (ECMWF). The purpose of this paper is to show that similar features are observable in fields of total column ozone from the Total Ozone Mapping Spectrometer (TOMS) satellite instrument. The observed features are episodic, have zonal (east-west) wavelengths of approximately 6,000-10,000 km and oscillate with periods of 5-10 days. The modes exhibit westward phase progression and eastward group velocity. A simple linear model is used to estimate the magnitude of total ozone perturbations induced by Rossby-gravity waves. The model is able to account for the magnitude of observed TOMS features as well as their asymmetry about the equator. The significance of finding Rossby-gravity waves in total ozone fields is that: (1) the report of similar features in ECMWF tropical fields is corroborated with an independent data set; and (2) the TOMS data set is demonstrated to possess surprising versatility and sensitivity to relatively smaller scale tropical phenomena.

Stanford, J. L.↗

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.↗

Estimating the urban bias of surface shelter temperatures using upper-air and satellite data. Part 2: Estimation of the urban bias

A methodology is presented for estimating the urban bias of surface shelter temperatures due to the effect of the urban heat island. 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, site-specific data to represent the local landscape, and satellite-derived data -- the normalized difference vegetation index (NDVI) and the Defense Meteorological Satellite Program (DMSP) nighttime brightness data -- to represent the urban and rural landscape. Local NDVI and DMSP values were calculated for each station using the mean NDVI and DMSP values from a 3 km x 3 km area centered over the given station. Regional NDVI and DMSP values were calculated to represent a typical rural value for each station using the mean NDVI and DMSP values from a 1 deg x 1 deg latitude-longitude area in which the given station was located. Models for the United States were then developed for monthly maximum, mean, and minimum temperatures using data from over 1000 stations in the U.S. Cooperative (COOP) Network and for monthly mean temperatures with data from over 1150 stations in the Global Historical Climate Network (GHCN). Local biases, or the differences between the model predictions using the observed NDVI and DMSP values, and the predictions using the background regional values were calculated and compared with the results of other research. The local or urban bias of U.S. temperatures, as derived from all U.S. stations (urban and rural) used in the models, averaged near 0.40 C for monthly minimum temperatures, near 0.25 C for monthly mean temperatures, and near 0.10 C for monthly maximum temperatures. The biases of monthly minimum temperatures for individual stations ranged from near -1.1 C for rural stations to 2.4 C for stations from the largest urban areas. The results of this study indicate minimal problems for global application once global NDVI and DMSP data become available.

Epperson, David L.↗

Variability of moisture in the upper troposphere as inferred from TOVS satellite observations and the ECMWF model analyses in 1989

Large-scale variability of moisture in the upper troposphere is examined using TIROS (Television Infrared Observation Satellite) Operational Vertical Sounder (TOVS) satellite observations and European Centre for Medium-Range Weather Forecasts (ECMWF) model analyses from 1989 in the latitude band from 40 deg N to 40 S. To compare these dissimilar datasets, upwelling radiances were computed for the 6-7 micrometer water vapor band from the ECMWF temperature and moisture analyses, and these computed radiances were compared to the corresponding TOVS satellite observations. The ECMWF-based radiances reproduce the general locations and seasonal cycle of the TOVS-observed moisture features, particularly after an improved convective parameterization scheme was adopted by ECMWF in May 1989. However, the ECMWF analysis scheme still results in much milder lateral moisture gradients and seasonal contrasts than indicated by the TOVS observations. Seasonally, the upper troposphere in each hemisphere dries in winter and moistens in summer, but there are regions in each hemisphere that run counter to this seasonal trend, apparently depending on continental- and monsoon-scale dynamics. Dynamically, the TOVS-observed regions of significant subtropical dryness are correlated with persistent subsidence indicated by ECMWF 300-mb vertical velocity analyses. The TOVS radiance observations indicate large variations in space and time of the upper-tropospheric moisture field, which are not fully captured by the ECMWF analyses.

Salathe, Eric P., Jr.↗

A regional simulation of marine boundary-layer clouds

A regional version of a two-layer model is used to simulate marine boundary-layer clouds over an area (20 deg - 42 deg and 117 - 145 deg W) of the eastern North Pacific. The large-scale conditions requried for the model are provided by European Centre for Medium-Range Weather Forecasts analyses and the Comprehensive Ocean-Atmosphere Data Set. The simulated cloud-top height and fractional cloudiness are compared with satellite data. The model predicts a realistic pattern of cloud-top height, although the heights appear to be overestimated in the southwest and underestimated along the eastern border of the domain. The overall simulated cloudiness is qualitatively comparable to the satellite derived, but the model predicts excessive cloudiness over cold ocean surfaces and insufficient cloudiness over an area centered at 30 deg N, 130 deg W. The possible causes for these discrepancies are discussed. Budget studies show that horizontal advection contributes substantially to the cloud-top height and cloud-layer moisture budgets. Thus, one may need to consider the effects of horizontal advection when validating steady-state solutions of a one-dimensional model. The simulated cloudiness is shown to be sensitive to large-scale subsidence, moisture above the clouds, and drizzle. An increase in the large-scale divergence by 2 x 10(exp -6)/sec decreases average cloudiness for the area by 30%. An increase in the moisture at 850 mb by 2 g/kg increases the cloudiness by 20%. The suppression of drizzle in the model increases the cloudiness by 34% and significantly changes its pattern.

Wang, Shouping↗

Space and time aliasing structure is monthly mean polar-orbiting satellite data

Monthly mean wind fields from the European Remote Sensing Satellite (ERS1) scatterometer are presented. A banded structure which resembles the satellite subtrack is clearly and consistently apparent in the isotachs as well as the u and v components of the routinely produced fields. The structure also appears in the means of data from other polar-orbiting satellites and instruments. An experiment is designed to trace the cause of the banded structure. The European Centre for Medium-Range Weather Forecast (ECMWF) gridded surface wind analyses are used as a control set. These analyses are also sampled with the ERS1 temporal-spatial samplig pattern to form a simulated scatterometer wind set. Both sets are used to create monthly averages. The banded structures appear in the monthly mean simulated data but do not appear in the control set. It is concluded that the source of the banded structure lies in the spatial and temporal sampling of the polar-orbiting satellite which results in undersampling. The problem involves multiple timescales and space scales, oversampling and under-sampling in space, aliasing in the time and space domains, and preferentially sampled variability. It is shown that commonly used spatial smoothers (or filters), while producing visually pleasing results, also significantly bias the true mean. A three-dimensional spatial-temporal interpolator is designed and used to determine the mean field. It is found to produce satisfactory monthly means from both simulated and real ERS1 data. The implications to climate studies involving polar-orbiting satellite data are discussed.

Zeng, Lixin↗

TOGA COARE Satellite data summaries available on the World Wide Web

Satellite data summary images and analysis plots from the Tropical Ocean Global Atmosphere Coupled Ocean-Atmosphere Response Experiment (TOGA COARE), which were initially prepared in the field at the Honiara Operations Center, are now available on the Internet via World Wide Web browsers such as Mosaic. These satellite data summaries consist of products derived from the Japanese Geosynchronous Meteorological Satellite IR data: a time-size series of the distribution of contiguous cold cloudiness areas, weekly percent high cloudiness (PHC) maps, and a five-month time-longitudinal diagram illustrating the zonal motion of large areas of cold cloudiness. The weekly PHC maps are overlaid with weekly mean 850-hPa wind calculated from the European Centre for Medium-Range Weather Forecasts (ECMWF) global analysis field and can be viewed as an animation loop. These satellite summaries provide an overview of spatial and temporal variabilities of the cloud population and a large-scale context for studies concerning specific processes of various components of TOGA COARE.

Chen, S. S.↗

A mesoscale modeling study of the atmospheric circulation of high southern latitudes

The meteorology of high southern latitudes during winter is simulated using a cloud-free version of The Pennsylvania State University-National Center for Atmospheric Research Mesoscale Model version 4 (MM4) with a 100-km horizontal resolution. Comparisons between idealized simulations of Antarctic with MM4 and with the mesoscale model of Parish and Waight reveal that both models produce similarly realistic velocity fields in the boundary layer. The latter model tends to produce slightly faster drainage winds over East Antarctica. The intensity of the katabatic winds produced by MM4 is sensitive to parameterizations of boundary layer fluxes. Two simulations performed with MM4 using analyses from the European Center for Medium-Range Weather Forecasts (ECMWF) for June 1988 as initial and boundary conditions. A simulation of the period from 000 UTC 2 June to 0000 UTC 8 June produces realistic synoptic phenomena including ridge development over East Antarctica, frontogenesis over the Amundsen Sea, and a katabatic surge over the Ross Ice Shelf. The simulated time-averaged fields for June 1988, particularly that of a 500-hPa height, are in good agreement with time-averaged fields analyzed by the ECMWF. The results of the simulations provide detailed features of the Antarctic winter boundary layer along the steeply sloping terrain. Highest boundary layer wind speeds averaged over the month-long simulation are approximately 20 m/s. The lack of latent heating in the simulations apparently results in some bias in the results. In particular, the cloud-free version of MM4 underpredicts the intensity of lows in the sea level pressure field.

Hines, K. M.↗

The use of stellite scatterometer winds to drive a primitive equation model of the Indian Ocean: The impact of bandlike sampling

The aim of this study is to evaluate the impact of the bandlike sampling of spaceborne scatterometers on the ability of scatterometer winds to successfully force the mean flow and seasonal cycle of an ocean model in the context of equatorial and tropical dynamics. The equatorial ocean is simulated with a four-layer, primitive equation, reduced gravity model of the Indian Ocean. The variable wind stress used in this study is derived from one year (1988) of 6-hour analyses of the 10-m wind vector over the Indian Ocean performed at the European Centre for Medium-Range Weather Forecasts (ECMWF). It is applied as a forcing at every grid point of the model to drive a reference circulation. Scatterometer winds are simulated from ECMWF winds, using the nominal configurations and orbital parameters of the European Remote Sensing 1 (ERS-1) and NASA Scatterometer (NSCAT) missions. The model is forced in real time under swaths with the raw scatterometer winds of ERS-1 and NSCAT, with a persistence condition (i.e., the wind is kept constsnt until the next passage of the satellite provides a new value). The circulation obtained for each of the scatterometer experiments is compared with the reference circulation. The seasonal circulation of the Indian Ocean with NSCAT winds is very similar to the reference. The perturbations introduced by the bandlike sampling and the persistance condition have an impact similar to that of a small uncorrelated noise added to the reference forcing. The persistence condition for ERS-1 does not give results which are as good as those obtained for NSCAT.

Barnier, Bernard↗

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.↗

Technical Report Series on Global Modeling and Data Assimilation. Volume 7: Proceedings of the Workshop on the GEOS-1 Five-year Assimilation

The primary objective of the three-day workshop on results from the Data Assimilation Office (DAO) five-year assimilation was to provide timely feedback from the data users concerning the strengths and weaknesses of version 1 of the Goddard Earth Observing System (GEOS-1) assimilated products. A second objective was to assess user satisfaction with the current methods of data access and retrieval. There were a total of 49 presentations, with about half (23) of the presentations from scientists from outside of Goddard. The first two days were devoted to applications of data: studies of the energy diagnostics, precipitation and diabatic heating, hydrological modeling and moisture transport, cloud forcing and validation, various aspects of intraseasonal, seasonal, and interannual variability, ocean wind stress applications, and validation of surface fluxes. The last day included talks from the National Meteorological Center (NMC), the National Center for Atmospheric Research (NCAR), the Center for Ocean-Land-Atmosphere Studies (COLA), the United States Navy, and the European Center for Medium Range Weather Forecasts (ECMWF).

Suarez, Max J.↗