SEARCH · Search NASA

Results for “general circulation model”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 415 records · Page 23

All sources NASA NTRS DOE OSTI

DOE OSTI2021

Retrospective Analysis and Bayesian Model Averaging of CMIP6 Precipitation in the Nile River Basin

The Nile River basin is one of the global hotspots vulnerable to climate change impacts because of a fast-growing population and geopolitical tensions. Previous studies demonstrated that general circulation models (GCMs) frequently show disagreement in the sign of change in annual precipitation projections. Here, we first evaluate the performance of 20 GCMs from phase six of the Coupled Model Intercomparison Project (CMIP6) benchmarked against a high-spatial-resolution precipitation dataset dating back to 1983 from Precipitation Estimation from Remotely Sensed Information Using Artificial Neural Networks–Climate Data Record (PERSIANN-CDR). Next, a Bayesian model averaging (BMA) approach is adopted to derive probability distributions of precipitation projections in the Nile basin. Retrospective analysis reveals that most GCMs exhibit considerable (up to 64% of mean annual precipitation) and spatially heterogenous bias in simulating annual precipitation. Moreover, it is shown that all GCMs underestimate interannual variability; thus, the ensemble range is underdispersive and is a poor indicator of uncertainty. The projected changes from the BMA model show that the value and sign of change vary considerably across the Nile basin. Specifically, it is found that projected changes in the two headwaters basins, namely, the Blue Nile and Upper White Nile, are 0.03% and -1.65%, respectively; both are statistically insignificant at α = 0.05. The uncertainty range estimated from the BMA model shows that the probability of a precipitation decrease is much higher in the Upper White Nile basin whereas projected change in the Blue Nile is highly uncertain both in magnitude and sign of change.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Convection-Permitting Simulations With the E3SM Global Atmosphere Model

This paper describes the first implementation of the Δx = 3.25 km version of the Energy Exascale Earth System Model (E3SM) global atmosphere model and its behavior in a 40-day prescribed-sea-surface-temperature simulation (January 20 through February 28, 2020). This simulation was performed as part of the DYnamics of the Atmospheric general circulation Modeled On Non-hydrostatic Domains (DYAMOND) Phase 2 model intercomparison. Effective resolution is found to be ~ 6x the horizontal dynamics grid resolution despite using a coarser grid for physical parameterizations. Despite this new model being in an immature and untuned state, moving to 3.25 km grid spacing solves several long-standing problems with the E3SM model. In particular, Amazon precipitation is much more realistic, the frequency of light and heavy precipitation is improved, agreement between the simulated and observed diurnal cycle of tropical precipitation is excellent, and the vertical structure of tropical convection and coastal stratocumulus look good. In addition, the new model is able to capture the frequency and structure of important weather events (e.g., tropical cyclones, extratropical cyclones including atmospheric rivers, and cold air outbreaks). Interestingly, this model does not get rid of the erroneous southern branch of the intertropical convergence zone nor the tendency for strongest convection to occur over the Maritime Continent rather than the West Pacific, both of which are classic climate model biases. Several other problems with the simulation are identified, underscoring the fact that this model is a work in progress.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS2005

Data Assimilation and Data Fusion for Planetary Atmospheres

The overarching goal of this Cooperative Agreement was to develop a model and procedures for the data assimilation of planetary spacecraft atmospheric observations. Data assimilation - in its application to weather analysis and prediction - is the process of finding an initial state of the meteorological variables (winds, temperatures, pressures, etc.) of an atmosphere, which, when propagated forward in time using a deterministic general circulation model, reproduces all of the available observations over that time to within the measurement and computational errors. With this definition, data assimilation is seen to be a natural extension of well-known least-squares minimization techniques. The primary complication arises from the scale of the problem: For the Martian atmosphere with the available nadir-viewing Thermal Emission Spectrometer data from Mars Global Surveyor, approximately 1,000,000 individual measurements of channel radiances (in the 15-micrometer region, where these radiances relate directly to the surface and atmospheric temperature) were made per day. A suitable general circulation model for dealing with this data set has on the order of 20,000 independent variables. After some spatial and temporal averaging of the data - which provides a necessary statistical estimate of the representativeness of the measurements, a crucial issue in data assimilation - the problem reduces in scale to the solution of approximately 50,000 equations for the 20,000 variables.

Houben, Howard↗

NASA NTRS2000

Observed Changes in the Gravity Field of Mars Due to Seasonal Atmospheric Processes

The atmosphere of Mars deposits approximately 30% of its mass at the winter pole as part of its seasonal cycle of CO2 exchange and sublimes it back to the atmosphere in the spring, thus creating an annual hemispheric cycle of mass re-distribution. Using X-band tracking data of the Mars Global Surveyor (MGS) spacecraft, we have detected the signature of changes in the low degree gravity field from March 1999 through August 2000, corresponding to about three-quarters of a Martian year. The observed variations show a general resemblance to predicted variations from general circulation models. Also observed are irregular changes that appear to be due to transient phenomena in the Martian atmosphere such as large dust storms that provide significant heat into the lower atmosphere, even in the polar regions. In addition, we have identified a change in the rate of rotation of Mars over this same period that appears somewhat smaller than anticipated, but in general agreement with general circulation models.

Smith, David E.↗

NASA NTRS1981

A stochastic-dynamical approach to the study of the natural variability of the climate

A method, suggested by Leith (1975), which employed stochastic-dynamic forecasts obtained from a general circulation model in such a way as to satisfy the definition of climatic noise, was used to validate assumptions accounting for the effects of external influences in estimating the climatic noise. Two assumptions were investigated: (1) that the weather fluctuations can be represented as a Markov process, and (2) that changing external conditions do not influence the atmosphere's statistical properties on short time scales. The general circulation model's simulation of the daily weather fluctuations was generated by performing integrations with prescribed climatological boundary conditions for random initial atmospheric states, with resulting dynamical forecasts providing an ensemble of simulated data for the autoregressive modeling of weather fluctuations. To estimate the climatic noise from the observational data (consisting of hourly values of sea level pressure and surface temperature at 54 U.S. stations for the month of January for the years 1949-1975) use of the short time-scale assumption is made. The simulated and observed data were found not to be consistent with either white noise or a Markov process of weather fluctuations. Good agreement was found between the results of the hypothetical testing of the simulated and the observed surface temperatures; and only partial support was found for the short time-scale assumption, i.e., for sea level pressure.

Straus, D. M.↗

NASA NTRS1993

Seasonal variation of surface and atmospheric cloud radiative forcing over the globe derived from satellite data

Global distributions of surface and atmospheric cloud radiative forcing parameters have been derived using parameterized radiation models with satellite meteorological data from the International Satellite Cloud Climatology Project, and directly measured top-of-atmosphere radiative fluxes from the Earth Radiation Budget Experiment. Specifically, shortwave, longwave, and total cloud forcing at the surface, and column-averaged values of longwave cloud forcing of the atmosphere were derived for the midseasonal months of April, July, and October 1985 and January 1986, covering a complete annual cycle. Seasonal variability is illustrated by comparing the results for July 1985 and January 1986, which represent the seasonal extremes. Surface shortwave cloud forcing is always negative, representing a cooling of the surface, with strongest cooling (-120 to -180 W/sq m) occurring over midlatitude storm tracks of the summer hemisphere. Surface longwave cloud forcing is always positive, representing a warming of the surface, with strongest warming (60 to 75 W/sq m) occurring over storm tracks of the winter hemisphere. Zonal averages show the entire summer hemisphere dominated by shortwave cooling, the middle and high latitudes of the winter hemisphere dominated by longwave warming, and a broad zone of transition in between. The globally averaged total cloud forcing amounts to a cooling throughout the year, ranging from a low of about -12 W/sq m for July 1985 to a high of about -25 W/sq m for January 1986. The longwave cloud forcing of the atmosphere shows a strong warming over deep convective regions in the tropics and a moderate cooling outside the tropics, amounting to a weak cooling (-2 to -5 W/sq m) in the global average. Comparisons of the results with general circulation model simulations show broad qualitative agreement regarding the locations of prominent warming and cooling regions. Quantitative comparisons, on the other hand, show significant differences between the magnitudes of warming and cooling in these regions. Most of the larger differences can be attributed to known deficiences of the general circulation model simulations. Comparisons with satellite-derived results of other investigators show much better agreement.

Gupta, Shashi K.↗

NASA NTRS2011

Mars-GRAM 2010: Improving the Precision of Mars-GRAM

It has been discovered during the Mars Science Laboratory (MSL) site selection process that the Mars Global Reference Atmospheric Model (Mars-GRAM) when used for sensitivity studies for Thermal Emission Spectrometer (TES) MapYear=0 and large optical depth values, such as tau=3, is less than realistic. Mars-GRAM's perturbation modeling capability is commonly used, in a Monte-Carlo mode, to perform high fidelity engineering end-to-end simulations for entry, descent, and landing (EDL). Mars-GRAM 2005 has been validated against Radio Science data, and both nadir and limb data from TES. Traditional Mars-GRAM options for representing the mean atmosphere along entry corridors include: (1) TES mapping year 0, with user-controlled dust optical depth and Mars-GRAM data interpolated from NASA Ames Mars General Circulation Model (MGCM) results driven by selected values of globally-uniform dust optical depth, or (2) TES mapping years 1 and 2, with Mars-GRAM data coming from MGCM results driven by observed TES dust optical depth. From the surface to 80 km altitude, Mars-GRAM is based on NASA Ames MGCM. Above 80 km, Mars-GRAM is based on the University of Michigan Mars Thermospheric General Circulation Model (MTGCM). MGCM results that were used for Mars-GRAM with MapYear=0 were from a MGCM run with a fixed value of tau=3 for the entire year at all locations. This choice of data has led to discrepancies that have become apparent during recent sensitivity studies for MapYear=0 and large optical depths. Unrealistic energy absorption by time-invariant atmospheric dust leads to an unrealistic thermal energy balance on the polar caps. The outcome is an inaccurate cycle of condensation/sublimation of the polar caps and, as a consequence, an inaccurate cycle of total atmospheric mass and global-average surface pressure. Under an assumption of unchanged temperature profile and hydrostatic equilibrium, a given percentage change in surface pressure would produce a corresponding percentage change in density at all altitudes. Consequently, the final result of a change in surface pressure is an imprecise atmospheric density at all altitudes.

Justh, H. L.↗

NASA NTRS2004

Lunar and Planetary Science XXXV: Special Session: Mars Climate Change

The session "Mars Climate Change" contained the following reports:Geological Evidence for Climate Change on Mars; A New Astronomical Solution for the Long Term Evolution of the Insolation Quantities of Mars; Interpreting Martian Paleoclimate with a Mars General Circulation Model; History and Progress of GCM Simulations on Recent Mars Climate Change; Northern and Southern Permafrost Regions on Mars with High Content of Water Ice: Similarities and Differences; Periods of Active Permafrost Layer Formation in the Recent Geological History of Mars; Microclimate Zones in the Dry Valleys of Antarctica: Implications for Landscape; Evolution and Climate Change on Mars; Geomorphic Evidence for Martian Ground Ice and Climate Change; Explaining the Mid-Latitude Ice Deposits with a General Circulation Model; Tharsis Montes Cold-based Glaciers: Observations and Constraints for Modeling and Preliminary Results; Ice Sheet Modeling: Terrestrial Background and Application to Arsia Mons Lobate Deposit, Mars; Enhanced Water-Equivalent Hydrogen on the Western Flanks of the Tharsis Montes and Olympus Mons: Remnant Subsurface Ice or Hydrate Minerals?; and New Age Mars.

Source record↗

NASA NTRS2004

Lunar and Planetary Science XXXV: Special Session: Mars Climate Change

The session "Mars Climate Change" included the following topics:Geological Evidence for Climate Change on Mars; A New Astronomical Solution for the Long Term Evolution of the Insolation Quantities of Mars; Interpreting Martian Paleoclimate with a Mars General Circulation Model; History and Progress of GCM Simulations on Recent Mars Climate Change; Northern and Southern Permafrost Regions on Mars with High Content of Water Ice: Similarities and Differences; Periods of Active Permafrost Layer Formation in the Recent Geological History of Mars; Microclimate Zones in the Dry Valleys of Antarctica: Implications for Landscape Evolution and Climate Change on Mars; Geomorphic Evidence for Martian Ground Ice and Climate Change; Explaining the Mid-Latitude Ice Deposits with a General Circulation Model; Tharsis Montes Cold-based Glaciers: Observations and Constraints for Modeling and Preliminary Results; Ice Sheet Modeling: Terrestrial Background and Application to Arsia Mons Lobate Deposit, Mars; Enhanced Water-Equivalent Hydrogen on the Western Flanks of the Tharsis Montes and Olympus Mons: Remnant Subsurface Ice or Hydrate Minerals?; and New Age Mars.

Source record↗

NASA NTRS1999

Satellite Data Assimilation

Satellite altimetry provides the only routine observation of a dynamic variable of the global ocean. For instance, measurements from the TOPEX/POSEIDON altimeter have provided dramatic basin-wide images of the space-time evolution of sea level associated with the 1997-98 El Nino event. The nature of large-scale ocean circulation is being studied that underlies such sea level changes measured by satellite altimetry. To first approximation, properties of large-scale (order 1000 km and larger) sea level variability is latitudinally dependent. An analysis using a general circulation model shows that sea level changes are largely due to wind-driven baroclinic (depth dependent) circulation in the tropics (within 20 degrees of the equator), but are primarily due to the expansion and contraction of near surface water in temperate latitudes (between 20 deg. and 40 deg.) forced by seasonal heating and cooling. In contrast, wind-driven barotropic (depth independent) circulation dominates sea level variability in high latitudes (40 deg. and higher) characterized by periods that are as short as a few days. The presence of such inhomogeneity and the significance of high-frequency, large-scale sea level changes had not been fully recognized prior to this study, and are summarized in Fukumori et al. Numerical models provide theoretical relationships among properties that can be inverted using observations so as to estimate the entire state of the ocean, including properties that are otherwise difficult to measure remotely. The process is data assimilation. An approximate Kalman filter and smoother have been devised to assimilate three years of TOPEX/POSEIDON sea level data into a global ocean general circulation model. The figure below demonstrates the skill of the assimilation, and shows the altimeter assimilated estimate being in closer agreement than the simulation is with independent in situ measurements of subsur ace temperature and velocity, consistent with formal uncertainty estimates. Results from this study have been summarized and submitted for publication. The study demonstrates the feasibility of global ocean data assimilation and illustrates applications in monitoring and understanding of processes controlling the evolution of the ocean. Additional information is contained in the original.

Fukumori, Ichiro↗

NASA NTRS2022

Scientific Data From Precipitation Driver Response Model Intercomparison Project

This data descriptor reports the main scientific values from General Circulation Models (GCMs) in the Precipitation Driver and Response Model Intercomparison Project (PDRMIP). The purpose of the GCM simulations has been to enhance the scientific understanding of how changes in greenhouse gases, aerosols, and incoming solar radiation perturb the Earth’s radiation balance and its climate response in terms of changes in temperature and precipitation. Here we provide global and annual mean results for a large set of coupled atmospheric-ocean GCM simulations and a description of how to easily extract files from the dataset. The simulations consist of single idealized perturbations to the climate system and have been shown to achieve important insight in complex climate simulations. We therefore expect this data set to be valuable and highly used to understand simulations from complex GCMs and Earth System Models for various phases of the Coupled Model Intercomparison Project

Precipitation Driver Response Model Intercompariso↗

NASA NTRS1999

Seasonal Distributions of Global Ocean Chlorophyll and Nutrients: Analysis with a Coupled Ocean General Circulation Biogeochemical, and Radiative Model

A coupled general ocean circulation, biogeochemical, and radiative model was constructed to evaluate and understand the nature of seasonal variability of chlorophyll and nutrients in the global oceans. The model is driven by climatological meteorological conditions, cloud cover, and sea surface temperature. Biogeochemical processes in the model are determined from the influences of circulation and turbulence dynamics, irradiance availability, and the interactions among three functional phytoplankton groups (diatoms, chorophytes, and picoplankton) and three nutrient groups (nitrate, ammonium, and silicate). Phytoplankton groups are initialized as homogeneous fields horizontally and vertically, and allowed to distribute themselves according to the prevailing conditions. Basin-scale model chlorophyll results are in very good agreement with CZCS pigments in virtually every global region. Seasonal variability observed in the CZCS is also well represented in the model. Synoptic scale (100-1000 km) comparisons of imagery are also in good conformance, although occasional departures are apparent. Agreement of nitrate distributions with in situ data is even better, including seasonal dynamics, except for the equatorial Atlantic. The good agreement of the model with satellite and in situ data sources indicates that the model dynamics realistically simulate phytoplankton and nutrient dynamics on synoptic scales. This is especially true given that initial conditions are homogenous chlorophyll fields. The success of the model in producing a reasonable representation of chlorophyll and nutrient distributions and seasonal variability in the global oceans is attributed to the application of a generalized, processes-driven approach as opposed to regional parameterization, and the existence of multiple phytoplankton groups with different physiological and physical properties. These factors enable the model to simultaneously represent the great diversity of physical, biological, chemical, and radiative environments encountered in the global oceans.

Gregg, Watson W.↗

DOE OSTI2021

Are 100 Ensemble Members Enough to Capture the Remote Atmospheric Response to +2°C Arctic Sea Ice Loss?

This study presents results from the Polar Amplification Multimodel Intercomparison Project (PAMIP) single-year time-slice experiments that aim to isolate the atmospheric response to Arctic sea ice loss at global warming levels of +2°C. Using two general circulation models (GCMs), the ensemble size is increased up to 300 ensemble members, beyond the recommended 100 members. After partitioning the response in groups of 100 ensemble members, the reproducibility of the results is evaluated, with a focus on the response of the midlatitude jet streams in the North Atlantic and North Pacific. Both atmosphere-only and coupled ocean–atmosphere PAMIP experiments are analyzed. Substantial differences in the midlatitude response are found among the different experiment subsets, suggesting that 100-member ensembles are still significantly influenced by internal variability, which can mislead conclusions. Despite an overall stronger response, the coupled ocean–atmosphere runs exhibit greater spread due to additional ENSO-related internal variability when the ocean is interactive. The lack of consistency in the response is true for anomalies that are statistically significant according to Student’s t and false discovery rate tests. This is problematic for the multimodel assessment of the response, as some of the spread may be attributed to different model sensitivities whereas it is due to internal variability. Here, we propose a method to overcome this consistency issue that allows for more robust conclusions when only 100 ensemble members are used.

58 GEOSCIENCES↗

NASA NTRS1994

Venus O2 visible and IR nightglow: Implications for lower thermosphere dynamics and chemistry

The National Center for Atmospheric Research thermospheric general circulation model for the Venus thermosphere is modified to examine two observed night airglow features, both of which serve as sensitive tracers of the thermospheric circulation. New O2 nightglow data from the Pioneer Venus Orbiter (PVO) star tracker (O2 Herzberg II at 400-800 nm) and ground-based telescopes (O2 IR at 1.27 microns) yield additional model constraints for estimating Venus winds over 100-130 km. Atomic oxygen, produced by dayside CO2 photolysis peaking near 110 km, and transported to the nightside by the global wind system, is partially destroyed through three-body recombination, yielding the O2 Herzberg II visible nightglow. This emission is very sensitive to horizontal winds at altitudes between 100 and 130 km. Other trace species catalytic reactions also contribute to the production of the very strong nightside infrared (1.27 microns) emission. This paper examines the dynamical and chemical implications of these new data using the Venus thermospheric general circulation model (VTGCM) as an analysis tool. Three-dimensional calculations are presented for both solar maximum and solar medium conditions, corresponding to early PVO (1979-1981) and PVO entry (mid-1992) time periods. Very distinct periods are identified in which zonal winds are alternately weak and strong in the Venus lower thermosphere. VTGCM sensitivity studies are conducted to assess the impacts of potential changes in thermospheric zonal and day-to-night winds, and eddy diffusion on the corresponding nightglow intensities. It appears that cyclostrophic balance extends above 80 km periodically, owing to a reversal of the upper mesosphere latitudinal temperature gradient, and thereby producing strong zonal winds and correspondingly modified O2 nightglow distributions that are observed.

Bougher, S. W.↗

DOE OSTI2021

Using TRMM Latent Heat as a Source to Estimate Convection Induced Gravity Wave Momentum Flux in the Lower Stratosphere

Due to a lack of observations, it is a big challenge to quantify the momentum flux transported from the troposphere into the lower stratosphere by gravity waves generated by convection. This limits our understanding of the dynamics of general circulation in the stratosphere and the realism of simulations of the Quasi Biennial Oscillation. In the past decade, some general circulation models have linked the momentum flux to the latent heating from subgrid scale convective precipitation with some success. However, there is still a large uncertainty in the sources of subgrid scale convection in these models. This study applies the links between the momentum flux and convective latent heating directly to the convective precipitation derived from 16 years of TRMM precipitation radar observations in the tropics and subtropics. The total and directional momentum flux at 100 hPa are derived at individual convective pixels by using the TRMM Spectrum Latent Heating product and large-scale wind profiles from the ERA-Interim dataset. For the first time, we are able to estimate the geographical distribution of momentum flux at 100 hPa in tropics and subtropics from the observed convective sources. The diurnal, seasonal, and interannual variations of the derived momentum flux are presented and discussed. These results could provide a reference to validate the dynamic coupling between the tropical troposphere and the stratosphere in the general circulation models.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

The Role of Radiative Interactions in Tropical Cyclone Development under Realistic Boundary Conditions

Abstract The impact of radiative interactions on tropical cyclone (TC) climatology is investigated using a global, TC-permitting general circulation model (GCM) with realistic boundary conditions. In this model, synoptic-scale radiative interactions are suppressed by overwriting the model-generated atmospheric radiative cooling rates with their monthly varying climatological values. When radiative interactions are suppressed, the global TC frequency is significantly reduced, indicating that radiative interactions are a critical component of TC development even in the presence of spatially varying boundary conditions. The reduced TC activity is primarily due to a decrease in the frequency of pre-TC synoptic disturbances (“seeds”), whereas the likelihood that the seeds undergo cyclogenesis is less affected. When radiative interactions are suppressed, TC genesis shifts toward coastal regions, whereas TC lysis locations stay almost unchanged; together the distance between genesis and lysis is shortened, reducing TC duration. In a warmer climate, the magnitude of TC reduction from suppressing radiative interactions is diminished due to the larger contribution from latent heat release with increased sea surface temperatures. These results highlight the importance of radiative interactions in modulating the frequency and duration of TCs.

Meteorology & Atmospheric Sciences↗

NASA NTRS1991

South Pacific convergence zone and global-scale circulations

The major goals of the research can be classified under three categories: (1) heat and moisture budget studies in the tropics; (2) the role of tropical forcing on subtropical wind maxima in the Southern Hemisphere; and (3) general circulation modeling of low latitude circulations. Research related to the first goal included budget estimates of precipitation, as well as a comparison to rainfall rates observed at numerous island stations. It also included a comparison between observed and satellite estimates of precipitable water. Results and methods used to achieve this goal are given. Included in the second research goal were diagnostic studies of the intraseasonal oscillation which is known to be an important component of organized tropical convection. Also, the role of extratropical forcing on subtropical wind maxima was examined through the use of Eliassen-Palm (EP) fluxes. Again, results and procedural details are discussed in section 3. To accomplish the third research goal, the GLA general circulation model (GCM) was ported to Purdue University and successfully used in a number of control and experimental runs to investigate the significance of tropical heating on subtropical and extratropical circulation features. As with the first two research goals, results and methodology are detailed. Documentation of publications, theses, conference preprint papers, and other scientific information relevant to the research results and significant accomplishments achieved under the contract are given. Finally, the conclusions of the sponsored research are given.

Vincent, Dayton G.↗

NASA NTRS2010

Forced and Unforced Variability of Twentieth Century North American Droughts and Pluvials

Research on the forcing of drought and pluvial events over North America is dominated by general circulation model experiments that often have operational limitations (e.g., computational expense, ability to simulate relevant processes, etc). We use a statistically based modeling approach to investigate sea surface temperature (SST) forcing of the twentieth century pluvial (1905-1917) and drought (1932-1939, 1948-1957, 1998-2002) events. A principal component (PC) analysis of Palmer Drought Severity Index (PDSI) from the North American Drought Atlas separates the drought variability into five leading modes accounting for 62% of the underlying variance. Over the full period spanning these events (1900-2005), the first three PCs significantly correlate with SSTs in the equatorial Pacific (PC 1), North Pacific (PC 2), and North Atlantic (PC 3), with spatial patterns (as defined by the empirical orthogonal functions) consistent with our understanding of North American drought responses to SST forcing. We use a large ensemble statistical modeling approach to determine how successfully we can reproduce these drought/pluvial events using these three modes of variability. Using Pacific forcing only (PCs 1-2), we are able to reproduce the 1948-1957 drought and 1905-1917 pluvial above a 95% random noise threshold in over 90% of the ensemble members; the addition of Atlantic forcing (PCs 1-2-3) provides only marginal improvement. For the 1998-2002 drought, Pacific forcing reproduces the drought above noise in over 65% of the ensemble members, with the addition of Atlantic forcing increasing the number passing to over 80%. The severity of the drought, however, is underestimated in the ensemble median, suggesting this drought intensity can only be achieved through internal variability or other processes. Pacific only forcing does a poor job of reproducing the 1932-1939 drought pattern in the ensemble median, and less than one third of ensemble members exceed the noise threshold (28%). Inclusion of Atlantic forcing improves the ensemble median drought pattern and nearly doubles the number of ensemble members passing the noise threshold (52%). Even with the inclusion of Atlantic forcing, the intensity of the simulated 1932-1939 drought is muted, and the drought itself extends too far into the southwest and southern Great Plains. To an even greater extent than the 1998-2002 drought, these results suggest much of the variance in the 1932-1939 drought is dependent on processes other than SST forcing. This study highlights the importance of internal noise and non SST processes for hydroclimatic variability over North America, complementing existing research using general circulation models.

Cook, Benjamin I.↗