bluuebunny/crossref_metadata_2025 huggingface dataset:
Seasonal variations between sampling and classical mean turbulent heat flux estimates in the eastern North Atlantic
Authors: Josey S. A. , Kent E. C. , Taylor P. K.
Published: July 2005

Abstract. The two commonly used statistical measures of the air-sea heat flux, the sampling and classical means, have been compared using hourly reports over a 7-year-period from a weather ship stationed in the NE Atlantic. The sampling mean is the average over all flux estimates in a given period, where individual flux estimates are determined from ship reports of meteorological variables using the well-known bulk formulae. The classical mean is the flux derived by substituting period-averaged values for each of the meteorological variables into the bulk formula (where the averaging period employed is the same as that over which the fluxes are to be determined). Monthly sampling and classical means are calculated for the latent and sensible heat fluxes. The monthly classical mean latent heat flux is found to overestimate the sampling mean by an amount which increases from 1–2 W m–2 in summer to 7 W m–2 in winter, on average, over the 7-year-period. In a given winter month, the excess may be as great as 15 W m–2, which represents about 10% of the latent heat flux. For the sensible heat flux, any seasonal variation between the two means is of the order of 1 W m–2 and is not significant compared to the interannual variation. The discrepancy between the two means for the latent heat flux is shown to arise primarily from a negative correlation between the wind speed and sea-air humidity difference, the effects of which are implicitly included in the sampling method but not in the classical. The influence of the dominant weather conditions on the sign and magnitude of this correlation are explored, and the large negative values that it takes in winter are found to depend on the typical track of the mid-latitude depressions with respect to the position sampled. In conclusion, it is suggested that sampling means should be employed where possible in future climatological studies.
Turbulent heat fluxes by profile and inertial dissipation methods: analysis of the atmospheric surface layer from shipboard measurements during the SOFIA/ASTEX and SEMAPHORE experiments
Authors: Dupuis H. , Weill A. , Katsaros K. , Taylor P. K.
Published: July 2005

Abstract. Heat flux estimates obtained using the inertial dissipation method, and the profile method applied to radiosonde soundings, are assessed with emphasis on the parameterization of the roughness lengths for temperature and specific humidity. Results from the inertial dissipation method show a decrease of the temperature and humidity roughness lengths for increasing neutral wind speed, in agreement with previous studies. The sensible heat flux estimates were obtained using the temperature estimated from the speed of sound determined by a sonic anemometer. This method seems very attractive for estimating heat fluxes over the ocean. However allowance must be made in the inertial dissipation method for non-neutral stratification. The SOFIA/ASTEX and SEMAPHORE results show that, in unstable stratification, a term due to the transport terms in the turbulent kinetic energy budget, has to be included in order to determine the friction velocity with better accuracy. Using the profile method with radiosonde data, the roughness length values showed large scatter. A reliable estimate of the temperature roughness length could not be obtained. The humidity roughness length values were compatible with those found using the inertial dissipation method.
The decay of radar echoes from meteors with particular reference to their use in the determination of temperature fluctuations near the mesopause
Authors: Jones W.
Published: July 2005

Abstract. The rate of decay of a radar echo from an ionised meteor train will be governed by the diffusion coefficient of the plasma and this in turn will depend on the temperature. Very recently the temperature fluctuations near the mesopause have been monitored by this means, by the recording of the decay times of underdense trains. The usual derivation of the precise expression relating the underdense echo decay time to the temperature contains two important assumptions, (i) that the train is created with a Gaussian ionisation profile, and (ii) that kinetic theory may be applied to calculate the diffusion coefficient. We investigate the effect of these assumptions, showing that the first assumption is unnecessary, an underdense backscatter echo decaying exponentially with a decay time equal to λ2/(32π2D), where λ is the wavelength and D the diffusion coefficient, independently of the initial distribution. However, the second assumption is shown to be incorrect, and whereas according to kinetic theory D∝T1/2/ρ, where T and ρ are the atmospheric temperature and density, the correct result is D∝Tρ. This leads to an appreciable correction to the results for the temperature fluctuations.
Some effects of a mean zonal thermocline gradient on planetary equatorial waves
Authors: Pontaud M. , Thual O.
Published: July 2005

Abstract. Planetary equatorial waves are studied with the shallow water equations in the presence of a mean zonal thermocline gradient. The interactions between this gradient and waves are represented by three non-linear terms in the equations: one in the wind-forcing formulation in the x-momentum equation, and two for the advection of mass and divergence of the velocity field in the continuity equation. When the mean gradient is imposed but small, these three (linearized) terms will perturb the behaviour of the equatorial waves. This paper gives a simple analytic treatment of this problem. The equatorial Kelvin mode is first solved with all three contributions, using a Wentzel-Kramers-Brillouin method. The Kelvin mode shows a spatial or/and temporal growth when the thermocline gradient is negative which is the usual situation in the equatorial Pacific ocean (deep thermocline in the west and shallow in the east). The more robust and efficient contribution comes from the advection term. The single effect of the advection of the mean zonal thermocline gradient is then studied for the Kelvin and planetary Rossby modes. The Kelvin mode remains unstable (damped), while the Rossby modes appear damped (unstable) for a negative (positive) thermocline gradient.
Electromagnetic proton cyclotron instability: heating of cool magnetospheric helium ions
Authors: Gary S. Peter , Yin Lin , Winske Dan
Published: July 2005

Abstract. The electromagnetic proton cyclotron anisotropy instability is excited if the hot proton temperature anisotropy, T⊥h/T\\mid\\mid h, is sufficiently large compared to unity, where the subscript h denotes the hot protons and the perpendicular and parallel symbols denote directions relative to the background magnetic field. This instability is important in the outer magnetosphere because it has been shown to lead to an upper bound on T⊥h/T\\mid\\mid h and to cool iron heating. Here one-dimensional initial-value hybrid simulations with spatial variations in the direction of the background magnetic field are used to study this instability in a homogeneous plasma model which represents three ionic constituents of the outer magnetosphere: hot anisotropic protons, cool, initially isotropic protons, and cool, initially isotropic singly ionized helium. These simulations show that the presence of a tenuous helium component does not significantly change the scalings of either the hot proton anisotropy upper bound or the heating of the cool protons. The simulations also show that the helium ion heating rate increases with β\\mid\\midh in contrast to the cool proton energization which decreases with this parameter. The prediction of this homogeneous plasma model, therefore, for cool ions subject to heating by the proton cyclotron instability is that the observed ratio of cool helium temperature to cool proton temperature should increase as β\\mid\\midh increases.
Excitation of VLF quasi-electrostatic oscillations in the ionospheric plasma
Authors: Lundin B. , Krafft C. , Matthieussent G. , Jiricek F. , Shmilauer J. , Triska P.
Published: July 2005

Abstract. A numerical solution of the dispersion equation for electromagnetic waves in a hot magnetized collisionless plasma has shown that, in a current-free ionospheric plasma, the distortion of the electron distribution function reproducing the downward flow of a thermal electron component and the compensating upward flow of the suprathermal electrons, which are responsible for the resulting heat flux, can destabilize quasi-electrostatic ion sound waves. The numerical analysis, performed with ion densities and electron temperature taken from the data recorded by the Interkosmos-24 (IK-24, Aktivny) satellite, is compared with a VLF spectrum registered at the same time on board. This spectrum shows a wide frequency band emission below the local ion plasma frequency. The direction of the electron heat flux inherent to the assumed model of VLF emission generation is discussed
Model studies of the influence of O2 photodissociation parameterizations in the Schumann-Runge bands on ozone related photolysis in the upper atmosphere
Authors: Koppers Gijs A. A. , Murtagh Donal P.
Published: July 2005

Abstract. A new parameterization for atmospheric transmission and O2 photodissociation in the Schumann-Runge band region has been developed and tested with a 1D radiative-photochemical model. The parameterization is based on the O2-column along the line of sight to the Sun and the local temperature. Line-by-line calculations have served as a benchmark for testing this method and several other, commonly used, parameterizations. The comparisons suggest that differences between the line-by-line calculations and currently accepted parameterizations can be reduced significantly by using the new method, particularly at large solar zenith angles. The production rate of O-atoms computed with this method shows less than 6% deviation compared to the line-by-line calculations at any altitude, all solar zenith angles and in all seasons. The largest errors are found toward the shorter wavelengths in the Schumann-Runge region at low altitudes. Transmittance is approximated to better than 4% at any altitude and/or solar zenith angle. The total O-production rate above 20 km is approximated to better than 2%. The new parameterization is easily implemented in existing photochemical models and in many cases it may simply replace the existing algorithm. The computational effort exceeds that of other parameterizations but in view of the total computation time needed for the actual calculation of the parameterized Schumann-Runge bands this should not lead to significant performance degeneration. The first 14 coefficients of the parameterization are included in this study. Both the complete sets of coefficients and a simple algorithm can be obtained by contacting the authors. A photochemical model study shows the largest effect of the parameterization method is on odd hydrogen concentrations. Subsequent interaction with an odd oxygen family causes differences in the ozone concentrations between the different parameterizations of more than 10% at selected altitudes. Although it is already established that deficiencies in the treatment of Schumann-Runge band absorption are unlikely to explain the current underestimation of ozone concentration at the stratopause in a variety of photochemical models, this study does show that the choice of parameterization has a large impact on the accuracy of the results at large solar zenith angles and in different seasons.
Effects of alternative cloud radiation parameterizations in a general circulation model
Authors: Lee Wan-Ho , Somerville Richard C. J.
Published: July 2005

Abstract. Using the National Center for Atmospheric Research (NCAR) general circulation model (CCM2), a suite of alternative cloud radiation parameterizations has been tested. Our methodology relies on perpetual July integrations driven by ±2 K sea surface temperature forcing. The tested parameterizations include relative humidity based clouds and versions of schemes involving a prognostic cloud water budget. We are especially interested in testing the effect of cloud optical thickness feedbacks on global climate sensitivity. All schemes exhibit negative cloud radiation feedbacks, i.e., cloud moderates the global warming. However, these negative net cloud radiation feedbacks consist of quite different shortwave and longwave components between a scheme with interactive cloud radiative properties and several schemes with specified cloud water paths. An increase in cloud water content in the warmer climate leads to optically thicker middle- and low-level clouds and in turn negative shortwave feedbacks for the interactive radiative scheme, while a decrease in cloud amount leads to a positive shortwave feedback for the other schemes. For the longwave feedbacks, a decrease in high effective cloudiness for the schemes without interactive radiative properties leads to a negative feedback, while no distinct changes in effective high cloudiness and the resulting feedback are exhibited for the scheme with interactive radiative properties. The resulting magnitude of negative net cloud radiation feed-back is largest for the scheme with interactive radiative properties. Even though the simulated values of cloud radiative forcing for the present climate using this method differ most from the observational data, the approach shows great promise for the future.
The use of interplanetary scintillation maps in the prediction of geomagnetic activity
Authors: Lucek E. A. , Clark T. D. G.
Published: July 2005

Abstract. Interplanetary scintillation (IPS), the twinkling of small angular diameter radio sources, arises from the interaction of the signal with small-scale plasma irregularities in the solar wind. The technique may be used to sense remotely the near-Earth heliosphere and has potential for tracking large-scale interplanetary disturbances from close to the Sun to the Earth. Such observations might be useful within routine geomagnetic forecasts, and we use data from the Mullard Radio Astronomy Observatory to test this suggestion. A forecast was based on the visual evaluation of each daily map. If an IPS event was observed then we proposed that any associated geomagnetic activity would occur either on that day, or during the following two days. We consider the success of these forecasts in predicting days when either an SSC/SI or an Ap value exceeding 30 were recorded. The identification of IPS events is necessarily subjective and so two observers compiled independent events lists, and the results were compared. Approximately half of the IPS events in each list were followed by a geomagnetic signature but comparison of the two lists showed that different days were being chosen. We also found that the forecasts had very high false alarm rates. Since IPS is sensitive to a volume we did not expect all events to be associated with a geomagnetic signature. However, the technique failed to forecast a large proportion of geomagnetic events and the association between IPS events and geomagnetic activity is not much better than would be expected by chance. Comparing the IPS forecasts with forecasts of Ap released by the Space Environment Services Center (SESC) we found that SESC correctly predicted a similar proportion of days when Ap\\geq30, but that the performance was significantly better than would be expected by chance, and had a much lower false alarm rate. We conclude that these IPS data cannot be used alone to produce reliable geomagnetic activity forecasts.
Substorm correlated absorption on a 3200 km trans-auroral HF propagation path
Authors: Milan S. E. , Jones T. B. , Lester M. , Warrington E. M. , Reeves G. D.
Published: July 2005

Abstract. A high-frequency transmitter located at Clyde River, NWT, Canada, and a receiver located near Boston, USA, provide a 3200 km trans-auroral, near-meridional propagation path over which the propagation characteristics have been measured. Out of the fourteen frequencies in the HF band sampled every hour for the duration of the experimental campaign (16 January–8 February 1989), the signal level measurements of 6.800 MHz transmissions were selected in order to determine the extent and occurrence of auroral absorption. The median level of auroral absorption along the path is found to increase with geomagnetic activity, quantified by the index Kp, with the increase being greater in the post-midnight sector than in the pre-midnight sector. This asymmetric behaviour is attributed to the precipitation of high energy electrons into the midnight and morning sector auroral D region. The measured diurnal variation in the median level of absorption is consistent with previous models describing the extent and magnitude of auroral absorption and electron precipitation. Individual substorms, identified from geosynchronous satellite data, are found to cause short-lived absorption events in the HF signal level of ~30 dB at 6.800 MHz. The occurrence of substorm correlated auroral absorption events is confined to the midnight and morning sectors, consistent with the location of the electron precipitation. The magnitude of absorption is related to the magnetotail stress during the substorm growth phase and the magnetotail relaxation during the substorm expansion phase onset. The absorption magnitude and the occurrence of substorms during the period of the campaign increase at times of high Kp , leading to an increase in median auroral absorption during disturbed periods.
A numerical model of the ionosphere, including the E-region above EISCAT
Authors: Diloy P.-Y. , Robineau A. , Lilensten J. , Blelly P.-L. , Fontanari J.
Published: July 2005

Abstract. It has been previously demonstrated that a two-ion (O+ and H+) 8-moment time-dependent fluid model was able to reproduce correctly the ionospheric structure in the altitude range probed by the EISCAT-VHF radar. In the present study, the model is extended down to the E-region where molecular ion chemistry (NO+ and O+2, essentially) prevails over transport; EISCAT-UHF observations confirmed previous theoretical predictions that during events of intense E×B induced convection drifts, molecular ions (mainly NO+) predominate over O+ ions up to altitudes of 300 km. In addition to this extension of the model down to the E-region, the ionization and heating resulting from both solar insolation and particle precipitation is now taken into account in a consistent manner through a complete kinetic transport code. The effects of E×B induced convection drifts on the E- and F-region are presented: the balance between O+ and NO+ ions is drastically affected; the electric field acts to deplete the O+ ion concentration. The [NO+]/[O+] transition altitude varies from 190 km to 320 km as the perpendicular electric field increases from 0 to 100 mV m-1. An interesting additional by-product of the model is that it also predicts the presence of a noticeable fraction of N+ ions in the topside ionosphere in good agreement with Retarding Ion Mass Spectrometer measurements onboard Dynamic Explorer.
Mechanisms of the electron density depletion in the SAR arc region
Authors: Pavlov A. V.
Published: July 2005

Abstract. This study compares the measurements of electron density and temperature and the integral airglow intensity at 630 nm in the SAR arc region and slightly south of this (obtained by the Isis 2 spacecraft during the 18 December 1971 magnetic storm), with the model results obtained using the time dependent one-dimensional mathematical model of the Earth\\'s ionosphere and plasmasphere. The explicit expression in the third Enskog approximation for the electron thermal conductivity coefficient in the multicomponent mixture of ionized gases and a simplified calculation method for this coefficient presents an opportunity to calculate more exactly the electron temperature and density and 630 nm emission within SAR arc region are used in the model. Collisions between N2 and hot thermal electrons in the SAR arc region produce vibrationally excited nitrogen molecules. It appears that the loss rate of O+(4S) due to reactions with the vibrationally excited nitrogen is enough to explain electron density depression by a factor of two at F-region heights and the topside ionosphere density variations within the SAR arc if the erosion of plasma within geomagnetic field tubes, during the main phase of the geomagnetic storm and subsequent filling of geomagnetic tubes during the recovery phase, are considered. To explain the disagreement by a factor 1.5 between the observed and modeled SAR arc electron densities an additional plasma drift velocity ~–30 m s–1 in the ion continuity equations is needed during the recovery phase. This additional plasma drift velocity is likely caused by the transition from convecting to corotating flux tubes on the equatorward wall of the trough. The electron densities and temperatures and 630 nm integral intensity at the SAR arc and slightly south of this region as measured for the 18 December 1971 magnetic storm were correctly described by the model without perpendicular electric fields. Within this model framework the effect of the perpendicular electric field ~100 mv m–1 with a duration ~1 h on the SAR arc electron density profiles was found to be large. However, this effect is small if ~1–2 h have passed after the electric field was set equal to zero.
Rocket observation of atomic oxygen and night airglow: Measurement of concentration with an improved resonance fluorescence technique
Authors: Kita K. , Imamura T. , Iwagami N. , Morrow W. H. , Ogawa T.
Published: July 2005

Abstract. An improved resonant fluorescence instrument for measuring atomic oxygen concentration was developed to avoid the Doppler effect and the aerodynamic shock effect due to the supersonic motion of a rocket. The shock effect is reduced by adopting a sharp wedge-shaped housing and by scanning of the detector field of view to change the distance between the scattering volume and the surface of the housing. The scanning enables us to determine absolute values of atomic oxygen concentration from relative variation of the scattered light signal due to the self-absorption. The instrument was calibrated in the laboratory, and the numerical simulation reproduced the calibration result. Using the instrument, the altitude profile of atomic oxygen concentration was observed by a rocket experiment at Uchinoura (31°N) on 28 January 1992. The data obtained from the rocket experiment were not perfectly free from the shock effect, but errors due to the effect were reduced by the data analysis procedure. The observed maximum concentration was 3.8× 1011 cm–3 at altitudes around 94 km. The systematic error is estimated to be less than ±0.7×1011 cm–3 and the relative random error is less than±0.07× 1011 cm–3at the same altitudes. The altitude profile of the OI 557.7-nm airglow was also observed in the same rocket experiment. The maximum volume emission rate was found to be 150 photons cm–3 s–1 at 94 km. The observed altitude profiles are compared with the MSIS model and other in situ observations.
First results from ground-based daytime optical investigation of the development of the equatorial ionization anomaly
Authors: Pallam Raju D. , Sridharan R. , Gurubaran S. , Raghavarao R.
Published: July 2005

Abstract. A meridional scanning OI 630.0-nm dayglow photometer was operated from Ahmedabad (17.2°N dip lat.) scanning a region towards the south in the upper atmosphere extending over ~5° in latitude from 10.2°N to 15.2°N dip latitude. From the spatial and temporal variabilities of the dayglow intensity in the scanning region we show for the first time, evidence for the passage of the crest of the equatorial ionization anomaly (EIA) in the daytime by means of a ground-based optical technique. The relationship between the daytime eastward electric field over the dip equator in the same longitude zone as inferred from the equatorial electrojet strength and the evolutionary pattern of EIA is clearly demonstrated. The latter as inferred from the dayglow measurements is shown to be consistent with our present understanding of the electrodynamical processes in the equatorial region. The present results reveal the potential of this ground-based optical technique for the investigation of ionospheric/thermospheric phenomena with unprecedented spatial and temporal resolution.
Multifractal scaling of the kinetic energy flux in solar wind turbulence
Authors: Marsch E. , Tu C.-Y. , Rosenbauer H.
Published: July 2005

Abstract. The geometrical and scaling properties of the energy flux of the turbulent kinetic energy in the solar wind have been studied. Using present experimental technology in solar wind measurements we cannot directly measure the real volumetric dissipation rate, ε(t), but are constrained to represent it by its surrogate the energy flux near the dissipation range at the proton gyro scale. There is evidence for the multifractal nature of the so defined dissipation field ε(t), a result derived from the scaling exponents of its statistical moments. The generalized dimension Dq has been determined and reveals that the dissipation field has a multifractal structure, which is not compatible with a scale-invariant cascade. The related multifractal spectrum α has been estimated for the first time for MHD turbulence in the solar wind. Its features resemble those obtained for turbulent fluids and other nonlinear multifractal systems. The generalized dimension Dq can for turbulence in high-speed streams be fitted well by the functional dependence of the p-model with a comparatively large parameter p1=0.87, indicating a strongly intermittent multifractal energy cascade. The experimental value for Dp/3 used in the scaling exponent s(p) of the velocity structure function gives an exponent that can describe some of the observations. The scaling exponent μ of the autocorrelation function of ε (t) has also been directly evaluated, being 0.37. Finally, the mean dissipation rate was determined, which could be used in solar wind heating models.
Seasonal variations in the occurrence of geomagnetic storms
Authors: Taylor John R. , Lester Mark , Yeoman Timothy K.
Published: July 2005

Abstract. Seasonal variations in the onset of magnetic storms are investigated. For the purposes of this study storms have been defined as events in which Dst falls below –50 nT for at least four consecutive hours. The storms have been classified as either storm sudden commencements (SSCs; storms initiated by a sudden commencement) or as storm gradual commencements (SGCs; all other storms). It is found that the semi-annual variation of magnetic activity is reflected in the occurrence statistics of SGC events only, indicative that the solar wind origin is different for SSCs and SGCs. It is suggested that the heliospheric latitude model of seasonal magnetic activity is relatively ineffective in modulating the previously observed seasonal variations in the occurrence of magnetic storms.
Observations of gravity waves in the upper and lower stratosphere by lidar and ozonesondes
Authors: Mitchell N. J. , McDonald A. J. , Reid S. J. , Price J. D.
Published: July 2005

Abstract. Lidar observations of Rayleigh backscatter and temperature profiles measured by ozonesondes have been used to investigate gravity waves in the upper and lower stratosphere, respectively, over Aberystwyth (52.4°N, 4.1°W). Both data sets have been used to investigate the vertical wavenumber spectrum of the wave field at high wavenumbers. Similar analytic techniques applied to each data set enable direct comparison of the spectra. The possibility of laminar structures generated by differential advection contaminating gravity-wave fields deduced from temperature and/or density measurement is discussed and the behaviour of the wave-field mean potential energy reveals a seasonal cycle throughout the stratosphere with a late winter maximum and a summer minimum.
Effect of atmospheric gases, surface albedo and cloud overlap on the absorbed solar radiation
Authors: Sinha Ashok
Published: July 2005

Abstract. Recent studies have provided new evidence that models may systematically underestimate cloud solar absorption compared to observations. This study extends previous work on this "absorption anomaly'' by using observational data together with solar radiative transfer parameterisations to calculate fs (the ratio of surface and top of the atmosphere net cloud forcings) and its latitudinal variation for a range of cloud types. Principally, it is found that (a) the zonal mean behaviour of fs varies substantially with cloud type, with the highest values obtained for low clouds; (b) gaseous absorption and scattering can radically alter the pattern of the variation of fs with latitude, but gaseous effects cannot in general raise fs to the level of around 1.5 as recently determined; (c) the importance of the gaseous contribution to the atmospheric ASR is such that whilst fs rises with surface albedo, the net cloud contribution to the atmospheric ASR falls; (d) the assumed form of the degree of cloud overlap in the model can substantially affect the cloud contribution to the atmospheric ASR whilst leaving the parameter fs largely unaffected; (e) even large uncertainties in the observed optical depths alone cannot account for discrepancies apparent between modelled and newly observed cloud solar absorption. It is concluded that the main source of the anomaly may derive from the considerable uncertainties regarding impure droplet microphysics rather than, or together with, uncertainties in macroscopic quantities. Further, variable surface albedos and gaseous effects may limit the use of contemporaneous satellite and ground-based measurements to infer the cloud solar absorption from the parameter fs.
Flux correction and the simulation of changing climate
Authors: Schneider Edwin K.
Published: July 2005

Abstract. A flaw is pointed out in the manner in which flux correction is currently applied to coupled atmosphere-ocean general circulation models. If a transient climate simulation were carried out using perfect initial data and a perfect model, then a perfect simulation would be made. However, if the model were flux corrected so that it is in equilibrium for current conditions, according to current practice, then errors in the simulation would grow initially to a finite amplitude and persist indefinitely. Larger errors would be produced by a simulation with the flux corrected model beginning from pre-industrial conditions than by a simulation beginning from current conditions. An example with a simple linear model is constructed to illustrate this point, and the relationship to the cold start problem is demonstrated. An optimal flux correction for the simple example is shown to be one which would eliminate the error in the current climate from a transient simulation begun sufficiently far in the past.
Generation and evolution of interplanetary slow shocks
Authors: Wu C.-C. , Wu S. T. , Dryer M.
Published: July 2005

Abstract. It is well known that most MHD shocks observed within 1 AU are MHD fast shocks. Only a very limited number of MHD slow shocks are observed within 1 AU. In order to understand why there are only a few MHD slow shocks observed within 1 AU, we use a one-dimensional, time-dependent MHD code with an adaptive grid to study the generation and evolution of interplanetary slow shocks (ISS) in the solar wind. Results show that a negative, nearly square-wave perturbation will generate a pair of slow shocks (a forward and a reverse slow shock). In addition, the forward and the reverse slow shocks can pass through each other without destroying their characteristics, but the propagating speeds for both shocks are decreased. A positive, square-wave perturbation will generate both slow and fast shocks. When a forward slow shock (FSS) propagates behind a forward fast shock (FFS), the former experiences a decreasing Mach number. In addition, the FSS always disappears within a distance of 150R⊙ (where R⊙ is one solar radius) from the Sun when there is a forward fast shock (with Mach number ≥1.7) propagating in front of the FSS. In all tests that we have performed, we have not discovered that the FSS (or reverse slow shock) evolves into a FFS (or reverse fast shock). Thus, we do not confirm the FSS-FFS evolution as suggested by Whang (1987).
Calculating the global mass exchange between stratosphere and troposphere
Authors: Grewe V. , Dameris M.
Published: July 2005

Abstract. Large-scale cross-tropopause mass fluxes are diagnosed globally from 1979 to 1989 for Northern Hemisphere winter conditions (December, January, and February). Results of different methods of approaches with regard to the definition of the tropopause and the way to calculate the mass fluxes are compared and discussed. The general pattern of the mass exchange from the tropopause into the stratosphere and vice versa agrees fairly well when using different methods, but the absolute values can differ up to 100%. An inspection of the temporal development of the mass fluxes for solstice conditions indicates a complex picture. Whereas a permanent significant downward flux from the stratosphere into the troposphere is detected for latitude regions nearly between 25°N and 40°N and between 30°S and 50°S (initiated by the poleward branches of the Hadley cells), a non-uniform behaviour is observed at higher latitude bands. Periods of strong mass exchange from the troposphere into the stratosphere are disrupted by periods of an opposite mass exchange. A comparison of the stratoshere-troposphere (ST) exchange with the exchange at higher altitudes through surfaces, quasi-parallel to the tropopause, excludes a general connection. Only a few strong upward directed ST mass exchange events have counterparts at higher altitudes. The composition of the stratosphere may be influenced directly by the ST exchange only in a thin layer above the tropopause.
Simulations of stratospheric sudden warmings in the Berlin troposphere-stratosphere-mesosphere GCM
Authors: Erlebach P. , Langematz U. , Pawson S.
Published: July 2005

Abstract. Stratospheric sudden warming events in the Northern Hemisphere of the Berlin TSM GCM are investigated. In about 50% of the simulated years (13 out of 28), major midwinter warmings occur. This agrees well with observations but, whereas real events tend to occur approximately every second season, those in the model are clustered, most of them occur in the period between years 15/16 and years 24/25. In most other years, minor warming events take place. The warming events are found earlier in the winter than in reality. Many of the observed characteristics of warming events are well captured by the model: pulses of wave activity propagate out of the troposphere; these transient events force the zonal-mean zonal wind in the stratosphere and coincide with increases of the temperature at the North Pole and cooling at low levels in the tropics; temperature changes of opposite sign are modelled at higher levels. Synoptically, the modelled stratosphere evolves quite realistically before the warmings: the cyclonic vortex is displaced from the Pole by an amplifying anticyclone. After minor warmings, the stratosphere remains too disturbed as the cyclonic centre does not return to the North Pole as quickly as in reality. In the aftermath of major warmings the cyclonic vortex is not fully eroded and the anticyclonic circulation does not develop properly over the Pole; furthermore, the wintertime circulation is not properly restored after the event.
Interannual variability of precipitable water
Authors: Kane R. P.
Published: July 2005

Abstract. The 12-month running means of the surface-to-500 mb precipitable water obtained from analysis of radiosonde data at seven selected locations showed three types of variability viz: (1) quasi-biennial oscillations; these were different in nature at different latitudes and also different from the QBO of the stratospheric tropical zonal winds; (2) decadal effects; these were prominent at middle and high latitudes and (3) linear trends; these were prominent at low latitudes, up trends in the Northern Hemisphere and downtrends in the Southern Hemisphere.
On the scattering behaviour of bullet-rosette and bullet-shaped ice crystals
Authors: Strauss B.
Published: July 2005

Abstract. The scattering behaviour of bullet-rosette and bullet-shaped ice particles is investigated using model calculations (ray tracing method) with special emphasis on the asymmetry factor g. Because the variability of the geometrical features of these particles is very large, some representative shapes are used in the calculations. The model is based on geometrical optics, and particles are assumed to be oriented randomly; a wavelength of 0.56 μm is considered; absorption is neglected. The scattering behaviour of bullet rosettes is compared to that of single branches out of the bullet rosette. It turns out that there are slight differences in the asymmetry factor values, depending on the lengths of the branches (∆g~0.02) and on the angles between the branches (∆g~0.01). Bullets show some special features in their phase functions due to the pyramid. The length of the particle influences the asymmetry factor (∆g~0.10), as does the shape of the pyramid (∆g~0.07). The influence of the pyramidal shape decreases with increasing particle length. Bullets were compared to hexagonally shaped columns. This was done for two columns, one as long as the columnar part of the bullet (length without pyramid), and one for a column as long as the bullet including the pyramid. Asymmetry factor values of bullets with a pyramidal angle of 28° deviate less than ∆g~0.01 from the range given by the two values of the columns.
The convection electrojet and the substorm electrojet
Authors: Kamide Y. , Nakamura R.
Published: July 2005

Abstract. Enhancements in the auroral electrojets associated with magnetospheric substorms result from those in either the electric field or the ionospheric conductivities, or both. Their relative importance varies significantly, even during a single substorm, depending on the location as well as on the substorm phases. It is predicted that different parts of the electrojets tend to respond in different ways to substorm activity. The unprecedented, unique opportunity for CLUSTER spacecraft observations of electric/magnetic fields and precipitating particles, combined with radar measurements of ionospheric quantities and with ground magnetometers, will provide us with crucial information regarding the physical nature of the separation between the "electric field-dominant'' and "conductivity-dominant'' auroral electrojets. This study also discusses the implications of these two auroral-electrojet components in terms of solar wind-magnetosphere-ionosphere interactions.
Dispersive O+ conics observed in the plasma-sheet boundary layer with CRRES/LOMICS during a magnetic storm
Authors: Wüest M. , Young D. T. , Thomsen M. F. , Barraclough B. L. , Singer H. J. , Anderson R. R.
Published: July 2005

Abstract. We present initial results from the Low-energy magnetospheric ion composition sensor (LOMICS) on the Combined release and radiation effects satellite (CRRES) together with electron, magnetic field, and electric field wave data. LOMICS measures all important magnetospheric ion species (H+, He++, He+, O++, O+) simultaneously in the energy range 60 eV to 45 keV, as well as their pitch-angle distributions, within the time resolution afforded by the spacecraft spin period of 30 s. During the geomagnetic storm of 9 July 1991, over a period of 42 min (0734 UT to 0816 UT) the LOMICS ion mass spectrometer observed an apparent O+ conic flowing away from the southern hemisphere with a bulk velocity that decreased exponentially with time from 300 km/s to 50 km/s, while its temperature also decreased exponentially from 700 to 5 eV. At the onset of the O+ conic, intense low-frequency electromagnetic wave activity and strong pitch-angle scattering were also observed. At the time of the observations the CRRES spacecraft was inbound at L~7.5 near dusk, magnetic local time (MLT), and at a magnetic latitude of –23°. Our analysis using several CRRES instruments suggests that the spacecraft was skimming along the plasma sheet boundary layer (PSBL) when the upward-flowing ion conic arrived. The conic appears to have evolved in time, both slowing and cooling, due to wave-particle interactions. We are unable to conclude whether the conic was causally associated with spatial structures of the PSBL or the central plasma sheet.
Decay of the Dst field of geomagnetic disturbance after substorm onset and its implication to storm-substorm relation
Authors: Iyemori T. , Rao D. R. K.
Published: July 2005

Abstract. In order to investigate the causal relationship between magnetic storms and substorms, variations of the mid-latitude geomagnetic indices, ASY (asymmetric part) and SYM (symmetric part), at substorm onsets are examined. Substorm onsets are defined by three different phenomena; (1) a rapid increase in the mid-latitude asymmetric-disturbance indices, ASY-D and ASY-H, with a shape of so-called `mid-latitude positive bay\\'; (2) a sharp decrease in the AL index; (3) an onset of Pi2 geomagnetic pulsation. The positive bays are selected using eye inspection and a pattern-matching technique. The 1-min-resolution SYM-H index, which is essentially the same as the hourly Dst index except in terms of the time resolution, does not show any statistically significant development after the onset of substorms; it tends to decay after the onset rather than to develop. It is suggested by a simple model calculation that the decay of the magnetospheric tail current after substorm onset is responsible for the decay of the Dst field. The relation between the IMF southward turning and the development of the Dst field is re-examined. The results support the idea that the geomagnetic storms and substorms are independent processes; that is, the ring-current development is not the result of the frequent occurrence of substorms, but that of enhanced convection caused by the large southward IMF. A substorm is the process of energy dissipation in the magnetosphere, and its contribution to the storm-time ring-current formation seems to be negligible. The decay of the Dst field after a substorm onset is explained by a magnetospheric energy theorem.
Numerical modelling of pollutant dispersion in sea breeze conditions
Authors: Camps J. , Massons J. , Soler M. R.
Published: July 2005

Abstract. The characteristics of air pollution in Tarragona (Spain) were investigated. Tarragona has an important petrochemical industry in a coastal region with a complex terrain. The numerical study was made in sea breeze conditions with a three-dimensional mesoscale model. Temporal and spatial variations of the wind fields have been used in the Eulerian equation for a non-reactive pollutant. The results of this study reveal the complexity of the dispersion patterns due to the combined effects of the sea breeze circulation and the orography. This work presents a comparison between the model output and the observed wind data by sodar and surface wind measurements. The evaluation shows that the model is capable of providing very realistic wind fields within this domain.
Predicting geomagnetic storms from solar-wind data using time-delay neural networks
Authors: Gleisner H. , Lundstedt H. , Wintoft P.
Published: July 2005

Abstract. We have used time-delay feed-forward neural networks to compute the geomagnetic-activity index Dst one hour ahead from a temporal sequence of solar-wind data. The input data include solar-wind density n, velocity V and the southward component Bz of the interplanetary magnetic field. Dst is not included in the input data. The networks implement an explicit functional relationship between the solar wind and the geomagnetic disturbance, including both direct and time-delayed non-linear relations. In this study we especially consider the influence of varying the temporal size of the input-data sequence. The networks are trained on data covering 6600 h, and tested on data covering 2100 h. It is found that the initial and main phases of geomagnetic storms are well predicted, almost independent of the length of the input-data sequence. However, to predict the recovery phase, we have to use up to 20 h of solar-wind input data. The recovery phase is mainly governed by the ring-current loss processes, and is very much dependent on the ring-current history, and thus also the solar-wind history. With due consideration of the time history when optimizing the networks, we can reproduce 84% of the Dst variance.
On the 18-day quasi-periodic oscillation in the ionosphere
Authors: Altadill D.
Published: July 2005

Abstract. The presence and persistence of an 18-day quasi-periodic oscillation in the ionospheric electron density variations were studied. The data of lower ionosphere (radio-wave absorption at equivalent frequency near 1 MHz), middle and upper ionosphere (critical frequencies f0E and f0F2) for the period 1970–1990 have been used in the analysis. Also, solar and geomagnetic activity data (the sunspot numbers Rz and solar radio flux F10.7 cm, and aN index respectively) were used to compare the time variations of the ionospheric with the solar and geomagnetic activity data. Periodogram, complex demodulation, auto- and cross-correlation analysis have been used. It was found that 18-day quasi-periodic oscillation exists and persists in the temporal variations of the ionospheric parameters under study with high level of correlation and mean period of 18–19 days. The time variation of the amplitude of the 18-day quasi-periodic oscillation in the ionosphere seems to be modulated by the long-term solar cycle variations. Such oscillations exist in some solar and geomagnetic parameters and in the planetary wave activity of the middle atmosphere. The high similarities in the amplitude modulation, long-term amplitude variation, period range between the oscillation of investigated parameters and the global activity of oscillation suggests a possible solar influence on the 18-day quasi-periodic oscillation in the ionosphere.
Shear flow instabilities in the Earth's magnetotail
Authors: Reddy R. V. , Lakhina G. S.
Published: July 2005

Abstract. Shear flow instability is studied in the Earth's magnetotail by treating plasma as compressible. A dispersion relation is derived from the linearized MHD equations using the oscillating boundary conditions at the inner central plasma sheet/outer central plasma sheet (OCPS) interface and OCPS/plasma-sheet boundary layer (PSBL) interface, whereas the surface-mode boundary condition is used at the PSBL/lobe interface. The growth rates and the real frequencies are obtained numerically for near-Earth (∣X∣~10–15 RE) and far-Earth (∣X∣~100 RE) magnetotail parameters. The periods and wavelengths of excited modes depend sensitively on the value of plasma-sheet half thickness, L, which is taken as L=5 RE for quiet time and L=1 RE for disturbed time. The plasma-sheet region is found to be stable for constant plasma flows unless MA3>1.25, where MA3 is the Alfvén Mach number in PSBL. For near-Earth magnetotail, the excited oscillations have periods of 2–20 min (quiet time) and 0.5–4 min (disturbed time) with typical transverse wavelengths of 2–30 RE and 0.5–6.5 RE, respectively; whereas for distant magnetotail, the analysis predicts the oscillation periods of ~8–80 min for quiet periods and 2–16 min for disturbed periods.