Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “spatial pattern”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 1,063 records · Page 59
All sourcesNASA NTRSDOE OSTI
NASA NTRS2002

Configuration and Intraseasonal Duration of Interannual Anomalies of the Great Plains Low-Level Jet

Despite the fact that the low-level jet of the southern Great Plains (the GPLLJ) of the U.S. is primarily a nocturnal phenomenon that virtually vanishes during the daylight hours, it is one of the most persistent and stable climatological features of the low-level continental flow during the warm-season months, May through August. We have used significant-level data to validate the skill of the GEOS-1 Data Assimilation System (DAS) in realistically detecting this jet and inferring its structure and evolution. We have then carried out a 15-year reanalysis with the GEOS-1 DAS to determine its climatology and mean diurnal cycle and to study its interannual variability. Interannual anomalies of the meridional flow associated with the GPLLJ are much smaller than the mean diurnal fluctuations, than random intraseasonal anomalies, and than the mean wind itself. There are three maxima of low-level meridional flow variance over the Great Plains and the Gulf of Mexico: a 1.2 m2 s-2 peak over the southeast Texas, to the east and south of the mean velocity peak, a 1.0 m2 s-2 peak over the western Gulf of Mexico, and a .8 m2 s-2 peak over the upper Great Plains (UGP), near the Nebraska/South Dakota border. Each of the three variance maxima corresponds to a spatially coherent, jet-like pattern of low-level flow interannual variability. There are also three dominant modes of interannual variability corresponding to the three variance maxima, but not in a simple one-to-one relationship. Cross-sectional profiles of mean southerly wind over Texas remain relatively stable and recognizable from year to year with only its eastward flank showing significant variability. This variability, however, exhibits a distinct, biennial oscillation during the first six to seven years of the reanalysis period and only then. This intermittent biennial oscillation (IBO, one of the three modes discussed in the previous paragraph) in the lowlevel flow is restricted to the region surrounding eastern Texas and is also evident in the NCEP/NCAR reanalysis data set from about 1978 to 1985 or 1986 and again from 1995 to 2000. It is evident as well in surface pressure in both the GEOS-1 and NCEP/NCAR sets. The interannual anomalies do not necessarily persist uniformly throughout an entire season, but can fluctuate from one part of the season to the next. To estimate the characteristic sub-seasonal time scales for coherence of these fluctuations, we have taken the weekly anomaly of low-level wind at each point of the domain from the climatological average for that given point and that given week of the season and computed the covariance of its fluctuations over all weeks and over all years with the weekly climatological anomaly of the meridional wind at each of the three reference points discussed above. The typical duration of a coherent interannual anomaly within a given warm season increases with decreasing latitude from 2 to 3 weeks over the UGP, to 6 to 7 weeks over eastern Texas. Coherence over the western Gulf of Mexico is intermediate between the two with a typical duration of 4 to 5 weeks. There appears to be evidence that the interannual anomalies over Texas the Gulf propagate to the UGP after a week and those over the Gulf propagate there after 2 to 3 weeks. There also appears to be some reverse propagation of interannual anomalies over the UGP to Texas and to the Gulf after a period of about one week. The interannual anomalies in southerly flow over eastern Texas seem to correlate well with interannual anomalies of surface temperature and (negative) ground wetness and over western Texas.

Helfand, H. M.↗
NASA NTRS2007

[Overview of the GSFC Effort to Build Microwave Detectors for CMB Polarization]

Filled arrays of bolometers are currently being employed for use in astronomy from the far-infrared through millimeter portion of the electromagnetic spectrum. Because of the large range of wavelengths for which such detectors are applicable, the size of a pixel (p) relative to the wavelength (lambda) will vary according to the specific application of a given available technology. As the pixel size becomes comparable to a wavelength of the absorbed radiation, correlations in the radiation field at the pixel can affect the imaging properties of the system independent of the correlations introduced by the front-end optics. We study the dependence of image fidelity and induced polarization on the size of the pixel by employing a formalism in which diffraction due to the pixel boundary is treated by propagating the second-order statistical correlations of the radiation field through a model optical system. We construct simulated polarized images of square pixels for various ratios of p lambda. For the limit in which few modes are supported by the pixel (p lambda lesssiml), we find that the diffraction due to the pixel edges is non-negligible and hence must be considered along with the telescope diffraction pattern in modeling the ultimate spatial resolution of an imaging system. For the case in which the pixel is over-moded (p lambda ggl), the geometric limit is approached as expected. This technique provides a quantitative approach to understand and optimize the imaging properties of planar detectors in the few-mode limit.

Chuss, David↗
NASA NTRS2008

Interferometric Rayleigh Scattering Measurement System

A method and apparatus for performing simultaneous multi-point measurements of multiple velocity components in a gas flow is described. Pulses of laser light are directed to a measurement region of unseeded gas to produce Rayleigh or Mie scattered light in a plurality of directions. The Rayleigh or Mie scattered light is collected from multiple directions and combined in a single collimated light beam. The Rayleigh or Mie scattered light is then mixed together with a reference laser light before it is passed through a single planar Fabry-Perot interferometer for spectral analysis. At the output of the interferometer, a high-sensitivity CCD camera images the interference fringe pattern. This pattern contains the spectral and spatial information from both the Rayleigh scattered light and the reference laser light. Interferogram processing software extracts and analyzes spectral profiles to determine the velocity components of the gas flow at multiple points in the measurement region. The Rayleigh light rejected by the interferometer is recirculated to increase the accuracy and the applicability of the method for measurements at high temperatures without requiring an increase in the laser energy.

Bivolaru, Daniel↗
NASA NTRS2017

High Sensitivity of Gross Primary Production in the Rocky Mountains to Summer Rain

n the catchments of the Rocky Mountains, peak snowpack is declining in response to warmer spring temperatures. To understand how this will influence terrestrial gross primary production (GPP), we compared precipitation data across the intermountain west with satellite retrievals of solar-induced fluorescence (SIF), a proxy for GPP. Annual precipitation patterns explained most of the spatial and temporal variability of SIF, but the slope of the response was dependent on site to site differences in the proportion of snowpack to summer rain. We separated the response of SIF to different seasonal precipitation amounts and found that SIF was approximately twice as sensitive to variations in summer rain than snowpack. The response of peak GPP to a secular decline in snowpack will likely be subtle, whereas a change in summer rain amount will have precipitous effects on GPP. The study suggests that the rain use efficiency of Rocky Mountain ecosystems is strongly dependent on precipitation form and timing.

solar-induced 15 Fuorescence (SIF)↗
NASA NTRS2021

Global irrigation contribution to wheat and maize yield

Irrigation is the largest sector of human water use and an important option for increasing crop production and reducing drought impacts. However, the potential for irrigation to contribute to global crop yields remains uncertain. Here, we quantify this contribution for wheat and maize at global scale by developing a Bayesian framework integrating empirical estimates and gridded global crop models on new maps of the relative difference between attainable rainfed and irrigated yield (ΔY). At global scale, ΔY is 34 ± 9% for wheat and 22 ± 13% for maize, with large spatial differences driven more by patterns of precipitation than that of evaporative demand. Comparing irrigation demands with renewable water supply, we find 30–47% of contemporary rainfed agriculture of wheat and maize cannot achieve yield gap closure utilizing current river discharge, unless more water diversion projects are set in place, putting into question the potential of irrigation to mitigate climate change impacts.

Irrigation↗
NASA NTRSDate not supplied

PM 2.5 Concentrations over Major Metropolitan Regions Inferred from Airborne High Spectral Resolution Lidar Measurements Using Machine Learning Regression

We use measurements of near-surface aerosol backscatter, extinction, and depolarization acquired by four NASA Langley Research Center airborne High Spectral Resolution Lidars (HSRLs) to develop a machine learning regression methodology to infer PM2.5 concentrations at the surface and aloft. These airborne HSRL measurements were acquired over major metropolitan regions in the United States and Asia during more than 170 flights since 2010. Hourly surface PM2.5 measurements from the EPA air quality system and similar networks in other countries acquired within 10 km and 15 minutes of these near-surface HSRL measurements are used to train models that compute PM2.5 concentrations from the HSRL measurements. We examine several regression methods and find that exponential Gaussian Process algorithms consistently give the best performance in terms of the lowest root-mean-square (RMS) errors and the highest correlations. Model performance varies significantly depending on various combinations of HSRL aerosol measurements (e.g., aerosol backscatter, extinction, depolarization, backscatter color ratios, lidar ratios, aerosol optical thickness) and retrievals (e.g., mixed layer height, aerosol type) used in the regressions. Models that use near-surface measurements of aerosol backscatter and aerosol intensive properties such as depolarization, backscatter color ratio, and lidar ratio typically give the best performance with RMS errors around 4 mg/m3 and correlation coefficients above 0.9. HSRL measurements were often acquired when the aircraft flew systematic “raster-scan” patterns for several hours over these cities. These flight patterns enabled measurements of the spatial, temporal, and vertical variabilities in the distributions of aerosol backscatter and aerosol intensive properties and allowed us to derive the corresponding variabilities in PM2.5 concentrations. We present examples of such variabilities over urban areas in the United States as well as Asia. We describe also how the distribution of surface PM2.5 varies with aerosol type and use these retrievals to examine model simulations of surface PM2.5 in these metropolitan regions. We also discuss how this methodology may be applied to measurements from satellite lidars such as CALIOP on CALIPSO and ATLID on EarthCARE.

lidar↗
DOE OSTI2023

A Bayesian model for multivariate discrete data using spatial and expert information with application to inferring building attributes

When modeling sparsely observed multivariate data, strong prior information elicited from experts can be used to bolster predictive accuracy and counteract sampling bias. Similarly, modeling autocorrelation in space can help make use of co-occurrence patterns present in many types of spatial data. To make use of both expert prior information and spatial structure, we propose a novel graphical model for a spatial Bayesian network developed specifically to address challenges in inferring the attributes of buildings from geographically sparse observational data. This model is implemented as the sum of a spatial multivariate Gaussian random field and a tabular conditional probability function in real-valued space prior to projection onto the probability simplex. This modeling form is especially suitable for the usage of prior information in the form of sets of atomic rules obtained from experts. To perform inference with missing data, we implement a Markov chain Monte Carlo scheme composed of alternating steps of Gibbs sampling of missing entries and Hamiltonian Monte Carlo for model parameters. A case study in building attribution is presented to highlight the advantages and limitations of this approach.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1999

Visualization of Subsurface Defects in Composites using a Focal Plane Array Infrared Camera

A technique for enhanced defect visualization in composites via transient thermography is presented in this paper. The effort targets automated defect map construction for multiple defects located in the observed area. Experimental data were collected on composite panels of different thickness with square inclusions and flat bottom holes of different depth and orientation. The time evolution of the thermal response and spatial thermal profiles are analyzed. The pattern generated by carbon fibers and the vignetting effect of the focal plane array camera make defect visualization difficult. An improvement of the defect visibility is made by the pulse phase technique and the spatial background treatment. The relationship between a size of a defect and its reconstructed image is analyzed as well. The image processing technique for noise reduction is discussed.

Plotnikov, Yuri A.↗
NASA NTRS2020

“It’s Raining Bits”: Patterns in Directional Precipitation Persistence Across the United States

The spatial and temporal ordering of precipitation occurrence impacts ecosystems, streamflow, and water availability. For example, both large-scale climate patterns and local landscapes drive weather events, and the typical speeds and directions of these events moving across a basin dictate the timing of flows at its outlet. We address the predictability of precipitation occurrence at a given location, based on the knowledge of past precipitation at surrounding locations. We identify ‘‘dominant directions of precipitation influence’’ across the continental United States based on a gridded daily dataset. Specifically, we apply information theory–based measures that characterize dominant directions and strengths of spatial and temporal precipitation dependencies. On a national average, this dominant direction agrees with the prevalent direction of weather movement from west to east across the country, but regional differences reflect topographic divides, precipitation gradients, and different climatic drivers of precipitation. Trends in these information relationships and their correlations with climate indices over the past 70 years also show seasonal and spatial divides. This study expands upon a framework of information-based predictability to answer questions about spatial connectivity in addition to temporal persistence. The methods presented here are generally useful to understand many aspects of weather and climate variability.

North America↗
NASA NTRS1973

Automated thematic mapping and change detection of ERTS-A images

The author has identified the following significant results. For the recognition of terrain types, spatial signatures are developed from the diffraction patterns of small areas of ERTS-1 images. This knowledge is exploited for the measurements of a small number of meaningful spatial features from the digital Fourier transforms of ERTS-1 image cells containing 32 x 32 picture elements. Using these spatial features and a heuristic algorithm, the terrain types in the vicinity of Phoenix, Arizona were recognized by the computer with a high accuracy. Then, the spatial features were combined with spectral features and using the maximum likelihood criterion the recognition accuracy of terrain types increased substantially. It was determined that the recognition accuracy with the maximum likelihood criterion depends on the statistics of the feature vectors. Nonlinear transformations of the feature vectors are required so that the terrain class statistics become approximately Gaussian. It was also determined that for a given geographic area the statistics of the classes remain invariable for a period of a month but vary substantially between seasons.

Gramenopoulos, N.↗
DOE OSTI2020

Pattern formation in a coupled membrane-bulk reaction-diffusion model for intracellular polarization and oscillations

Reaction-diffusion systems have been widely used to study spatio-temporal phenomena in cell biology, such as cell polarization. Coupled bulk-surface models naturally include compartmentalization of cytosolic and membrane-bound polarity molecules. Here we study the distribution of the polarity protein Cdc42 in a mass-conserved membrane-bulk model, and explore the effects of diffusion and spatial dimensionality on spatio-temporal pattern formation. We first analyze a one-dimensional (1-D) model for Cdc42 oscillations in fission yeast, consisting of two diffusion equations in the bulk domain coupled to nonlinear ODEs for binding kinetics at each end of the cell. In 1-D, our analysis reveals the existence of symmetric and asymmetric steady states, as well as anti-phase relaxation oscillations typical of slow-fast systems. We then extend our analysis to a two-dimensional (2-D) model with circular bulk geometry, for which species can either diffuse inside the cell or become bound to the membrane and undergo a nonlinear reaction-diffusion process. We also consider a nonlocal system of PDEs approximating the dynamics of the 2-D membrane-bulk model in the limit of fast bulk diffusion. In all three model variants we find that mass conservation selects perturbations of spatial modes that simply redistribute mass. In 1-D, only anti-phase oscillations between the two ends of the cell can occur, and in-phase oscillations are excluded. In higher dimensions, no radially symmetric oscillations are observed. Instead, the only instabilities are symmetry-breaking, either corresponding to stationary Turing instabilities, leading to the formation of stationary patterns, or to oscillatory Turing instabilities, leading to traveling and standing waves. Codimension-two Bogdanov—Takens bifurcations occur when the two distinct instabilities coincide, causing traveling waves to slow down and to eventually become stationary patterns. Our work clarifies the effect of geometry and dimensionality on behaviors observed in ma.ss-conserved cell polarity models.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1981

Interactive initialization of heat flux parameters for numerical models using satellite temperature measurements

Day/night Heat Capacity Mapping Mission image pairs over Kansas and Indiana were used to examine the spatial variation of moisture availability on the mesoscale. Of particular concern was whether patterns of moisture availability and thermal inertia, as well as surface heat fluxes, respond to significant spatial variations in the rainfall pattern. A temperature analysis for the Indiana case is presented.

Carlson, T. N.↗
NASA NTRS1998

Multi-Scale Fractal Analysis of Image Texture and Pattern

Fractals embody important ideas of self-similarity, in which the spatial behavior or appearance of a system is largely independent of scale. Self-similarity is defined as a property of curves or surfaces where each part is indistinguishable from the whole, or where the form of the curve or surface is invariant with respect to scale. An ideal fractal (or monofractal) curve or surface has a constant dimension over all scales, although it may not be an integer value. This is in contrast to Euclidean or topological dimensions, where discrete one, two, and three dimensions describe curves, planes, and volumes. Theoretically, if the digital numbers of a remotely sensed image resemble an ideal fractal surface, then due to the self-similarity property, the fractal dimension of the image will not vary with scale and resolution. However, most geographical phenomena are not strictly self-similar at all scales, but they can often be modeled by a stochastic fractal in which the scaling and self-similarity properties of the fractal have inexact patterns that can be described by statistics. Stochastic fractal sets relax the monofractal self-similarity assumption and measure many scales and resolutions in order to represent the varying form of a phenomenon as a function of local variables across space. In image interpretation, pattern is defined as the overall spatial form of related features, and the repetition of certain forms is a characteristic pattern found in many cultural objects and some natural features. Texture is the visual impression of coarseness or smoothness caused by the variability or uniformity of image tone or color. A potential use of fractals concerns the analysis of image texture. In these situations it is commonly observed that the degree of roughness or inexactness in an image or surface is a function of scale and not of experimental technique. The fractal dimension of remote sensing data could yield quantitative insight on the spatial complexity and information content contained within these data. A software package known as the Image Characterization and Modeling System (ICAMS) was used to explore how fractal dimension is related to surface texture and pattern. The ICAMS software was verified using simulated images of ideal fractal surfaces with specified dimensions. The fractal dimension for areas of homogeneous land cover in the vicinity of Huntsville, Alabama was measured to investigate the relationship between texture and resolution for different land covers.

Emerson, Charles W.↗
DOE OSTI2021

Multi-atom quasiparticle scattering interference for superconductor energy-gap symmetry determination

Complete theoretical understanding of the most complex superconductors requires a detailed knowledge of the symmetry of the superconducting energy-gap Δ$^α_k$, for all momenta k on the Fermi surface of every band α. While there are a variety of techniques for determining |Δ$^α_k$|, no general method existed to measure the signed values of Δ$^α_k$. Recently, however, a technique based on phase-resolved visualization of superconducting quasiparticle interference (QPI) patterns, centered on a single non-magnetic impurity atom, was introduced. In principle, energy-resolved and phase-resolved Fourier analysis of these images identifies wavevectors connecting all k-space regions where Δ$^α_k$ has the same or opposite sign. But use of a single isolated impurity atom, from whose precise location the spatial phase of the scattering interference pattern must be measured, is technically difficult. Here we introduce a generalization of this approach for use with multiple impurity atoms, and demonstrate its validity by comparing the Δ$^α_k$ it generates to the Δ$^α_k$ determined from single-atom scattering in FeSe where s ± energy-gap symmetry is established. Finally, to exemplify utility, we use the multi-atom technique on LiFeAs and find scattering interference between the hole-like and electron-like pockets as predicted for Δ$^α_k$ of opposite sign.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1998

Visually induced adaptation in three-dimensional organization of primate vestibuloocular reflex

The adaptive plasticity of the spatial organization of the vestibuloocular reflex (VOR) has been investigated in intact and canal-plugged primates using 2-h exposure to conflicting visual (optokinetic, OKN) and vestibular rotational stimuli about mutually orthogonal axes (generating torsional VOR + vertical OKN, torsional VOR + horizontal OKN, vertical VOR + horizontal OKN, and horizontal VOR + vertical OKN). Adaptation protocols with 0.5-Hz (+/-18 degrees ) head movements about either an earth-vertical or an earth-horizontal axis induced orthogonal response components as high as 40-70% of those required for ideal adaptation. Orthogonal response gains were highest at the adapting frequency with phase leads present at lower and phase lags present at higher frequencies. Furthermore, the time course of adaptation, as well as orthogonal response dynamics were similar and relatively independent of the particular visual/vestibular stimulus combination. Low-frequency (0. 05 Hz, vestibular stimulus: +/-60 degrees ; optokinetic stimulus: +/-180 degrees ) adaptation protocols with head movements about an earth-vertical axis induced smaller orthogonal response components that did not exceed 20-40% of the head velocity stimulus (i.e., approximately 10% of that required for ideal adaptation). At the same frequency, adaptation with head movements about an earth-horizontal axis generated large orthogonal responses that reached values as high as 100-120% of head velocity after 2 h of adaptation (i.e., approximately 40% of ideal adaptation gains). The particular spatial and temporal response characteristics after low-frequency, earth-horizontal axis adaptation in both intact and canal-plugged animals strongly suggests that the orienting (and perhaps translational) but not inertial (velocity storage) components of the primate otolith-ocular system exhibit spatial adaptability. Due to the particular nested arrangement of the visual and vestibular stimuli, the optic flow pattern exhibited a significant component about the third spatial axis (i.e., orthogonal to the axes of rotation of the head and visual surround) at twice the oscillation frequency. Accordingly, the adapted VOR was characterized consistently by a third response component (orthogonal to both the axes of head and optokinetic drum rotation) at twice the oscillation frequency after earth-horizontal but not after earth-vertical axis 0.05-Hz adaptation. This suggests that the otolith-ocular (but not the semicircular canal-ocular) system can adaptively change its spatial organization at frequencies different from those of the head movement.

Non-NASA Center↗
NASA NTRS2017

Classifying Multi-Model Wheat Yield Impact Response Surfaces Showing Sensitivity to Temperature and Precipitation Change

Crop growth simulation models can differ greatly in their treatment of key processes and hence in their response to environmental conditions. Here, we used an ensemble of 26 process-based wheat models applied at sites across a European transect to compare their sensitivity to changes in temperature (minus 2 to plus 9 degrees Centigrade) and precipitation (minus 50 to plus 50 percent). Model results were analysed by plotting them as impact response surfaces (IRSs), classifying the IRS patterns of individual model simulations, describing these classes and analysing factors that may explain the major differences in model responses. The model ensemble was used to simulate yields of winter and spring wheat at four sites in Finland, Germany and Spain. Results were plotted as IRSs that show changes in yields relative to the baseline with respect to temperature and precipitation. IRSs of 30-year means and selected extreme years were classified using two approaches describing their pattern. The expert diagnostic approach (EDA) combines two aspects of IRS patterns: location of the maximum yield (nine classes) and strength of the yield response with respect to climate (four classes), resulting in a total of 36 combined classes defined using criteria pre-specified by experts. The statistical diagnostic approach (SDA) groups IRSs by comparing their pattern and magnitude, without attempting to interpret these features. It applies a hierarchical clustering method, grouping response patterns using a distance metric that combines the spatial correlation and Euclidian distance between IRS pairs. The two approaches were used to investigate whether different patterns of yield response could be related to different properties of the crop models, specifically their genealogy, calibration and process description. Although no single model property across a large model ensemble was found to explain the integrated yield response to temperature and precipitation perturbations, the application of the EDA and SDA approaches revealed their capability to distinguish: (i) stronger yield responses to precipitation for winter wheat than spring wheat; (ii) differing strengths of response to climate changes for years with anomalous weather conditions compared to period-average conditions; (iii) the influence of site conditions on yield patterns; (iv) similarities in IRS patterns among models with related genealogy; (v) similarities in IRS patterns for models with simpler process descriptions of root growth and water uptake compared to those with more complex descriptions; and (vi) a closer correspondence of IRS patterns in models using partitioning schemes to represent yield formation than in those using a harvest index. Such results can inform future crop modelling studies that seek to exploit the diversity of multi-model ensembles, by distinguishing ensemble members that span a wide range of responses as well as those that display implausible behaviour or strong mutual similarities.

Classification; Climate change; Crop model; Ensemb↗
DOE OSTI2024

Quantifying Dissolution Dynamics in Porous Media Using a Spatial Flow Focusing Profile

Abstract The diverse range of patterns in porous media formed by dissolution processes depends on the relative magnitude of flow, transport, and chemical reactions at pore surfaces. However, distinguishing between regimes often relies solely on qualitative, visual comparisons of emergent structures. Here, we propose a quantitative measure capable of identifying different regimes using the concept of the spatial flow focusing profile, which segments the medium into cross sections along the flow direction to calculate the flow focusing index for each section. We employ this measure in numerical simulations of a dissolving porous medium using a pore network model. We obtain a morphological phase diagram of dissolution patterns, which we characterize using the flow focusing profile. In particular, we demonstrate that analyzing the temporal changes in the profile allows one to quantitatively distinguish between wormholing and channeling. The transition between them is shown to be affected by the heterogeneity of the system.

58 GEOSCIENCES↗
NASA NTRS2015

Relative Contributions of Mean-State Shifts and ENSO-Driven Variability to Precipitation Changes in a Warming Climate

The El Nino-Southern Oscillation (ENSO) is an important driver of regional hydroclimate variability through far-reaching teleconnections. This study uses simulations performed with Coupled General Circulation Models (CGCMs) to investigate how regional precipitation in the 21st century may be affected by changes in both ENSO-driven precipitation variability and slowly-evolving mean rainfall. First, a dominant, time-invariant pattern of canonical ENSO variability (cENSO) is identified in observed SST data. Next, the fidelity with which 33 state-of-the-art CGCMs represent the spatial structure and temporal variability of this pattern (as well as its associated precipitation responses) is evaluated in simulations of 20th century climate change. Possible changes in both the temporal variability of this pattern and its associated precipitation teleconnections are investigated in 21st century climate projections. Models with better representation of the observed structure of the cENSO pattern produce winter rainfall teleconnection patterns that are in better accord with 20th century observations and more stationary during the 21st century. Finally, the model-predicted 21st century rainfall response to cENSO is decomposed into the sum of three terms: 1) the 21st century change in the mean state of precipitation; 2) the historical precipitation response to the cENSO pattern; and 3) a future enhancement in the rainfall response to cENSO, which amplifies rainfall extremes. By examining the three terms jointly, this conceptual framework allows the identification of regions likely to experience future rainfall anomalies that are without precedent in the current climate.

Celine J W Bonfils↗