Relation of PC index to the geomagnetic storm Dst variation

Relationships between the polar cap magnetic activity index PC and the magnetic storm Dst index have been studied for the magnetic storms with duration more 12 h and peak value Dst<?30 nT and, observed in 1998–2002 and 2004–2005. Along with PC index the geoeffective interplanetary electric field Em was also examined. It has been found that all examined storms, lying in range from ?30 to ?373 nT, started when the PC index and, correspondingly, the Em field firmly exceeded the threshold >2 mV/m. In particular, the “anomalous” magnetic storm on January 21–22, 2005 occurring under conditions of northward IMF B_Z (Du et al., 2008) is usual phenomena fitted well with the threshold restriction owing to the large IMF By component input. The maximal storm depression (the peak value of Dst) is linearly related to the quantities Em and PC, averaged for the time interval from the storm beginning to the storm maximum. The correlation between Dst and PC is more steady and larger than correlation between Dst and Em, the latter being dependent on Em value (effect of “Dst saturation”). The moment of the firm descent of the Em and PC quantities below the threshold level ～2 mV/m is indicative of the depression damping and transition to the recovery phase. The results are consistent with the similar peculiarities revealed for substorms development (Troshichev and Janzhura, 2009) and support the conclusion that the PC index is a reliable proxy characterizing the solar wind energy having been entered into the magnetosphere.     

Does the polar cap disappear under an extended strong northward IMF?

The suggestion that the polar cap can completely disappear under certain northward IMF conditions is still controversial. We know that the size of the polar cap is strongly controlled by the interplanetary magnetic field (IMF). Under a southward IMF, the polar cap is usually large and filled with weak diffuse polar rain electrons. The polar cap shrinks under a northward IMF. Here we use the global auroral images and coincident particle measurements on May 15, 2005 to show that the discrete arcs (due to precipitation of both electrons and ions) expanded from the dayside oval to the nightside oval and filled the whole polar ionosphere after a long (8 h) and strong (～5–30 nT) northward IMF B_z, The observations suggested that the polar cap disappeared under a closed magnetosphere.     

Spectacular low- and mid-latitude electrical fields and neutral winds during a superstorm

In November 2004, a major magnetic storm occurred, a lengthy portion of which was recorded by the Upper Atmospheric Radar Chain. On the 9th and 10th, the Jicamarca Radar detected the highest magnitude penetrating electric fields (±3 mV/m) and vertical drifts (±120 m/s) ever seen at this premiere facility. These large and variable drifts were highly correlated with the interplanetary magnetic and electric fields and created a double F layer on the dayside and unusual TEC behavior throughout the low-latitude zone. These solar wind-induced drifts both suppressed and generated irregularities at the magnetic equator at different times. Large-scale thermospheric disturbances were generated by high-latitude heating and tracked through the middle- to low-latitude zones where both parallel and perpendicular plasma drifts created major ionospheric changes. The auroral oval was located at a magnetic L shell of about three for many hours.     

Long-term thermospheric neutral wind observations over the northern polar cap

We study the solar dependence of the thermospheric dynamics based on more than 20 years Fabry–Perot interferometer O 6300 Å emission observation of polar cap thermospheric wind from three stations: Thule (76.53°N, 68.73°W, MLAT 86N), Eureka (80.06°N, 86.4°W, MLAT 89N), and Resolute (74.72°N, 94.98°W, MLAT 84N) in combination with the National Center for Atmospheric Research Thermosphere Ionosphere Electrodynamics General Circulation Model (NCAR-TIEGCM). All three stations showed a dominant diurnal oscillation in both the meridional and zonal components, which is a manifestation of anti-sunward thermospheric wind in the polar cap. The three-station observations and the TIEGCM simulation exhibit varying degree of correlations between the anti-sunward thermospheric wind and solar F10.7 index. The diurnal oscillation is stronger at Eureka (∼150 m/s) than that at Resolute (∼100 m/s) according to both observations and TIEGCM simulation. The semidiurnal oscillation is stronger at Resolute (∼20 m/s) than that at Eureka based (∼10 m/s) on data and model results. These results are consistent with a two-cell convection pattern in the polar cap thermospheric winds. The Thule results are less consistent between the model and observations. The simulated meridional wind diurnal and semidiurnal oscillations are stronger than those observed.     

Response of the inner and outer magnetosphere to solar wind density fluctuations during the recovery phase of a moderate magnetic storm

We examine the geomagnetic field and space plasma disturbances developing simultaneously in the solar wind, in the inner and outer magnetosphere, and on the ground from 0730 to 2030 UT on April 11, 1997 during the recovery phase of a moderate magnetic storm. The fluctuations of the solar wind density, H-component of the geomagnetic field, and power of Pc1–2 (0.1–5 Hz) waves at middle and low latitudes evolve nearly simultaneously. These fluctuations also match very well with variations of density and flux of the magnetospheric plasma at the geosynchronous orbit, and of the geomagnetic field at the geosynchronous orbit and northern polar cap. The time delay between the occurrence of disturbances in different magnetosphere regions matches the time of fast mode propagation. These disturbances are accompanied by the generation of Pc1–2 waves at mid- and high-latitude observatories in nearly the same frequency range. A scenario of the evolution of wave phenomena in different magnetospheric domains is proposed.     

Polar cap particle precipitation and aurora: Review and commentary

Polar rain has a beautiful set of symmetry properties, individually established, but not previously discussed collectively, which can be organized by a single unifying principle. The key polar rain properties are favored hemisphere (controlled by the interplanetary magnetic field B_x), dawn/dusk gradient (IMF B_y), merging rate (IMF B_z or more generally dΦ_MP/dt), nightside/dayside gradient, and seasonal effect. We argue that all five properties involve variants on a single theme: the further downstream a field line exits the magnetosphere (or less directly points toward the solar wind electron heat flux), the weaker the polar rain. This effect is the result of the requirements of charge quasi-neutrality, and because the ion thermal velocity declines and the tailward ion bulk flow velocity rises moving down tail from the frontside magnetopause.Polar cap arcs (or more properly, high-latitude sun-aligned arcs) are largely complementary to the polar rain, occurring most frequently when the dayside merging rate is low, and thus when polar rain is weak. Sun-aligned arcs are often considered as originating either in the polar rain or the expansion of the plasma sheet into the polar cap. In fact three quite distinct types of sun-aligned high-latitude arcs exist, two common, and one rare. One type of arc occurs as intensifications of the polar rain, and is common, but weak, typically <0.1 ergs/cm² s, and lacks associated ion precipitation. A second category of Sun-aligned arcs with energy flux >0.1 ergs/cm² s usually occurs adjacent to the auroral oval, and includes ion precipitation. The plasma regime of these common, and at times intense, arcs is often distinct from the oval which they abut. Convection alone does not specify the open/closed nature of these arcs, because multiple narrow convection reversals are common around such arcs, and the arcs themselves can be embedded within flows that are either sunward or anti-sunward. These observational facts do not neatly fit into either a plasma sheet origin or a polar rain origin (e.g., the necessity to abut the auroral oval, and the presence of ions does not fit the properties of polar rain, which can in any event be nearly absent for northward interplanetary magnetic field). One theory is that such arcs are associated with merging tailward of the cusp. Both of these common types of sun-aligned arcs fade within about 30 min of a southward IMF turning.The third, and rarest, category of sun-aligned arcs are intense, well detached from the auroral oval, contain plasma sheet origin ion precipitation as well as electrons, and persist for hours after a southward turning. These intense detached sun-aligned arcs can rapidly cross the polar cap, sometimes multiple times. Most events discussed in the literature as “theta-aurora” do not fit into this category (for example, although they may appear detached in images, they abut the oval in particle data, and do not have the persistence of detached events under southward IMF turnings). It is possible that no single theory can account for all three types of sun-aligned arcs.Solar energetic particle (SEP) events are at times used to demarcate polar cap open/closed boundaries. Although this works at times, examples exist where this method fails (e.g., very quiet conditions for which SEP reaches below L=4), and the method should be used with caution. Finally, it is shown that, although it is rare, the polar cap can at times completely close.     

Detection of large scale TIDs associated with the dayside cusp using SuperDARN data

Variations in the dayside ionosphere parameters as a result of a large-scale acoustic gravity wave (LS AGW) were studied for the 17 February 1998 substorm using the super dual auroral radar network (SuperDARN) measurements. This event was characterised by a sharp rise in the AE index with a maximum of ～900 nT. The source of the disturbance responsible for the LS AGW appears to have been located within the plasma convection throat and in the dayside cusp region. The location of the source was obtained from studies of a number of datasets including high-latitude convection maps, data from 4 DMSP satellites and networks of ground-based magnetometers. The propagation of the LS AGWs caused quasi-periodic variations in the skip distance (with an amplitude up to 220–260 km) of the ground backscatter measured by up to 6 SuperDARN radars, including Goose Bay and Kapuskasing, resulting in two large-scale travelling ionospheric disturbances (LS TIDs). The LS TIDs had wave periods of 1.5 and 2 h, a velocity of ～400 m/s for both, and wavelengths of 2200 and 2900 km, respectively. These quasi-periodic variations were also present in the peak electron density and height of the F2 layer measured by the Goose Bay ionosonde. The numerical simulation of the inverse problem show good agreement between Goose Bay radar and Goose Bay ionosonde measurements. But these LS TIDs would be difficult to deduce from the ground based ionospheric station data alone, because h_mF2 variations were 10–40 km only and f_OF2 variations between 10% and 20%. The results demonstrate how important SuperDARN radars can be, and that this is a more powerful technique than routine ground-based sounding for studies of weak quasi-periodic variations in the dayside subauroral ionosphere related to LS AGW.     

Hybrid model for nowcasting and forecasting the K index

A new empirical model nowcasting and predicting a proxy to the geomagnetic K index is developed, which is based on the combined use of solar wind parameters and ground-based magnetic data. The present approach implements the previously developed solar wind-based MAK model, calibrating its values with magnetogram-derived K index. The new model is named as Hybrid Dourbes K (HDK) model. The HDK nowcast model provides the quantity Kdf, obtained by solar wind-based Ksw and corrected with a combination of differences between several past values of Kd and Ksw. The model error of the nowcast Kdf is found to be 0.38 KU, or nearly twice less than that of the MAK model. Kdf has a good predictability. Prediction made by weighted extrapolation 6 h ahead carries an error of 1.0 KU, while for the first 1 h the error is 0.58 KU only.     

Comparison of the Hill–Siscoe polar cap potential theory with the Weimer and AMIE models

The magnetic storm on November 2004 was characterized by a high solar wind pressure and thus offers a unique opportunity to test the Hill–Siscoe formula (H–S) for the polar cap potential (PCP). To estimate the polar cap potential, we use the Weimer Statistical Convection Model (WCM), and the Assimilative Mapping of Ionospheric Electrodynamics Model (AMIE), based on ingestion of a number of data sets. H–S is in excellent agreement with WCM, and with AMIE during times when DMSP is used in the latter. The implication is that the AMIE conductivity model yields conductivities that are too high by a factor of 2–3. Both H–S and WCM display saturation effects, although WCM is more severe. The two methods track well until an IEF of about 20 mV/m occurs, where H–S continues to increase while WCM levels off. Even at high electric field values, the pressure increases the denominator of the H–S formula by 60%, keeping the potential lower than its saturation value. There are several H–S points above 250 kV, even up to 400 kV, that are not found in WCM and occur right after a rapid transition from B_z north to south. For B_z north, we find evidence for a saturation effect on the PCP at large IEF, little effect as a function of solar wind velocity, and an increase of the PCP with increasing pressure. This seems to rule out viscous interaction but may involve geometric changes in the high-altitude polar cusp that affect recombination there for B_z north.     

Thunderstorm-associated responses in the vertical motion of the mid-latitude F-region ionosphere

Mid-latitude Digisonde Doppler velocities, auroral electrojet (AE) indices and cloud-to-ground (CG) lightning strokes during August 2003–2004 were used to study the perturbations in the F-region vertical drift associated with terrestrial thunderstorms. A superposed epoch analysis (SEA) showed that the F-region vertical drifts V_z had a net descent of ～0.6 m s^?1 peaking ～3 h after lightning. Stronger downward perturbations of up to ～0.9 m s^?1 were observed in the afternoon on the day prior to lightning days. The perturbations were less significant on the day after and insignificant during the remaining intervals up to 144 h on either side of the lightning. The stronger responses on the day before are consistent with causality because the lightning times were merely proxies for the physical mechanisms involved. The actual causes are unclear, but we discuss the possible roles of lightning-induced ionisation enhancements, intense electric fields penetrating upward from electrified clouds, and atmospheric gravity waves (AGWs) radiated from thunderstorms or from the accompanying tropospheric fronts. There is no doubt that the behaviour of the mid-latitude F-region is controlled by the thermospheric winds and the solar wind-magnetosphere electrical generators, but our results suggest that electrified clouds also account for a significant, albeit relatively small component of the ionospheric variability.     

Variations in the geomagnetic field and auroras during the main phase of a large magnetic storm of November 20, 2003

The substorm characteristics during the main phase of a large magnetic storm of November 20, 2003, are studied based on the data of TV observations of auroras and auroral absorption at Tixie Bay station and at the global network of magnetic stations. The contribution of auroral particles, responsible for the emission of discrete auroras, has been estimated based on an analysis of the spatial-time variations in the auroral luminosity intensity. This contribution accounted for ～40% of the total luminous flux, which is approximately twice as large as was previously observed in substorm disturbances. Responses of the solar wind and IMF parameters in substorms and variations in the magnetic indices, characterizing geomagnetic activity in the northern polar cap and ring current (PCN, ASY-H and SYM-H), have been detected. The spatial-time distribution of the equivalent ionospheric currents has been constructed, and the total value of these currents along the meridian has been determined based on the [Popov et al., 2001] method and using the IMAGE magnetic data. It has been obtained that the maximal total equivalent ionospheric current in the premidnight sector (～2000 MLT) leads the minimal value of the SYM-H index by ～1.5 h.     

PFISR observations of strong azimuthal flow bursts in the ionosphere and their relation to nightside aurora

Flow bursts within the ionosphere are the ionospheric signatures of flow bursts in the plasma sheet and have been associated with poleward boundary intensifications (PBIs). Some PBIs extend equatorward from the polar cap boundary, where they can be roughly divided into north–south-aligned and east–west-aligned structures. In this paper, we present two flow burst events observed by the new Poker Flat Advanced Modular Incoherent Scatter Radar (PFISR) in the pre-midnight auroral zone on 28 April 2007, one towards the west and the other towards the east. In both cases, enhanced flows lasted for about 8–10 min with peak velocities exceeding 1500 m/s. The concurrently measured electron density showed that the flow bursts occurred in low conductivity regions. However, near the poleward (equatorward) edge of the westward (eastward) flow burst, strong electron density enhancements were observed in the E region, indicating the presence of discrete auroral arcs. Auroral images from the Polar spacecraft were available at the time of the eastward flow burst and they indicate that this burst was associated with an east–west-aligned auroral structure that connected at later MLT to a north–south structure. In addition, simultaneous precipitating particle energy spectrum measured by the the Defense Meteorological Satellites Program (DMSP) F13 satellite reveals that this auroral structure resulted from mono-energetic electron precipitation associated with a significant field-aligned potential drop. These observations show direct evidence of the relationship between flow bursts, field-aligned currents and auroral intensifications, and suggest that eastward/westward flow bursts are associated with east–west-oriented PBI structures that have extended well within the plasma sheet. This is in contrast to the equatorward-directed flow that has been previously inferred for PBIs near the polar cap boundary and for north–south auroral structures. This paper illustrates the use of the PFISR radar for studying the magnetosphere–ionosphere coupling of flow bursts.     

Propagating plasma patch and associated neutral stream flows

Plasma patches are regions of enhanced ionization that are created in the dayside cusp or equatorward of the cusp in the sunlit hemisphere during northward interplanetary magnetic field. After formation, and a change to a southward interplanetary magnetic field, they drift across the polar cap with the prevailing convection speed. As a plasma patch propagates, charge exchange reactions occur, which lead to the production of both ion and neutral particles throughout the patch. In the region directly above the patch, an upward jet of H⁺ and O⁺ forms. This ion jet, in turn, acts to produce an upward flux of neutral H and O stream particles because of charge exchange reactions between the ion jet and the background neutral atmosphere. A three-dimensional, time-dependent model of the ion and neutral polar winds was used in order to study the evolution of the neutral stream particles that are produced in a ‘representative’ propagating plasma patch, with the anticipation that the neutral stream particles produced by the ion jet would display a distinct signature. However, the outflow of neutral H atoms above a patch is only slightly visible in the simulation due to a continuous outflow flux of H (∼10⁹ cm⁻² s⁻¹) across the entire polar cap. On the other hand, the upward flux of neutral O from the patch is more dependent on both the state of the ionosphere and the amount of heating, with increased upward fluxes over areas where the heating is high. Typically, the upward neutral O streams are predominantly located in the pre-midnight auroral oval.     

Global distribution and climatological features of the 5–6-day planetary waves seen in the SABER/TIMED temperatures (2002–2007)

The paper is focused on the global spatial structure, seasonal and interannual variability of the ～5-day Rossby (W1) and ～6-day Kelvin (E1) waves derived from the SABER/TIMED temperature measurements for 6 full years (January 2002–December 2007). The latitude structure of the ～5-day W1 wave is related to the gravest symmetric wave number 1 Rossby wave. The vertical structure of the ～5-day Rossby wave amplitude consists of double-peaked maxima centred at ～80–90 km and ～105–110 km. This wave has a vertically propagating phase structure from the stratosphere up to 120 km altitude with a mean vertical wavelength of ～50–60 km. The ～6-day E1 wave is an equatorially trapped wave symmetric about the equator and located between 20°N and 20°S. Its seasonal behaviour indicates some equinoctial and June solstice amplifications, while the vertical phase structure indicates that this is a vertically propagating wave between 20–100 km altitudes with a mean vertical wavelength of ～25 km.     

Numerical modeling of the light ion trough and heat balance of the topside ionosphere in quiet geomagnetic conditions

The results from the numerical calculations of the global distribution of topside ionospheric parameters such as H⁺ ions and ion and electron temperatures up to 1500 km height are presented for equinoctial conditions at solar minimum. Calculations are carried out using the Global Self-consistent Model of Thermosphere, Ionosphere and Protonosphere (GSM TIP) developed in WD IZMIRAN, and using a new calculation block for electric fields due to dynamo and of magnetospheric origin. A comparison of two sets of calculations of magnetospheric convection electric field for a given potential difference is carried out, one through polar caps and other through field aligned currents of first zone. It is shown that the distribution of the electric potential obtained through field aligned currents of first zone is more self-consistent than that through polar caps. The light ion trough in H⁺ ions is deeper and occupies larger region for the potential difference through polar cap. For a given potential difference through field aligned current, at 1500 km, the maximum ion temperature is 150 K higher, minimum ion temperature is 200 K lower and maximum electron temperature is 100 K higher than those obtained for the same potential difference through polar caps. It is concluded that for modeling the electric field of magnetospheric origin, it is necessary to use the potential difference through field aligned current of first zone instead of through polar caps.     

Simultaneous observations of small multi-scale structures in an auroral arc

Auroral arcs can develop small-scale distortions known as vortex streets or curls. Other common and somewhat larger spatially periodic distortions are auroral folds. In this event study we present simultaneous wide and narrow field imager observations of a third kind of structuring, on even smaller spatial scales. Boundary undulations, or “ruffs”, have been observed to form on the edge of an auroral arc and they occur superimposed on curls, folds or at times of auroral shear. The undulations typically have wavelengths of less than 3 km and amplitudes of less than 800 m. They are observed to move on the edge of the arc, with velocities of about 11 km/s. These observations, with multi-scale deformations, reveal a much more intricate structuring of auroral arcs than previously found.     

High resolution OSL dating back to MIS 5e in the central Sea of Okhotsk

Marine sediments contain important archives of past ocean and climate changes, but at high latitudes the absence of carbonate has prevented the construction of accurate chronological models. We have begun a study to (1) determine the accuracy of luminescence ages in deep-sea marine sediments, e.g. by comparison with marine oxygen isotope stratigraphy where possible, (2) describe changes in sedimentation rate through time, and (3) test whether it is possible to date back to marine isotope stage 5e (MIS 5e). We show here that optical dating of fine grains of quartz from the central Sea of Okhotsk is able to provide an accurate and precise chronology for the reconstruction of the palaeoceanic and palaeoclimatic environment at our site. The upper 6.5 m of the 18.42 m long core MR0604-PC07A is believed, based on its magnetic susceptibility and the oxygen isotope (δ¹⁸O) records to contain the last ～150 ka. Forty OSL samples were taken from this upper part of the core. The single-aliquot regenerative-dose (SAR) procedure is used for equivalent dose (D_e) determination. The luminescence characteristics of fine-grained quartz (4–11 μm) extracted from the core are described. The OSL signal is dominated by the fast component and a dose recovery test shows that we can accurately measure a known dose given in the laboratory prior to any heat treatment. Dose rates were determined using high-resolution gamma spectrometry, and vary between 0.4 and 1.6 Gy/ka. The OSL ages from this section lie between ～140 ka and ～15 ka and are in very good agreement with the δ¹⁸O stratigraphy up to MIS 5e. A clear change in sedimentation rate is identified: between ～139 and 110 ka, the sedimentation rate was ～0.09 m/ka, but then from ～110 to 15 ka, the sedimentation rate decreases to a constant value of ～0.04 m/ka. Our data confirm that OSL dating using widely distributed fine-grain quartz has great potential for dating deep-sea sediments. Because luminescence methods use clastic materials, they do not depend on the presence of biogenic carbonate. As a result it is now likely that we can establish a chronology in regions of the ocean that were previously undatable.     

Astrochronology of the late Jurassic Kimmeridge Clay (Dorset,England) and implications for Earth system processes

The Late Jurassic Kimmeridge Clay Formation (KCF) is an economically important, organic-rich source rock of Kimmeridgian–Early Tithonian age. The main rock types of the KCF in Dorset, UK, include grey to black laminated shale, marl, coccolithic limestone, and dolostone, which occur with an obvious cyclicity at astronomical timescales. In this study, we examine two high-resolution borehole records (Swanworth Quarry 1 and Metherhills 1) obtained as part of a Rapid Global Geological Events (RGGE) sediment drilling project. Datasets examined were total organic carbon (TOC), and borehole wall microconductivity by Formation Microscanner (FMS). Our intent is to assess the rhythmicity of the KCF with respect to the astronomical timescale, and to discuss the results with respect to other key Late Jurassic geological processes. Power spectra of the untuned data reveal a hierarchy of cycles throughout the KCF with ～ 167 m, ～ 40 m, 9.1 m, 3.8 m and 1.6 m wavelengths. Tuning the ～ 40 m cycles to the 405-kyr eccentricity cycle shows the presence of all the astronomical parameters: eccentricity, obliquity, and precession index. In particular, ～ 100-kyr and 405-kyr eccentricity cycles are strongly expressed in both records. The 405-kyr eccentricity cycle corresponds to relative sea-level changes inferred from sequence stratigraphy. Intervals with elevated TOC are associated with strong obliquity forcing. The 405-kyr-tuned duration of the lower KCF (Kimmeridgian Stage) is 3.47 Myr, and the upper KCF (early part of the Tithonian Stage, elegans to fittoni ammonite zones) is 3.32 Myr. Two other chronologies test the consistency of this age model by tuning ～ 8–10 m cycles to 100-kyr (short eccentricity), and ～ 3–5 m cycles to 36-kyr (Jurassic obliquity). The ‘obliquity-tuned’ chronology resolves an accumulation history for the KCF with a variation that strongly resembles that of Earth's orbital eccentricity predicted for 147.2 Ma to 153.8 Ma. There is evidence for significant non-deposition (up to 1 million years) in the lowermost KCF (baylei–mutabilis zones), which would indicate a Kimmeridgian/Oxfordian boundary age of 154.8 Ma. This absolute calibration allows assignment of precise numerical ages to zonal boundaries, sequence surfaces, and polarity chrons of the lower M-sequence.     

Crustal velocity structure in the southern edge of the Central Alborz (Iran)

Inversion of local earthquake travel times and joint inversion of receiver functions and Rayleigh wave group velocity measurements were used to derive a simple model for the velocity crustal structure beneath the southern edge of the Central Alborz (Iran), including the seismically active area around the megacity of Tehran. The P and S travel times from 115 well-located earthquakes recorded by a dense local seismic network, operated from June to November 2006, were inverted to determine a 1D velocity model of the upper crust. The limited range of earthquake depths (between 2 km and 26 km) prevents us determining any velocity interfaces deeper than 25 km. The velocity of the lower crust and the depth of the Moho were found by joint inversion of receiver functions and Rayleigh wave group velocity data. The resulting P-wave velocity model comprises an upper crust with 3 km and 4 km thick sedimentary layers with P wave velocities (V_p) of ～5.4 and ～5.8 km s^?1, respectively, above 9 km and 8 km thick layers of upper crystalline crust (V_p ～6.1 and ～6.25 km s^?1 respectively). The lower crystalline crust is ～34 km thick (V_p ～ 6.40 km s^?1). The total crustal thickness beneath this part of the Central Alborz is 58 ± 2 km.     

The polar wind: Recent observations

The polar wind is an ambipolar outflow of thermal plasma from the high-latitude ionosphere to the magnetosphere, and it primarily consists of H⁺, He⁺ and O⁺ ions and electrons. Statistical and episodic studies based primarily on ion composition observations on the ISIS-2, DE-1, Akebono and Polar satellites over the past four decades have confirmed the existence of the polar wind. These observations spanned the altitude range from 1000 to ∼50,500 km, and revealed several important features in the polar wind that are unexpected from “classical” polar wind theories. These include the day–night asymmetry in polar wind velocity, which is 1.5–2.0 times larger on the dayside; appreciable O⁺ flow at high altitudes, where the velocity at 5000–10,000 km is of 1–4 km/s; and significant electron temperature anisotropy in the sunlit polar wind, in which the upward-to-downward electron temperature ratio is 1.5–2. These features are attributable to a number of “non-classical” polar wind ion acceleration mechanisms resulting from strong ionospheric convection, enhanced electron and ion temperatures, and escaping atmospheric photoelectrons. The observed polar wind has an averaged ion temperature of ∼0.2–0.3 eV, and a rate of ion velocity increase with altitude that correlates strongly with electron temperature and is greatest at low altitudes (<4000 km for H⁺). The rate of velocity increase below 4000 km is larger at solar minimum than at solar maximum. Above 4000 km, the reverse is the case. This suggests that the dominant polar wind ion acceleration process may be different at low and high altitudes, respectively. At a given altitude, the polar wind velocity is highly variable, and is on average largest for H⁺ and smallest for O⁺. Near solar maximum, H⁺, He⁺, and O⁺ ions typically reach a velocity of 1 km/s near 2000, 3000, and 6000 km, respectively, and velocities of 12, 7, and 4 km/s, respectively, at 10,000 km altitude. Near solar minimum, the velocity of all three species is smaller at high altitudes. Observationally it is not always possible to unambiguously separate an energized “non-polar-wind” ion such as a low-energy “cleft ion fountain” ion that has convected into a polar wind flux tube from an energized “polar-wind” ion that is accelerated locally by “non-classical” polar-wind ion acceleration mechanisms. Significant questions remain on the relative contribution between the cleft ion fountain, auroral bulk upflow, and the topside polar-cap ionosphere to the O⁺ polar wind population at high altitudes, the effect of positive spacecraft charging on the lowest-energy component of the H⁺ polar wind population, and the relative importance of the various classical and non-classical ion acceleration mechanisms. These questions pose several challenges in future polar wind observations: These include measurement of the lowest-energy component in the presence of positive spacecraft potential, definitive determination and if possible active control of the spacecraft potential, definitive discrimination between polar wind and other inter-mixed thermal ion populations, measurement of the three-dimensional ion drift velocity vector and the parallel and perpendicular ion temperatures or the detailed three-dimensional velocity distribution function, and resolution of He⁺ and other minor ion species in the polar wind population.