一种TEC地图构建方法与全球电子总含量的气候学研究

Global positioning system (GPS) can not only provide precise service for navigation and timing, but also be used to investigate the ionospheric variation. From the GPS observations, we can obtain total electron content (TEC), so-called GPS TEC, which is used to characterize the ionospheric structure. This thesis mainly concerns about GPS TEC data processing and ionospheric climatological analysis as follows. Firstly, develop an algorithm for high-resolution global ionospheric TEC mapping. According to current algorithms in global TEC mapping, we propose a practical way to calibrate the original GPS TEC with the existing GIM results. We also finish global/local TEC mapping by model fitting with the processed GPS TEC data; in practice, we apply it into the local TEC mapping in Southeast of China and obtain some initial results. Next, suggest a new method to calculate equivalent ionospheric global electron content (GEC). We calculate such an equivalent GEC with the TEC data along the geographic longitude 120°E. With the climatological analysis, we can see that GEC climatological variation is mainly composed of three factors: solar cycle, annual and semiannual variations. Solar cycle variation is dominant among them, which indicates the most prominent influence; both annual and semiannual variations play a secondary role and are modulated by solar activity. We construct an empirical GEC model driven by solar activity and seasonal factors on the basis of partial correlation analysis. Generally speaking, our researches not only show that GPS is advantageous in now-casting ionospheric TEC as an important observation, but also show that GEC may become a new index to describe the solar influence on the global ionosphere since the great correlation between GEC and solar activity factor indicates the close relationship between the ionosphere and solar activity.

计算层状介质中电磁场的层矩阵法

As an important branch of electrical prospecting method, the artificial source frequency domain electromagnetism method has received more and more attention. But when conducts the fundamental research, people often isolated study some concrete method, so the research results of one method are very difficult to apply to another method directly. This article will possess the artificial source frequency domain EM method to an 1D model simply. It is stratified medium model, with an electric or magnetic source in or outside of it. Then take the horizontal electric dipole source as an example to introduce how to computing the EM field in stratified medium. Because layer matrix is the key of establishing equations, so we call it the layer-matrix method. The key of layer-matrix method is establishing equations by using layer matrixes in wavenumber(kx, ky, z) domain, then obtains the electromagnetic field value of wavenumber domain. After Fourier transform, we can get electromagnetic field of any position in spatial domain. The layer matrix technique theoretically can calculate electromagnetic field of any position for any source, is suitable for many kinds of electromagnetic method. After introduction of the layer matrix method, this article has done some CSAMT, MCSEM and Wireless Electro-Magnetic Method (WEM) modeling with layer matrix method separately. In CSAMT modeling, we get electromagnetic field dissemination characteristics considering wave number of the air, and obtain three-dimensional distribution characteristics of the electromagnetic field. In MCSEM modeling, we get electromagnetic field dissemination characteristics with and without considering the airwave, and obtain three-dimensional distribution characteristics of electromagnetic field. In WEM modeling, we get electromagnetic field’s difference between considering the ionosphere and not considering it, and recognize the ionosphere’s influence of electromagnetic field. With the layer matrix technique, we have got some new understandings of EM dissemination rules of different situations. All analysis results indicate that the layer-matrix technique is credible and effective, and are worthy of further thorough research and development.

太阳耀斑和日食期间电离层行为的统计和模拟研究

The ionosphere is the ionized component of the Earth's upper atmosphere. Solar EUV radiation is the source of ionospheric ionization. Thus the ionosphere is affected strongly by the variations in solar radiation. Solar flares and solar eclipses can induce remarkable short time changes in solar radiation: the solar radiation would increase suddenly during solar flares and decrease significantly during solar eclipses. Solar flare and eclipse events not only affect directly the photochemical processes, but also affect the dynamic processes, and even affect the neutral atmosphere, which is strongly coupled with the ionosphere. The study on the ionospheric response to solar flares and eclipses can advance our knowledge on the ionosphere and its photochemical and dynamic processes and help us to evaluate the ionospheric parameters (such as ion loss coefficients). In addition, the study on the ionospheric responses to solar flares and eclipses is an important part of the ionospheric space weather, which can provide guides for space weather monitoring. This thesis devotes to the study on the ionospheric responses to solar flares and solar eclipses. I have developed two models to simulate the variations of solar EUV radiation during solar flares and solar eclipses, and involved in developing a 2D mid- and low-latitude ionospheric model. On the basis of some observed data and the ionospheric model, I study the temporal and spatial variations of the ionosphere during solar flares and eclipses, and investigate the influences of solar activity, solar zenith angle, neutral gas density, and magnetic dip angle on the ionospheric responses to solar flares and solar eclipses. The main points of my works and results are summarized as follows. 1. The ionospheric response to the X17.2 solar flare on October 28, 2003 was modeled via using a one-dimension theoretical ionospheric model. The simulated variation of TEC is in accordance with the observations, though there are some differences in the amplitude of the variation. Then I carried out a series of simulations to explore the local time and seasonal dependences of the ionospheric responses to solar flares. These calculations show that the ionospheric responses are largely related with the solar zenith angle (SZA). During the daytime (small SZA), most of the increases in electron density occur at altitudes below 300 km with a peak at around 115 km; whereas around sunrise and sunset (SZA>90°), the strongest ionospheric responses occur at much higher altitudes. The TEC increases slower at sunrise than at sunset, which is caused by the difference in the evolution of SZA at sunrise and sunset: SZA decreases with time at sunrise and increase with time at sunset. The ionospheric response is largest in summer and smallest in winter, which is also related to the seasonal difference of SZA. 2. Based on the observations from the ionosondes in Europe and the ionospheric model, I investigated the differences of the ionosphere responses to solar eclipses between the E-layer and F1-layer. Both the observation and simulation show that the decrease in foF1 due to the solar eclipses is larger than that in foE. This effect is due to that the F1 region locates at the transition height between the atomic ion layer and the molecular ion layer. With the revised model of solar radiation during solar flares, our model calculates the radiations from both the inside and outside of photosphere. Large discrepancy can be found between the observations and the calculations with an unrevised model, while the calculations with the revised model consist with the observations. 3. I also explore the effects of the F2-layer height, local time, solar cycle, and magnetic dip angle on the ionospheric responses to solar eclipses via using an ionospheric model and study on the solar zenith angle and the dip dependences by analyzing the data derived from 23 ionosonde stations during seven eclipse events. Both the measured and simulated results show that these factors have significant effect on the ionospheric response. The larger F2-layer height causes the smaller decrease in foF2, which is because that the electron density response decreases with height. The larger dip results in the smaller eclipse effect on the F2 layer, because the larger dip would cause the more diffusion from the top ionosphere which can make up for the plasma loss. The foF2 response is largest at midday and decreases with the increasing SZA. The foF2 response is larger at high solar activity than at low solar activity. The simulated results show that the local time and solar activity discrepancy of the eclipse effect mainly attribute to the difference of the background neutral gas density. 4. I carried out a statistical study on the latitudinal dependence of the ionospheric response to solar eclipses and modeled this latitudinal dependence by the ionospheric model. Both the observations and simulations show that the foF2 and TEC responses have the same latitudinal dependence: the eclipse effects on foF2 and TEC are smaller at low latitudes than at middle latitudes; at the middle latitudes (>40°), the eclipse effect decreases with increasing latitude. In addition, the simulated results show the change in electron temperature at the heights of above 300 km of low latitudes is much smaller than that at the same heights of middle latitudes. This is due to the smaller decrease in photoelectron production rate at its conjugate low heights. 5. By analyzing the observed data during the October 3, 2005 solar eclipse, I find some significant disturbances in the conjugate region of the eclipse region, including a decrease in Te, an increase in foF2 and TEC, and an uprising in hmF2. I also simulated the ionosphere behavior during this eclipse using a mid-low latitude ionospheric model. The simulations reproduce the measured ionospheric disturbances mentioned above in the conjugated hemisphere. The simulations show that the great loss of arriving photoelectron heat from the eclipse region is the principal driving source for the disturbances in the conjugate hemisphere.

电离层垂直TEC映射函数的实验观测与统计特性

When used in the determining the total electron content (TEC), which may be the most important ionospheric parameter, the worldwide GPS observation brings a revolutionary change in the ionospheric science. There are three steps in the data processing to retrieve GPS TEC: (1) to estimate slant TEC from the measurements of GPS signals; (2) to map the slant TEC into vertical; and (3) to interpolate the vertical TEC into grid points. In this scientific dissertation we focus our attention on the second step, the mapping theory and method to convert slant TEC into vertical. This is conventionally done by multiplying on the slant TEC a mapping function which is usually determined by certain models of electron density profile. Study of the vertical TEC mapping function is of significance in GPS TEC measurement. This paper first reviews briefly the three steps in GPS TEC mapping process. Then we compare the vertical TEC mapping function which were respectively calculated from the electron density profiles of the ionospheric model and retrieved from the observation of worldwide GPS TEC. We also perform the statistical analysis on the observational mapping functions. The main works and results are as follows: 1. We calculated the vertical TEC mapping functions for both SLM and Chapman models, and discussed the modulation of the ionosphere height to the mapping functions. We use two simple models, single layer model (SLM) and Chapman models, of the ionospheric electron density profiles to calculate the vertical TEC mapping function. In the case of the SLM, we discuss the control of the ionospheric altitude, i.e., the layer height hipp, to the mapping function. We find that the mapping function decreases rapidly as hipp increases. For the Chapman model we study also the control mapping function by both ionospheric altitude indicated by the peak electron density height hmF2, and the scale height, H, which present the thickness of the ionosphere. It is also found that the mapping function decreases rapidly as hmF2 increases. and it also decreases as H increases. 2. Then we estimate the mapping functions from the GPS observations and compare them with those calculated from the electron density models. We first, proposed a new method to estimate the mapping functions from GPS TEC data. This method is then used to retrieve the observational mapping function from both the slant TEC (TECS) provided by International GPS Service (IGS）and vertical TEC provide by JPL Global Ionospheric Maps (GIMs). Then we compare the observational mapping function with those calculated from the electron density models, SLM and Chapman. We find that the values of the observational mapping functions are much smaller than that from the model mapping functions, when the zenith angle is large enough. We attribute this to the effect of the plasmasphere which is above about 1000 km. 3. We statistically analyze the observational mapping functions and reveal their climatological changes. Observational mapping functions during 1999-2007 are used in our statistics. The main results are as follows. (1) The observational mapping functions decrease obviously with the decrement of the solar activity which is represented by the F10.7 index; (2) In annual variations of the observational mapping functions, the semiannual component is found at low-latitudes, and the remarkable seasonal variations at mid- and high-latitudes. (3) The diurnal variation of the observational mapping functions is that they are large in daytime and small at night, they become extremely small in the early morning before sunrise. (4) The observational mapping functions change with latitudes that they are smaller at lower latitudes and larger at higher. All of the above variations of the observational mapping functions are explained by the existence of the plasmasphere, which changes more slowly with time and more rapidly with latitude than the ionosphere does . In summary, our study on the vertical TEC mapping function imply that the ionosphere height has a modulative effect on the mapping function. We first propose the concept of the 'observational mapping functions' , and provide a new method to calculate them. This is important in improving the TEC mapping. It may also possible to retrieving the plasmaspheric information from GPS observations.

Millstone Hill电离层F层等效风场的气候学研究

Neutral winds and electric fields in the ionospheric F layer play important roles in the variations of the ionosphere, and also affect the thermospheric circulation via the close coupling between the ionosphere and the thermosphere. By now, the neutral winds and electric drifts are generally observed with ground-based Fabry-Perot interferometers (FPI) and incoherent scatter radars (ISR), rockets, and satellite-borne instrument. Based on the servo theory, the ionospheric equivalent winds, which include the information of both the neutral winds and electric fields, can be derived from these characteristic parameters observed by ionosondes. This indirect derivation has potential values in climatological researches and space weather forecast. With the data set of the incoherent scatter radar observations at Millstone Hill, USA, from 1976 to 2006, we statistically analyzed the climatological variations of the vertical component of the equivalent winds (VEWs) over Millstone Hill, which are derived from the ionospheric key parameters (the peak electron number density and peak height of the F2 layer, NmF2 and hmF2) on the basis of the servo theory, Liu's method, and measurements from the ion line-of-sight velocity as well. The main results of this analysis are summarized as follows: (1) The values of VEWs over Millstone Hill during nighttime are stronger than in the daytime, and the upward drift dominates most of the day. In 1993, Hagan found that the component of the neutral winds in the magnetic meridion in daytime is weaker than during nighttime under both solar maximum and minimum conditions; he also found that the equatorward winds dominate most of the day. Both results suggest that the thermosphere in Millstone Hill is modulated by the aurorally driven high-latitude circulation cell; that is, during geomagnetic quiet periods, the average auroral activity is strong enough to drive thermospheric circulation equatorward for most of the day at Millstone Hill. Moreover, since ion drag is the strongest during daytime when F region densities are enhanced by photoionization, the wind speeds are smaller during the daytime than in the nighttime. (2) There is equinoctial symmetry in VEWs at Millstone Hill. The amplitudes and phases of VEWs in spring are quite similar to those in autumn. In contrast, the nighttime upward drift in winter is weaker than in summer and the difference becomes more significant with increasing solar activity. This solstice asymmetry indicates that, the aurorally driven circulation in the northern hemisphere at Millstone Hill latitude is weaker in winter due to arctic darkness, because the subsolar point is in the southern hemisphere. (3) The comparison of the VEWs derived from three methods, i.e., the servo theory, Liu's method, and the ISR ion line-of-sight velocity measurements, indicates that the amplitudes and main phase tendencies of these VEWs accord well with each other during nighttime hours. However, the case in the daytime is relatively worse. This daytime discrepancy can be explained in terms of the effects of photochemical processes and the choices of the servo constants. A larger servo constant gives a stronger plasma drift in daytime. Therefore, this study tells how important to choose a suitable constant for deriving VEWs at Millstone Hill.

通过第一性原理计算研究内地核的结构与成分

The determination of the composition and structure of the Earth’s inner core has long been the major subject in the study of the Earth’s deep interior. It’s widely believed that the Earth’s core is formed by iron with a fraction of nickel. However, light elements must exist in the inner core because the earth core is less dense than pure iron-nickel alloy (~2-3% in the solid inner core and ~6-7% in the liquid outer core). The questions are what and how much light element is there in the iron-nickel alloy. Besides the composition, the crystal structure of the iron with or without light element is also not well known. According to the seismological observations, the sound waves propagate 3-4% faster along the spin axis than in the equatorial plane. That means the inner core is anisotropic. The densest structure of iron-nickel alloy should be h.c.p structure under the very high pressures. However, the h,c,p structure does not propagate waves anisotropic ally. Then what is the structure of the iron-nickel alloy or the iron-nickle-light element alloy. In this study, we tried to predict the composition and the structure of the inner core through ab initio calculation of the Gibbs free energy, which is a function of internal energy, density and entropy. We conclude that the h.c.p structure is more stable than the b.c.c structure under high pressure and 0 K, but with the increase of temperature, the free energy of the b.c.c structure is decreasing much faster than the h.c.p structure caused by the vibration of the atomics, so the b.c.c structure is more stable at high temperatures. With the addition of light elements (S or Si or both), the free energy of b.c.c. decreases even faster, about 3at% of Si not only explains why the inner core is about 2-3 % lighter than the iron-nickle alloy, but also reasons why the inner core is anisotropic, since the b.c.c. structure becomes more stable than the h.c.p structure at 5500-6000K and b.c.c. is anisotropic in propagating seismic waves. Therefore, we infer that the inner core of the earth is formed by b.c.c iron and a fraction of nickel plus ~3at.% Si, with a temperature higher than 5500K, which is consistent with the studies from other approaches.

中国南方热带地区田洋玛珥湖24万年以来的气候变化及湖泊演变

A continuous long (224m) and high-resolution core TY2 was recovered from paleo-maar-lake Tianyang, tropical South China. Based on the diatom records of the upper 130-m core, this paper focuses on exploring climate change and the lake evolution history in tropical South China during the past 240ka. The most typical and unique characteristics of the diatom assemblages is that, Aulacoseira granulata was dominant or absolutely dominant species (80-90%) during most parts of the 130-m core, while Cyclotella stelligera var. tenuis and Fragilaria construens var. venter were subdominant species in only limited parts of the lower and upper core, respectively. Time scale is always the biggest problem for the study of TY2 core, so although diatom is seldom used for establishing time scale, here we attempt this by correlating the diatom-reconstructed temperature sequence with the time scale of ODP core 806B from Equatorial Western Pacific. Verified by the few most reliable ages from TY2 core and the parallel core TYl, a rather reasonable and reliable time scale was established. 01S 7/6 falls at the depth of 100m (ca. 194kaBP), OIS 6/5 at 75m (ca. 132kaBP), OIS 5/4 at 46m (ca. 75kaBP), OIS 4-3 at 35m (ca. 60kaBP). Qualitative and quantitative environmental reconstructions are made on the basis of diatom assemblage ecotype and EDDI dataset. Correlation of diatom-reconstructed temperature and moisture changes of Core TY2 with pollen-reconstructed temperature and rainfall sequence of Core TYl proves that the results are quite consistent in most periods. Thus the reconstruction results from diatom are quite reliable, and probably have a much higher resolution than pollen results. Combined with lithological and magnetic susceptibility variations, the diatom analysis reveals that, the general climate in tropical South China during the past 240ka was warm and wet. On the time scale of glacial-interglacial, warm and wet, cool and dry are not always synchronous. It was relatively warm-wet during the penultimate interglacial, cool-dry during the penultimate glacial, warm-dry during the last interglacial, and cooler-drier during the last glacial. In contrast, on the time scale of subglacial-subinterglacial scales, warm and dry, cool and wet corresponds very obviously, showing very clear 21-23 ka precession cycle. Analysis also shows that, the water of Tianyang paleo-maar-lake was generally warm, turbulent, turbid, meso-trophic, slightly alkaline, low conductivity and fresh during the past 240 kaBP, with small variations in some parts. Tianyang paleolake experienced shallow to semi-deep lake in OIS7d, open shallow lake in OIS7c-OIS5b, shallow coastal lake in OIS5a-OIS4c, swamp in OIS4b, and then completely dried up in OIS3c. The lake evolution was mainly controlled by temperature and precipitation changes in tropical China. While temperature and precipitation changes were probably controlled by the migrations of monsoon rainband and the evaporation rate, which was in turn controlled by the evolution of East Asian monsoon. Therefore, when the summer monsoon was strongest the climate was warm-dry, when stronger the climate was warm-wet; when the winter monsoon was strongest the climate was cool-dry, stronger cool-wet. This mechanism caused the warm-dry sub interglacial and cool-wet subglacial climate in the tropical South China.

岱海岩芯沉积的粒度组成与湖区全新世降水变化历史

Sediment cores DH99a and DH99b recovered in the central part of Daihai Lake in north-central China were analysed at 2- to 4-crn intervals for grain-size distribution. Grain-size distributions of the lake sediments are inferred to be a proxy for past changes in East Asian monsoon precipitation, such that greater silt-size percentage and higher median grain size reflect increased monsoonal precipitation rates. The grain-size record of Daihai Lake sediments spanning the last ca 11,000 yr indicates that the monsoonal precipitation in the lake region can be divided into three stages: the Early, Middle and Late Holocene. During the Early Holocene before ca 7900 cal yr BP, the median grain size (Md) and the silt-fraction content were relatively low and constant, suggesting relatively low precipitation over the lake region. The Middle Holocene between ca 7900 and 3100 cal yr BP was marked by intensified and highly variable monsoonal precipitation, as indicated by high and variable Md values and silt contents of the lake sediments. During this period, average precipitation rate gradually increased from ca 7900 to 6900 cal yr BP, displayed intense oscillations between ca 6900 and 4400 cal yr BP, and exhibited a decreasing trend while fluctuating from ca 4400 to 3100 cal yr BP Although generally high during the Middle Holocene, both the Md and the silt content assumed distinctly low values at the short intervals of ca 6500-6400, 6000-5900, 5700-5600, 4400-4200 cal yr BP, implying that monsoonal precipitation might have been significantly reduced during these intervals. During the Late Holocene since ca 3100 cal yr BP, grain-size values suggest that precipitation decreased. However, during the Late Holocene, relatively higher Md values and silt contents occurring between ca 1700 to 1000 cal yr BP may denote an intensification of hydrological cycles in the lake area. Changes in the East Asian monsoonal precipitation were not only directly linked with the changing seasonality of solar insolation resulting from progressive changes in the Earth's orbital parameters, but also may have been closely related to variations in the temperature and size of the Western Pacific Warm Pool, in the intensity of the El Nino-Southern Oscillation, and in the path and strength of the North Equatorial Current in the western Pacific.

喀斯特地质灾害的形成条件及其环境背景——以贵州省为例

贵州喀斯特环境地质灾害类型多样、分布广泛、突发性强,其发生、发育过程不但具有明显的季节性,且具有同步叠加、连续发生及交替出现等特点;地质灾害的危害程度、暴发频率也在逐年增大。其形成过程主要受岩性、地貌、灾害性降水和人类经济活动等条件影响。地球内动力地质作用决定了喀斯特地质灾害的环境背景,如广布的碳酸盐岩和陡峻的地貌特征等,特别是挽近期的构造运动对近、现代地质灾害具有明显的控制用;人为活动改变了岩(土)体的平衡状态而导致环境地质灾害发生。

Modeling subauroral polarization streams equatorward of the plasmapause footprints

Pint?r, B.; Thom, S. D.; Balthazor, R.; Vo, H.; Bailey, G. J., Modeling subauroral polarization streams equatorward of the plasmapause footprints, Journal of Geophysical Research, Volume 111, Issue A10, CiteID A10306 RAE2008

Evidence for the tongue of ionization under northward interplanetary magnetic field conditions

Pryse, Sian; Middleton, H. R.; Kersley, L.; Bust, G. S., 'Evidence for the tongue of ionization under northward interplanetary magnetic field conditions', Journal of Geophysical Research (2005) 110(A7) pp.A07301 RAE2008

Large-scale structure of the fast solar wind

Breen, Andrew; Bisi, M.M.; Fallows, R.A.; Habbal, S.R., (2007) 'Large-scale structure of the fast solar wind', Journal of Geophysical Research 112(A6) pp.A06101 RAE2008

The dayside high-latitude trough under quiet geomagnetic conditions: Radio tomography and the CTIP model

Pryse, Sian; Dewis, K.L.; Middleton, H.R.; Balthazor, R.L., (2005) 'The dayside high-latitude trough under quiet geomagnetic conditions: Radio tomography and the CTIP model', Annales Geophysicae 23(4) pp.1199-1206 RAE2008

Contrasting the efficiency of radiation belt losses caused by ducted and nonducted whistler-mode waves from ground-based transmitters

It has long been recognized that whistler-mode waves can be trapped in plasmaspheric whistler ducts which guide the waves. For nonguided cases these waves are said to be "nonducted", which is dominant for L < 1.6. Wave-particle interactions are affected by the wave being ducted or nonducted. In the field-aligned ducted case, first-order cyclotron resonance is dominant, whereas nonducted interactions open up a much wider range of energies through equatorial and off-equatorial resonance. There is conflicting information as to whether the most significant particle loss processes are driven by ducted or nonducted waves. In this study we use loss cone observations from the DEMETER and POES low-altitude satellites to focus on electron losses driven by powerful VLF communications transmitters. Both satellites confirm that there are well-defined enhancements in the flux of electrons in the drift loss cone due to ducted transmissions from the powerful transmitter with call sign NWC. Typically, ∼80% of DEMETER nighttime orbits to the east of NWC show electron flux enhancements in the drift loss cone, spanning a L range consistent with first-order cyclotron theory, and inconsistent with nonducted resonances. In contrast, ∼1% or less of nonducted transmissions originate from NPM-generated electron flux enhancements. While the waves originating from these two transmitters have been predicted to lead to similar levels of pitch angle scattering, we find that the enhancements from NPM are at least 50 times smaller than those from NWC. This suggests that lower-latitude, nonducted VLF waves are much less effective in driving radiation belt pitch angle scattering. Copyright 2010 by the American Geophysical Union.

Intensification of summer rainfall variability in the southeastern United States during recent decades

The variability of summer precipitation in the southeastern United States is examined in this study using 60-yr (1948-2007) rainfall data. The Southeast summer rainfalls exhibited higher interannual variability with more intense summer droughts and anomalous wetness in the recent 30 years (1978-2007) than in the prior 30 years (1948-77). Such intensification of summer rainfall variability was consistent with a decrease of light (0.1-1 mm day-1) and medium (1-10 mm day-1) rainfall events during extremely dry summers and an increase of heavy (.10 mm day-1) rainfall events in extremely wet summers. Changes in rainfall variability were also accompanied by a southward shift of the region of maximum zonal wind variability at the jet stream level in the latter period. The covariability between the Southeast summer precipitation and sea surface temperatures (SSTs) is also analyzed using the singular value decomposition (SVD) method. It is shown that the increase of Southeast summer precipitation variability is primarily associated with a higher SST variability across the equatorial Atlantic and also SST warming in the Atlantic. © 2010 American Meteorological Society.