Search NASA⌕ Search

SEARCH · Search NASA

Results for “empirical orthogonal functions”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
NASA NTRS2008

Mechanisms of Interannual Variations of the Meridional Overturning Circulation of the North Atlantic Ocean

The authors investigate the nature of the interannual variability of the meridional overturning circulation (MOC) of the North Atlantic Ocean using an Estimating the Circulation and Climate of the Ocean (ECCO) assimilation product for the period of 1993-2003. The time series of the first empirical orthogonal function of the MOC is found to be correlated with the North Atlantic Oscillation (NAO) index, while the associated circulation anomalies correspond to cells extending over the full ocean depth. Model sensitivity experiments suggest that the wind is responsible for most of this interannual variability, at least south of 40(deg)N. A dynamical decomposition of the meridional streamfunction allows a further look into the mechanisms. In particular, the contributions associated with 1) the Ekman flow and its depth-independent compensation, 2) the vertical shear flow, and 3) the barotropic gyre flowing over zonally varying topography are examined. Ekman processes are found to dominate the shorter time scales (1.5-3 yr), while for longer time scales (3-10 yr) the MOC variations associated with vertical shear flow are of greater importance. The latter is primarily caused by heaving of the pycnocline in the western subtropics associated with the stronger wind forcing. Finally, how these changes in the MOC affect the meridional heat transport (MHT) is examined. It is found that overall, Ekman processes explain a larger part of interannual variability (3-10 yr) for MHT (57%) than for the MOC (33%).

meridional heat transport↗
NASA NTRS2007

TRMM-TMI Satellite Observed Soil Moisture and Vegetation Density (1998-2005) Show Strong Connection with El Nino in Eastern Australia

Spatiotemporal patterns in soil moisture and vegetation water content across mainland Australia were investigated from 1998 through 2005, using TRMMITMI passive microwave observations. The Empirical Orthogonal Function technique was used to extract dominant spatial and temporal patterns in retrieved estimates of moisture content for the top 1-cm of soil (theta) and vegetation moisture content (via optical depth tau). The dominant temporal theta and tau patterns were strongly correlated to El Nino/Southern Oscillation (ENSO) in spring (3 = 0.90), and to a progressively lesser extent autumn, summer and winter. The Indian Ocean Dipole (IOD) index also explained part of the variation in spring 8 and z. Cluster analysis suggested that the regions most affected by ENS0 are mainly located in eastern Australia. The results suggest that the drought conditions experienced in eastern Australia since 2000 an clearly expressed in these satellite observations have a strong connection with ENSO patterns.

Liu, Yi↗
NASA NTRS2010

Warm Season Subseasonal Variability and Climate Extremes in the Northern Hemisphere: The Role of Stationary Rossby Waves

This study examines the nature of boreal summer subseasonal atmospheric variability based on the new NASA Modern-Era Retrospective analysis for Research and Applications (MERRA) for the period 1979-2010. An analysis of the June, July and August subseasonal 250hPa v-wind anomalies shows distinct Rossby wave-like structures that appear to be guided by the mean jets. On monthly subseasonal time scales, the leading waves (the first 10 rotated empirical orthogonal functions or REOFs of the 250hPa v-wind) explain about 50% of the Northern Hemisphere vwind variability, and account for more than 30% (60%) of the precipitation (surface temperature) variability over a number of regions of the northern middle and high latitudes, including the U.S. northern Great Plains, parts of Canada, Europe, and Russia. The first REOF in particular, consists of a Rossby wave that extends across northern Eurasia where it is a dominant contributor to monthly surface temperature and precipitation variability, and played an important role in the 2003 European and 2010 Russian heat waves. While primarily subseasonal in nature, the Rossby waves can at times have a substantial seasonal mean component. This is exemplified by REOF 4 which played a major role in the development of the most intense anomalies of the U.S. 1988 drought (during June) and the 1993 flooding (during July), though differed in the latter event by also making an important contribution to the seasonal mean anomalies. A stationary wave model (SWM) is used to reproduce some of the basic features of the observed waves and provide insight into the nature of the forcing. In particular, the responses to a set of idealized forcing functions are used to map the optimal forcing patterns of the leading waves. Also, experiments to reproduce the observed waves with the SWM using MERRA-based estimates of the forcing indicate that the wave forcing is dominated by sub-monthly vorticity transients.

Schubert, Siegfried↗
NASA NTRS2011

Hyperspectral Remote Sensing of Atmosphere and Surface Properties

Atmospheric Infrared Sounder (AIRS), Infrared Atmospheric Sounding Interferometer (IASI), and Cross-track Infrared Sounder (CrIS) are all hyper-spectral satellite sensors with thousands of spectral channels. Top of atmospheric radiance spectra measured by these sensors contain high information content on atmospheric, cloud, and surface properties. Exploring high information content contained in these high spectral resolution spectra is a challenging task due to computation e ort involved in modeling thousands of spectral channels. Usually, only very small fractions (4{10 percent) of the available channels are included in physical retrieval systems or numerical weather forecast (NWP) satellite data assimilations. We will describe a method of simultaneously retrieving atmospheric temperature, moisture, cloud, and surface properties using all available spectral channels without sacrificing computational speed. The essence of the method is to convert channel radiance spectra into super-channels by an Empirical Orthogonal Function (EOF) transformation. Because the EOFs are orthogonal to each other, about 100 super-channels are adequate to capture the information content of the radiance spectra. A Principal Component-based Radiative Transfer Model (PCRTM) developed at NASA Langley Research Center is used to calculate both the super-channel magnitudes and derivatives with respect to atmospheric profiles and other properties. There is no need to perform EOF transformations to convert super channels back to spectral space at each iteration step for a one-dimensional variational retrieval or a NWP data assimilation system. The PCRTM forward model is also capable of calculating radiative contributions due to multiple-layer clouds. The multiple scattering effects of the clouds are efficiently parameterized. A physical retrieval algorithm then performs an inversion of atmospheric, cloud, and surface properties in super channel domain directly therefore both reducing the computational need and preserving the information content of the IASI measurements. The inversion algorithm is based on a non-linear Levenberg-Marquardt method with climatology covariance matrices and a priori information as constraints. One advantage of this approach is that it uses all information content from the hyper-spectral data so that the retrieval is less sensitive to instrument noise and eliminates the need for selecting a sub-set of the channels.

Liu, Xu↗
NASA NTRS2011

Annual Cycle of Surface Longwave Radiation

The annual cycles of upward and downward longwave fluxes at the Earth s surface are investigated by use of the NASA/GEWEX Surface Radiation Budget Data Set. Because of the immense difference between the heat capacity of land and ocean, the surface of Earth is partitioned into these two categories. Principal component analysis is used to quantify the annual cycles. Over land, the first principal component describes over 95% of the variance of the annual cycle of the upward and downward longwave fluxes. Over ocean the first term describes more than 87% of these annual cycles. Empirical orthogonal functions show the corresponding geographical distributions of these cycles. Phase plane diagrams of the annual cycles of upward longwave fluxes as a function of net shortwave flux show the thermal inertia of land and ocean.

Mlynczak, Pamela E.↗
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.↗
NASA NTRS2011

Experiments in Reconstructing Twentieth-Century Sea Levels

One approach to reconstructing historical sea level from the relatively sparse tide-gauge network is to employ Empirical Orthogonal Functions (EOFs) as interpolatory spatial basis functions. The EOFs are determined from independent global data, generally sea-surface heights from either satellite altimetry or a numerical ocean model. The problem is revisited here for sea level since 1900. A new approach to handling the tide-gauge datum problem by direct solution offers possible advantages over the method of integrating sea-level differences, with the potential of eventually adjusting datums into the global terrestrial reference frame. The resulting time series of global mean sea levels appears fairly insensitive to the adopted set of EOFs. In contrast, charts of regional sea level anomalies and trends are very sensitive to the adopted set of EOFs, especially for the sparser network of gauges in the early 20th century. The reconstructions appear especially suspect before 1950 in the tropical Pacific. While this limits some applications of the sea-level reconstructions, the sensitivity does appear adequately captured by formal uncertainties. All our solutions show regional trends over the past five decades to be fairly uniform throughout the global ocean, in contrast to trends observed over the shorter altimeter era. Consistent with several previous estimates, the global sea-level rise since 1900 is 1.70 +/- 0.26 mm/yr. The global trend since 1995 exceeds 3 mm/yr which is consistent with altimeter measurements, but this large trend was possibly also reached between 1935 and 1950.

Ray, Richard D.↗
NASA NTRS2005

Observed Land Impacts on Clouds, Water Vapor, and Rainfall at Continental Scales

How do the continents affect large-scale hydrological cycles? How important can one continent be to the climate system? To address these questions, five years of National Aeronautics and Space Administration (NASA) Terra Moderate Resolution Imaging Spectroradiometer (MODIS) observations, Tropical Rainfall Measuring Mission (TRMM) observations, and the Global Precipitation Climatology Project (GPCP) global precipitation analysis, were used to assess the land impacts on clouds, ~ainfall, and water vapor at continental scales. At these scales, Empirical Orthogonal Function (EOF) and continentally averaged analyses illustrate that continents as integrated regions enhance the seasonality of atmospheric and surface hydrological parameters. Specifically, the continents of Eurasia and North America enhance the seasonality of cloud optical thickness, cirrus fraction, rainfall, and water vapor. Over land, both liquid water and ice cloud effective radii are smaller than over oceans, primarily because land has more aerosol particles. In addition, different continents have similar impacts on hydrological variables in terms of seasonality, but differ in magnitude. For example, in winter, North America and Eurasia increase cloud optical thickness to 17.5 and 16, respectively, while in summer, Eurasia has much smaller cloud optical thicknesses than North An1erica. Such different land impacts are determined by each continent's geographical condition, land cover, and land use. These new understandings help further address the land-ocean contrasts on global climate, help validate global climate model simulated land-atmosphere interactions, and shed light on interpreting climate change over land.

Jin, Menglin↗
NASA NTRS2012

A Technique Using Principal Component Analysis to Compare Seasonal Cycles of Earth Radiation from CERES and Model Computations

A method for quantitatively comparing the seasonal cycles of two global data sets is presented. The seasonal cycles of absorbed solar radiation (ASR) and outgoing longwave radiation (OLR) have been computed from an eight-year data set from the Clouds and Earth's Radiant Energy System (CERES) scanning radiometers and from a model data set produced by the NASA Goddard Space Flight Center Global Modeling and Assimilation Office. To compare the seasonal cycles from these two data sets, principal component (PC) analysis is used, where the PCs express the time variations and the corresponding empirical orthogonal functions (EOFs) describe the geographic variations. Ocean has a long thermal response time compared to land, so land and ocean are separated for the analysis. The root-mean square values for the seasonal cycles of ASR and OLR are extremely close for the two data sets. The first three PCs are quite close, showing that the time responses and magnitudes over the globe are very similar. The agreement between the two sets of PCs is quantified by computing the matrix of inner products of the two sets. For ASR over land, the first PCs of CERES and the model agree to better than 99.9%. The EOF maps are similar for most of the globe, but differ in a few places, and the agreement of the EOF maps is likewise quantified. Maps of differences between the annual cycles show regions of agreement and disagreement.

Smith, G. Louis↗
NASA NTRS2013

Application of Spectral Analysis Techniques in the Intercomparison of Aerosol Data: 1. an EOF Approach to the Spatial-Temporal Variability of Aerosol Optical Depth Using Multiple Remote Sensing Data Sets

Many remote sensing techniques and passive sensors have been developed to measure global aerosol properties. While instantaneous comparisons between pixel-level data often reveal quantitative differences, here we use Empirical Orthogonal Function (EOF) analysis, also known as Principal Component Analysis, to demonstrate that satellite-derived aerosol optical depth (AOD) data sets exhibit essentially the same spatial and temporal variability and are thus suitable for large-scale studies. Analysis results show that the first four EOF modes of AOD account for the bulk of the variance and agree well across the four data sets used in this study (i.e., Aqua MODIS, Terra MODIS, MISR, and SeaWiFS). Only SeaWiFS data over land have slightly different EOF patterns. Globally, the first two EOF modes show annual cycles and are mainly related to Sahara dust in the northern hemisphere and biomass burning in the southern hemisphere, respectively. After removing the mean seasonal cycle from the data, major aerosol sources, including biomass burning in South America and dust in West Africa, are revealed in the dominant modes due to the different interannual variability of aerosol emissions. The enhancement of biomass burning associated with El Niño over Indonesia and central South America is also captured with the EOF technique.

variability↗
NASA NTRS2013

Northern North Atlantic Sea Surface Height and Ocean Heat Content Variability

The evolution of nearly 20 years of altimetric sea surface height (SSH) is investigated to understand its association with decadal to multidecadal variability of the North Atlantic heat content. Altimetric SSH is dominated by an increase of about 14 cm in the Labrador and Irminger seas from 1993 to 2011, while the opposite has occurred over the Gulf Stream region over the same time period. During the altimeter period the observed 0-700 m ocean heat content (OHC) in the subpolar gyre mirrors the increased SSH by its dominantly positive trend. Over a longer period, 1955-2011, fluctuations in the subpolar OHC reflect Atlantic multidecadal variability (AMV) and can be attributed to advection driven by the wind stress ''gyre mode'' bringing more subtropical waters into the subpolar gyre. The extended subpolar warming evident in SSH and OHC during the altimeter period represents transition of the AMV from cold to warm phase. In addition to the dominant trend, the first empirical orthogonal function SSH time series shows an abrupt change 2009-2010 reaching a new minimum in 2010. The change coincides with the change in the meridional overturning circulation at 26.5N as observed by the RAPID (Rapid Climate Change) project, and with extreme behavior of the wind stress gyre mode and of atmospheric blocking. While the general relationship between northern warming and Atlantic meridional overturning circulation (AMOC) volume transport remains undetermined, the meridional heat and salt transport carried by AMOC's arteries are rich with decade-to-century timescale variability.

Sea-surface↗
NASA NTRS2014

The East Atlantic - West Russia Teleconnection in the North Atlantic: Climate Impact and Relation to Rossby Wave Propagation

Large-scale winter teleconnection of the East Atlantic - West Russia (EA-WR) over the Atlantic and surrounding regions is examined in order to quantify its impacts on temperature and precipitation and identify the physical mechanisms responsible for its existence. A rotated empirical orthogonal function (REOF) analysis of the upper-tropospheric monthly height field captures successfully the EA-WR pattern and its interannual variation, with the North Atlantic Oscillation as the first mode. EA-WRs climate impact extends from eastern North America to Eurasia. The positive (negative) EA-WR produces positive (negative) temperature anomalies over the eastern US, western Europe and Russia east of Caspian Sea, with negative (positive) anomalies over eastern Canada, eastern Europe including Ural Mountains and the Middle East. These anomalies are largely explained by lower-tropospheric temperature advections. Positive (negative) precipitation anomalies are found over the mid-latitude Atlantic and central Russia around 60E, where lower-level cyclonic (anticyclonic) circulation anomaly is dominant. The eastern Canada and the western Europe are characterized by negative (positive) precipitation anomalies.The EA-WR is found to be closely associated with Rossby wave propagation. Wave activity fluxes show that it is strongly tied to large-scale stationary waves. Furthermore, a stationary wave model (SWM) forced with vorticity transients in the mid-latitude Atlantic (approximately 40N) or diabatic heat source over the subtropical Atlantic near the Caribbean Sea produces well-organized EA-WR-like wave patterns, respectively. Sensitivity tests with the SWM indicate improvement in the simulation of the EA-WR when the mean state is modified to have a positive NAO component that enhances upper-level westerlies between 40-60N.

EA/WR↗
NASA NTRS2013

North Pacific Decadal Variability in the GEOS-5 Atmosphere-Ocean Model

This study examines the mechanisms of the Pacific decadal oscillation (PDO) in the GEOS-5 general circulation model. The model simulates a realistic PDO pattern that is resolved as the first empirical orthogonal function (EOF) of winter sea surface temperature (SST). The simulated PDO is primarily forced by Aleutian low through Ekman transport and surface fluxes, and shows a red spectrum without any preferred periodicity. This differs from the observations, which indicate a greater role of El Nino-Southern Oscillation (ENSO) forcing, and likely reflects the too short time scale of the simulated ENSO. The geostrophic transport in response to the Aleutian low is limited to the Kuroshio-Oyashio Extension, and is unlikely the main controlling factor in this model, although it reinforces the Ekman-induced SST anomalies. The delay between the Aleutian low and the PDO is relatively short (1 year) suggesting that the fast Ekman response (rather than Rossby wave propagation) sets the SST pattern immediately following an Aleutian low fluctuation. The atmospheric feedback (response to the SST) is only about 25 of the forcing and never evolves into an Aleutian low completely, instead projecting onto the North Pacific Oscillation (NPO), a meridional dipole in sea level pressure (SLP). The lack of preferred periodicity and weak atmospheric response bothindicate a coupled oscillation is an unlikely mechanism for the PDO in this model. In agreement with recent studies, the NPO is correlated with the North Pacific Gyre Oscillation (NPGO), which is another leading EOF of the North Pacific SST. A possible connection between the PDO and the NPGO is discussed.

Pacific Decadal Variability↗
NASA NTRS2013

Trends and Variations of Ocean Surface Latent Heat Flux: Results from GSSTF2c Data Set

Trends and variations of Goddard Satellite-based Surface Turbulent Fluxes (GSSTF) version 2c (GSSTF2c) latent heat flux (LHF) are examined. This version of LHF takes account of the correction in Earth incidence angle. The trend of global mean LHF for GSSTF2c is much reduced relative to GSSTF version 2b Set 1 and Set 2 for the same period 1988-2008. Temporal increase of GSSTF2c LHF in the two decades is 11.0%, in which 3.1%, 5.8%, and 2.1% are attributed to the increase in wind, the increase in sea surface saturated air humidity, and the decrease in near-surface air humidity, respectively. The first empirical orthogonal function of LHF is a conventional El Nino Southern Oscillation (ENSO) mode. However, the trends in LHF are independent of conventional ENSO phenomena. After removing ENSO signal, the pattern of LHF trends is primarily determined by the pattern of air-sea humidity difference trends.

global mean↗
NASA NTRS2016

North Pacific Decadal Variability: Insights from a Biennial ENSO Environment

This study examines the mechanisms of the Pacific decadal oscillation (PDO) in the NASA GEOS-5 general circulation model (GCM). Similar to several other state-of-the-art GCMs, the El Niño-Southern Oscillation (ENSO) simulated by the GEOS-5 has a strong biennial periodicity. Since this is a model bias that precludes a strong role of ENSO, it provides a unique environment to assess the other leading mechanisms of North Pacific decadal variability. Despite the biennial ENSO periodicity, the model simulates a realistic PDO pattern in the North Pacific that is resolved as the first empirical orthogonal function (EOF) of winter mean sea surface temperature (SST). The spectrum of the PDO indicates no preferred periodicity. The SST anomalies associated with the PDO, particularly its basin wide structure, are primarily forced by the Aleutian low through Ekman transport. The slow geostrophic transport in association with the meridional adjustment of the subtropical gyre is limited to a narrow region in the Kuroshio-Oyashio extension, north of 40degN. The atmosphere's response to the PDO, while weak, projects onto the North Pacific Oscillation (NPO), a meridional dipole in sea level pressure. Both the lack of preferred periodicity and the weak atmospheric response indicate an air-sea coupled oscillation is an unlikely mechanism in this model. In agreement with recent studies, the NPO is correlated with the North Pacific Gyre Oscillation, which is another leading EOF of North Pacific SST variability. The results emphasize the role of atmospheric variability in the North Pacific SST modes, thereby bringing into question the potential for their predictability.

GCM↗
NASA NTRS2017

Climate Controls on Tree Growth in the Western Mediterranean

The first large-scale network of tree-ring chronologies from the western Mediterranean (WM; 32 deg N-43 deg N, 10 deg W-17 deg E) is described and analyzed to identify the seasonal climatic signal in indices of annual ring width. Correlation and rotated empirical orthogonal function analyses are applied to 85 tree-ring series and corresponding gridded climate data to assess the climate signal embedded in the network. Chronologies range in length from 80 to 1129 years. Monthly correlations and partial correlations show overall positive associations for Pinus halepensis (PIHA) and Cedrus atlantica (CDAT) with winter (December-February) and spring (March-May) precipitation across this network. In both seasons, the precipitation correlation with PIHA is stronger, while CDAT chronologies tend to be longer. A combination of positive correlations between growth and winter-summer precipitation and negative partial correlations with growing season temperatures suggests that chronologies in at least part of the network reflect soil moisture and the integrated effects of precipitation and evapotranspiration signal. The range of climate response observed across this network reflects a combination of both species and geographic influences. Western Moroccan chronologies have the strongest association with the North Atlantic Oscillation.

chronology development↗
NASA NTRS2016

Degree of Ice Particle Surface Roughness Inferred from Polarimetric Observations

The degree of surface roughness of ice particles within thick, cold ice clouds is inferred from multidirectional, multi-spectral satellite polarimetric observations over oceans, assuming a column-aggregate particle habit. An improved roughness inference scheme is employed that provides a more noise-resilient roughness estimate than the conventional best-fit approach. The improvements include the introduction of a quantitative roughness parameter based on empirical orthogonal function analysis and proper treatment of polarization due to atmospheric scattering above clouds. A global 1-month data sample supports the use of a severely roughened ice habit to simulate the polarized reflectivity associated with ice clouds over ocean. The density distribution of the roughness parameter inferred from the global 1- month data sample and further analyses of a few case studies demonstrate the significant variability of ice cloud single-scattering properties. However, the present theoretical results do not agree with observations in the tropics. In the extra-tropics, the roughness parameter is inferred but 74% of the sample is out of the expected parameter range. Potential improvements are discussed to enhance the depiction of the natural variability on a global scale.

surface↗
NASA NTRS2017

An Investigation of Dust Storms Observed with the Mars Color Imager

Daily global imaging by the Mars Color Imager (MARCI) continues the record of the Mars Orbiter Camera (MOC) and has allowed creation of a long-duration record of Martian dust storms. We observe dust storms over the first two Mars years of the MARCI record, including tracking individual storms over multiple sols, as well as tracking the growth and recession of the seasonal polar caps. Using the combined 6 Mars year record of textured dust storms (storms with visible textures on the observed dust cloud tops), we study the relationship between textured dust storm activity and meteorology (as simulated by the MarsWRF general circulation model) and surface properties. We find that textured dust storms preferentially occur in places and seasons with above average surface wind stress. Textured dust storm occurrence also has a modest linear anti-correlation with surface albedo (0.43) and topography (0.40). Lastly, we perform an empirical orthogonal function (EOF) analysis on the distribution of occurrence of textured dust storms and find that over 50 of the variance in textured dust storm activity can be explained by two EOF modes. We associate the first EOF mode with cap-edge storms just before Ls = 180deg and the second EOF mode with flushing dust storms that occur from Ls = 180-210deg and again near Ls = 320deg.

Guzewich, Scott D.↗