Long-term relationships between sunspots,Ca-plages and the ionosphere

Smith and King (1981) recently found a relationship for Slough values of foE and foE2 using relative sunspot number (R) and area of photospheric faculae (A_F) from Greenwich publications. Since Greenwich stopped publishing the area of photospheric faculae, it is suggested in this paper that the area of Ca-plages (A_P) may be substituted for A_F for calculating foE and foF2. The values of foF2 calculated using the relationship given by Smith and King, with A_P and A_F separately, are found to have a good correlation (r=0.988). Similar equations showing the relationship between the foF2 values, sunspot numbers and Ca-plage area have been obtained for a few Indian stations. It has been found that in the regression equation, while the coefficients of the sunspot number and the Ca-plage area do not change significantly from one station to another, the constant term shows a latitude variation similar to the geomagnetic anomaly.     

A comparison of foF2 obtained from a time-dependent ionospheric model with Argentine Islands' data for quiet conditions

In this study a comparison is made of the Utah State University Time-Dependent Ionospheric Model (TDIM) and an ionosonde data set from Argentine Islands. This study is unique in that the Argentine Islands data set of foF2 spans complete diurnal, seasonal and solar cycle conditions for low geomagnetic activity. The TDIM reproduces these foF2 variations extremely well. Although the observed winter and summer solstice foF2 diurnal curves have opposite phases, they are readily modelled. At equinox where a sharp transition occurs from winter to summer, or vice versa, the monthly average is complicated by this feature and hence the TDIM does not reproduce the diurnal fine structure.The neutral wind induced vertical plasma drift is the only free parameter in this study. All the other inputs are fixed for the specific solar, seasonal and diurnal conditions. A vertical plasma drift variation is presented; although simplistic, it couples the geographic and geomagnetic frames. With additional information such as hmF2, it would be possible to deduce a unique vertically induced drift pattern.     

Comparison of foF2 characteristics at similar low latitudes but different longitudes in the South American region

A comparison of the foF2 diurnal variation characteristics at Huancayo (12°S, 75°W, dip angle +2°)and Fortaleza(4°S, 38°W, dip angle −2° in 1976) showed in general larger midday foF2 at Fortaleza. In contrast, a comparison of foF2 at Tucuman (27°S, 65°W, dip angle − 22°) and Sào Paulo (24°S, 47°W, dip angle −22° in IGY period) showed smaller midday foF2 at São Paulo. Thus, a weaker fountain effect in the eastern part of the South American region is indicated.     

foF2 hysteresis variations and the semi-annual geomagnetic wave

Seasonal and solar cycle variations of the foF2 hysteresis magnitude are investigated. Data for the noon foF2 monthly medians for Slough (51.48°N, 0.57°W), the monthly means for the sunspot numbers, and for the geomagnetic activity index aa(N) for the northern hemisphere for the period 1933–1986, covering solar cycle from 17 to 21, are used. It is found that: (1) the greatest negative amplitudes of the foF2 hysteresis variation are near the equinoxes, and (2) the solar cycle average noon foF2 hysteresis magnitude is linearly correlated with the solar cycle average semi-annual geomagnetic amplitude of the aa-index. These results support the hypothesis that the foF2 hysteresis is due to the geomagnetic activity variation during the sunspot cycle.     

Observed and calculated F2 peak heights and derived meridional winds at mid-latitudes over a full solar cycle

The peak height of the F2 layer, hmF2, has been calculated using the ‘servo’ model of Rishbeth et al. [(1978), J. atmos. terr. Phys. 40, 767], combined with the hedin et al. [(1988), J. geophys. Res. 93, 9959] neutral wind model. The results are compared with observed values at noon and midnight derived from ionosonde measurements at two mid-latitude stations, Boulder and Wallops Island, over a full solar cycle. The reduced height of the F2 layer, zmF2, is also computed for the same period using the observed hmF2 values and the MSIS-86 model. Day-night, seasonal, and solar cycle variations in zmF2 are attributed to neutral composition changes and winds. Anomalously low values of hmF2 and zmF2 during summer both at solar minimum and during the solar cycle maximum in magnetic activity may be associated with increases in the molecular to atomic ion concentration ratio. Under these circumstances the F2 peak may lie significantly below the O⁺ peak height calculated by the servo model. Neutral meridional winds at Wallops Island are derived from the servo model using the observed hmF2 values and the calculated O⁺ ‘balance height’. It is shown that if the anomalously low hmF2 values are used, unrealistically large poleward winds are derived, which are inconsistent with both theory and observations made using other techniques. For most conditions the F2 peak is clearly an O⁺ peak, and daily mean winds at hmF2 derived from the servo model are consistent with the hedin et al. (1988) wind model. Unexpectedly, the results do not show an abrupt transition in the thermospheric circulation at the equinoxes. Diurnal curves of the servo model winds reveal a larger day-night difference at solar minimum than at solar maximum.     

The influence of neutral temperatures and winds on the F-Iayer height

Observations of neutral winds and temperatures obtained using a FabryPerot interferometer at Beveridge (37°28′S, 145°6′E) have been combined with h'F measurements from ionosondes at Canberra (35°21′S, 149°10′E) and Hobart (42°54′S, 147° 12′E). Data from 16 nights have been used to study the response (height change) of the F2-layer to changes in neutral wind and temperature. The observations have been compared with the ‘servo’ model of Rishbeth. It is found that the ‘night stationary level’ of the F2-layer depends on temperature, with the height changing by (13 ± 6) km per 100K. This agrees well with the prediction of the ‘servo’ model. There is reasonable overall agreement between the observations and the model predictions for the change in height produced by a given meridional wind. However, there is considerable scatter in the individual comparisons due to the approximations used to apply the theory to the observations. In particular, the effect of electric fields on the F2-layer height has been ignored.     

Comparison of total electron content calculated using the IRI with observations in China

Values of total electron content (TEC) calculated using the International Reference Ionosphere (IRI-86 and IRI-90) are compared with the observations at Xinxiang based on the Faraday rotation measurement. It is found that the IRI gives acceptable values with respect to the observations during low solar activity. Generally the IRI-90 is better than the IRI-86 and the URSI coefficients are better than the CCIR coefficients in the calculation of TEC. Making use of the foF2 and M(3000)F₂ calculated using the Asia Oceania Region F₂-layer mapping (AOR) instead of using the CCIR or the URSI coefficients, the IRI gives more accurate TEC values. In October-April during high solar activity, however, the IRI obviously underestimates TEC in the daytime, which could be due to an improper topside electron density profile.     

Seasonal variations of day-time ionisation flows inferred from a comparison of calculated and observed NmF2

This paper reports on a comparison of calculated and observed monthly mean day-time ionospheric F2-peak density (NmF2) at a chain of stations from Japan to Australia for both solar minimum (1976) and solar maximum (1980). Nm values are calculated using the MSIS model for the observed peak heights (hmF2) and a simplified version of the continuity equation for day-time equilibrium conditions. The observed NmF2 values are always higher than the calculated ones in winter. This implies that a substantial downward flow of ionisation from above into the winter ionosphere is induced by the strongly poleward winter neutral wind which drives the ionisation down the field lines, lowering the peak height hmF2. In summer, winds are smaller, and the fluxes are more upward in comparison to winter. The seasonal variation of the ionisation fluxes and neutral winds are estimated for solar minimum, and compared with results of detailed calculations.     

Mesopause temperatures and integrated band brightnesses calculated from airglow OH emissions recorded at Maynooth (53.2°N, 6.4°W) during 1993

Spectra of the hydroxyl emissions in the wavelength range 1.0–1.6 μm, which originate at mesopause altitudes, have been obtained, using a Fourier transform spectrometer at Maynooth (53.2°N, 6.4°W), on all suitable nights during the period January–December 1993. Rotational temperatures and integrated band brightnesses have been calculated from the spectra of the OH(3, 1) and (4, 2) vibration-rotation bands. The mean annual temperatures calculated over all measurements were T(3, 1)=200±19 K and T(4, 2)=206±19 K, where the uncertainty represents the standard deviation on the measurements. Harmonic analysis of the nightly averaged temperature values revealed an amplitude of 27 ± 1 K and a phase of 95 ± 2 days in the annual variation of the (3, l) band at our latitude. The semiannual component was found to have an amplitude of 7 ± 1 K and a phase of −51 ± 9 days for this band. Results for the (4, 2) band were identical in both amplitude and phase for the annual component, while the semiannual component gave an amplitude of 8 ± 1 K and a phase of − 43 ± 7 days. These results are compared with data recorded by the SME satellite, and with the predictions of the MSISE-90 model for a station at 53° latitude. Temperatures predicted by the MSISE-90 model for Maynooth are consistently below the values obtained in this study by 15–20 K. Excellent agreement is observed between the absolute value of temperature, in the case of the SME satellite, and in the amplitude and phase of the annual variation predicted by MSISE-90. The phase of the semiannual component observed in our data deviates somewhat from the −99 ± 1 days predicted by MSISE-90.The annual mean brightness of the OH (3, 1) band was found to be 75 ± 18 kR, while that of the (4, 2) band was 106±26 kR. Diurnal variations generally showed a steady decrease from dusk to dawn, apart from a brief period in June and July. Monthly average values of band brightness have been calculated for each band and are compared with the predictions of a recent photochemical model (Le Texier et al., 1987). The model shows some elements of agreement with our observations, particularly a pair of maxima near the equinoxes, but it does not predict the broad winter maximum observed in both bands at this latitude.     

Some results on mean tidal structure and day-to-day variability over Arecibo

A long series of monthly 36 h observations, plus a series of observations at the same local time each day over two months at Arecibo, have been analyzed for tidal structure and variability. Some of the results are as follows. (1) A diurnal tide with vertical structure similar to that of the S_1,1 mode dominates the wind field up to heights of the order of 110 km over Arecibo. A semidiurnal oscillation dominates above that height. For non-winter conditions the semidiurnal oscillation in the 80–200 km region closely resembles the S_2,2 mode, though there is the possibility of contributions from higher order modes in the 100–120 km region. (2) Larger semidiurnal amplitudes are observed in the lower thermosphere for winter conditions. The data appears roughly consistent with Bernard's (1979) hypothesis that S_2,2, S_2,4 and S_2,5 modes are thermally excited, with the S_2,4 and S_2,5 modes out of phase in the meteor region in summer and in phase in winter. (3) The day-to-day variability of the tides is at least half the amplitude of the mean oscillations. The maximum wind variability was observed to occur in the 100–110 km region where the diurnal tide is strongly dissipated. (4) The day-to-day deviations in the wind and temperature oscillations from a long term mean at one local time tend to be wave-like structures which are generally correlated from day to day. The structures tend to move upwards, i.e., appear at a later local time, from day to day.     

The variation in δ13C values in bone collagen for two wild herbivore populations: Implications for palaeodiet studies

The δ¹³C values for un-defatted bone gelatin of individual members of wild populations of mule deer (Odocoileus hemionus hemionus) and field vole (Microtus townsendii) have been measured to determine the variability within each population. Means and standard deviations for the mule deer sample (N=55) and vole sample (N=22) are −20·6±0·5‰ and −23·2±0·6‰, respectively. The distributions of δ¹³C values for each of these populations are normal and we observed no correlation in δ¹³C values with age or sex.     

Equatorial penetration of magnetic disturbance effects in the thermosphere and ionosphere

The neutral dynamic and electrodynamic coupling between high and low latitudes, and the mutual interactions between these two processes, are investigated. For 22 March 1979, when a sudden increase in magnetic activity occurred, we have analyzed the following experimental data: (a) neutral densities and cross-track neutral winds as a function of latitude (0°–80°) near 200 km from a satellite-borne accelerometer; (b) hourly mean H-component magnetic data from the Huancayo Observatory (0.72°S, 4.78°E; dipole geomagnetic coordinates) magnetometer; and (c) hourly mean foF2 measurements from the ionosonde at Huancayo. Comparisons are also made with a self-consistent thermosphere-ionosphere general circulation model and with observationally-based empirical models of winds and density.In concert with the increase in magnetic activity to K_p levels of 5–7, a nighttime (2230 LT) westward intensification of the neutral wind approaching 400 ± 100 ms⁻¹ occurred near the magnetic equator on 22 March 1979, accompanied by a 35% increase in neutral mass density. About 2 h after each of two substorm commencements associated with periods of southward IMF, ∼100γ and ∼200γ reductions in the daytime Huancayo H-component (corrected for ring current effects) are interpreted in terms of ∼0.5 and ∼1.0 mVm⁻¹ westward perturbation electric fields, respectively. An intervening 2-hour period of northward IMF preceded a positive equatorial magnetic perturbation of about 200γ. Time scales for field variations are a few hours, suggesting that processes other than Alfven shielding are involved. Variations in f₀F2 (∼ ± 1.0 MHz) over Huancayo are consistent with the inferred electric fields and magnetic variations. Similar equatorial perturbations are found through examination of other magnetic disturbances during 1979.     

Observations of long delayed echoes

Routine observations of ducted magnetospheric echoes have been made at 1.91 MHz using a ground-based transmitter with a peak power of 17 kW and a pulse length of 150 μs. The positively identified echoes were recorded during periods of moderate solar activity, with delay times of approximately 260 ms, corresponding to L values near three.     

Ionospheric plasma convection in the southern hemisphere

The first ionospheric plasma convection maps ordered by the y- and z-components of the IMF using only data from the southern hemisphere are presented. These patterns are determined from line-of-sight velocity measurements of the Polar Anglo-American Conjugate Experiment (PACE) located at Halley, Antarctica, with the majority of the observations coming from 65°–75° magnetic latitude. For IMF Bz positive and negative conditions, the observed plasma motions are consistent with a standard two cell pattern. For the periods from dusk through midnight to dawn, flow speeds are at least twice as large for Bz negative component compared with Bz positive. The observations about noon are significantly different from each other. For Bz positive, little ordered plasma motion is observed. For Bz negative, there are large anti-sunward flows the orientation of which is ordered by IMF By. These By orientated flows are consistent with theoretical predictions, and are anti-symmetric to those reported from the northern hemisphere. The two most significant differences from previous observations are that the convection reversal in the late morning sector for By negative conditions occurs at about a 4° lower latitude than the Heppner and Maynard (1987) model. This may be due to a seasonal bias in the PACE dataset. Also, the separatrix between eastward and westward flow near midnight has a very different shape dependent upon the orientation of IMF By. For positive By conditions, the separatrix is observed at progressively lower latitudes at later local times, but for By negative conditions, the separatrix appears at increasingly higher latitudes at later times.     

A servo-model of the night-time F2-layer above St. Santin

A modified form of the ionospheric servo-model is used to describe the night-time F2-layer above St. Santin. Data taken by the incoherent scatter radar on nine nights in 1974–1977 were used to determine the height profiles of electron density, electron and ion temperature and electric field. The servo-model was then used to compute the theoretical variation through the night of the height of the F2 peak and the field-aligned plasma velocity, using gas concentrations and horizontal pressure gradients derived from the MSIS79 atmospheric model. On magnetically quiet nights these calculated values agreed closely with the observations. On disturbed nights, however, the calculations and observations began to diverge an hour or so after the onset of a substorm. The divergence could be explained by an enhanced southward wind.     

Windspeed Maps for Tasmania and Their Use in Wind Energy Assessment

Laplacian smoothing spline functions have been used in a detailed spatial analysis of monthly and annual mean windspeed values in Tasmania. The techniques used here have been described in detail by Hutchinson et aL and the present paper significantly extends their results for Tasmania. The windspeed maps are discussed in detail with particular emphasis on their use for the assessment of Tasmania's wind energy potential     

Comparisons between EISCAT observations and model calculations of the high latitude ionosphere

Calculations using a numerical model of the convection dominated high latitude ionosphere are compared with observations made by EISCAT as part of the UK-POLAR Special Programme. The data used were for 24–25 October 1984, which was characterized by an unusually steady IMF, with Bz < 0 and By > 0; in the calculations it was assumed that a steady IMF implies steady convection conditions. Using the electric field models of Heppner and Maynard (1983) appropriate to By > 0 and precipitation data taken from Spiroet al. (1982), we calculated the velocities and electron densities appropriate to the EISCAT observations. Many of the general features of the velocity data were reproduced by the model. In particular, the phasing of the change from eastward to westward flow in the vicinity of the Harang discontinuity, flows near the dayside throat and a region of slow flow at higher latitudes near dusk were well reproduced. In the afternoon sector modelled velocity values were significantly less than those observed. Electron density calculations showed good agreement with EISCAT observations near the F-peak, but compared poorly with observations near 211 km. In both cases, the greatest disagreement occurred in the early part of the observations, where the convection pattern was poorly known and showed some evidence of long term temporal change. Possible causes for the disagreement between observations and calculations are discussed and shown to raise interesting and, as yet, unresolved questions concerning the interpretation of the data. For the data set used, the late afternoon dip in electron density observed near the F-peak and interpreted as the signature of the mid-latitude trough is well reproduced by the calculations. Calculations indicate that it does not arise from long residence times of plasma on the nightside, but is the signature of a gap between two major ionization sources, viz. photoionization and particle precipitation.     

Recent advances in magnetospheric substorm research

More than two decades of magnetospheric exploration have led to a reasonably clear morphological picture of geomagnetic substorms, which is often summarized in terms of the near-Earth neutral line (NENL) model of substorms. Although this qualitative theory is quite comprehensive and explains a great many observations, it is hard pressed to explain both recent observations of consistently earthward flow within 19 R_E and also the prompt onset of magnetic turbulence at 8 R_E at the time of substorm onset. Other theories have recently been proposed which tend to be more quantitative, but which explain a more limited number of substorm observations. The challenge seems to be to understand the essential physics of these various quantitative theories and integrate them into a larger structure such as is provided by the near-Earth neutral line model.     

Wave form and spectral distribution of the electromagnetic field of lightning currents

In a resonant wave guide model of lightning currents, two impulse type standing waves can exist: aperiodic waves of the Bruce-Golde form (type 1), and damped oscillations (type 2). The electromagnetic waves generated by these two types of lightning currents are calculated for various distances and compared with observations. It is shown that the measured wave forms of sferics at distances smaller than about 300 km of return strokes (R strokes) as well as of intracloud strokes (K strokes) are generated mainly by type 2 lightning currents. The channel parameters like channel length and channel diameter derived from the observed sferics are 19 km and 4.6 cm, respectively, for the average R stroke, and 4 km and 1.6 cm for typical K strokes. The large values of the lengths probably correspond to the real lengths of oblique and tortuous channels, rather than to their vertical elevation. The finite electric conductivity of the earth modifies the high frequency component of the wave forms. With decreasing conductivity and/or increasing distance, the rise times of the radiation component to its first maximum increases and the maximum amplitude decreases.Typical rise times for R-strokes are about 3 μs consistent with the observations if the electric conductivity of the earth is of the order of 10⁻³ S/m. The spectral functions of the wave forms are also calculated. The spectral amplitude of the average type 2 R-stroke has its maximum near 4 kHz, and that of the type 2 K-stroke maximizes near 35 kHz. Within the high frequency region at frequencies greater than about 300 kHz. the spectral amplitudes decay proportional to the reciprocal third power of the frequency. The radiation component in the far field contains 7% (16%) of the total electromagnetic energy of type 2 R (K) strokes.     

Multi‐Domain User‐Generated Content Based Model to Enrich Road Network Data for Multi‐Criteria Route Planning

By utilizing today's web‐based technologies, people can act as sensors and share their perceptions, emotions and observations in a variety of data forms, such as images, videos, texts, Global Positioning System (GPS) trajectories and maps. These forms are collectively called user‐generated content (UGC). These data are in different domains and have a multi‐modality nature. Although recent efforts have probed the acquisition of local knowledge by using single‐domain UGC data in specific applications, such efforts have not thus far presented a model considering multi‐domain UGC specifically to enrich road network data. This article aims at presenting such a model wherein, with the help of each data domain of UGC, one aspect of people knowledge about the road segment is obtained. These different aspects of knowledge are integrated using a Skyline operator to support multi‐criteria route finding. We name this model ERSBU (enriching road segments based on UGC). In ERSBU, road segments are basic spatial units, and their subjective properties have been extracted by using available UGC. The scenic score for each road segment was computed by using geo‐tagged Panoramio photos. The accessibility level of each road segment to different facilities was calculated based on data captured from Wikimapia and OpenStreetMap. Moreover, for measuring the movement popularity of each road segment, Wikiloc and Everytrail GPS trajectories were utilized. For the implementation of the ERSBU model, Tehran region 6 was considered the case study area. The Evaluation of the results proved that road segments that achieved a high score based on knowledge extracted from UGC also mostly gained top scores by analyzing traditional maps. ERSBU allows users to accomplish more‐qualitative path finding by considering the multi‐view characteristics of road segments.