Bibliography - William F Stern
- Delworth, Thomas L., Anthony J Broccoli, Anthony Rosati, Ronald J Stouffer, Ventakramani Balaji, J A Beesley, W F Cooke, Keith W Dixon, John Dunne, Krista A Dunne, J W Durachta, Kirsten L Findell, Paul Ginoux, Anand Gnanadesikan, C Tony Gordon, Stephen Griffies, Rich Gudgel, Matthew J Harrison, Isaac Held, Richard S Hemler, Larry Horowitz, Stephen A Klein, Thomas R Knutson, P J Kushner, A R Langenhorst, H C Lee, Shian-Jiann Lin, Jian Lu, S Malyshev, P C D Milly, V Ramaswamy, J L Russell, M Daniel Schwarzkopf, Elena Shevliakova, Joseph J Sirutis, Michael J Spelman, William F Stern, Michael Winton, Andrew T Wittenberg, Bruce Wyman, Fanrong Zeng, and Rong Zhang, 2006: GFDL's CM2 Global Coupled Climate Models. Part I: Formulation and Simulation Characteristics. Journal of Climate, 19(5), doi:10.1175/JCLI3629.1.
[ Abstract ]The formulation and simulation characteristics of two new global coupled climate models developed at NOAA's Geophysical Fluid Dynamics Laboratory (GFDL) are described. The models were designed to simulate atmospheric and oceanic climate and variability from the diurnal time scale through multicentury climate change, given our computational constraints. In particular, an important goal was to use the same model for both experimental seasonal to interannual forecasting and the study of multicentury global climate change, and this goal has been achieved.
Two versions of the coupled model are described, called CM2.0 and CM2.1. The versions differ primarily in the dynamical core used in the atmospheric component, along with the cloud tuning and some details of the land and ocean components. For both coupled models, the resolution of the land and atmospheric components is 2° latitude × 2.5° longitude; the atmospheric model has 24 vertical levels. The ocean resolution is 1° in latitude and longitude, with meridional resolution equatorward of 30° becoming progressively finer, such that the meridional resolution is 1/3° at the equator. There are 50 vertical levels in the ocean, with 22 evenly spaced levels within the top 220 m. The ocean component has poles over North America and Eurasia to avoid polar filtering. Neither coupled model employs flux adjustments.
The control simulations have stable, realistic climates when integrated over multiple centuries. Both models have simulations of ENSO that are substantially improved relative to previous GFDL coupled models. The CM2.0 model has been further evaluated as an ENSO forecast model and has good skill (CM2.1 has not been evaluated as an ENSO forecast model). Generally reduced temperature and salinity biases exist in CM2.1 relative to CM2.0. These reductions are associated with 1) improved simulations of surface wind stress in CM2.1 and associated changes in oceanic gyre circulations; 2) changes in cloud tuning and the land model, both of which act to increase the net surface shortwave radiation in CM2.1, thereby reducing an overall cold bias present in CM2.0; and 3) a reduction of ocean lateral viscosity in the extratropics in CM2.1, which reduces sea ice biases in the North Atlantic.
Both models have been used to conduct a suite of climate change simulations for the 2007 Intergovernmental Panel on Climate Change (IPCC) assessment report and are able to simulate the main features of the observed warming of the twentieth century. The climate sensitivities of the CM2.0 and CM2.1 models are 2.9 and 3.4 K, respectively. These sensitivities are defined by coupling the atmospheric components of CM2.0 and CM2.1 to a slab ocean model and allowing the model to come into equilibrium with a doubling of atmospheric CO2. The output from a suite of integrations conducted with these models is freely available online (see http://nomads.gfdl.noaa.gov/).
Manuscript received 8 December 2004, in final form 18 March 2005
- Jones, C, D E Waliser, K M Lau, and William F Stern, 2004: Global occurrences of extreme precipitation and the Madden-Julian Oscillation: Observations and predictability. Journal of Climate, 17(23), 4575-4589.
[ Abstract PDF ]This study investigates 1) the eastward propagation of the Madden-Julian oscillation (MJO) and global occurrences of extreme precipitation, 2) the degree to which a general circulation model with a relatively realistic representation of the MJO simulates its influence on extremes, and 3) a possible modulation of the MJO on potential predictability of extreme precipitation events. The observational analysis shows increased frequency of extremes during active MJO phases in many locations. On a global scale, extreme events during active MJO periods are about 40% higher than in quiescent phases of the oscillation in locations of statistically significant signals.
A 10-yr National Aeronautics and Space Administration (NASA) Goddard Laboratory for the Atmospheres (GLA) GCM simulation with fixed climatological SSTs is used to generate a control run and predictability experiments. Overall, the GLA model has a realistic representation of extremes in tropical convective regions associated with the MJO, although some shortcomings also seem to be present. The GLA model shows a robust signal in the frequency of extremes in the North Pacific and on the west coast of North America, which somewhat agrees with observational studies. The analysis of predictability experiments indicates higher success in the prediction of extremes during an active MJO than in quiescent situations. Overall, the predictability experiments indicate the mean number of correct forecasts of extremes during active MJO periods to be nearly twice the correct number of extremes during quiescent phases of the oscillation in locations of statistically significant signals.
- Jones, C, D E Waliser, K M Lau, and William F Stern, 2004: The Madden-Julian Oscillation and its impact on Northern Hemisphere weather predictability. Monthly Weather Review, 132(6), 1462-1471.
[ Abstract PDF ]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 uses numerical model experiments to investigate the influence of the MJO activity on weather predictability in the midlatitudes of the Northern Hemisphere. The National Aeronautics and Space Administration (NASA) Goddard Laboratory for the Atmospheres (GLA) general circulation model was used in a 10-yr simulation with fixed climatological SSTs to generate a control dataset 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 and standardized root-mean-square errors in the midlatitudes of the Northern Hemisphere (20°–60°N) were computed to assess predictability characteristics. The analyses of 500-hPa geopotential height, 200-hPa streamfunction, and 850-hPa zonal wind component systematically show larger predictability (~2–3 days) during periods of active MJO as opposed to quiescent episodes of the oscillation. While further studies are necessary to investigate possible model sensitivity, the results shown here highlight the importance of the MJO in modulating weather variability and show the importance of improving the representation of the MJO in operational numerical weather forecast models.
- Zheng, Y, D E Waliser, William F Stern, and C Jones, 2004: The role of coupled sea surface temperatures in the simulation of the tropical intraseasonal oscillation. Journal of Climate, 17(21), 4109-4134.
[ Abstract PDF ]This study compares the tropical intraseasonal oscillation (TISO) variability in the Geophysical Fluid Dynamics Laboratory (GFDL) coupled general circulation model (CGCM) and the stand-alone atmospheric general circulation model (AGCM). For the AGCM simulation, the sea surface temperatures (SSTs) were specified using those from the CGCM simulation. This was done so that any differences in the TISO that emerged from the two simulations could be attributed to the coupling process and not to a difference in the mean background state. The comparison focused on analysis of the rainfall, 200-mb velocity potential, and 850-mb zonal wind data from the two simulations, for both summer and winter periods, and included comparisons to analogous diagnostics using NCEP–NCAR reanalysis and Climate Prediction Center (CPC) Merged Analysis of Precipitation (CMAP) rainfall data.
The results of the analysis showed three principal differences in the TISO variability between the coupled and uncoupled simulations. The first was that the CGCM showed an improvement in the spatial variability associated with the TISO mode, particularly for boreal summer. Specifically, the AGCM exhibited almost no TISO variability in the Indian Ocean during boreal summer—a common shortcoming among AGCMs. The CGCM, on the other hand, did show a considerable enhancement in TISO variability in this region for this season. The second was that the wavenumber -frequency spectra of the AGCM exhibited an unrealistic peak in variability at low wavenumbers (1-3, depending on the variable) and abouit 3 cycles yr-1 (cpy). This unrealistic peak of variability was absent in the CGCM, which otherwise tended to show good agreement with the observations. The third difference was that the AGCM showed a less realistic phase lag between the TISO-related convection and SST anomalies. In particular, the CGCM exhibited a near-quadrature relation between precipitation and SST anomalies, which is consistent with observations, while the phase lag was reduced in the AGCM by about 1.5 pentads (~ 1 week). The implications of the above results, including those for the notions of "perfect SST" and "two tier" experiments, are discussed, as are the caveats associated with the study's modeling framework and analysis.
- Waliser, D E., K Jin, I-S Kang, William F Stern, S D Schubert, M L C Wu, K M Lau, M-I Lee, V Krishnamurthy, A Kitoh, G A Meehl, V Y Galin, V Satyan, S K Mandke, G Wu, Y Liu, and C-K Park, 2003: AGCM simulations of intraseasonal variability associated with the Asian summer monsoon. Climate Dynamics, 21(5-6), 423-446.
[ Abstract PDF ]The intraseasonal variability associated with the Asian summer monsoon as simulated by a number of atmospheric general circulation models (AGCMs) are analyzed and assessed against observations. The model data comes from the Monsoon GCM Intercomparison project initiated by the CLIVAR/Asian-Australian Monsoon Panel. Ten GCM groups, i.e., the Center for Ocean-Land-Atmosphere Studies (COLA), Institute of Numerical Mathematics (DNM), Goddard Space Flight Center (GSFC), Geophysical Fluid Dynamics Laboratory (GFDL), Institute of Atmospheric Physics (IAP), Indian Institute of Tropical Meteorology (IITM), Meteorological Research Institute (MRI), National Center for Atmospheric Research (NCAR), Seoul National University (SNU), and the State University of New York (SUNY), participated in the intraseasonal component of the project. Each performed a set of 10 ensemble simulations for 1 September 1996-31 August 1998 using the same observed weekly SST values but with different initial conditions. The focus is on the spatial and seasonal variations associated with intraseasonal variability (ISV) of rainfall, the structure of each model''s principal mode of spatial-temporal variation of rainfall [i.e. their depiction of the Intraseasonal Oscillation (ISO)], the teleconnection patterns associated with each model''s ISO, and the implications of the models'' ISV on seasonal monsoon predictability. The results show that several of the models exhibit ISV levels at or above that found in observations with spatial patterns of ISV that resemble the observed pattern. This includes a number of rather detailed features, including the relative distribution of variability between ocean and land regions. In terms of the area-averaged variance, it is found that the fidelity of a model to represent NH summer versus winter ISV appears to be strongly linked. In addition, most models'' ISO patterns do exhibit some form of northeastward propagation. However, the model ISO patterns are typically less coherent, lack sufficient eastward propagation, and have smaller zonal and meridional spatial scales than the observed patterns, and are often limited to one side or the other of the maritime continent. The most pervasive and problematic feature of the models'' depiction of ISV and/or their ISO patterns is the overall lack of variability in the equatorial Indian Ocean. In some cases, this characteristic appears to result from some models forming double convergence zones about the equator rather than one region of strong convergence on the equator. This shortcoming results in a poor representation of the local rainfall pattern and also significantly influences the models'' representations of the global-scale teleconnection patterns associated with the ISO. Finally, analysis of the model ensemble shows a positive relationship between the strength of a model''s ISV of rainfall and its intra-ensemble variability of seasonal monsoon rainfall. The implications of this latter relation are discussed in the context of seasonal monsoon predictability.
- Waliser, D E., K M Lau, William F Stern, and C Jones, 2003: Potential predictability and the Madden-Julian Oscillation. Bulletin of the American Meteorological Society, 84(1), 33-50.
[ Abstract PDF ]Ensembles of twin-predictability experiments suggest that successful dynamical forecasts of the Madden-Julian oscillation may offer one avenue for bridging the gap between medium to long-range weather forecasting and short-term climate prediction.
- Waliser, D E., William F Stern, S D Schubert, and K M Lau, 2003: Dynamic predictability of intraseasonal variability associated with the Asian summer monsoon. Quarterly Journal of the Royal Meteorological Society, 129(594), 2897-2925.
[ Abstract PDF ]The objective of this study is to estimate the limit of dynamical predictability for the tropical intraseasonal oscillation (ISO) associated with the Asian summer monsoon. Ensembles of 'twin' predictability experiments were carried out with the National Aeronautics and Space Administration Goddard Laboratory for the Atmospheres atmospheric general circulation model using specified annual cycle sea surface temperatures. Initial conditions were taken from a 10-year control simulation during periods of strong ISO activity identified via extended empirical orthogonal function analysis of 30-90 day band-passed tropical rainfall. From this analysis, 21 cases were chosen for each of four distinct phases of the ISO cycle, making 84 cases in total. Two different sets of small random perturbations were added to these 84 initial states. Simulations were then performed for 90 days from each of these 168 perturbed initial conditions. A measure of potential predictability was constructed based on a ratio of the signal associated with the ISO, in terms of rainfall or 200 hPa velocity potential (VP200), and the mean-square difference between sets of twin forecasts. This ratio indicates that the limit of predictability for this model's ISO extends out to about 25 days for VP200 and to about 15 days for rainfall. The predictability measure shows modest dependence on the strength of the ISO events, with stronger events having a greater limit of predictability by a few days. It also exhibits a dependence on the phase of the ISO, with greater predictability for the convective phase at short (~> 5 days) lead times and for the suppressed phase at longer (~< 15 days) lead times. The implications of these results as well as their associated model and analysis caveats are discussed.
- Kang, I-S, K Jin, K M Lau, J Shukla, V Krishnamurthy, S D Schubert, D E Waliser, William F Stern, and V Satyan, et al., 2002: Intercomparison of atmospheric GCM simulated anomalies associated with the 1997/98 El Niño. Journal of Climate, 15(19), 2791-2805.
[ Abstract PDF ]The atmospheric anomalies for the 1997/98 El Niño–Southern Oscillation (ENSO) period have been analyzed and intercompared using the data simulated by the atmospheric general circulation models (GCMs) of 11 groups participating in the Monsoon GCM Intercomparison Project initiated by the Climate Variability and Prediction Program (CLIVAR)/Asian–Australian Monsoon Panel. Each participating GCM group performed a set of 10 ensemble simulations for 1 September 1996–31 August 1998 using the same sea surface temperature (SST) conditions but with different initial conditions. The present study presents an overview of the intercomparison project and the results of an intercomparison of the global atmospheric anomalies during the 1997/98 El Niño period. Particularly, the focus is on the tropical precipitation anomalies over the monsoon–ENSO region and the upper-tropospheric circulation anomalies in the Pacific–North American (PNA) region.
- Kang, I-S, K Jin, B Wang, K M Lau, J Shukla, V Krishnamurthy, S D Schubert, D E Waliser, and William F Stern, et al., 2002: Intercomparison of the climatological variations of Asian summer monsoon precipitation simulated by 10 GCMs. Climate Dynamics, 19(5-6), 383-395.
[ Abstract PDF ]We assess the overall performance of state-of-the-art atmospheric GCMs in simulating the climatological variations of summer monsoon rainfall over the Asian-Western Pacific region and the systematic errors that are common to a group of GCMs. The GCM data utilized are obtained from 10 GCM groups participating in the CLIVAR Monsoon GCM Intercomparison Project. The model composite shows that the overall spatial pattern of summer monsoon rainfall is similar to the observed, although the western Pacific rainfall is relatively weak. For the simulated precipitation over the western Pacific, the models can be classified into two categories. The first category of models simulates the precipitation more confined to the equatorial region and weaker precipitation in the subtropical western Pacific compared to the observed. The second category of models simulates large precipitation in the subtropical western Pacific but the region is shifted to the north by 5-10°. None of the models realistically reproduce the observed Mei-yu rain band in the region from the East China Sea to the mid Pacific. Most of the models produce a rain band along the continental side of East Asia. The climatological variations of simulated summer rainfall are examined in terms of their amplitude and their principal EOF modes. All models simulate larger amplitudes of the climatological seasonal variation of Indian summer monsoon than the observed, though most models simulate smaller amplitudes in the western Pacific. The ten model composite produces four leading EOF modes over the Asian-western Pacific region, which are remarkably similar to the observed counterparts. The first and second eigenmodes, respectively, represent the smoothed seasonal march of broad-scale monsoon and the onsets of the Indian and East Asian summer monsoon. The third and fourth modes relate to the climatological intraseasonal oscillation (CISO). In contrast to the model composite, several models fail to reproduce the first principal mode, and most models do not reproduce the observed modes higher than the second. The CISO of precipitation is also examined over the Indian monsoon and the East Asia-western Pacific monsoon regions separately.
- Anderson, Jeffrey L., H van den Dool, A Barnston, W Chen, William F Stern, and Jeff J Ploshay, 1999: Present-day capabilities of numerical and statistical models for atmospheric extratropical seasonal simulation and prediction. Bulletin of the American Meteorological Society, 80(7), 1349-1361.
[ Abstract PDF ]A statistical model and extended ensemble integrations of two atmospheric general circulation models (GCMs) are used to simulate the extratropical atmospheric response to forcing by observed SSTs for the years 1980 through 1988. The simulations are compared to observations using the anomaly correlation and root-mean-square error of the 700-hPa height field over a region encompassing the extratropical North Pacific Ocean and most of North America. On average, the statistical model is found to produce considerably better simulations than either numerical model, even when simple statistical corrections are used to remove systematic errors from the numerical model simulations. In the mean, the simulation skill is low, but there are some individual seasons for which all three models produce simulations with good skill.
An approximate upper bound to the simulation skill that could be expected from a GCM ensemble, if the model's response to SST forcing is assumed to be perfect, is computed. This perfect model predictability allows one to make some rough extrapolations about the skill that could be expected if one could greatly improve the mean response of the GCMs without significantly impacting the variance of the ensemble. These perfect model predictability skills are better than the statistical model simulations during the summer, but for the winter, present-day statistical forecasts already have skill that is as high as the upper bound for the GCMs. Simultaneous improvements to the GCM mean response and reduction in the GCM ensemble variance would be required for these GCMs to do significantly better than the statistical model in winter. This does not preclude the possibility that, as is presently the case, a statistical blend of GCM and statistical predictions could produce a simulation better than either alone.
Because of the primitive state of coupled ocean-atmosphere GCMs, the vast majority of seasonal predictions currently produced by GCMs are performed using a two-tiered approach in which SSTs are first predicted and then used to force an atmospheric model; this motivates the examination of the simulation problem. However, it is straightforward to use the statistical model to produce true forecasts by changing its predictors from simultaneous to precursor SSTs. An examination of the decrease in skill of the statistical model when changed from simulation to prediction mode is extrapolated to draw conclusions about the skill to be expected from good coupled GCM predictions.
- Vitart, F, Jeffrey L Anderson, and William F Stern, 1999: Impact of large-scale circulation on tropical storm frequency, intensity, and location, simulated by an ensemble of GCM integrations. Journal of Climate, 12(11), 3237-3254.
[ Abstract PDF ]Tropical storms simulated by a nine-member ensemble of GCM integrations forced by observed SSTs have been tracked by an objective procedure for the period 1980-88. Statistics on tropical storm frequency, intensity, and first location have been produced. Statistical tools such as the chi-square and the Kolmogorov-Smirnov test indicate that there is significant potential predictability of interannual variability of simulated tropical storm frequency, intensity, and first location over most of the ocean basins. The only common point between the nine members of the ensemble is the SST forcing. This implies that SSTs play a fundamental role in model tropical storm frequency, intensity, and first location interannual variability. Although the interannual variability of tropical storm statistics is clearly affected by SST forcing in the GCM, there is also a considerable amount of noise related to internal variability of the model. An ensemble of atmospheric model simulations allows one to filter this noise and gain a better understanding of the mechanisms leading to interannual tropical storm variability.
An EOF analysis of local SSTs over each ocean basin and a combined EOF analysis of vertical wind shear, 850-mb vorticity, and 200-mb vorticity have been performed. Over some ocean basins such as the western North Atlantic, the interannual frequency of simulated tropical storms is highly correlated to the first combined EOF, but it is not significantly correlated to the first EOF of local SSTs. This suggests that over these basins the SSTs have an impact on the simulated tropical storm statistics from a remote area through the large-scale circulation as in observations. Simulated and observed tropical storm statistics have been compared. The interannual variability of simulated tropical storm statistics is consistent with observations over the ocean basins where the model simulates a realistic interannual variability of the large-scale circulation.
- Anderson, Jeffrey L., H van den Dool, A Barnston, W Chen, William F Stern, and Jeff J Ploshay, 1998: Capabilities of dynamical and statistical methods for atmospheric extratropical seasonal prediction In Proceedings of the Twenty-Second Annual Climate Diagnostics and Prediction Workshop, Springfield, VA, NTIS, 46-49.
- D'Andrea, F, S Tibaldi, M Blackburn, G J Boer, M Déqué, M R Dix, B Dugas, L Ferranti, and William F Stern, et al., 1998: Northern Hemisphere atmospheric blocking as simulated by 15 atmospheric general circulation models in the period 1979-1988. Climate Dynamics, 14(6), 385-407.
[ Abstract PDF ]As a part of the Atmospheric Model Intercomparison Project (AMIP), the behaviour of 15 general circulation models has been analyzed in order to diagnose and compare the ability of the different models in simulating Northern Hemisphere midlatitude atmospheric blocking. In accordance with the established AMIP procedure, the 10-year model integrations were performed using prescribed, time-evolving monthly mean observed SSTs spanning the period January 1979 - December 1988. Atmospheric observational data (ECMWF analyses) over the same period have been also used to verify the model results. The models involved in this comparison represent a wide spectrum of model complexity, with different horizontal and vertical resolution, numerical techniques and physical parameterizations, and exhibit large differences in blocking behaviour. Nevertheless, a few common features can be found, such as the general tendency to underestimate both blocking frequency and the average duration of blocks. The problem of the possible relationship between model blocking and model systematic errors has also been assessed, although without resorting to ad-hoc numerical experimentation it is impossible to relate with certainty particular model deficiencies in representing blocking to precise parts of the model formulation.
- Stern, William F., and Anthony Rosati, 1998: Issues of orographic adjustment of SSTs in a spectral/coupled GCM In Research Activities in Atmospheric and Oceanic Modelling, Report No. 27, WMO/TD-No. 865, Geneva, Switzerland, World Meteorological Organization, 6.21-6.22.
- Vitart, F, Jeffrey L Anderson, and William F Stern, 1998: Evaluation of the skill of an ensemble of GCM integrations in simulating seasonal tropical storm frequency, intensity and location In Proceedings of the Twenty-Second Annual Climate Diagnostics and Prediction Workshop, Springfield, VA, NTIS, 38-41.
- Vitart, F, Jeffrey L Anderson, and William F Stern, 1998: Simulation of the internally variability of tropical storm frequency, intensity and location in an ensemble of GCM integrations In Research Activities in Atmospheric and Oceanic Modelling, WMO/TD No. 865, Geneva, Switzerland, World Meteorological Organization, 6.28-6.29.
- Yang, X-Q, Jeffrey L Anderson, and William F Stern, 1998: Reproducible forced modes in AGCM ensemble integrations and potential predictability of atmospheric seasonal variations in the extratropics In Proceedings of the Twenty-Second Annual Climate Diagnostics and Prediction Workshop, Springfield, VA, NTIS, 50-53.
- Yang, X-Q, Jeffrey L Anderson, and William F Stern, 1998: Reproducible forced modes in AGCM ensemble integrations and potential predictability of atmospheric seasonal variations in the extratropics. Journal of Climate, 11(11), 2942-2959.
[ Abstract PDF ]An approach to assess the potential predictability of the extratropical atmospheric seasonal variations in an ensemble of atmospheric general circulation model (AGCM) integrations has been proposed in this study by isolating reproducible forced modes and examining their contributions to the local ensemble mean. The analyses are based on the monthly mean output of an eight-member ensemble of 10-yr Atmospheric Model Intercomparison Project integrations with a T42L18 AGCM.
An EOF decomposition applied to the ensemble anomalies shows that there exist some forced modes that are less affected by the internal process and thus appear to be highly reproducible. By reconstructing the ensemble in terms of the more reproducible forced modes and by developing a quantitative measure, the potential predictability index (PPI), which combines the reproducibility with the local variance contribution, the local ensemble mean over some selective geographic areas in the extratropics was shown to result primarily from reproducible forced modes rather than internal chaotic fluctuations. Over those regions the ensemble mean is potentially predictable. Extratropical potentially predictable regions are found mainly over North America and part of the Asian monsoon regions. Interestingly, the potential predictability over some preferred areas such as Indian monsoon areas and central Africa occasionally results primarily from non-ENSO-related boundary forcing, although ENSO forcing generally dominates over most of the preferred areas.
The quantitative analysis of the extratropical potential predictability with PPI has shown that the preferred geographic areas have obvious seasonality. For the 850-hPa temperature, for example, potentially predictable regions during spring and winter are confined to Alaska, northwest Canada, and the southeast United States, the traditional PNA region, while during summer and fall they are favored over the middle part of North America. It has also been shown that the boreal summer season (June-August) possesses the largest potentially predictable area, which seems to indicate that it is a favored season for the extratropical potential predictability. On the contrary, boreal winter (December-February) appears to have a minimum area of extratropical potential predictability.
The results have been compared with the more traditional statistical tests for potential predictability and with observations from the National Centers for Environmental Prediction reanalysis, which indicates that the PPI analysis proposed here is successful in revealing extratropical potential predictability determined by the external forcing.
- Stern, William F., 1997: Investigations into the role of tropical intraseasonal oscillations in seasonal NWP In CAS/JSC Working Group on Numerical Experimentation - Research Activities in Atmospheric & Oceanic Modelling, WMO/ICSU/IOC World Climate Research Programme, Report No. 25, WMO/TD-No. 792, Geneva, Switzerland, World Meteorological Organization, 8.59-8.60.
- Vitart, F, Jeffrey L Anderson, and William F Stern, 1997: Simulation of interannual variability of tropical storm frequency in an ensemble of GCM integrations. Journal of Climate, 10(4), 745-760.
[ Abstract PDF ]The present study examines the simulation of the number of tropical storms produced in GCM integrations with a prescribed SST. A 9-member ensemble of 10-yr integrations (1979-88) of a T42 atmospheric model forced by observed SSTs has been produced; each ensemble member differs only in the initial atmospheric conditions. An objective procedure for tracking model-generated tropical storms is applied to this ensemble during the last 9 yrs of the integrations (1980-88). The seasonal and monthly variations of tropical storm numbers are compared with observations for each ocean basin.
Statistical tools such as the Chi-square test, the F test, and the t test are applied to the ensemble number of tropical storms, leading to the conclusion that the potential predictability is particularly strong over the western North Pacific and the eastern North Pacific, and to a lesser extent over the western North Atlantic. A set of tools including the joint probability distribution and the ranked probability score are used to evaluate the simulation skill of this ensemble simulation. The simulation skill over the western North Atlantic basin appears to be exceptionally high, particularly during years of strong potential predictability.
- Anderson, Jeffrey L., and William F Stern, 1996: Evaluating the potential predictive utility of ensemble forecasts. Journal of Climate, 9(2), 260-269.
[ Abstract PDF ]A method is presented for determining when an ensemble of model forecasts has the potential to provide some useful information. An ensemble forecast of a particular scalar quantity is said to have potential predictive utility when the ensemble forecast distribution is significantly different from an appropriate climatological distribution. Here, the potential predictive utility is measured using Kuiper's statistical test for comparing two discrete distributions. More traditional measures of the potential usefulness of an ensemble forecast based on ensemble mean or variance discard possibly valuable information by making implicit assumptions about the distributions being compared
Application of the potential predictive utility to long integrations of an atmospheric general circulation model in a boundary value problem (an ensemble of Atmospheric Model Intercomparison Project integrations) reveals a number of features about the response of a GCM to observed sea surface temperatures. In particular, the ensemble of forecasts is found to have potential predictive utility over large geographic areas for a number of atmospheric fields during strong El Niño-Southern Oscillation anomalous events. Unfortunately, there are only limited areas of potential predictive utility for near-surface fields and precipitation outside the regions of the tropical oceans. Nevertheless, the method presented here can identify all areas where the GCM ensemble may provide useful information, whereas methods that make assumptions about the distribution of the ensemble forecast variables may not be able to do so.
- D'Andrea, F, and William F Stern, et al., 1996: Northern Hemisphere Atmospheric Blocking as Simulated by 15 Atmospheric General Circulation Models in the Period 1979-1988 (Results from an AMIP Diagnostic Subproject, WCRP-96. WMO/TD-No. 784, WMO, Geneva, Switzerland: World Meteorological Organization, 25 pp.
[ Abstract ]As a part of the Atmospheric Model Intercomparison Project (AMIP), the behavior of 15 General Circulation Models has been analyzed in order to diagnose and compare the ability of the different models in simulating midlatitude atmospheric blocking. In accordance to the established AMIP procedure, the 10-year model integrations were performed using prescribed, time evolving monthly mean observed SSTs spanning the period January 1979- December 1988. Atmospheric observational data (ECMWF analysis) over the same period have been also used to verify the models results.
The models involved in this comparison represent a wide spectrum of model complexity, with different horizontal and vertical resolution, numerical techniques and physical parameterizations, and exhibit large differences in blocking behavior. Nevertheless, a few common features can be found, such as the general tendency to underestimate both blocking frequency and the average duration of blocks.
The relation between model blocking and systematic errors has also been assessed, although without resorting to ad-hoc numerical experimentation it is impossible to relate with certainty particular model deficiencies in representing blocking to precise parts of the model formulation.
- Stern, William F., and Jeffrey L Anderson, 1996: Interannual variability of tropical intraseasonal oscillations in the GFDL/DERF GCM inferred from an ensemble of AMIP integrations In 11th Conference on Numerical Weather Prediction, Boston, MA, American Meteorological Society, 15-16.
- Stern, William F., and C Tony Gordon, 1996: Specifying and modeling snow cover in multi-year GCM integrations In Research Activities in Atmospheric and Oceanic Modelling, CAS/JSC Working Group on Numerical Experimentation, Report No. 23 WMO/TD No. 734, World Meteorological Organization, 4.43-4.44.
- Vitart, F, Jeffrey L Anderson, and William F Stern, 1996: Potential predictability of tropical storms in an ensemble of forecasts In Proceedings of the 20th Annual Climate Diagnostics Workshop, U.S. Dept. of Commerce/NOAA/NWS, 263-266.
- Vitart, F, Jeffrey L Anderson, and William F Stern, 1996: Potential predictability of tropical storms in ensemble GCM simulations In Research Activities in Atmospheric and Oceanic Modelling, CAS/JSC Working Group on Numerical Experimentation, Report No. 23 WMO/TD No. 734, World Meteorological Organization, 6.32.
- Anderson, Jeffrey L., and William F Stern, 1995: A method of evaluating the predictive ability of ensemble forecasts In Proceedings of the 19th Annual Climte Diagnostics Workshop, Springfield, VA, NTIS, 472-475.
- Gleckler, P J., Kikuro Miyakoda, and William F Stern, et al., 1995: Cloud-radiative effects on implied oceanic energy transports as simulated by atmospheric general circulation models. Geophysical Research Letters, 22(7), 791-794.
[ Abstract PDF ]This paper summarizes the ocean surface net energy flux simulated by fifteen atmospheric general circulation models constrained by realistically-varying sea surface temperatures and sea ice as part of the Atmospheric Model Intercomparison Project. In general, the simulated energy fluxes are within the very large observational uncertainties. However, the annual mean oceanic meridional heat transport that would be required to balance the simulated surface fluxes is shown to be critically sensitive to the radiative effects of clouds, to the extent that even the sign of the Southern Hemisphere ocean heat transport can be affected by the errors in simulated cloud-radiation interactions. It is suggested that improved treatment of cloud radiative effects should help in the development of coupled atmosphere-ocean general circulation models.
- Miyakoda, Kikuro, Joseph J Sirutis, Anthony Rosati, C Tony Gordon, Rich Gudgel, William F Stern, Jeffrey L Anderson, and A Navarra, 1995: Atmospheric parameterizations in coupled air-sea models used for forecasts of ENSO In Proceedings of the International Scientific Conference on the Tropical Ocean Global Atmosphere (TOGA) Programme, WCRP-91, WMO/TD No. 717, Geneva, Switzerland, World Meteorological Organization, 802-806.
[ Abstract ]In order to investigate the feasibility of seasonal forecasts, a prediction system is developed. Here the main theme is the study of atmospheric physics parameterization for coupled air-sea modeling. The oceanic GCM uses 1 degree global grid with a finer resolution in the equatorial belt. The atmospheric GCM has the spectral T30 representation, which includes all of the usual physics parameterizations. Using a first version of the model (Coupled Model I) and a set of appropriate initial conditions, the capability of El Niño and La Niña forecasting with a 13 month lead time was tested, resulting in successful forecasts of the 1982/83 and 1988/89 events (Rosati et al., 1995b). However, longer runs of this system have revealed a sizable systematic error in simulations with a tendency to cool most of the world ocean, particularly the western tropical Pacific, and also without an adequate annual cycle of the SST in the eastern tropical Pacific.
In order to improve some of these features, particularly the ENSO phenomena, various versions of the atmospheric parameterizations and mountain representation are incorporated into the atmospheric GCM, and the model simulations are examined. The experiments are divided into two steps: one is with the uncoupled atmospheric model, and the other is with the coupled model. In the first step, five year simulations are carried out with the observed SST prescribed, and the results are compared with observations, which enables one to make the critical validation of the model. The second step is to couple the atmospheric and oceanic models, and integrate them from a January 1982 initial condition for 7 years, and also for another initial condition, i.e., January, 1988 for 13 months.
Compared with the boundary forced simulation, the coupling process introduces more degree of freedom, with increase of the sensitivity as well as the complexity considerably. In particular, the El Niño simulation is sensitive to any change of physics. For this reason, the objective of the simulation is focused only on the equatorial Pacific process and secondly the Indian monsoon, as opposed to the overall improvement of the general circulation. In other words, the approach is close to that of mechanistic modeling with specific targets rather than that of a GCM with broader objectives. The research is proceeding in two directions. One is: investigating the model's sensitivity for El Niño and La Niña processes to variation in a coupling parameter. The second is: after a number of trial-and-error experiments on various combinations of the parameterizations, the second atmospheric model, i.e., Model II, is selected. It is shown that Coupled Model II performs substantially better in some aspects but worse in other aspects than Coupled Model I. The improvement is found in the SST: warming occurs not only over the equatorial Pacific but also over the whole globe. The SST increase is achieved by the strong effect of the cumulus convection. On the other hand, some deficiencies remain the same in both models, i.e., the large positive errors of the SST in the eastern oceans, the lack of an annual cycle of the SST in the eastern equatorial Pacific, and the failure in forecast of the second El Niño. In summary, the prediction of the Southern Oscillation has been achieved by the two models for a full first cycle but not for the second cycle .
- Stern, William F., 1995: Tropical intraseasonal variability in the GFDL/DERF GCM during AMIP In Proceedings of the First International AMIP Scientific Conference, WCRP-92, WMO/TD No. 273, Geneva, Switzerland, World Meteorological Organization, 137-142.
- Stern, William F., and Kikuro Miyakoda, 1995: Feasibility of seasonal forecasts inferred from multiple GCM simulations. Journal of Climate, 8(5), 1071-1085.
[ Abstract PDF ]Assuming that SST provides the major lower boundary forcing for the atmosphere, observed SSTs are prescribed for an ensemble of atmospheric general circulation model (GCM) simulations. The ensemble consists of 9 "decadal" runs with different initial conditions chosen between 1 January 1979 and 1 January 1981 and integrated about 10 years. The main objective is to explore the feasibility of seasonal forecasts using GCMs. The extent to which the individual members of the ensemble reproduce the solutions of each other (i.e., reproducibility) may be taken as an indication of potential predictability. In addition, the ability of a particular GCM to produce realistic solutions, when compared with observations, must also be addressed as part of the predictability problem.
A measure of reproducibility may be assessed from the spread among ensemble members. A normalized spread index, can be defined at any point in space and time, as the variability of the ensemble normalized by the climatological seasonal variability. In the time mean it is found that the reproducibility is significantly below unity for certain regions. Low values of the spread index are seen generally in the Tropics, whereas the extratropics does not exhibit a high degree of reproducibility. However, if one examines plots in time of seasonal mean for the U.S. region, for example, it is found that for certain periods this index is much less than unity, perhaps implying "occasional potential predictability." In this regard, time series of ensemble mean soil moisture and precipitation over the United States are compared with corresponding observations. This study reveals some marginal skill in simulating periods of drought and excessive wetness over the United States during the 1980s (i.e., the droughts of 1981 and 1988, and the excessive wetness during the 1982/83 El Niño). In addition, by focusing on regions of better time-averaged reproducibility - that is, the southeast United States and northeast Brazil - a clearer indication of a relationship between good reproducibility and seasonal predictability seems to emerge.
- Stern, William F., and Kikuro Miyakoda, 1995: Interannual variability and reproducibility from multiple GCM simulations In Proceedings of the 19th Annual Climate Diagnostics Workshop, Springfield, VA, NTIS, 92-95.
- Ploshay, Jeff J., William F Stern, and Kikuro Miyakoda, 1992: FGGE reanalysis at GFDL. Monthly Weather Review, 120(9), 2083-2108.
[ Abstract PDF ]The reanalysis of FGGE [First GARP (Global Atmospheric Research Program) Global Experiment] data for 128 days during two special observing periods has been performed, using an improved data-assimilation system and the revised FGGE level II dataset. The data-assimilation scheme features forward continuous (in time) data injection in both the original and the new systems. However, the major revisions in the new system include a better first guess and a more efficient dynamical balancing for the assimilation of observed data. The results of the implementation of this system are assessed by intercomparisons among the new FGGE analysis of other institutions such as ECMWF (European Centre for Medium-Range Weather Forecasts) and NMC (National Meteorological Center, Washington, D.C.), and also the original GFDL (Geophysical Fluid Dynamics Laboratory) analysis. The quality of the new GFDL analysis is now comparable to those of the other two institutions. However, the moisture analysis appears to be appreciably different, suggesting that the cumulus convection parameterizations and the boundary-layer moisture fluxes in the models are responsible for this discrepancy.
A detailed investigation of the results has been carried out by comparing the analyses with radiosonde observations. This verification reveals that temperature and wind differences have been reduced considerably from the original to the new GFDL analysis; they are now competitive with those of ECMWF and NMC, while with regard to the geopotential height, differences of the GFDL reanalysis are larger than the original GFDL as well as the ECMWF and the NMC. A comparative study is also made with UCLA analyses over Asia in connection with the Indian monsoon. The results indicate that the qualities of both analyses are comparable. The capability of representing Madden-Julian oscillations in the reanalysis and in the ECMWF and old GFDL analysis is investigated by comparing with satellite observations. It is revealed that these oscillations are successfully reproduced by the new analysis; however, the agreement with the satellite data is not quite satisfactory. The utilization of satellite-observed wind (satobs) and aircraft data (aireps) in the data assimilation needs particular care. It appears that the quality control of these data in the GFDL reanalysis is too restrictive; in other words, the toss-out criterion of wind data is too small. A consequence of the failure to accept some single-level data turns out to be a fairly large discrepancy in representing the maximum wind speed in the analysis. It is also discussed that the current forward continuous-injection scheme is not adequate to obtain diabatic quantities for the archive.
- Stern, William F., and Jeff J Ploshay, 1992: A scheme for continuous data assimilation. Monthly Weather Review, 120(7), 1417-1432.
[ Abstract PDF ]Major revisions to the Geophysical Fluid Dynamics Laboratory's (GFDL) continuous data assimilation system have been implemented and tested. Shortcomings noted during the original processing of data from FGGE [First GARP (Global Atmospheric Research Program) Global Experiment ] served as the basis for these improvements. This new system has been used to reanalyze the two FGGE special observing periods. The main focus here will be on assessing the changes to the assimilation system using comparisons of rerun test results with results from the original FGGE processing.
The key new features in the current system include: a reduction in the assimilation cycle from 12 to 6 h; the use of a 6-h forecast first guess for the OI (optimum-interpolation analysis) as opposed to the previous use of persistence as a first guess; an extension of the OI search range from 250 to 500 km with an increase in the maximum number of observations used per analysis point from 8 to 12; the introduction of incremental linear normal-mode initialization, eliminating the periodic nonlinear normal-mode initialization; and an increase in the horizontal resolution of the assimilating model from 30 waves to 42 waves, rhomboidally truncated. Tests of the new system show a significant reduction in the level of noise, improved consistency between mass and momentum analyses, and a better fit of the analyses to observations. In addition, the new system has demonstrated a greater ability to resolve rapidly moving and deepening transient features, with an indication of less rejection of surface pressure data.
In addition to the quantities archived during the original FGGE data processing, components of diabatic heating from the assimilating model have also been archived. They should be used with caution to the extent that they reflect model bias and spinup in addition to real features of the general circulation.
- Stern, William F., and Kikuro Miyakoda, 1988: Systematic errors in GFDL's extended range prediction spectral GCM In Workshop on Systematic Errors in Models of the Atmosphere, CAS/JSC Working Group on Numerical Experimentation, Report No. 12, WMO/TD No. 273, World Meteorological Organization, 78-85.
- Stern, William F., R T Pierrehumbert, Joseph J Sirutis, Jeff J Ploshay, and Kikuro Miyakoda, 1986: Recent developments in the GFDL extended-range forecasting system In Short- and Medium-Range Numerical Weather Prediction, Collection of papers presented at the WMO/IUGG NWP Symposium, World Meteorological Organization, 359-363.
[ Abstract ]An assessment is made of the areas of focus for improving extended-range forecasting. Two topics currently being researched involve the reduction of systematic error by improving a GCM's accuracy and the refinement of the transition between the data assimilation phase and the forecasting phase.
Subgrid-scale orographic parameterizations have been the subject of recent model improvement activities. Results are shown for an envelope orography with an N48L9 gridpoint model and using a mountain gravity wave drag scheme with an R42L18 spectral model. In both cases there is an encouraging reduction in the systematic errors.
Proper initialization of tropical features, i.e., 40-50 day waves, may be crucial for extended-range predictions in the extra-tropics as well as the tropics. Using a continuous data assimilation scheme the 40-50 day oscillations in the tropics appear to be well maintained from the assimilation to the forecast phase. However, the assimilation system underestimates precipitation and evaporation rates.
- Gordon, C T., William F Stern, and R D Hovanec, 1985: A simple scheme for generating two layers of radiatively constrained effective clouds in GCMs. Journal of Geophysical Research, 90(D6), 10,563-10,585.
[ Abstract PDF ]GCM-dependent, radiatively constrained cloud amount fields could be preferable to currently available observed fields for calculating radiative fluxes in GCM's used in long-range weather forecast studies. Motivated by this premise, we formulate an economical effective cloud algorithm for GCM's called "SATCLD," which utilizes compact, readily accessible analyses of observed satellite-derived radiative flux data. We then examine the plausibility of preliminary effective cloud fields. Analysis of SATCLD and other cloud fields and associated radiative fluxes (diagnosed by our GCM's cloud radiation model from observed atmospheric temperature and water vapor data) also provides some insight into biases in our GCM's cloud radiation model and surface albedo field. "SATCLD" generates effective low and high cloud amounts at each GCM grid point by minimizing the sum of the squares of the local residual (i.e., model-diagnosed minus observed) shortwave and longwave radiative fluxes. The preliminary SATCLD effective cloud amount fields seem plausible in many respects, based upon comparison with the satellite-derived 3DNEPH and surface-based SFCOBS analyses. In the tropics, the SATCLD effective high cloud amount is rather well correlated with 3DNEPH, while systematic differences in low cloud amount are evident in July. Off the west coasts of Central and South America and southern Africa, the SATCLD effective low cloud resembles SFCOBS in July. At mid-latitudes, the strongest similarities are between SATCLD versus 3DNEPH high cloud amount in July and SATCLD versus SFCOBS low cloud amount over the oceans. The SATCLD analysis is ill conditioned in the polar night region. Limitations of the present scheme as well as deficiencies in our GCM's cloud-radiation model and surface albedo fields and in the archived satellite data are discussed. Suggestions are made for reducing discrepancies between effective versus real clouds without sacrificing consistency between GCM-diagnosed versus observed radiative fluxes.
- Gordon, C T., R D Hovanec, and William F Stern, 1984: Analyses of monthly mean cloudiness and their influence upon model-diagnosed radiative fluxes. Journal of Geophysical Research, 89(D3), 4713-4738.
[ Abstract PDF ]Two different monthly mean analyses of low, middle, and high cloud amounts for January 1977 and July 1979 are compared: 3DNEPH is a condensed version (northern hemisphere only) of the Air Force 3D-Neph analysis, which incorporates satellite data plus surface observations of clouds and auxilliary meteorological data. SFCOBS is objectively analyzed from surface observations of clouds. The SFCOBS and 3DNEPH analyses of low cloud amounts agree qualitatively in the winter extratropics. The 3DNEPH ITCZ is much more sharply defined than the SFCOBS. The sensitivity of radiative fluxes to 3DNEPH, SFCOBS, and zonal mean 3DNEPH clouds is then evaluated. The fluxes are diagnosed by a cloud-radiation model utilizing "observed" monthly mean temperature and water vapor fields and are verified against satellite data. The outgoing longwave radiative flux clearly verifies best for 3DNEPH clouds and worst for zonal mean 3DNEPH clouds in the tropics. It is predominately controlled by surface temperature in the winter extratropics. Generally speaking, the shortwave fluxes do not verify as well as the longwave fluxes. Also, outside of the winter extratropics, the net radiative fluxes correlate poorly with observation. Biases in the zonal mean long and shortwave fluxes can be reduced by adjusting other cloud-related parameters. Based upon the above results, it may be worthwhile to construct a monthly mean cloud climatology from a condensed version of the 3D-Neph. However, alternative strategies should also be explored, such as the development of cloud analysis schemes that constrain the model-diagnosed net radiative flux to be consistent with observation.
- Gordon, C T., and William F Stern, 1984: Medium range prediction by a GFDL global spectral model: results for three winter cases and sensitivity to dissipation. Monthly Weather Review, 112(2), 217-245.
[ Abstract PDF ]A preliminary evaluation is made of the medium range predictive capability of a GFDL global spectral model of the atmosphere, based upon three winter blocking cases. Analogous forecasts by a GFDL global grid point model provide a background standard of comparison. The spectral model is rhomboidally truncated at wavenumber 30, has 9 sigma levels, incorporates sub-grid scale physical processes commonly associated with general circulation models and employs semi-implicit time differencing. The grid point model has somewhat finer horizontal resolution and fairly similar sub-grid scale physical processes, and employs explicit time differencing. The spectral model is up to 6 times more economical. The level of forecast skill for the 5 to 15 day range is generally less than practically useful and is more case-dependent than spectral versus grid point model-dependent. In the most successful case, i.e., 16 January 1979, an observed Atlantic blocking ridge is simulated quite well, especially by the spectral model. The predicted Atlantic ridges tend to retard approaching upstream transient disturbances. A zonal bias of the midlatitude circulation, which develops in all three spectral and grid point model predictions is most pronounced in the spectral model forecast from 1 January 1977. Results of a diffusion sensitivity experiment and other evidence suggest that insufficient frictional dissipation may have enhanced the zonal bias of the above forecast. The bias diminishes, consistent with a redistribution of spectral kinetic energy among zonal wavenumbers 0, 1 and 2, if a static stability-dependent parameterization of vertical mixing or stronger Δ4 horizontal diffusion are used. Also, the predicted-enstrophy spectrum at midlatitudes steepens, given the stronger Δ4 horizontal diffusion.
- Miyakoda, Kikuro, C Tony Gordon, R Caverly, William F Stern, Joseph J Sirutis, and W Bourke, 1983: Simulation of a blocking event in January 1977. Monthly Weather Review, 111(4), 846-869.
[ Abstract PDF ]January 1977 was a month noted for its extraordinary weather over North America. The winter was dominated by two persistent large amplitude ridges positioned over the west coast of North America and the Icelandic region of the Atlantic Ocean. A very intense trough reached deep into the eastern United States and caused one of the coldest Januaries on record. One-month integrations of various GCMs were conducted in order to test their ability to simulate this blocking event. Reasonably high resolution finite difference and spectral models available at GFDL were used. Each GCM was integrated from three different analyses of the initial conditions. For some models, a fairly accurate forecast was obtained and considerable skill was recognized in the simulation of the 30-day evolution in terms of the 5-day or 10-day mean flow fields, including the period of record breaking coldness over the eastern United States. The main conclusion is that proper treatment of the subgrid-scale processes as well as sufficient spatial resolution are essential for the simulations of this phenomenon as an initial value problem. Weak zonal wind poleward of about 40 degrees N and upstream of the blocking ridge appears to be crucial for the successful simulation of the sustained blocking ridge.
- Gordon, C T., and William F Stern, 1982: A description of the GFDL global spectral model. Monthly Weather Review, 110(7), 625-644.
[ Abstract PDF ]A multi-level, global, spectral transform model of the atmosphere, based upon spherical harmonics, has been developed at GFDL. The basic model has nine sigma levels in the vertical and rhomboidal spectral truncation at wavenumber 30. However, finer spectral or vertical resolution versions are available as well. The model's efficient semi-implicit time differencing scheme does not appear to adversely affect medium range predictions. The model has physical processes commonly associated with grid point GCM's. Two unique features are a linearized virtual temperature correction and an optional, spectrally-computed non-linear horizontal diffusion scheme. A parameterization of vertical mixing based upon the turbulent closure method is also optional.
The GFDL spectral model has been widely utilized at GFDL for extended range weather prediction experiments. In addition, it has been adapted and applied to climate studies, four-dimensional data assimilation experiments and even to the atmosphere of Venus. These applications are briefly reviewed.
- Gordon, C T., and William F Stern, 1982: Sensitivity of GCM-diagnosed radiative fluxes to "observed" monthly mean distributions of cloudiness In Proceedings of the Sixth Annual Climate Diagnostics Workshop, Washington, DC, NOAA, 268-278.
Direct link to page: http://www.gfdl.noaa.gov/bibliography/resultstest.php?author=1098