Pulsating of the generalized ion and neutral polar winds

A three-dimensional, time-dependent model of the ion and neutral polar winds was used to study their dynamic evolution during the May 4, 1998 magnetic storm. The simulation tracked the dynamics of five species (O⁺, H⁺, H_s, O_s, and electrons) and covered a 9-h period. During the storm, D_st decreased to −210 nT, Ap reached 300, and K_p was elevated. The IMF B_z component was southward at the start of the storm and for several hours thereafter and then turned northward. However, the magnetospheric energy input to the ionosphere exhibited a 50-min oscillation, with the plasma convection and particle precipitation patterns expanding and contracting in a periodic manner. As a consequence, the ion and neutral polar winds pulsated with an approximate 50-min period. The H⁺ and O⁺ ions displayed cyclic upflows and downflows in the topside ionosphere as well as a highly structured spatial distribution that varied with time. The vertical flux of the neutral H_s atoms was upward at the top of the ionosphere, but the magnitude varied in a cyclic manner in response to the oscillating stormtime energy input. The vertical flux of neutral O_s atoms was downward at the top of the ionosphere and varied significantly with the stormtime energy input. For H⁺, O⁺, and H_s, the maximum total (integrated) vertical flux during the storm was upward at the top of the ionosphere, with values of 8–9×10²⁵ particles/s for H⁺, 2–4×10²⁶ particles/s for O⁺, and 2–3×10²⁷ particles/s for H_s. The corresponding total vertical O_s flux was predominately downward, with only localized areas with positive fluxes.     

Quiet-condition hmF2, NmF2, and B0 variations at Jicamarca and comparison with IRI-2001during solar maximum

We use the measurements of the Jicamarca digisonde to examine the variations in F2 layer peak electron density (NmF2), its height (hmF2), and the F2 layer thickness parameter (B₀) near the dip equator. The hourly ionograms during geomagnetic quiet-conditions for a 12-month period close to the maximum solar activity, April 1999–March 2000, are used to calculate the monthly averages of these parameters, for each month. The averages are compared with the International Reference Ionosphere (IRI)-2001 model values. The results show that the higher hmF2 values during daytime, associated with the upward velocity, are mainly responsible for the greater values of NmF2 and B₀; while the nighttime lower hmF2, related to the downward velocity, are responsible for the smaller NmF2 and B₀. For daytime, hmF2 and NmF2 are correlated with the solar activity in the equinoctial and summer months. The hmF2 and B₀ peaks at sunset with an associated sharp decrease in NmF2 are presented in the equinoctial and summer months, but not in the winter months. Comparison of the measured hmF2 values with the International Radio Consultative Committee (CCIR) maps used in IRI-2001 (IRI-CCIR) reveals an IRI overestimate in hmF2 during daytime. The most significant discrepancy is that the IRI-CCIR does not model the post-sunset peak in hmF2. For the NmF2 comparison, the values obtained from both the CCIR and URSI maps are generally close to the observed values. For the B₀ comparison, the highest discrepancy between the observation and the Gulyaeva option (IRI-Gulyaeva) is the location of the annual maximum for the daytime values, also the winter daytime predictions are too low. Additionally, the significant negative difference between the observation and the B₀-table option (IRI-B₀-table) provides a slightly better prediction, except for 0400–1000 LT when the model significantly overestimates. The post-sunset peak in B₀ at some months is predicted by neither the IRI-Gulyaeva nor the IRI-B₀-table options.     

Time variations of geomagnetic activity indices Kp and Ap: an update

Kp and Ap indices covering the period 1932 to 1995 are analysed in a fashion similar to that attempted by Bartels for the 1932–1961 epoch to examine the time variations in their characteristics. Modern analysis techniques on the extended data base are used for further insight. The relative frequencies of occurrence of Kp with different magnitudes and the seasonal and solar cycle dependences are seen to be remarkably consistent despite the addition of 35 years of observations. Many of the earlier features seen in the indices and special intervals are shown to be replicated in the present analysis. Time variations in the occurrence of prolonged periods of geomagnetic calm or of enhanced activity are presented and their relation to solar activity highlighted. It is shown that in the declining phase the occurrence frequencies of Kp = 4–5 (consecutively over 4 intervals) can be used as a precursor for the maximum sunspot number to be expected in the next cycle. The semiannual variation in geomagnetic activity is reexamined utilising not only the Ap index but also the occurrence frequencies of Kp index with different magnitudes. Lack of dependence of the amplitude of semiannual variation on sunspot number is emphasised. Singular spectrum analysis of the mean monthly Ap index shows some distinct periodic components. The temporal evolution of 44 month, 21 month and 16 month oscillations are examined and it is postulated that while QBO and the 16 month oscillations could be attributed to solar wind and IMF oscillations with analogous periodicity, the 44 month variation is associated with a similar periodicity in recurrent high speed stream caused by sector boundary passage. It is reconfirmed that there could have been only one epoch around 1940 when solar wind speed could have exhibited a 1.3-year periodicity comparable to that seen during the post-1986 period.     

The sunspot cycle, the QBO, and the total ozone over Northeastern Europe: a connection through the dynamics of stratospheric circulation

The interaction between the factors of the quasi-biennial oscillation (QBO) and the 11-year solar cycle is considered as an separate factor influencing the interannual January–March variations of total ozone over Northeastern Europe. Linear correlation analysis and the running correlation method are used to examine possible connections between ozone and solar activity at simultaneous moment the QBO phase. Statistically significant correlations between the variations of total ozone in February and, partially, in March, and the sunspot numbers during the different phases of QBO are found. The running correlation method between the ozone and the equatorial zonal wind demonstrates a clear modulation of 11-y solar signal for February and March. Modulation is clearer if the QBO phases are defined at the level of 50 hPa rather than at 30 hPa. The same statistical analyses are conducted also for possible connections between the index of stratospheric circulation C₁ and sunspot numbers considering the QBO phase. Statistically significant connections are found for February. The running correlations between the index C₁ and the equatorial zonal wind show the clear modulation of 11-y solar signal for February and March. Based on the obtained correlations between the interannual variations of ozone and index C₁, it may be concluded that a connection between solar cycle – QBO – ozone occurs through the dynamics of stratospheric circulation.     

Observations of the response time of high-latitude ionospheric convection to variations in the interplanetary magnetic field using EISCAT and IMP-8 data

We have combined ∼300 h of tristatic measurements of the field-perpendicular F region ionospheric flow measured overhead at Tromsø by the EISCAT UHF radar, with simultaneous IMP-8 measurements of the solar wind and interplanetary magnetic field (IMF) upstream of the Earth’s magnetosphere, in order to examine the response time of the ionospheric flow to changes in the north-south component of the IMF (B_z). In calculating the flow response delay, the time taken by field changes observed by the spacecraft to first effect the ionosphere has been carefully estimated and subtracted from the response time. Two analysis methods have been employed. In the first, the flow data were divided into 2 h-intervals of magnetic local time (MLT) and cross-correlated with the “half-wave rectifier” function V²B_s, where V is the solar wind speed, and B_s is equal to IMF B_z if the latter is negative, and is zero otherwise. Response delays, determined from the time lag of the peak value of the cross-correlation coefficient, were computed versus MLT for both the east-west and north-south components of flow. The combined data set suggests minimum delays at ∼1400 MLT, with increased response times on the nightside. For the 12-h sector centred on 1400 MLT, the weighted average response delay was found to be 1.3 ± 0.8 min, while for the 12-h sector centred on 0200 MLT the weighted average delay was found to increase to 8.8 ± 1.7 min. In the second method we first inspected the IMF data for sharp and enduring (at least ∼5 min) changes in polarity of the north-south component, and then examined concurrent EISCAT flow data to determine the onset time of the corresponding enhancement or decay of the flow. For the case in which the flow response was timed from whichever of the flow components responded first, minimum response delays were again found at ∼1400 MLT, with average delays of 4.8 ± 0.5 min for the 12-h sector centred on 1400 MLT, increasing to 9.2 ± 0.8 min on the nightside. The response delay is thus found to be reasonably small at all local times, but typically ∼6 min longer on the nightside compared with the dayside. In order to make an estimate of the ionospheric information propagation speed implied by these results, we have fitted a simple theoretical curve to the delay data which assumes that information concerning the excitation and decay of flow propagates with constant speed away from some point on the equatorward edge of the dayside open-closed field line boundary, taken to lie at 77° magnetic latitude. For the combined cross-correlation results the best-fit epicentre of information propagation was found to be at 1400 MLT, with an information propagation phase speed of 9.0 km s⁻¹. For the combined event analysis, the best-fit epicentre was also found to be located at 1400 MLT, with a phase speed of 6.8 km s⁻¹.     

Scattering attenuation and intrinsic absorption using uniform and depth dependent model – Application to full seismogram envelope recorded in Northern Chile

Two seismic wave attenuation factors, scatteringattenuation Q _s ^-1 and intrinsicabsorption Q _i ^-1 are measured using theMultiple Lapse Time Window (MLTW) analysis method forthree different frequency bands, 1–2, 2–4, and 4–8 Hz.Data from 54 temporally deployed seismic stationslocated in northern Chile are used. This methodcompares time integrated seismic wave energies withsynthetic coda wave envelopes for a multiple isotropicscattering model. In the present analysis, the waveenergy is assumed to decay with distance in proportionto1/GSF·exp(- (Q _s ^-1+Q _i ^-1)· r/v), where r, and v are the propagationdistance, angular frequency and S wave velocity,respectively, and GSF is the geometricalspreading factor. When spatial uniformity of Q _s ^-1, Q _i ^-1 and v isassumed, i.e. GSF = 4r ², theestimates of the reciprocal of the extinction length,L _e ^-1 (= (Q _s ^-1+Q _i ^-1)·/v), are 0.017,0.012 and 0.010 km^-1, and those of the seismicalbedo, B ₀ (= Q _s ^-1/ (Q _s ^-1+Q _i ^-1)), are 0.48, 0.40and 0.34 for 1–2, 2–4 and 4–8 Hz, respectively, whichindicates that scattering attenuation is comparable toor smaller than intrinsic absorption. When we assumea depth dependent velocity structure, we also findthat scattering attenuation is comparable to orsmaller than intrinsic absorption. However, since thequantitative estimates of scattering attenuationdepend on the assumed velocity structure (strength ofvelocity discontinuity and/or Moho depth), it isimportant to consider differences in velocitystructure models when comparing attenuation estimates.     

Energy mode distribution: An analysis of the ratio of anti-Stokes to Stokes amplitudes generated by a pair of counterpropagating Langmuir waves

Results from plasma wave experiments in spacecrafts give support to nonlinear interactions involving Langmuir, electromagnetic, and ion-acoustic waves in association with type III solar radio bursts. Starting from a general form of Zakharov equation (Zakharov, V.E., 1985. Collapse and self-focusing of Langmuir waves. Hand-book of Plasma Physics Cap.2, 81–121) the equations for electric fields and density fluctuations (density gratings) induced by a pair of counterpropagating Langmuir waves are obtained. We consider the coupling of four triplets. Each two triplets have in common the Langmuir pump wave (forward or backward wave) and a pair of independent density gratings. We numerically solve the dispersion relation for the system, extending the work of (Alves, M.V., Chian, A.C.L., Moraes, M.A.E., Abalde, J.R., Rizzato, F.B., 2002. A theory of the fundamental plasma emission of type- III solar radio bursts. Astronomy and Astrophysics 390, 351–357). The ratio of anti-Stokes (AS) (ω₀+ω) to Stokes (S) (ω₀-ω^*) electromagnetic mode amplitudes is obtained as a function of the pump wave frequency, wave number, and energy. We notice that the simultaneous excitation of AS and S distinguishable modes, i.e., with Re{ω}=ω_r≠0, only occurs when the ratio between the pump wave amplitudes, r is ≠1 and the pump wave vector k₀ is , W₀ being the forward pump wave energy. We also observe that the S mode always receives more energy.     

Total Solar Irradiance Observations

The record of total solar irradiance (TSI) during the past 35?years has overlapping observations from space which can be merged to a composite, and three are available, namely the PMOD, the ACRIM and the IRMB composites. There are important differences between them, which are discussed in detail in order to find the best representation of solar variability during the last three cycles, for the following discussions of solar irradiance variability. Moreover, the absolute value of TSI from TIM on SORCE is 1,361?Wm^?2, substantially lower than the value 1,365?Wm^?2, which was observed by the classical radiometers. New results from specific experiments are now available, which are discussed in order to define the value to be used in, e.g., climate models. The most important issue regarding the recent TSI records is the low value observed during the minimum in 2009, which is 25% of a typical cycle amplitude lower than the value in 1996. The validity of this low value has been confirmed by comparing all existing TSI observations during cycle 23. On the other hand, activity indices, such as the sunspot number, the 10.7-cm radio flux (F10.7), the CaII and MgII indices and also the Ly-α irradiance or the frequency changes in low-order p modes, show a much smaller decreases relative to their respective typical cycle amplitude. It is most likely that an increasing contrast of the facular and network elements with decreasing magnetic field is responsible for this discrepancy. The value of TSI at minima is correlated with the open magnetic field of the Sun, B _R,?at minima. Using B _R at minima, interpolated linearly in between as a fourth component of a proxy model based on the photometric sunspot index and on the MgII index improves the explanation of the variance of TSI over the full period of the last three solar cycles to 84.7%. Results from other models are also discussed.     

Erosion of the inner magnetosphere during geomagnetic storms

Using the empirical magnetic field model dependent on the Dst index and solar wind dynamic pressure, we calculated the behaviour of the contour B = B_s in the equatorial plane of the magnetosphere where B_s is the magnetic field in the subsolar point at the magnetopause. The inner domain of the magnetosphere outlined by this contour contains the bulk of geomag-netically trapped particles. During quiet time the boundary of the inner magnetosphere passes at the distance ∼10R_E at noon and at ∼7R_E at midnight. During very intense storms this distance can be reduced to 4–5 R_E for all MLT. The calculation results agree well with the satellite measurements of the magneto-pause location during storms. The ionospheric projection of the B = B_s contour calculated with the Euler potential technique is close to the equatorward edge of the auroral oval.     

Solar wind velocity distribution on the heliospheric current sheet during Carrington rotations 1787-1795

The solar wind velocity distribution in the heliosphere is best represented using a v-map, where velocity contours are plotted in heliographic latitude-longitude coordinates. It has already been established that low-speed regions of the solar wind on the source surface correspond to the maximum bright regions of the K-corona and the neutral line of the coronal magnetic field. In this analysis, v-maps on the source surface for Carrington rotations (CRs) 1787-1795, during 1987, have been prepared using the interplanetary scintillation measurements at Research Institute of Atmospherics (RIA), Nagoya Univ., Japan. These v-maps were then used to study the time evolution of the low-speed (\leq450 km s⁻¹) belt of the solar wind and to deduce the distribution of solar wind velocity on the heliospheric current sheet. The low-speed belt of the solar wind on the source surface was found to change from one CR to the next, implying a time evolution. Instead of a slow and systematic evolution, the pattern of distribution of solar wind changed dramatically at one particular solar rotation (CR 1792) and the distributions for the succeeding rotations were similar to this pattern. The low-speed region, in most cases, was found to be close to the solar equator and almost parallel to it. However, during some solar rotations, they were found to be organised in certain longitudes, leaving regions with longitudinal width greater than 30^○ free of low-speed solar wind, i.e. these regions were occupied by solar wind with velocities greater than 450 km s⁻¹. It is also noted from this study that the low-speed belt, in general, followed the neutral line of the coronal magnetic field, except in certain cases. The solar wind velocity on the heliospheric current sheet (HCS) varied in the range 300–585 km s⁻¹ during the period of study, and the pattern of velocity distribution varied from rotation to rotation.     

A straightforward estimation of the maximum sunspot number for cycle 23

Using the annual number of geomagnetically quiet days (aa < 20 γ) for the year after the solar minimum, this precursor method predicts that the maximum sunspot number for cycle 23 will be 140 + 32, indicating that cycle 23 will be similar to cycles 21 and 22.     

Vertical velocities associated with gravity waves measured in the mesosphere and lower thermosphere with the EISCAT VHF radar

The EISCAT VHF radar (69.4°N, 19.1°E) has been used to record vertical winds at mesopause heights on a total of 31 days between June 1990 and January 1993. The data reveal a motion field dominated by quasi-monochromatic gravity waves with representative apparent periods of 30–40 min, amplitudes of up to 2.5 m s^–1 and large vertical wavelength. In some instances waves appear to be ducted. Vertical profiles of the vertical-velocity variance display a variety of forms, with little indication of systematic wave growth with height. Daily mean variance profiles evaluated for consecutive days of recording show that the general shape of the variance profiles persists over several days. The mean variance evaluated over a 10 km height range has values from 1.2 m²s^–2 to 6.5 m²s^–2 and suggests a semi-annual seasonal cycle with equinoctial minima and solsticial maxima. The mean vertical wavenumber spectrum evaluated at heights up to 86 km has a slope (spectral index) of -1.36 ± 0.2, consistent with observations at lower heights but disagreeing with the predictions of a number of saturation theories advanced to explain gravity-wave spectra. The spectral slopes evaluated for individual days have a range of values, and steeper slopes are observed in summer than in winter. The spectra also appear to be generally steeper on days with lower mean vertical-velocity variance.     

Quasi-biennial and quasi-triennial oscillations in geomagnetic activity indices

Data for geomagnetic activity index aa for 1868–1994 were subjected to spectral analysis for 12 intervals each of 11 consecutive years. In each interval, QBO and QTO (quasi-biennial and quasi-triennial oscillations) were observed at ∼ 2.00, 2.15, 2.40, 2.70 y and ∼ 3.20, 3.40 y, but not all in all intervals. These fluctuations are absent near (2–3 y before and after) the sunspot minima and are present only as 2 or 3 peaks in aa indices, one near or before the sunspot maximum and the other (one or two, generally the larger ones) in the declining phase of the sunspot cycle. Comparison with the solar wind (1965 onwards) showed a fairly good match, indicating that the aa variations were mostly due to similar variations in the solar wind, which must have their origin in solar physical processes. A few aa variations did not match with solar wind. When compared with terrestrial phenomena, no match was found with stratospheric low-latitude zonal wind QBO; but some QTO in aa matched QTO in ENSO (El Nino/ Southern Oscillation). This may or may not be a chance coincidence and needs further exploration.     

Characteristics of VHF radiowave scintillations over a solar cycle (1983–1993) at a low-latitude station: Waltair (17.7°N, 83.3°E)

The characteristics of VHF radiowave scintillations at 244 MHz (FLEETSAT) during a complete solar cycle (1983–93) at a low-latitude station, Waltair (17.7°N, 83.3°E), are presented. The occurrence of night-time scintillations shows equinoctial maxima and summer minima in all the epochs of solar activity, and follows the solar activity. The daytime scintillation occurrence is negatively correlated with the solar activity and shows maximum occurrence during the summer months in a period of low solar activity. The occurrence of night-time scintillations is inhibited during disturbed days of high solar activity and enhanced during low solar activity.     

Ecosystem water use efficiency in an irrigated cropland in the North China Plain

The eddy covariance technique and the cuvette method were used to investigate water use efficiency in an irrigated winter wheat (Triticum asetivum L.)/summer maize (Zea mays L.) rotation system in the North China Plain. The results show that ecosystem water use efficiency (WUE_e) changed diurnally and seasonally. Daily maximal WUE_e appeared in the morning. WUE_e generally peaked in late April in wheat field and in late July/early August in maize field. From 2003 to 2006, seasonal mean WUE_e was 6.7–7.4 mg CO₂ g⁻¹ H₂O for wheat and 8.4–12.1 mg CO₂ g⁻¹ H₂O for maize. WUE_e was much lower than canopy water use efficiency (WUE_c) under small leaf area index (LAI) but very close to WUE_c under large LAI. With the increase in LAI, WUE_e enlarged rapidly under low LAI but slowly when LAI was higher than one. WUE_e was greater on the cloudy days than on the sunny days. Under the same solar radiation, WUE_e was higher in the morning than in the afternoon. The ratio of internal to ambient CO₂ partial pressure (C_i/C_a) decreased significantly with the increase in photosynthetically active radiation (PAR) when PAR was lower than the critical values (around 500 and 1000 μmol m⁻² s⁻¹ for wheat and maize, respectively). Beyond critical PAR, C_i/C_a was approximately constant at 0.69 for wheat and 0.42 for maize. Therefore, when LAI and solar radiation was large enough, WUE_e has negative correlation with vapor pressure deficit in both of irrigated wheat and maize fields.     

A self-consistent estimate of O+ + N2– rate coefficient and total EUV solar flux with λ < 1050 Å using EISCAT observations

There are differences between existing models of solar EUV with < 1050 Å and between laboratory measurements of the O⁺ + N₂ – reaction rate coefficient, both parameters being crucial for the F2-region modeling. Therefore, indirect aeronomic estimates of these parameters may be useful for qualifying the existing EUV models and the laboratory measured O⁺ + N₂ – rate coefficient. A modified self-consistent method for daytime F2-region modeling developed by Mikhailov and Schlegel was applied to EISCAT observations (32 quiet summer and equinoctial days) to estimate the set of main aeronomic parameters. Three laboratory measured temperature dependencies for the O⁺ + N₂ – rate coefficient were used in our calculations to find self-consistent factors both for this rate coefficient and for the solar EUV flux model from Nusinov. Independent of the rate coefficient used, the calculated values group around the temperature dependence recently measured by Hierl et al. in the 850–1400 K temperature range. Therefore, this rate coefficient may be considered as the most preferable and is recommended for aeronomic calculations. The calculated EUV flux shows a somewhat steeper dependence on solar activity than both, the Nusinov and the EUVAC models predict. In practice both EUV models may be recommended for the F2-region electron density calculations with the total EUV flux shifted by ±25% for the EUVAC and Nusinov models, correspondingly.     

Short-period rayleigh waves recorded in the western carpathians during dss

Summary The periods, the maximum absolute displacement amplitudes and the maximum particle velocities of the surface waves, propagating in the weathered layer are investigated. Dependences of the parameters under discussion on the distance r (km) between sites and shot points are expressed for distances from 3.6 to 38.6 km by the functions: T(s)=0.40r ^0.30 , A _max (m)==502.73r ^–1.93 and v _max (mm s –1)=7.95 r –2.22.     

A report on long-term trends and variabilities in middle atmospheric temperature over Gadanki (13.5°N, 79.2°E)

The present study reports long-term variabilities and trends in the middle atmospheric temperature (March 1998–2008) derived from Rayleigh backscattered signals received by the Nd:YAG lidar system at Gadanki (13.5°N, 79.2°E). The monthly mean temperature compositely averaged for the years 1998–2008 shows maximum temperature of 270 K in the months of March–April and September at altitudes between 45 and 55 km. The altitude profile of trend coefficients estimated from the 10 years of temperature observations using regression analysis shows that there exists cooling at the rate with 1σ uncertainty of 0.12±0.1 K/year in the lower stratospheric altitudes (35–42 km) and 0.2±0.08 K/year at altitudes near 55–60 km. The trend is nearly zero (no significant cooling or warming) at altitudes 40–55 km. The regression analysis reveals the significant ENSO response in the lower stratosphere (1 K/SOI) and also in mesosphere (0.6 K/SOI). The solar cycle response shows negative maxima of 1.5 K/100F10.7 units at altitudes 36 km, 41 km and 1 K/100F10.7 units at 57 km. The response is positive at mesospheric altitude near 67 km (1.3 K/100F10.7 units). The amplitudes and phases of semiannual, annual and quasi-biennial oscillations are estimated using least squares method. The semiannual oscillation shows larger amplitudes at altitudes near 35, 45, 62 and 74 km whereas the annual oscillation peaks at 70 km. The quasi-biennial oscillations show larger amplitudes below 35 km and above 70 km. The phase profiles of semiannual and annual oscillations show downward propagation.     

Effect of the IMF B y component on the ionospheric flow overhead at EISCAT: observations and theory

We have analysed a database of 300 h of tristatic ionospheric velocity measurements obtained overhead at Tromsø (66.3° magnetic latitude) by the EISCAT UHF radar system, for the presence of flow effects associated with the y-component of the IMF. Since it is already known that the flow depends upon IMF B_z, a least-squares multivariate analysis has been used to determine the flow dependence on both IMF B_y and B_z simultaneously. It is found that significant flow variations with IMF B_y occur, predominantly in the midnight sector (2100/0300 MLT), but also pre-dusk (1600/1700 MLT), which are directed eastward for IMF B_y positive and westward for IMF B_y negative. The flows are of magnitude 20/30 m s^–1 nT^–1 in the midnight sector, and smaller, 10/20 m s^–1 nT^–1, pre-dusk, and are thus associated with significant changes of flow of order a few hundred m s^–1 over the usual range of IMF B_y of about ±5 nT. At other local times the IMF B_y-related perturbation flows are much smaller, less than 5 m s^–1 nT^–1, and consistent with zero within the uncertainty estimates. We have investigated whether these IMF B_y-dependent flows can be accounted for quantitatively by a theoretical model in which the equatorial flow in the inner magnetosphere is independent of IMF B_y, but where distortions of the magnetospheric magnetic field associated with a penetrating component of the IMF B_y field changes the mapping of the field to the ionosphere, and hence the ionospheric flow. We find that the principal flow perturbation produced by this effect is an east-west flow whose sense is determined by the north-south component of the unperturbed flow. Perturbations in the north-south flow are typically smaller by more than an order of magnitude, and generally negligible in terms of observations. Using equatorial flows which are determined from EISCAT data for zero IMF B_y, to which the corotation flow has been added, the theory predicts the presence of zonal perturbation flows which are generally directed eastward in the Northern Hemisphere for IMF B_y positive and westward for IMF B_y negative at all local times. However, although the day and night effects are therefore similar in principle, the model perturbation flows are much larger on the nightside than on the dayside, as observed, due to the day-night asymmetry in the unperturbed magnetospheric magnetic field. Overall, the model results are found to account well for the observed IMF B_y-related flow perturbations in the midnight sector, in terms of the sense and direction of the flow, the local time of their occurrence, as well as the magnitude of the flows (provided the magnetic model employed is not too distorted from dipolar form). At other local times the model predicts much smaller IMF B_y-related flow perturbations, and thus does not account for the effects observed in the pre-dusk sector.     

A fracture mechanics study of subcritical tensile cracking of quartz in wet environments

Load relaxation and cross-head displacement rate-change experiments have been used to establish log₁₀ stress intensity factor (K) versus log₁₀ crack velocity (v) diagrams for double torsion specimens, of synthetic quartz cracked on thea plane in liquid water and moist air.For crack propagation normal toz and normal tor at 20°C,K _Ic (the critical stress intensity factor) was found to be 0.852±0.045 MN·m^–3/2 and 1.002±0.048 MN·m^–3/2, respectively.Subcritical crack growth at velocities from 10^–3 m·s^–1 to 10^–9 m·s^–1 at temperatures from 20°C to 80°C is believed to be facilitated by chemical reaction between the siloxane bonds of the quartz and the water or water vapour of the environment (stress corrosion). The slopes, of isotherms in theK-v diagrams are dependent upon crystallographic orientation. The isotherms have a slope of 12±0.6 for cracking normal tor and 19.9±1.7 for cracking normal toz. The activation enthalpy for crack propagation in the former orientation in liquid water at temperatures from 20°C to 80°C is 52.5±3.8 kJ·mole^–1.A discussion is presented of the characteristics of theK-v diagrams for quartz.