Showing papers in "Annales Geophysicae in 2006"

Detection efficiency of the VLF World-Wide Lightning Location Network (WWLLN): initial case study

Abstract: . An experimental Very Low Frequency (VLF) World-Wide Lightning Location Network (WWLLN) has been developed through collaborations with research institutions across the world, providing global real-time locations of lightning discharges. As of April 2006, the network included 25 stations providing coverage for much of the Earth. In this paper we examine the detection efficiency of the WWLLN by comparing the locations from this network with lightning location data purchased from a commercial lightning location network operating in New Zealand. Our analysis confirms that WWLLN favours high peak current return stroke lightning discharges, and that discharges with larger currents are observed by more stations across the global network. We then construct a first principles detection efficiency model to describe the WWLLN, combining calibration information for each station with theoretical modelling to describe the expected amplitudes of the VLF sferics observed by the network. This detection efficiency model allows the prediction of the global variation in WWLLN lightning detection, and an estimate of the minimum CG return stroke peak current required to trigger the network. There are strong spatial variations across the globe, primarily due to station density and sensitivity. The WWLLN is currently best suited to study the occurrence and impacts of high peak-current lightning. For example, in 2005 about 12% of the global elve-producing lightning will have been located by the network. Since the lightning-EMP which produce elves has a high mean rate (210 per minute) it has the potential to significantly influence the ionosphere on regional scales.

Temporal and spatial variations in TEC using simultaneous measurements from the Indian GPS network of receivers during the low solar activity period of 2004-2005

Abstract: . With the recent increase in the satellite-based navigation applications, the ionospheric total electron content (TEC) and the L-band scintillation measurements have gained significant importance. In this paper we present the temporal and spatial variations in TEC derived from the simultaneous and continuous measurements made, for the first time, using the Indian GPS network of 18 receivers located from the equator to the northern crest of the equatorial ionization anomaly (EIA) region and beyond, covering a geomagnetic latitude range of 1° S to 24° N, using a 16-month period of data for the low sunspot activity (LSSA) years of March 2004 to June 2005. The diurnal variation in TEC at the EIA region shows its steep increase and reaches its maximum value between 13:00 and 16:00 LT, while at the equator the peak is broad and occurs around 16:00 LT. A short-lived day minimum occurs between 05:00 to 06:00 LT at all the stations from the equator to the EIA crest region. Beyond the crest region the day maximum values decrease with the increase in latitude, while the day minimum in TEC is flat during most of the nighttime hours, i.e. from 22:00 to 06:00 LT, a feature similar to that observed in the mid-latitudes. Further, the diurnal variation in TEC show a minimum to maximum variation of about 5 to 50 TEC units, respectively, at the equator and about 5 to 90 TEC units at the EIA crest region, which correspond to range delay variations of about 1 to 8 m at the equator to about 1 to 15 m at the crest region, at the GPS L1 frequency of 1.575 GHz. The day-to-day variability is also significant at all the stations, particularly during the daytime hours, with maximum variations at the EIA crest regions. Further, similar variations are also noticed in the corresponding equatorial electrojet (EEJ) strength, which is known to be one of the major contributors for the observed day-to-day variability in TEC. The seasonal variation in TEC maximizes during the equinox months followed by winter and is minimum during the summer months, a feature similar to that observed in the integrated equatorial electrojet (IEEJ) strength for the corresponding seasons. In the Indian sector, the EIA crest is found to occur in the latitude zone of 15° to 25° N geographic latitudes (5° to 15° N geomagnetic latitudes). The EIA also maximizes during equinoxes followed by winter and is not significant in the summer months in the LSSA period, 2004–2005. These studies also reveal that both the location of the EIA crest and its peak value in TEC are linearly related to the IEEJ strength and increase with the increase in IEEJ.

Local structure of the magnetotail current sheet: 2001 Cluster observations

Abstract: . Thirty rapid crossings of the magnetotail current sheet by the Cluster spacecraft during July-October 2001 at a geocentric distance of 19 RE are examined in detail to address the structure of the current sheet. We use four-point magnetic field measurements to estimate electric current density; the current sheet spatial scale is estimated by integration of the translation velocity calculated from the magnetic field temporal and spatial derivatives. The local normal-related coordinate system for each case is defined by the combining Minimum Variance Analysis (MVA) and the curlometer technique. Numerical parameters characterizing the plasma sheet conditions for these crossings are provided to facilitate future comparisons with theoretical models. Three types of current sheet distributions are distinguished: center-peaked (type I), bifurcated (type II) and asymmetric (type III) sheets. Comparison to plasma parameter distributions show that practically all cases display non-Harris-type behavior, i.e. interior current peaks are embedded into a thicker plasma sheet. The asymmetric sheets with an off-equatorial current density peak most likely have a transient nature. The ion contribution to the electric current rarely agrees with the current computed using the curlometer technique, indicating that either the electron contribution to the current is strong and variable, or the current density is spatially or temporally structured.

A summary of WIND magnetic clouds for years 1995-2003: model-fitted parameters, associated errors and classifications

Abstract: . Interplanetary magnetic clouds (MCs) have been identified for the first 8.6 years of the WIND mission, and their magnetic field structures have been parameter-fitted by a static, force free, cylindrically-symmetric model (Lepping et al., 1990) with various levels of success. This paper summarizes various aspects of the results of the model fitting by providing: seven estimated model fit-parameter values for each of the 82 MCs found, their objectively determined quality estimates, closest approach vectors (in two coordinate frames), fit-parameter errors for the cases of acceptable quality (50 cases, or 61%), axial magnetic fluxes, axial current densities, and total axial current - as well as some examples of MC profiles for various conditions and "categories" for each case (e.g. Bz: N→S or S→N, etc.). MC quality is estimated from a quantitative consideration of a large set of parameters, such as the chi-squared of the model fit, degree of asymmetry of the B profile, and a comparison of two means of estimating radius. This set of MCs was initially identified by visual inspection of relevant field and plasma data. Each resulting MC candidate is then tested through the use of the MC parameter model, for various adjusted durations to determine the best fit, which helps to refine the boundary-times. The resulting MC set is called Set 1. Another, larger, set (Set 2) of MCs is identified through an automated program whose criteria are based on general MC plasma and field characteristics at 1AU determined through past experience. Set 1 is almost fully contained within Set 2, whose frequency of occurrence better matches that of the sunspot cycle than Set 1. The difference-set (Set 2-Set 1) is referred to as the magnetic cloud-like (MCL) set, whose members do not very well represent good flux ropes through modeling. We present a discussion of how a MC's front boundary is specifically identified in terms of multi-parameter considerations (i.e. any one or more of: increase in B, directional discontinuity, magnetic hole in B, drop in proton plasma beta, B-fluctuation level change, proton temperature drop, etc.), as well as through the application of the flux rope model. Also presented are examples of unusual MCs, as well as some commonly occurring relationships, such as the existence and frequency (approx. 1/2 the time) of upstream interplanetary shocks, and less frequent internal shocks.

Flux Transfer Events: 1. generation mechanism for strong southward IMF

Abstract: . We use a global numerical model of the interaction of the solar wind and the interplanetary magnetic field with Earth's magnetosphere to study the formation process of Flux Transfer Events (FTEs) during strong southward IMF. We find that: (i) The model produces essentially all observational features expected for FTEs, in particular the bipolar signature of the magnetic field BN component, the correct polarity, duration, and intermittency of that bipolar signature, strong core fields and enhanced core pressure, and flow enhancements; (ii) FTEs only develop for large dipole tilt whereas in the case of no dipole tilt steady magnetic reconnection occurs at the dayside magnetopause; (iii) the basic process by which FTEs are produced is the sequential generation of new X-lines which makes dayside reconnection inherently time dependent and leads to a modified form of dual or multiple X-line reconnection; (iv) the FTE generation process in this model is not dependent on specific assumptions about microscopic processes; (v) the average period of FTEs can be explained by simple geometric arguments involving magnetosheath convection; (vi) FTEs do not develop in the model if the numerical resolution is too coarse leading to too much numerical diffusion; and (vii) FTEs for nearly southward IMF and large dipole tilt, i.e., near solstice, should only develop in the winter hemisphere, which provides a testable prediction of seasonal modulation. The semiannual modulation of intermittent FTE reconnection versus steady reconnection is also expected to modulate magnetospheric and ionospheric convection and may thus contribute to the semiannual variation of geomagnetic activity.

Cluster observations of reconnection due to the Kelvin-Helmholtz instability at the dawnside magnetospheric flank

Abstract: . On 3 July 2001, the four Cluster satellites traversed along the dawnside magnetospheric flank and observed large variations in all plasma parameters. The estimated magnetopause boundary normals were oscillating in the z-direction and the normal component of the magnetic field showed systematic 2–3 min bipolar variations for 1 h when the IMF had a small positive bz-component and a Parker-spiral orientation in the x,y-plane. Brief 33 s intervals with excellent deHoffman Teller frames were observed satisfying the Walen relation. Detailed comparisons with 2-D MHD simulations indicate that Cluster encountered rotational discontinuities generated by Kelvin-Helmholtz instability. We estimate a wave length of 6 RE and a wave vector with a significant z-component.

Why is there more ionosphere in January than in July? The annual asymmetry in the F2-layer

Abstract: . Adding together the northern and southern hemisphere values for pairs of stations, the combined peak electron density NmF2 is greater in December-January than in June–July. The same applies to the total height-integrated electron content. This "F2-layer annual asymmetry" between northern and southern solstices is typically 30%, and thus greatly exceeds the 7% asymmetry in ion production due to the annual variation of Sun-Earth distance. Though it was noticed in ionospheric data almost seventy years ago, the asymmetry is still unexplained. Using ionosonde data and also values derived from the International Reference Ionosphere, we show that the asymmetry exists at noon and at midnight, at all latitudes from equatorial to sub-auroral, and tends to be greater at solar minimum than solar maximum. We find a similar asymmetry in neutral composition in the MSIS model of the thermosphere. Numerical computations with the Coupled Thermosphere-Ionosphere-Plasmasphere (CTIP) model give a much smaller annual asymmetry in electron density and neutral composition than is observed. Including mesospheric tides in the model makes little difference. After considering possible explanations, which do not account for the asymmetry, we are left with the conclusion that dynamical influences of the lower atmosphere (below about 30 km), not included in our computations, are the most likely cause of the asymmetry.

On the validity of the ionospheric pierce point (IPP) altitude of 350 km in the Indian equatorial and low-latitude sector

Abstract: . The GPS data provides an effective way to estimate the total electron content (TEC) from the differential time delay of L1 and L2 transmissions from the GPS. The spacing of the constellation of GPS satellites in orbits are such that a minimum of four GPS satellites are observed at any given point in time from any location on the ground. Since these satellites are in different parts of the sky and the electron content in the ionosphere varies both spatially and temporally, the ionospheric pierce point (IPP) altitude or the assumed altitude of the centroid of mass of the ionosphere plays an important role in converting the vertical TEC from the measured slant TEC and vice versa. In this paper efforts are made to examine the validity of the IPP altitude of 350 km in the Indian zone comprising of the ever-changing and dynamic ionosphere from the equator to the ionization anomaly crest region and beyond, using the simultaneous ionosonde data from four different locations in India. From this data it is found that the peak electron density height (hpF2) varies from about 275 to 575 km at the equatorial region, and varies marginally from 300 to 350 km at and beyond the anomaly crest regions. Determination of the effective altitude of the IPP employing the inverse method suggested by Birch et al. (2002) did not yield any consistent altitude in particular for low elevation angles, but varied from a few hundred to one thousand kilometers and beyond in the Indian region. However, the vertical TEC computed from the measured GPS slant TEC for different IPP altitudes ranging from 250 to 750 km in the Indian region has revealed that the TEC does not change significantly with the IPP altitude, as long as the elevation angle of the satellite is greater than 50 degrees. However, in the case of satellites with lower elevation angles (

Survey of large-amplitude flapping motions in the midtail current sheet

Abstract: . We surveyed fast current sheet crossings (flapping motions) over the distance range 10–30 RE in the magnetotail covered by the Geotail spacecraft. Since the local tilts of these dynamic sheets are large and variable in these events, we compare three different methods of evaluating current sheet normals using 4-s/c Cluster data and define the success criteria for the single-spacecraft-based method (MVA) to obtain the reliable results. Then, after identifying more than ~1100 fast CS crossings over a 3-year period of Geotail observations in 1997–1999, we address their parameters, spatial distribution and activity dependence. We confirm that over the entire distance covered and LT bins, fast crossings have considerable tilts in the YZ plane (from estimated MVA normals) which show a preferential appearance of one (YZ kink-like) mode that is responsible for these severe current sheet perturbations. Their occurrence is highly inhomogeneous; it sharply increases with radial distance and has a peak in the tail center (with some duskward shift), resembling the occurrence of the BBFs, although there is no one-to-one local correspondence between these two phenomena. The crossing durations typically spread around 1 min and decrease significantly where the high-speed flows are registered. Based on an AE index superposed epoch study, the flapping motions prefer to appear during the substorm expansion phase, although a considerable number of events without any electrojet and auroral activity were also observed. We also present statistical distributions of other parameters and briefly discuss what could be possible mechanisms to generate the flapping motions.

Thermal, atmospheric and ionospheric anomalies around the time of the Colima M7.8 earthquake of 21 January 2003

Abstract: . The paper examines the possible relationship of anomalous variations of different atmospheric and ionospheric parameters observed around the time of a strong earthquake (Mw 7.8) which occurred in Mexico (state of Colima) on 21 January 2003. These variations are interpreted within the framework of the developed model of the Lithosphere-Atmosphere-Ionosphere coupling. The main attention is focused on the processes in the near ground layer of the atmosphere involving the ionization of air by radon, the water molecules' attachment to the formed ions, and the corresponding changes in the latent heat. Model considerations are supported by experimental measurements showing the local diminution of air humidity one week prior to the earthquake, accompanied by the anomalous thermal infrared (TIR) signals and surface latent heat flux (SLHF) and anomalous variations of the total electron content (TEC) registered over the epicenter of the impending earthquake three days prior to the main earthquake event. Statistical processing of the data of the GPS receivers network, together with various other atmospheric parameters demonstrate the possibility of an early warning of an impending strong earthquake.

The structure of flux transfer events recovered from Cluster data

Abstract: The structure and formation mechanism of a total of five Flux Transfer Events (FTEs), encountered on the equatorward side of the northern cusp by the Cluster spacecraft, with separation of ~5000 km, are studied by applying the Grad-Shafranov (GS) reconstruction technique to the events The technique generates a magnetic field/plasma map of the FTE cross section, using combined magnetic field and plasma data from all four spacecraft, under the assumption that the structure is two-dimensional (2-D) and time-independent The reconstructed FTEs consist of one or more magnetic flux ropes embedded in the magnetopause, suggesting that multiple X-line reconnection was involved in generating the observed FTEs The dimension of the flux ropes in the direction normal to the magnetopause ranges from about 2000 km to more than 1 RE The orientation of the flux rope axis can be determined through optimization of the GS map, the result being consistent with those from various single-spacecraft methods Thanks to this, the unambiguous presence of a strong core field is confirmed, providing evidence for component merging The amount of magnetic flux contained within each flux rope is calculated from the map and, by dividing it by the time interval between the preceding FTE and the one reconstructed, a lower limit of the reconnection electric field during the creation of the flux rope can be estimated; the estimated value ranges from ~011 to ~026 mV m-1, with an average of 019 mV m-1 This can be translated to the reconnection rate of 0038 to 0074, with an average of 0056 Based on the success of the 2-D model in recovering the observed FTEs, the length of the X-lines is estimated to be at least a few RE

Applications of Kalman filters based on non-linear functions to numerical weather predictions

Abstract: . This paper investigates the use of non-linear functions in classical Kalman filter algorithms on the improvement of regional weather forecasts. The main aim is the implementation of non linear polynomial mappings in a usual linear Kalman filter in order to simulate better non linear problems in numerical weather prediction. In addition, the optimal order of the polynomials applied for such a filter is identified. This work is based on observations and corresponding numerical weather predictions of two meteorological parameters characterized by essential differences in their evolution in time, namely, air temperature and wind speed. It is shown that in both cases, a polynomial of low order is adequate for eliminating any systematic error, while higher order functions lead to instabilities in the filtered results having, at the same time, trivial contribution to the sensitivity of the filter. It is further demonstrated that the filter is independent of the time period and the geographic location of application.

A statistical analysis of the location and width of Saturn's southern auroras

Abstract: . A selection of twenty-two Hubble Space Telescope images of Saturn's ultraviolet auroras obtained during 1997–2004 has been analysed to determine the median location and width of the auroral oval, and their variability. Limitations of coverage restrict the analysis to the southern hemisphere, and to local times from the post-midnight sector to just past dusk, via dawn and noon. It is found that the overall median location of the poleward and equatorward boundaries of the oval with respect to the southern pole are at ~14° and ~16° co-latitude, respectively, with a median latitudinal width of ~2°. These median values vary only modestly with local time around the oval, though the poleward boundary moves closer to the pole near noon (~12.5°) such that the oval is wider in that sector (median width ~3.5°) than it is at both dawn and dusk (~1.5°). It is also shown that the position of the auroral boundaries at Saturn are extremely variable, the poleward boundary being located between 2° and 20° co-latitude, and the equatorward boundary between 6° and 23°, this variability contrasting sharply with the essentially fixed location of the main oval at Jupiter. Comparison with Voyager plasma angular velocity data mapped magnetically from the equatorial magnetosphere into the southern ionosphere indicates that the dayside aurora lie poleward of the main upward-directed field-aligned current region associated with corotation enforcement, which maps to ~20°–24° co-latitude, while agreeing reasonably with the position of the open-closed field line boundary based on estimates of the open flux in Saturn's tail, located between ~11° and ~15°. In this case, the variability in location can be understood in terms of changes in the open flux present in the system, the changes implied by the Saturn data then matching those observed at Earth as fractions of the total planetary flux. We infer that the broad (few degrees) diffuse auroral emissions and sub-corotating auroral patches observed in the dayside sector at Saturn result from precipitation from hot plasma sub-corotating in the outer magnetosphere in a layer a few Saturn radii wide adjacent to the magnetopause, probably having been injected either by Dungey-cycle or Vasyliunas-cycle dynamics on the nightside.

From pre-storm activity to magnetic storms: A transition described in terms of fractal dynamics

Abstract: . We show that distinct changes in scaling parameters of the Dst index time series occur as an intense magnetic storm approaches, revealing a gradual reduction in complexity. The remarkable acceleration of energy release – manifested in the increase in susceptibility – couples to the transition from anti-persistent (negative feedback) to persistent (positive feedback) behavior and indicates that the occurence of an intense magnetic storm is imminent. The main driver of the Dst index, the VBSouth electric field component, does not reveal a similar transition to persistency prior to the storm. This indicates that while the magnetosphere is mostly driven by the solar wind the critical feature of persistency in the magnetosphere is the result of a combination of solar wind and internal magnetospheric activity rather than solar wind variations alone. Our results suggest that the development of an intense magnetic storm can be studied in terms of "intermittent criticality" that is of a more general character than the classical self-organized criticality phenomena, implying the predictability of the magnetosphere.

A statistical study on the effect of earthquakes on the ionosphere, based on the subionospheric LF propagation data in Japan

Abstract: . A superimposed epoch analysis has been undertaken, in order to find the correlation of the ionospheric perturbations with seismic activity. We take the wave path from the Japanese LF transmitter (frequency=40 kHz) and an observing station of Kochi (wave path length of 770 km), and a much longer period (of five years) than before, is considered. This subionospheric LF propagation can be called "an integrated measurement" in the sense that any earthquakes in the LF sensitive area just around the great-circle path can influence the observed LF signals, so that we define the "effective magnitude" (Meff) by integrating the total energy from different earthquakes in the sensitive area on a current day and by converting it back into magnitude. A superimposed epoch analysis for the effective magnitude greater than 6.0 has yielded that the ionosphere is definitely disturbed in terms of both amplitude and dispersion, and that these perturbations tend to take place prior to an earthquake. The statistical z-test has also been performed, which has indicated that the amplitude is definitely depleted 2–6 days before the earthquake day and also that the dispersion is very much enhanced during the same period. This statistical study has given strong support to the existence of seismo-ionospheric perturbations for high seismic activity.

Enhanced gravity-wave activity and interhemispheric coupling during the MaCWAVE/MIDAS northern summer program 2002

Abstract: . We present new sensitivity experiments that link observed anomalies of the mesosphere and lower thermosphere at high latitudes during the MaCWAVE/MIDAS summer program 2002 to enhanced planetary Rossby-wave activity in the austral winter troposphere. We employ the same general concept of a GCM having simplified representations of radiative and latent heating as in a previous study by Becker et al. (2004). In the present version, however, the model includes no gravity wave (GW) parameterization. Instead we employ a high vertical and a moderate horizontal resolution in order to describe GW effects explicitly. This is supported by advanced, nonlinear momentum diffusion schemes that allow for a self-consistent generation of inertia and mid-frequency GWs in the lower atmosphere, their vertical propagation into the mesosphere and lower thermosphere, and their subsequent dissipation which is induced by prescribed horizontal and vertical mixing lengths as functions of height. The main anomalies in northern summer 2002 consist of higher temperatures than usual above 82 km, an anomalous eastward mean zonal wind between 70 and 90 km, an altered meridional flow, enhanced turbulent dissipation below 80 km, and enhanced temperature variations associated with GWs. These signals are all reasonably described by differences between two long-integration perpetual model runs, one with normal July conditions, and another run with modified latent heating in the tropics and Southern Hemisphere to mimic conditions that correspond to the unusual austral winter 2002. The model response to the enhanced winter hemisphere Rossby-wave activity has resulted in both an interhemispheric coupling through a downward shift of the GW-driven branch of the residual circulation and an increased GW activity at high summer latitudes. Thus a quantitative explanation of the dynamical state of the northern mesosphere and lower thermosphere during June-August 2002 requires an enhanced Lorenz energy cycle and correspondingly enhanced GW sources in the troposphere, which in the model show up in both hemispheres.

Study of spatial and temporal characteristics of L-band scintillations over the Indian low-latitude region and their possible effects on GPS navigation

Abstract: . The scintillation data (S4-index) at the L-band frequency of 1.575GHz, recorded from a total of 18 GPS receivers installed at different locations in India under the GAGAN project, have provided us with a unique opportunity, for the first time in the Indian region, to make a simultaneous study of spatio-temporal and intensity characteristics of the trans-ionospheric scintillations during the 18-month, low sunspot activity (LSSA) period from January 2004 to July 2005. During this period, the occurrence of scintillations is found to be maximum around the pre-midnight hours of equinox months, with very little activity during the post-midnight hours. No significant scintillation activity is observed during the summer and winter months of the period of observation. The intensity (S4 index) of the scintillation activity is stronger around the equatorial ionization anomaly (EIA) region in the geographic latitude range of 15° to 25° N in the Indian region. These scintillations are often accompanied by the TEC depletions with durations ranging from 5 to 25 min and magnitudes from 5 to 15 TEC units which affect the positional accuracy of the GPS by 1 to 3 m. Further, during the intense scintillation events (S4>0.45≈10 dB), the GPS receiver is found to lose its lock for a short duration of 1 to 4 min, increasing the error bounds effecting the integrity of the SBAS operation. During the present period of study, a total of 395 loss of lock events are observed in the Indian EIA region; this number is likely to increase during the high sunspot activity (HSSA) period, creating more adverse conditions for the trans-ionospheric communications and the GPS-based navigation systems.

Relation between substorm characteristics and rapid temporal variations of the ground magnetic field

Abstract: . Auroral substorms are one of the major causes of large geomagnetically induced currents (GIC) in technological systems. This study deals with different phases of the auroral substorm concerning their severity from the GIC viewpoint. Our database consists of 833 substorms observed by the IMAGE magnetometer network in 1997 (around sunspot minimum) and 1999 (rising phase of the sunspot cycle), divided into two classes according to the Dst index: non-storm (Dst>-40 nT, 696 events) and storm-time ones (Dst

The anomalous warming of summer 2003 in the surface layer of the Central Ligurian Sea (Western Mediterranean)

Abstract: Meteorological and sea temperature data from the ODAS Italia 1 buoy (Ligurian Sea, Western Mediterranean) are used to study the anomalous warming of summer 2003 at sea. The event was related to the record heat wave that interested much of Europe from June to September of that year. The data show that the anomalous warming was prevalently confined to within a few meters below the sea surface. On the contrary, the temperatures in the underlying layers were lower than usual. The limited vertical propagation of heat is ascribed to the high temperature difference that arose between the surface and the deeper layers due to protracted calm weather conditions. The degree of penetration of heat deduced from the observations is consistent with that computed on the basis of an energetic argument, wherein the wind constitutes the sole supply of kinetic energy, while the heating is viewed as the source of potential energy that must be "subtracted" by mixing. The results support the hypothesis that the scanty energy from the wind is mainly responsible for the development of the temperature anomaly at the sea surface.

A global climatology for equatorial plasma bubbles in the topside ionosphere

Abstract: . We have developed a global climatology of equatorial plasma bubble (EPB) occurrence based on evening sector plasma density measurements from polar-orbiting Defense Meteorological Satellite Program (DMSP) spacecraft during 1989-2004. EPBs are irregular plasma density depletions in the post-sunset ionosphere that degrade communication and navigation signals. More than 14400 EPBs were identified in ~134000 DMSP orbits. DMSP observations basically agree with Tsunoda's (1985) hypothesis that EPB rates peak when the terminator is aligned with the Earth's magnetic field, but there are also unpredicted offsets in many longitude sectors. We present an updated climatology for the full database from 1989-2004 along with new plots for specific phases of the solar cycle: maximum 1989-1992 and 1999-2002, minimum 1994-1997, and transition years 1993, 1998, and 2003. As expected, there are significant differences between the climatologies for solar maximum and minimum and between the two solar maximum phases as well. We also compare DMSP F12, F14, F15, and F16 observations at slightly different local times during 2000-2004 to examine local time effects on EPB rates. The global climatologies developed using the DMSP EPB database provide an environmental context for the long-range prediction tools under development for the Communication/Navigation Outage Forecasting System (C/NOFS) mission.

Oscillatory magnetic flux tube slippage in the plasma sheet

Abstract: . Cluster observations in the magnetotail revealed an abundance of strongly inclined current sheets. We determine the magnetic configuration of a particular subset of such phenomena: a series of sheet crossings, having significantly differing inclinations and occurring during quiet conditions. These wave-like variations appear to propagate azimuthally and their magnetic amplitude and magnetic gradient (current density) inside the sheet are proportional to their steepness (degree of inlcination). In spite of significant normal direction changes between neighboring crossings up to 150°, the magnetic field direction inside the neutral sheet remains almost constant. The wavelengths and spatial amplitudes are of the order of 2–5 RE. These observations are interpreted as crossings of a quasi-periodic dynamical structure produced by almost vertical slippage motion of the neighboring magnetic flux tubes in the high-β plasma sheet, rather than large-scale flapping of a stationary structure.

Theoretical aspects of the onset of Indian summer monsoon from perturbed orography simulations in a GCM

Abstract: A theory is proposed to determine the onset of the Indian Summer Monsoon (ISM) in an Atmospheric General Circulation Model (AGCM). The onset of ISM is delayed substantially in the absence of global orography. The impact of orography over different parts of the Earth on the onset of ISM has also been investigated using five additional perturbed simulations. The large difference in the date of onset of ISM in these simulations has been explained by a new theory based on the Surface Moist Static Energy (SMSE) and vertical velocity at the mid-troposphere. It is found that onset occurs only after SMSE crosses a threshold value and the large-scale vertical motion in the middle troposphere becomes upward. This study shows that both dynamics and thermodynamics play profound roles in the onset of the monsoon.

Simulated solar wind plasma interaction with the Martian exosphere: influence of the solar EUV flux on the bow shock and the magnetic pile-up boundary

Abstract: The solar wind plasma interaction with the Martian exosphere is investigated by means of 3-D multi-species hybrid simulations. The influence of the solar EUV flux on the bow shock and the magnetic pile-up boundary is examined by comparing two simulations describing the two extreme states of the solar cycle. The hybrid formalism allows a kinetic description of each ions species and a fluid description of electrons. The ionization processes (photoionization, electron impact and charge exchange) are included self-consistently in the model where the production rate is computed locally, separately for each ionization act and for each neutral species. The results of simulations are in a reasonable agreement with the observations made by Phobos 2 and Mars Global Surveyor spacecraft. The position of the bow shock and the magnetic pile-up boundary is weakly dependent of the solar EUV flux. The motional electric field creates strong asymmetries for the two plasma boundaries.

Stratospheric and mesospheric temperature variations for the quasi-biennial and semiannual (QBO and SAO) oscillations based on measurements from SABER (TIMED) and MLS (UARS)

Abstract: . We present the zonal mean temperature variations for the quasi-biennial oscillation (QBO) and the semiannual oscillation (SAO) based on data from SABER on the TIMED spacecraft (years 2002 to 2004) and from MLS on the UARS mission (1992 to 1994). The SABER measurements provide the rare opportunity to analyze data from one instrument over a wide altitude range (15 to 95 km), while MLS data were taken in the 16 to 55 km altitude range a decade earlier. The results are presented for latitudes from 48° S to 48° N. New results are obtained for the QBO, especially in the upper stratosphere and mesosphere, and at mid-latitudes. At Equatorial latitudes, the QBO amplitudes show local peaks, albeit small, that occur at different altitudes. From about 20 to 40 km, and within about 15° of the Equator, the amplitudes can approach 3.5° K for the stratospheric QBO (SQBO). For the mesospheric QBO (MQBO), we find peaks near 70 km, with temperature amplitudes reaching 3.5° K, and near 85 km, the amplitudes approach 2.5° K. Morphologically, the amplitude and phase variations derived from the SABER and MLS measurements are in qualitative agreement. As a function of latitude, the QBO amplitudes tend to peak at the Equator but then increase again pole-ward of about 15° to 20°. The phase progression with altitude varies more gradually at the Equator than at mid-latitudes. Many of the SAO results presented are also new, in part because measurements were not previously available or were more limited in nature. At lower altitudes near 45 km, within about 15° of the Equator, the temperature amplitudes for the stratospheric SAO (SSAO) reveal a local maximum of about 5° K. At higher altitudes close to the Equator, our results show separate peaks of about 7° K near 75 and 90 km for the mesospheric SAO (MSAO). In the SAO results, significant inter-annual differences are evident, with the amplitudes being largest in 2002 relative to 2003 and 2004. As in the case for the QBO, the SAO temperature amplitudes go through minima away from the Equator, and then increase towards mid latitudes, especially at altitudes above 55 km. We compare our findings with previously published empirical results, and with corresponding results from the numerical spectral model (NSM). Although not a focus of this study, we also show results for the inter-annual variations (which appear to be generated at least in part by the QBO) of the migrating diurnal tide. In the upper mesosphere, their amplitudes can approach 20° K, and they are derived jointly with the zonal-mean components.

Electron velocity distribution and lion roars in the magnetosheath

Abstract: Whistler waves which are termed "lion roars" in the magnetosheath are studied using data obtained by the Spectrum Analyser (SA) of the Spatio-Temporal Analysis of Field Fluctuations (STAFF) experiment aboard Cluster. Kinetic theory is then employed to obtain the theoretical expression for the whistler wave with electron temperature anisotropy which is believed to trigger lion roars in the magnetosheath. This allows us to compare theory and data. This paper for the first time studies the details of the electron velocity distribution function as measured by the Plasma Electron And Current Experiment (PEACE) in order to investigate the underlying causes for the different types of lion roars found in the data. Our results show that while some instances of lion roars could be locally generated, the source of others must be more remote regions of the magnetosheath.

Electric field measurements on Cluster : comparing the double-probe and electron drift techniques

Abstract: The four Cluster satellites each carry two instruments designed for measuring the electric field: a double-probe instrument (EFW) and an electron drift instrument (EDI). We compare data from the two instruments in a representative sample of plasma regions. The complementary merits and weaknesses of the two techniques are illustrated. EDI operations are confined to regions of magnetic fields above 30 nT and where wave activity and keV electron fluxes are not too high, while EFW can provide data everywhere, and can go far higher in sampling frequency than EDI. On the other hand, the EDI technique is immune to variations in the low energy plasma, while EFW sometimes detects significant nongeophysical electric fields, particularly in regions with drifting plasma, with ion energy (in eV) below the spacecraft potential (in volts). We show that the polar cap is a particularly intricate region for the double-probe technique, where large nongeophysical fields regularly contaminate EFW measurments of the DC electric field. We present a model explaining this in terms of enhanced cold plasma wake effects appearing when the ion flow energy is higher than the thermal energy but below the spacecraft potential multiplied by the ion charge. We suggest that these conditions, which are typical of the polar wind and occur sporadically in other regions containing a significant low energy ion population, cause a large cold plasma wake behind the spacecraft, resulting in spurious electric fields in EFW data. This interpretation is supported by an analysis of the direction of the spurious electric field, and by showing that use of active potential control alleviates the situation.

Geomagnetic D st index forecast based on IMF data only

Abstract: . In the past years several operational Dst forecasting algorithms, based on both IMF and solar wind plasma parameters, have been developed and used. We describe an Artificial Neural Network (ANN) algorithm which calculates the Dst index on the basis of IMF data only and discuss its performance for several individual storms. Moreover, we briefly comment on the physical grounds which allow the Dst forecasting based on IMF only.

Spectral width of SuperDARN echoes: measurement, use and physical interpretation

Abstract: . The Doppler velocity and spectral width are two important parameters derived from coherent scatter radar systems. The Super Dual Auroral Radar Network (SuperDARN) is capable of monitoring most of the high latitude region where different boundaries of the magnetosphere map to the ionosphere. In the past, the spectral width, calculated from SuperDARN data, has been used to identify the ionosphere footprints of various magnetosphere boundaries. In this paper we examine the way the spectral width is presently estimated from the radar data and describe several recommendations for improving the algorithm. Using the improved algorithm, we show that typical spectral width values reported in the literature are most probably overestimated. The physical interpretation of the cause of various magnitudes of the spectral width is explored in terms of the diffusion and dynamics of ionospheric plasma irregularities.

Measurements of spectral snow albedo at Neumayer, Antarctica

Abstract: . Spectral albedo in high resolution, from 290 to 1050 nm, has been measured at Neumayer, Antarctica, (70°39' S, 8°15' W) during the austral summer 2003/2004. At 500 nm, the spectral albedo nearly reaches unity, with slightly lower values below and above 500 nm. Above 600 nm, the spectral albedo decreases to values between 0.45 and 0.75 at 1000 nm. For one cloudless case an albedo up to 1.01 at 500 nm could be determined. This can be explained by the larger directional component of the snow reflectivity for direct incidence, combined with a slightly mislevelled sensor and the snow surface not being perfectly horizontal. A possible explanation for an observed decline in albedo is an increase in snow grain size. The theoretically predicted increase in albedo with increasing solar zenith angle (SZA) could not be observed. This is explained by the small range of SZA during albedo measurements, combined with the effect of changing snow conditions outweighing the effect of changing SZA. The measured spectral albedo serves as input for radiative transfer models, describing radiation conditions in Antarctica.

Plasma environment of Titan: a 3-D hybrid simulation study

Abstract: . Titan possesses a dense atmosphere, consisting mainly of molecular nitrogen. Titan's orbit is located within the Saturnian magnetosphere most of the time, where the corotating plasma flow is super-Alfvenic, yet subsonic and submagnetosonic. Since Titan does not possess a significant intrinsic magnetic field, the incident plasma interacts directly with the atmosphere and ionosphere. Due to the characteristic length scales of the interaction region being comparable to the ion gyroradii in the vicinity of Titan, magnetohydrodynamic models can only offer a rough description of Titan's interaction with the corotating magnetospheric plasma flow. For this reason, Titan's plasma environment has been studied by using a 3-D hybrid simulation code, treating the electrons as a massless, charge-neutralizing fluid, whereas a completely kinetic approach is used to cover ion dynamics. The calculations are performed on a curvilinear simulation grid which is adapted to the spherical geometry of the obstacle. In the model, Titan's dayside ionosphere is mainly generated by solar UV radiation; hence, the local ion production rate depends on the solar zenith angle. Because the Titan interaction features the possibility of having the densest ionosphere located on a face not aligned with the ram flow of the magnetospheric plasma, a variety of different scenarios can be studied. The simulations show the formation of a strong magnetic draping pattern and an extended pick-up region, being highly asymmetric with respect to the direction of the convective electric field. In general, the mechanism giving rise to these structures exhibits similarities to the interaction of the ionospheres of Mars and Venus with the supersonic solar wind. The simulation results are in agreement with data from recent Cassini flybys.