The equatorial electrojet and the day-to-day variability of Sq

The day-to-day variability of Sq(H) at equatorial stations has been studied, using the correlation between daily ranges for the IQSY at a large number of pairs of stations. It is found that the daily ranges at stations under the equatorial electrojet are significantly less well correlated with those at other equatorial stations than is the case for pairs of equatorial stations not involving an electrojet station, for comparable distances of separation of the pairs of stations.     

Identifying bilateral pairs of deer (Odocoileus sp.) bones: how symmetrical is symmetrical enough?

Identifying bilaterally paired bones in zooarchaeological collections rests on the assumption that left and right skeletal elements from the same organism will be symmetrical. Bilaterally paired bones in an individual are, however, typically asymmetrical to some degree, demanding that the question “How symmetrical is symmetrical enough to identify a bilateral pair?” be answered with control data. The degree of symmetry chosen, or tolerance, will influence both how many true pairs in an archaeological collection are not identified (type I error) and how many false pairs are identified (type II error). Bivariate measures of 60 pairs of astragali and 48 pairs of distal humeri of deer (Odocoileus virginianus and O. hemionus) indicate that both sorts of error are frequent even when a conservative level of tolerance (≤average asymmetry) is used. Simulation of an archaeological collection indicates that as sample size increases, frequencies of both kinds of error increase. Application of the matching criteria and tolerance level to an archaeological collection underscores that the analytical requirements of identifying paired bones are steep.     

Kinetic analysis of thermoluminescence glow curves in feldspar: evidence for a continuous distribution of energies

The thermoluminescence (TL) glow curves from feldspars have been the subject of numerous studies, because of their importance in luminescence dating and dosimetry. This paper presents new experimental TL glow curves in a plagioclase feldspar, measured using the T _max-T _stop technique of glow curve analysis. Kinetic analysis of the experimental results is carried out for a freshly irradiated sample, as well as for a sample which has undergone optical treatment using infrared light for 100 s at 50°C. Application of the initial rise method of analysis indicates that the TL signals from both samples can be characterized by a continuous distribution of energy levels. By subtracting the TL glow curves measured at successive T _stop values, a series of TL glow curves is obtained which are analyzed using the empirical general order kinetics. It is found that all TL glow curves obtained by this subtractive procedure can be described accurately by the same general order parameter b ～1.7. In a second attempt to analyze the same TL glow curves and possibly extract information about the underlying luminescence process, the shape of TL glow curves is analyzed using a recently proposed physical kinetic model which describes localized electronic recombination in donor-acceptor pairs. Within this model, recombination is assumed to take place via the excited state of the donor, and nearest-neighbor recombinations take place within a random distribution of centers. This recent model has been used recently to describe successfully several types of luminescence signals. This paper shows that it is possible to obtain good fits to the experimental data using either one of these two approaches.     

Hemp in ancient rope and fabric from the Christmas Cave in Israel: talmudic background and DNA sequence identification

The “Christmas Cave”, a cave in the Qidron Valley near the Dead Sea and Qumran, has yielded a complex collection of plant-derived rope and fabric artifacts. Using polymerase chain reaction (PCR) to amplify DNA of the samples, we estimated the sizes and determined restriction patterns and base sequences of chloroplast genes, primarily rbcL (gene for the large subunit of ribulose bisphosphate carboxylase). DNA was successfully extracted from all samples, but was limited to sizes of approximately 200–300 base pairs. As expected, the DNA extracted from the samples was identified as coming primarily from flax (Linum usitatissamum L.), but two samples had a significant fraction, and all samples had at least a trace, of hemp (Cannabis sativa L.) DNA. Artifacts from the Christmas Cave were thought to date from Roman times, but it was thought possible that some could be much older. Accelerator mass spectrometry (AMS)-based ¹⁴C dating confirmed that the samples contained representatives from both the Roman and Chalcolithic periods. This paper provides a synthesis of DNA, isotope, and literary analysis to illuminate textile history at the Christmas Cave site.     

Seasonal dating and the valve-pairing technique in shell-midden analysis

Three-dimensional excavation of 23 artifact fragments and 4382 shells was carried out in the shell layer of an abandoned dwelling pit. Principal component analysis using three locational variables was used for drawing the central axes and the regression plane of the shell layer. Valve-pairing carried out on the clam Meretrix lusoria indicated that 380 pairs could be established among a total of 2089 unbroken valves. Preservation of shells was estimated at 70% by comparing the number of shells found in the shell layer with the clam population size which was in turn calculated by a probability formula using paired valves among samples. In addition to these data, the distance separating pairs gave some idea of the movement of shells. Seasonal dating by daily growth line counting showed that the shell layer could be subdivided into at least three main seasonal layers and that the period of accumulation was approximately 500 days.     

Riming electrification mechanism for charge generation within a thundercloud

A time dependent model for the rate of growth of the electric field within a thundercloud by the process of graupel formation (riming electrification) has been presented. The parameters of the model are: p₀, the precipitation intensity; q, the charge acquired by the graupel in each collision it makes with an ice crystal and 〈p〉, the charge transfer efficiency. Sets of values for the field growth and the maximum field have been obtained by varying these parameters.It is found that the estimated field for the initiation of a lightning flash (3.4 × 10⁵Vm⁻¹) within a time of about 1200 s can be achieved, taking a reasonable value of q = −1.67 × 10⁻¹⁴ C only for P₀ ≳ 5.56 × 10⁻⁶ms⁻¹(20 mmh⁻¹) and 〈p〉 ≳ 0.5. The maximum attainable electric field, E_max, itself is not sensitive to the value of q, within a reasonable range of variation in it (for a given p₀), but its rate of growth is: it grows faster, if g is larger.     

In-situ studies of equatorial spread-F over SHAR—steep gradients in the bottomside F-region and transitional wavelength results

RH-560 rockets instrumented with Langmuir probes were launched from SHAR, India (dip 11°N) for in-situ studies of electron density irregularities associated with equatorial spread-F (ESF) when the F-region plasma was drifting down and strong range spread-F was observed with an ionosonde at SHAR. A high variability was observed in the steepness of the base of the F-region. The bases were found to be steeper during the periods when the F-region plasma was drifting down. On one of the flights irregularities were observed in the region around 280 km where the gradients in electron density were downwards, indicating that the gradient drift instability is the main mechanism for their generation. Assuming a power law of the type P_k∝ kⁿ for irregularities of transitional scale (20–200 m), it was found that the spectral index n ranges between −1.5 and −4.6, when the mean integrated spectral power P_T of the irregularities in the above scale size range varied from −45 to −12 db. A relationship between n and P_T was observed and can be represented by a Gaussian function using the above expression; the altitude variation of n normalized for a P_T value of −10 db showed that the nature of spectral index remains the same between 230 km and the apogee of the rocket. This is at variance with the observations of Kelley et al. [(1982), J. geophys. Res. 87, 1575] that 280 km is the threshold altitude for the steep drift wave type of spectra to a shallower spectra.     

Comparison of incoherent scatter observations of electron density,and electron and ion temperature at Millstone Hill with the International Reference Ionosphere

Millstone Hill incoherent scatter (IS) observations of electron density (N_e, electron temperature (T_e) and ion temperature (T_i) are compared with the International Reference Ionosphere (IRI-86) for both noon and midnight, for summer, equinox and winter, at both solar maximum (1979–1980) and solar minimum (1985–1986). The largest difference inN_e is found in the topside, where values of N_e given by IRI-86 are generally larger than those obtained from IS measurements, by a factor which increases with increasing height, and which has a mean value near two at 600 km. Apart from the bottom of the profile, which is tied to the CIRA neutral temperature, the IRI-86 T_e model has no solar cycle variation. However, the IS measurements during the summer reveal larger T_e at solar maximum than at solar minimum. At other seasons higher T_e at solar maximum occurs only during the daytime at the greater heights. Nighttime T_e is shown by the IS radar to be generally larger in winter than in summer, an effect not included in the IRI. This is apparently due to photoelectron heating during winter from the sunlit ionosphere conjugate to Millstone Hill. The day-night difference in T_i given by IRI-86 above 600km is not as large in the IS measurements.     

THE OPTIMAL LOCATIONS OF MULTIPLE BRIDGES,CONNECTING FACILITIES,OR PRODUCT VARIETIES*

ABSTRACT This paper considers the optimal locations of two or more facilities, and the optimal number of facilities, when trips are made in pairs. The results are the same as standard models of spatial competition when there is perfect matching, but not when there is random matching. The first interpretation is bridges across a river, with residential locations on one side matched perfectly or randomly to jobs on the other side. The second interpretation is connecting facilities, such as tennis courts or restaurants where pairs of consumers meet. The third interpretation is product differentiation, with husbands and wives jointly choosing from among varieties.     

Effective electron collision frequency measurements in the E- and F-regions

The collision frequency v in the ionosphere has often been determined by measuring differences in the amplitude and group path of two closely spaced signals reflected in the region of high group retardation. In this paper we describe a method of measuring v using a CW double-side-band modulated signal reflected obliquely in the ionosphere. This allows v to be determined on a continuous basis and it is found that the value of v obtained is 1–5 × 10⁴ s⁻¹ for the E-region and ~ 10³ s⁻¹ for the F-region. It is shown that measurements made just after sunset, when the E-region is still present, are more representative of E-region values than F-region.     

Three solar cycles of day-time southern hemisphere Es activity

The hourly values of the ionosonde parameters f_oE_s and f_bE_s for times near mid-day have been examined over three solar cycles 1947–82. The daily parameters were studied for two southern hemisphere stations: Christchurch (temperate zone) (−43°.6′ geographic, −48°.1′ geomagnetic) and Rarotonga (subtropical) (−21°.2′ geographic, −20°.7′ geomagnetic). The data represent the longest such analyzed sporadi-E record in the literature. The seasonal variations of per cent occurrence of f₀E_s ⩾ 5 MHz and f_bE_s ⩾ 4 MHz for both stations show no sunspot dependence, strong activity in local summer and, in the case of Christchurch, minor enhancements in local winter. Long-term f_bE_s occurrence is closely associated with Zurich sunspot number R_z (correlation coefficients 0.93 for Christchurch and 0.76 for Rarotonga) with all seasons, showing the same in-phase R_z control for both stations. Long-term f_oE_s occurrence for both stations and for all seasons exhibits no dependence on R_z. There is evidence of a 6 year period in seasonal f_oE_s activity for Rarotonga. The data show a large and unexplained systematic decrease in f_oE_s occurrence for both stations and for all seasons by a factor of about 3 from 1947 to 1982.     

Fractal analysis of gravity wave spectra in the middle atmosphere

Wind fluctuations in the middle atmosphere behave like colored noise processes. They have a continuum of scales without dominant features and a power spectrum density (PSD) that often decays with frequency ƒ as ƒ^−β. Spectral index β is generally obtained through least-square fit to PSD estimated by Fourier methods. Graphs of colored noise have fractal plane-filling properties depending on β. An efficient method for finding β using the fractal dimension (D), based on analysis of 1/ƒ noise in galactic X-ray luminosities by McHardy I. and Czerny B., (1987, Nature325, 696), is described. An empirical relation is found between D and β and its validity is confirmed in limiting cases. Then D is obtained from power-law dependence of a length metric L(μ) on scale μ. The method is applied to middle-atmospheric velocity data from the Poker Flat radar in Alaska. Variations of D follow those in β, from an earlier analysis by Bemraet al., (1986, Handbook for MAP20, 216), but show an offset of 0.1–0.2 even after corrections for outliers, gaps, and additive noise. Usefulness of this method for screening data as an aid to spectral analysis is examined.     

Using the STARE double pulse system to obtain the mean velocity and velocity of long-lived irregularities

The STARE and SABRE auroral radars use double pulses to measure one value of the autocorrelation function (ACF) to obtain the Doppler velocity of the irregularities in the E-region auroral zone. If the pulse separation of each double pulse is τ, then the measured velocity tends towards the mean velocity V_DP as τ → 0, and towards the velocity of the long-lived irregularities V_LL as τ → ∞. A practical implementation using these results is presented, which allows measurement of these two physically important velocities from just four pulse spacings, rather than the 11 or 15 currently used to obtain the full velocity spectrum in the STARE and SABRE radars, and thus gains greater precision or better time resolution. The two line-of-sight velocities, V_DP and V_LL, are in general not equal, with V_LL greater in magnitude than V_DP (sometimes much greater). Comparison of V_DP and V_LL for the two STARE stations and one of the SABRE stations shows interesting differences, and a number of statistically significant features. V_LL probably measures the velocity of primary irregularities whenever these exist with the appropriate k vector. Both velocities ought to be measured routinely, and also for comparison with other techniques.     

Effects of the ionospheric conductivity on the longitudinal structure of the F-region equatorial anomaly

A longitudinal comparison has been undertaken of the vertical ion velocities in terms of ion mobilities and the eastward electrojet currents in terms of the electric conductivities overhead the magnetic equator at African and West Asian regions to search for longitudinal differences in the developments of the equatorial anamoly. It is shown that the vertical ion velocities and electrojet currents in both regions are different. It is proposed that longitudinal differences in (k²_1i+k²_2i)/k_2i and K_1i/K_2i ratios give rise to a different vertical ion velocities in these regions. This is likely to cause different equatorial F-layer plasma fountain magnitudes resulting in different developments of the equatorial anomaly.     

The D e (T,t) plot: A straightforward self-diagnose tool for post-IR IRSL dating procedures

This study presents a new self-diagnose method for the recently developed post-IR infrared stimulated luminescence (pIRIR) dating protocols. This criterion studies the dependence of equivalent dose (D _e) on measurement-temperature (T) and time (t), by applying the D _e (t) analysis to the IRLS and pIRIR signals measured under different temperatures, and combines these D _e (t) plots into one, so-called the D _e (T, t) plot. The pattern of the D _e (T, t) plot is shown to be affected by anomalous fading, partial bleaching and non-bleachable signal. A D _e plateau can be achieved in the D _e (T, t) plot only when the effects of these factors are insignificant. Therefore, this plot can be used as a self-diagnose tool for the validity of pIRIR results. The D _e (T, t) analysis has been applied to four recently developed pIRIR protocols, using aeolian samples with different ages. The results show that this self-diagnose tool can be applied to different pIRIR protocols for validating the pIRIR dating results and evaluating the pIRIR measurement conditions.     

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.     

Spread-Es structure producing apparent small scale structure in the F-region

When transmitting on 5.8 MHz the Bribie Island HF radar array synthesizes a beam that is 2.5 wide. The beam can be steered rapidly across the sky or left to dwell in any direction to observe the fading rates of echoes within a small cone of angles. With the beam held stationary, the time scale associated with deep fading of F-region echoes is usually more than 5 min. This is consistent with the focusing and defocusing effects caused by the passage of ever-present medium-scale travelling ionospheric disturbances (TIDs). On occasion the time scale for deep fading is much shorter, of the order of tens of seconds or less, and this is thought to be due to the interference of many echoes from within the beam of the radar. It is shown that the echoes are not due to scatter from fine structure in the F-region, but rather due to the creation of multiple F-region paths with differing phase lengths by small, refracting irregularities in underlying, transparent spread sporadic-E, (Spread-E_s). The natural drift of the Spread-E_s causes the phase paths of the different echoes to change in different ways causing the interference.Two methods are used to investigate the rapidly fading F-region signals. Doppler sorting of the refracted F-region signal does not resolve echoes in angle of arrival suggesting that many echoes exist within a Fresnel zone [Whitehead and Monro (1975), J. atmos. terr. Phys. 37, 1427]. Statistical analysis of F-region amplitude data indicates that when the range spread in E_s is severe on ionograms, then a modified Rayleigh distribution caused by the combination of 10 or so echoes is most appropriate. Using knowledge of the refracting process the scale of E_s structure is deduced from these results. Both methods find a Spread-E_s irregularity size of the order of 1 km or less. It is proposed that the Rayleigh type F-region signals seen by Jacobsonet al. [(1991b), J. atmos. terr. Phys. 53, 63] are F-region signals refracted by spread-E_s.     

Electron temperature and electron density in the F-region of the ionosphere. I. Observed relationship

Ionospheric data from three incoherent scatter stations over the height range 225–450 km were studied for all daylight hours over a wide range of solar conditions. The relationship between electron temperature T_e, electron density Nand solar flux at 10.7 cm wavelength S_10.7 was expressed as T_e = A−B·(N−5 × 10¹¹) + C·(S_10.7−750), where N is in units of m⁻³ and S_10.7 in kJy.This provided a very satisfactory expression for all data taken at Malvern and St. Santin between 0800 and 1600 LT. For data taken at Arecibo, however, the linearity broke down at low electron densities. The data from all three stations were therefore divided into two sets according to electron density and reexamined.ForN < 5 × 10¹¹ m⁻³ B increased steadily with height and decreased steadily with latitude.For N > 5 × 10¹¹ m⁻³ B did not appear to vary with height, with season or with latitude. C was approximately constant for all sets of data.The different mechanisms involved in the heat balance of the electron population are discussed and a qualitative explanation for the relationship is proposed.     

Modifications in the nighttime low-latitude ionosphere after southward turning of the IMF

Sharp decreases in Φ_oF2 are found to occur frequently in the nighttime low-latitude ionosphere after southward turning of the IMF B_z component, especially under isolated B_z turnings, i.e. when the IMF has been northward for at least 6 h before its turning. These decreases occur simultaneously (within a 1-h time interval) with the B_z turning. The effect is observed both when a substorm or a magnetic storm begins after B_z has turned southward, and when a noticeable substorm does not occur. The effect is more pronounced after midnight and a maximum at 03 LT. Short-term (with scale times of about 1 h) variations of Φ_oF2 and h_mF2 for B_z southward turning are analysed using a large amount of ground-based and topside sounding data. The decreases in Φ_oF2 are shown to occur at first over the magnetic equator and then, during the second hour after the turning, at the crests of the equatorial anomaly. The ionosphere returns to its undisturbed state, on average, in 4–5 h (if other disturbing agents do not arise). These decreases are suggested to be caused by modifications in the electric field in the low-latitude ionosphere associated with B_z southward turning.