A statistical study of underestimates of wind speeds by VHF radar

Comparisons are made between horizontal wind measurements carried out using a VHF-radar system at Aberystwyth (52.4°N, 4.1°W) and radiosondes launched from Aberporth, some 50 km to the southwest. The radar wind results are derived from Doppler wind measurements at zenith angles of 6° in two orthogonal planes and in the vertical direction. Measurements on a total of 398 days over a 2-year period are considered, but the major part of the study involves a statistical analysis of data collected during 75 radiosonde flights selected to minimise the spatial separation of the two sets of measurements. Whereas good agreement is found between the two sets of wind direction, radar-derived wind speeds show underestimates of 4–6% compared with radiosonde values over the height range 4–14 km. Studies of the characteristics of this discrepancy in wind speeds have concentrated on its directional dependence, the effects of the spatial separation of the two sets of measurements, and the influence of any uncertainty in the radar measurements of vertical velocities. The aspect sensitivity of radar echoes has previously been suggested as a cause of underestimates of wind speeds by VHF radar. The present statistical treatment and case-studies show that an appropriate correction can be applied using estimates of the effective radar beam angle derived from a comparison of echo powers at zenith angles of 4.2° and 8.5°.     

Temperature and ozone variations in the stratosphere

Examined are temperature and ozone variations in the Northern Hemisphere stratosphere during the period 1958–77, as estimated from radiosondes rocketsondes, ozonesondes, and Umkehr measurements. The temperature variation in the low tropical stratosphere is a combination of the variation associated with the quasi-biennial oscillation, and a variation nearly out of phase with the pronounced 3-yearly temperature oscillation (Southern Oscillation) present in the tropical troposphere since 1963. Based on radiosonde and rocketsonde data, the quasibiennial temperature oscillation can be traced as high as the stratopause, the phase varying with both height and latitude. However, the rocketsonde-derived temperature decrease of several degrees Celsius in the 25–55 km layer of the Western Hemisphere between 1969 (sunspot maximum) and 1976 (sunspot minimum) is not apparent in high-level radiosonde data, so that caution is advised with respect to a possible solar-terrestrial relation.There has been a strong quasi-biennial oscillation in ozone in the 8–16 km layer of the north polar region, with ozone minimum near the time of quasi-biennial west wind maximum at a height of 20 km in the tropics. A quasi-biennial oscillation in ozone (of similar phase) is also apparent from both ozonesonde data and Umkehr measurements in 8–16 and 16–24 km layers of north temperate latitudes, but not higher up. Both measurement techniques also suggest a slight overall ozone decrease in the same layers between 1969 and 1976, but no overall ozone change in the 24–32 km layer. Umkehr measurements indicate a significant 6–8% increase in ozone amount in all stratospheric layers between 1964 and 1970, and in 1977 the ozone amount in the 32–46 km layer was still 4% above average despite the predicted depletion due to fluorocarbon emissions. The decrease in ozone in the 32–46 km, layer of mid latitudes following the volcanic eruptions of Agung and Fuego is believed to be mostly fictitious and due to the bias introduced into the Umkehr technique by stratospheric aerosols of volcanic origin. Above-average water vapor amounts in the low stratosphere at Washington, DC, appear closely related to warm tropospheric temperatures in the tropics, presumably reflecting variations in strength of the Hadley circulation.     

Mid-latitude summer response of the middle atmosphere to short-term solar UV changes

Temperature and wind data obtained with Rayleigh lidar since 1979 and Russian rockets since 1964 are analyzed to deduce the summer response of the middle atmosphere to short-term solar UV changes. The equivalent width of the 1083 nm He I line is used as a proxy to monitor the short-term UV flux changes. Spectral analyses are performed on 108-day windows to extract the 27-day component from temperature, wind and solar data sets. Linear regressions between these spectral harmonics show some significant correlations around 45 km at mid-latitudes. For large 27-day solar cycles, amplitudes of 2 K and 6 m s⁻¹ are calculated for temperature data series over the south of France (44°N), and on wind data series over Volgograd (49°N), respectively. Cross-spectrum analyses have indicated correlations between these atmospheric parameters and the solar proxy with a phase lag of less than 2 days. These statistically correlative results, which provide good qualitative agreement with numerical simulations, are both obtained at mid-latitude. However, the observed amplitudes are larger than expected, with numerical models suggesting that dynamical processes such as equatorial or gravity waves may be responsible.     

Intercalibration of HRDI and WINDII wind measurements

The High Resolution Doppler Imager (HRDI) and the Wind Imaging Interferometer (WINDII) instruments, which are both on the Upper Atmosphere Research Satellite, measure winds by sensing the Doppler shift in atmospheric emission features. Because the two observation sets are frequently nearly coincident in space and time, each provides a very effective validation test of the other. Discrepancies due to geophysical differences should be much smaller than for comparisons with other techniques (radars, rockets, etc.), and the very large sizes of the coincident data sets provide excellent statistics for the study. Issues that have been examined include relative systematic offsets and the wind magnitudes obtained with the two systems. A significant zero wind position difference of 6 m s^–1 is identified for the zonal component, and it appears that this arises from an absolute perturbation in WINDII winds of -4 m s^–1 and in HRDI of +2 m s^–1. Altitude offsets appear to be relatively small, and do not exceed 1 km. In addition, no evidence is found for the existence of a systematic wind speed bias between HRDI and WINDII. However, considerable day-to-day variability is found in the quality of the agreement, and RMS differences are surprisingly large, typically in the range of 20-30 m s^–1.     

Seismic experiments on Showa-Shinzan lava dome using firework shots

Seismic experiments were conducted on Showa-Shinzan, a parasitic lava dome of volcano Usu, Hokkaido, which was formed during 1943–1945 activity. Since we found that firework shots fired on the ground can effectively produce seismic waves, we placed many seismometers on and around the dome during the summer festivals in 1984 and 1985. The internal structure had been previously studied using a prospecting technique employing dynamite blasts in 1954. The measured interval velocity across the dome in 1984 ranges 1.8–2.2 km/s drastically low compared to the results (3.0–4.0 km/s) in 1954; in addition, the velocity is 0.3–0.5 km/s higher than that in the surrounding area. The variation of the observed first arrival amplitudes can be explained by geometrical spreading in the high velocity lava dome. These observations show a marked change in the internal physical state of the dome corresponding to a drop in the measured highest temperature at fumaroles on the dome from 800°C in 1947 to 310°C in 1986.     

High time resolution monitoring of tropospheric temperature with a Radio Acoustic Sounding System (RASS)

We have observed the time-height variation of the temperature field in the upper troposphere using a Radio Acoustic Sounding System (RASS) which consists of the MU radar and a high-power acoustic transmitter. The fast beam steerability of the MU radar has made it possible to measure temperature profiles in a fairly wide height range in the upper troposphere (5–11 km), even under intense wind conditions. Observations were continued for about 32 hr on 24–26 December, 1986 with a time-height resolution of 30 min and 150 m. During the observation period, the tropospheric jet was so intense that the acoustic wavefronts were severely distorted. Using wind velocity profiles observed by the MU radar we have numerically estimated the propagation of acoustic wavefronts, and further determined favorable pointing directions for the MU radar to receive significant backscattering from refractive index fluctuations produced by the acoustic waves. Conventional radiosonde soundings were carried out every 6 hr, which showed a temperature decrease of 4 K/day in the upper troposphere during the observation period. Temperature profiles taken by RASS agree well with the radiosonde results.     

Thermal Structure of the Ionian Slab

We propose a thermal model of the subducting Ionian microplate. The slab sinks in an isothermal mantle, and for the boundary conditions we take into account the relation between the maximum depth of seismicity and the thermal parameter L_th of the slab, which is a product of the age of the subducted lithosphere and the vertical component of the convergence rate. The surface heat-flux dataset of the Ionian Sea is reviewed, and a convective geotherm is calculated in its undeformed part for a surface heat flux of 42 mW m^–2, an adiabatic gradient of 0.6 mK m^–1, a mantle kinematic viscosity of 10¹⁷ m² s^–1 and an asthenosphere potential temperature of 1300°C. The calculated temperature-depth distribution compared to the mantle melting temperature indicates the decoupling limit between lithosphere and asthenosphere occurs at a depth of 105 km and a temperature of 1260°C. A 70–km thick mechanical boundary layer is found. By considering that the maximum depth of the seismic events within the slab is 600 km, a L_th of 4725 km is inferred. For a subduction rate equal to the spreading rate, the corresponding assimilation and cooling times of the microplate are about 7 and 90 Myr, respectively. The thermal model assumes that the mantle flow above the slab is parallel and equal to the subducting plate velocity of 6 cm yr^–1, and ignores the heat conduction down the slab dip. The critical temperature, above which the subduced lithosphere cannot sustain the stress necessary to produce seismicity, is determined from the thermal conditions governing the rheology of the plate. The minimum potential temperature at the depth of the deepest earthquake in the slab is 730°C.     

New terrestrial heat flow values from Hungary

Summary Temperature and conductivity measurements show, that in the Southern part of Transdanubia (the part of Hungary which lies Westwards from Danube) the heat flow is about 2–2.4·10^–6 cal/cm² sec. Eastward from the Danube, in the Hungarian Plain estimates are even higher, and vary between 2.3·10^–6 and 2.8·10^–6 cgs. The gradient of temperature is everywhere quite high, 5.0 resp. 5.8·10^–4 deg. C/cm on the average. Thus, at a depth of 1000 m, the virgin rock temperature is about 60–70 deg. C, at 2000 m about 110–130 deg. C.     

Electric field effects on ionospheric and thermospheric parameters above the EISCAT station for summer conditions

Numerical calculations of the thermospheric and ionospheric parameters above EISCAT are presented for quiet geomagnetic conditions in summer. The Global Self-consistent Model of the Thermosphere, Ionosphere and Protonosphere (GSM TIP) was used. The numerical results were obtained both with a self-consistent calculation of the electric fields of magnetospheric and dynamo-action origin and with the magnetospheric electric fields only. It was found that the dynamo-electric field has some effect on the ionospheric convection pattern during quiet geomagnetic conditions. It has a marked effect mainly on the zonal neutral wind component above EISCAT (±20m/s at 140 km altitude). We have studied the effects of various field-aligned current (FAC) distributions on thermosphere/ionosphere parameters and we show that a qualitative agreement can be obtained with region-I and -II FAC zones at 75° and 65° geomagnetic latitude, respectively. The maximum FAC intensities have been assumed at 03–21 MLT for both regions with peak values of 2.5 × 10^–7 Am^–2 (region I) and 1.25 × 10^–7 A m^–2 (region II). These results are in agreement with statistical potential distribution and FAC models constructed by use of EISCAT data. The lack of decreased electron density in the night-time sector as observed by the EISCAT radar was found to be due to the spatial distribution of ionospheric convection resulting from electric fields of magnetospheric origin.     

Temperature modelling along the Taratashskiy profile crossing the ural mountains

Summary The temperature-depth distribution was calculated to a depth of 70 km along the 520 km long Taratashskiy refraction profile crossing the Ural Mts., approximately along latitude 56°N. The steady-state model was solved numerically using the finite-difference method, the vertical distribution of heat production was derived from the observed seismic velocities. It was shown that at the Moho boundary, the mantle heat flow varied between 10 and 25 mWm^–2, and the Moho temperature amounted to 300–550°C for the two versions studied.     

Changes of Ozone Content at Middle Latitudes During Forbush Decreases in Cosmic Rays

The variations of total ozone at Alma-Ata (43°N, 76 °E) and ozone profiles obtained by balloon sounding at Tateno (36°N, 140°E), Wallops Island (38°N, 75°W) and Cagliari (39°N, 9°E) in the periods of Forbush decreases (FD) in galactic cosmic rays have been analysed. A decrease of total ozone was observed in the initial stage of the FD and an increase 10–11 days later. The average total deviations calculated using the superposed epoch method for 9 FD events are equal to 30 D. U. in the positive and to –18 D. U. in the negative phase. The changes of average ozone profiles, associated with 26 FD events, are more significant in the lower stratosphere and upper troposphere. The decrease of the partial ozone pressure at a height of 12–15 km is about 30 mb. These vertical variations of ozone coincide with the average changes of the respective temperature profiles. A cooling, on the average, of 3°C was observed at 12–15 km, and a heating of 4°C below this level.     

Drifter observations of the Hebrides slope current and nearby circulation patterns

The mean flow at and around the Hebrides and Shetland Shelf slope is measured with ARGOS tracked drifters. Forty-two drifters drogued at 50 m were deployed in three circles over the Hebrides slope at 56.15°N in two releases, one on 5th December, 1995 and the second on 5–9th May, 1996. The circles span a distance of some 20 km from water depths of 200 m to 1200 m. Drifters are initially advected poleward along-slope by the Hebrides slope current at between 0.05 and 0.70 m s^–1 in a laterally constrained (25–50 km wide) jet-like flow. Drifters released in winter remained in the slope current for over 2000 km whilst summer drifters were lost from the slope current beyond the Wyville-Thomson Ridge, a major topographic feature at 60°N. Dispersion from the slope region into deeper waters occurs at bathymetric irregularities, particularly at the Anton Dohrn Seamount close to which the slope current is found to bifurcate, both in summer and winter, and at the Wyville-Thomson Ridge where drifters move into the Faeroe Shetland Channel. Dispersion onto the continental shelf occurs sporadically along the Hebrides slope. The initial dispersion around the Hebrides slope is remarkably sensitive to initial position, most of the drifters released in shallower water moving onto the shelf, whilst those in 1000 m or more are mostly carried away from the slope into deeper water near the Anton Dohrn Seamount. The dispersion coefficients estimated in directions parallel and normal to the local direction of the 500 m contour, approximately the position of the slope current core, are approximately 8.8 × 10³ m² s^–1 and 0.36 × 10³ m² s^–1, respectively, during winter, and 11.4 × 10³ m² s^–1 and 0.36 x 10³ m² s^–1, respectively, during summer. At the slope there is a minimum in across-slope mean velocity, Reynolds stress, and across-slope eddy correlations. The mean across-slope velocity associated with mass flux is about 4 × 10^–3 m s^–1 shelfward across the shelf break during winter and 2 × 10^–3 m s^–1 during summer. The drifters also sampled local patterns of circulation, and indicate that the source of water for the seasonal Fair Isle and East Shetland currents are the same, and drawn from Atlantic overflows at the Hebrides shelf.     

Seasonal variabilities of low-latitude mesospheric winds

Observations of mesospheric winds over a period of four years with the partial reflection radar at Tirunelveli (8.7°N, 77.8°E), India, are presented in this study. The emphasis is on describing seasonal variabilities in mean zonal and meridional winds in the altitude region 70–98 km. The meridional winds exhibit overall transequatorial flow associated with differential heating in the Northern and Southern Hemispheres. At lower altitudes (70–80 km) the mean zonal winds reveal easterly flow during summer and westerly flow during winter, as expected from a circulation driven by solar forcing. In the higher altitude regime (80–98 km) and at all altitudes during equinox periods, the mean zonal flow is subjected to the semi-annual oscillation (SAO). The interannual variability detected in the occurrence of SAO over Tirunelveli has also been observed in the data sets obtained from the recent UARS satellite mission. Harmonic analysis results over a period of two years indicate the presence of long-period oscillations in the mean zonal wind at specific harmonic periods near 240, 150 and 120 days. Results presented in this study are discussed in the context of current understanding of equatorial wave propagation.     

Subsurface Temperature Changes Due to the Crustal Magmatic Activity – Numerical Simulation

Geothermal aspects of the hypothesis, relating the earthquake swarms in the West Bohemia/Vogtland seismoactive region to magmatic activity, are addressed. A simple 1-D geothermal model of the crust was used to assess the upper limit of the subsurface heating caused by magma intrusion at the assumed focal depth of 9 km. We simulated the process by solving the transient heat conduction equation numerically, considering the heat of magma crystallization to be gradually released in the temperature interval 1100°C to 900°C. The temperature field prior to the intrusion was in steady-state with a surface temperature of 10°C and heat flow of 80 mWm ^–2 , the temperature at the 9 km depth was 270°C. The results suggest that the temperature and heat flow in the uppermost 1 km of the crust begin to grow 100 ka after the intrusion emplacement only, and that the amplitudes of the changes for the realistic lateral extent (a few kilometres) of the intrusion are very small. It was also found that the rate of magma solidification depends strongly on the thickness of the intrusion. It takes about 100 years for a 50 m thick sill to cool down from 1100°C to 600°C, which value represents the lower limit of the solidus temperature. The same cooling takes only 60 days if the sill is 2 m thick. If the nature of the strongly reflected boundaries, interpreted from the January 1997 Nový Kostel seismograms, is connected with the fresh emplacement of magma, the calculated cooling rates have a predictive potential for the temporal changes of the waveforms.     

First observation of quasi-2-day oscillations in ionospheric plasma frequency at fixed heights

The existence and development of the quasi-2-day oscillations in the plasma frequency variations of the F region at northern middle latitudes are investigated. A new approach to study the quasi-2-day oscillations is presented, using a methodology that allows us to do such a study at fixed heights. The hourly values of plasma frequency at fixed heights, from 170 km to 220 km at 10 km step, obtained at the Observatori de lEbre station (40.8°N, 0.5°E) during 1995 are used for analysis. It is found that quasi-2-day oscillations exist and persisted in the ionospheric plasma frequency variations over the entire year 1995 for all altitudes investigated. The dominant period of oscillation ranges from 42 to 56 h. The amplitude of oscillation is from 0.1 MHz to 1 MHz. The activity of the quasi-2-day oscillation is better expressed during the summer half year when several enhancements, about 15–30 days in duration, were observed. The largest enhancements of the oscillation occurred during early June, July and early August; i. e., near and after the summer solstice when the 2-day wave in the middle neutral atmosphere typically displays its largest activity in the Northern Hemisphere. The results obtained may help us understand better the possible influencing mechanisms between the 2-day wave in the middle neutral atmosphere and the ionospheric quasi-2-day oscillations.     

Fabry-Perot interferometer measurements of neutral winds and F2 layer variations at the magnetic equator

This letter presents some night-time observations of neutral wind variations at F2 layer levels near the dip equator, measured by the Fabry-Perot interferometer set up in 1994 at Korhogo (Ivory Coast, geographic latitude 9.25°N, longitude 355°E, dip latitude -2.5°). Our instrument uses the 630 nm (O¹D) line to determine radial Doppler velocities of the oxygen atoms between 200 and 400 km altitude. First results for November 1994 to March 1995 reveal persistent eastward flows, and frequent intervals of southward winds of larger than 50 ms^–1 velocity. Compared with the simultaneous ionospheric patterns deduced from the three West African equatorial ionosondes at Korhogo, Ouagadougou (Burkina-Faso, dip latitude +1.5°) and Dakar (Sénégal, dip latitude + 5°), they illustrate various impacts of the thermospheric winds on F2 layer density: (1) on the mesoscale evolution (a few 10³ km and a few 100 minutes scales) and (2) on local fluctuations (hundreds of km and tens of minutes characteristic times). We report on these fluctuations and discuss the opportunity to improve the time-resolution of the Fabry-Perot interferometer at Korhogo.     

First observation of mesospheric wind shear as high as 330 m s−1 km−1

Mesospheric wind profiles with an altitude resolution of 25 m have been obtained by means of radar tracking of foil chaff clouds. Such experiments were performed during winter 1990 at Biscarrosse, France (44^°N, 1^°W). On one flight, a wind shear as high as 330 m s⁻¹ km⁻¹ at 87.4 km and a region of dynamical instability between 86 and 88 km was measured. This wind shear is believed to be the largest value ever measured in the mesosphere. The region of dynamical instability results from a superposition of two wave motions, and is found to link well with enhanced turbulence and small-scale wave activity.     

Stratospheric warming effects on the tropical mesospheric temperature field

Temperature observations at 20–90 km height and 5–15°N during the winter of 1992–1993, 1993–1994 and 2003–2004, from the Wind Imaging Interferometer (WINDII) and Microwave Limb Sounder (MLS) experiments on the Upper Atmosphere Research Satellite (UARS) satellite and the Sounding the Atmosphere using Broadband Emission Radiometry (SABER) experiment on the Thermosphere, Ionosphere, Mesosphere Energetics and Dynamics (TIMED) satellite are analyzed together with MF radar winds and UK Meteorological Office (UKMO) assimilated fields. Mesospheric cooling is observed at the time of stratospheric warming at the tropics correlative with stratospheric warming events at middle and high latitudes. Planetary waves m=1 with periods of 4–5, 6–8, 10 and 12–18 days are found to dominate the period. Westward 7- and 16-day waves at the tropics appear enhanced by stationary planetary waves during sudden stratospheric warming events.     

Temperature of an active lava channel from spectral measurements,Kilauea Volcano,Hawaii

A narrow band spectroradiometer was used to determine the characteristic temperatures of a very active channeled lava flow for the phase 50 eruption of Pu'u 'O'o on the East Rift Zone of Kilauea Volcano, Hawaii. During the twilight of 19 February 1992, 14 spectra of this activity were acquired over a 51 minute interval [18.29 to 19.20 Hawaiian Standard Time (HST)], from which the thermal distribution of energy of two 18 m² areas, one near the center and one near the margin of the flow, may be investigated. A twocomponent thermal mixing model applied to the data taken of the center of the channel gave, in the most powerful instance (1.8x10⁵ W/m²), a crust temperature of 940° C, a hot component temperature of 1120°C and a hot radiating area of 60% of the total area. A simultaneous spectrum acquired near the channeled flow margin yielded a crust temperature of 586° C and a hot area of only 1.2% of the total area radiating at 1130° C. Average radiant flux densities recorded for the center of the lava channel (1.3x10⁵ W/m² average) are much greater than previous measurements of lava lakes (4.9x10³ W/m²) or recently emplaced lava flows (maximum of 7.2x10⁴ W/m²). The energetic nature of this eruption is shown by satellite measurements made at 02.33 HST on 22 February 1992 by the Advanced Very High Resolution Radiometer in Band 2 (0.72–1.10 m). These show the utility of using existing satellites with moderate resolution (1 km x 1 km pixels) and high temporal coverage (eight overpasses each day for Hawaii) as potential thermal alarms for rapidly assessing the hazard potential of large volcanic eruptions.     

The effect of water on stress relaxation of faulted and unfaulted sandstone

A series of stress relaxation experiments have been carried out on faulted and intact Tennessee sandstone to explore the influence of pore water on strength at different strain rates. Temperatures employed were 20, 300 and 400°C, effective confining pressure was 1.5 kb and strain rates as low as 10^–10 sec^–1 were achieved. Most samples were prefaulted at 2.5 kb confining pressure and room temperature. This is thought to have secured a reproducible initial microstructure.The strength of the dry rock was almost totally insensitive to strain rate in the range 10^–4 to 10^–10 sec^–1. In contrast, the strength of the wet rock decreased rapidly with strain rate at rates less than 10^–6 sec^–1. Brittle fracture of the quartz grains which constitute this rock is the most characteristic mode of failure under the test conditions used.The experimental data are discussed in terms of the possible deformation rate controlling processes, and it is suggested that in the wet experiments at intermediate to high strain rates (10^–7 to 10^–4 sec^–1) the observed deformation rate is controlled by the kinetics of water assisted stress corrosion, whilst deformation at low strain rates (ca. 10^–9 sec^–1) is controlled by a pressure solution process.The results have implications for the rheology of fault rocks at depths of perhaps 10 to 15 km in sialic crust.