Validation of HRDI MLT winds with meteor radars

A validation study of the mesospheric and lower-thermospheric (MLT) wind velocities measured by the High-Resolution Doppler Imager (HRDI) on board the Upper-Atmosphere Research Satellite (UARS) has been carried out, comparing with observations by meteor radars located at Shigaraki, Japan and Jakarta, Indonesia. The accuracy of the HRDI winds relative to the meteor radars is obtained by a series of simultaneous wind measurements at the time of UARS overpasses. Statistical tests on the difference in the wind vectors observed by HRDI and the meteor radars are applied to determine whether the wind speed has been overestimated by HRDI (or underestimated by the MF radars) as previously noticed in HRDI vs. MF radar comparisons. The techniques employed are the conventional t-test applied to the mean values of the paired wind vector components as well as wind speeds, and two nonparametric tests suitable for testing the paired wind speed. The square-root transformation has been applied before the Mests of the wind speed in order to fit the wind-speed distribution function to the normal distribution. The overall results show little evidence of overestimation by HRDI (underestimation by meteor radars) of wind velocities in the MLT region. Some exceptions are noticed, however, at the altitudes around 88 km, where statistical differences occasionally reach a level of significance of 0.01. The validation is extended to estimate the precision of the wind velocities by both HRDI and meteor radars. In the procedure, the structure function defined by the mean square difference of the observed anomalies is applied in the vertical direction for the profile data. This method assumes the isotropy and the homogeneity of variance for the physical quantity and the homogeneity of variance for the observational errors. The estimated precision is about 6m s^– for the Shigaraki meteor radar, 15 m s^–1 for the Jakarta meteor radar, and 20 m s^–1 for HRDI at 90-km altitude. These values can be used ot confirm the statistical significance of the wind field obtained by averaging the observed winds.     

Mean vertical wind in the mesosphere-lower thermosphere region (80–120 km) deduced from the WINDII observations on board UARS

The WINDII interferometer placed on board the Upper Atmosphere Research Satellite measures temperature and wind from the O(¹S) green-line emission in the Earths mesosphere and lower thermosphere. It is a remote-sensing instrument providing the horizontal wind components. In this study, the vertical winds are derived using the continuity equation. Mean wind annually averaged at equinoxes and solstices is shown. Ascendance and subsidence to the order of 1–2 cm s^–1 present a seasonal occurrence at the equator and tropics. Zonal Coriolis acceleration and adiabatic heating and cooling rate associated to the mean meridional and vertical circulations are evaluated. The line emission rate measured together with the horizontal wind shows structures in altitude and latitude correlated with the meridional and vertical wind patterns. The effect of wind advection is discussed.     

东亚MLT区域平均纬向风再评估——WINDII测量分析结果

本文利用1991年11月至1997年8月期间美国WINDII/UARS获得的风场测量数据对东亚上空纬向风进行考察. 研究结果给出了位于120°E 子午圈中90～120 km之间平均纬向风的典型结构及其季节特征,与在武汉开展流星雷达探测结果进行比较的结果说明卫星测量分析结果在对季节特征的描述方面与地基测量有相当好的一致性;较好的一致性还表现在与过去从HRDI/UARS数据中得到的月平均纬向风. 这些说明卫星探测结果有相当好的代表性. 与国际标准大气CIRA-86月平均纬向风开展比较的结果显示,从100 km高度开始这两种卫星数据分析结果都与CIRA-86结果表现出严重偏离,例如在赤道和低纬度地区某些高度,CIRA-86纬向风在全年的大部分时段中表现出与卫星数据分析结果风向不一致. 分析结果还显示WINDII纬向风和HRDI纬向风分析结果之间表现出一个幅度约20 m·s^-1的系统偏差,考虑到本文分析过程中采用了通过归并36天测量数据来消除周日变化影响的方案,同时参考其他研究工作中对MLT纬向风周日潮幅度的描述,两种卫星数据分析结果之间的系统偏差可能部分来自大气潮汐的影响.     

Thermal tides and studies to tune the mechanistic tidal model using UARS observations

Monthly simulations of the thermal diurnal and semidiurnal tides are compared to High-Resolution Doppler Imager (HRDI) and Wind Imaging Interferometer (WINDII) wind and temperature measurements on the Upper-Atmosphere Research Satellite (UARS). There is encouraging agreement between the observations and the linear global mechanistic tidal model results both for the diurnal and semidiurnal components in the equatorial and mid-latitude regions. This gives us the confidence to outline the first steps of an assimilative analysis/interpretation for tides, dissipation, and mean flow using a combination of model results and the global measurements from HRDI and WINDII. The sensitivity of the proposed technique to the initial guess employed to obtain a best fit to the data by tuning model parameters is discussed for the January and March 1993 cases, when the WINDII day and night measurements of the meridional winds between 90 and 110 km are used along with the daytime HRDI measurements. Several examples for the derivation of the tidal variables and decomposition of the measured winds into tidal and mean flow components using this approach are compared with previous tidal estimates and modeling results for the migrating tides. The seasonal cycle of the derived diurnal tidal amplitudes are discussed and compared with radar observation between 80 and 100 km and 40° S and 40°N.     

Numerical simulation of a mesolow over a Gulf Stream filament

A three-dimensional mesoscale numerical model is used to investigate mesoscale circulation over a Gulf Stream filament. Two numerical experiments are performed with different initial uniform ambient wind speeds (U=0.1 m s^–1, 3.5 m s^–1 and 7 m s^–1) for a typical winter day. It is found that for both low and moderate winds, a closed mesoscale circulation forms over the Gulf Stream filament. When the Gulf Stream filament was removed, the model did not predict a mesoscale circulation. The modeled circulation over the filament is in agreement with the observations, suggesting that the atmospheric circulations over the filaments may be an important mechanism in the U.S. East Coast cyclogenesis.     

The summertime 12-h wind oscillation with zonal wavenumber s = 1 in the lower thermosphere over the South Pole

Meteor radar measurements of winds near 95 km in four azimuth directions from the geographic South Pole are analyzed to reveal characteristics of the 12-h oscillation with zonal wavenumber one (s = 1). The wind measurements are confined to the periods from 19 January 1995 through 26 January 1996 and from 21 November 1996 through 27 January 1997. The 12-h s = 1 oscillation is found to be a predominantly summertime phenomenon, and is replaced in winter by a spectrum of oscillations with periods between 6 and 11.5 h. Both summers are characterized by minimum amplitudes (5–10 ms^–1) during early January and maxima (15–20 ms^–1) in November and late January. For 10-day means of the 12-h oscillation, smooth evolutions of phase of order 4–6 h occur during the course of the summer. In addition, there is considerable day-to-day variability (±5–10 ms^–1 in amplitude) with distinct periods (i.e., 5 days and 8 days) which suggests modulation by planetary-scale disturbances. A comparison of climatological data from Scott Base, Molodezhnaya, and Mawson stations suggests that the 12-h oscillation near 78°S is s = 1, but that at 68°S there is probably a mixture between s = 1 and other zonal wavenumber oscillations (most probably s = 2). The mechanism responsible for the existence of the 12-h s = 1 oscillation has not yet been identified. Possible origins discussed herein include in situ excitation, nonlinear interaction between the migrating semidiurnal tide and a stationary s = 1 feature, and thermal excitation in the troposphere.     

Observations of the spring–neap modulation of the gravitational circulation in a partially mixed estuary

The long-term variability of the non-tidal circulation in Southampton Water, a partially mixed estuary, was investigated using 71-day acoustic Doppler current profiler (ADCP) time series. The data show evidence that the spring–neap tidal variability of the turbulent mixing modulates the strength of the non-tidal residual circulation, with subtidal neap tide surface flows reaching 0.12 m s^–1 compared to <0.05 m s^–1 at spring tides. The amplitude of the neap-tide events in this non-tidal circulation is shown to be related to a critical value of the tidal currents, illustrating the strong dependence on tidal mixing. The results suggest that the dominant mechanism for generating these neap-tide circulation events is the baroclinic forcing of the horizontal density gradient, rather than barotropic forcing associated with ebb-induced periodic stratification. While tidal turbulence is thought to be the dominant control on this gravitational circulation, there is evidence of the additional effect of wind-driven mixing, including the effects of wind fetch and possibly wave development with along-estuary winds being more efficient at mixing the estuary than across-estuary winds. Rapid changes in atmospheric pressure also coincided with fluctuations in the gravitational circulation. The observed subtidal flows are shown to be capable of rapidly flushing buoyant material out of the estuary and into the coastal sea at neap tides.Responsible Editor: Iris Grabemann     

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.     

Lengths and hazards from channel-fed lava flows on Mauna Loa, Hawai‘i, determined from thermal and downslope modeling with FLOWGO

Using the FLOWGO thermo-rheological model we have determined cooling-limited lengths of channel-fed (i.e. a) lava flows from Mauna Loa. We set up the program to run autonomously, starting lava flows from every 4th line and sample in a 30-m spatial-resolution SRTM DEM within regions corresponding to the NE and SW rift zones and the N flank of the volcano. We consider that each model run represents an effective effusion rate, which for an actual flow coincides with it reaching 90% of its total length. We ran the model at effective effusion rates ranging from 1 to 1,000 m³ s^–1, and determined the cooling-limited channel length for each. Keeping in mind that most flows extend 1–2 km beyond the end of their well-developed channels and that our results are non-probabilistic in that they give all potential vent sites an equal likelihood to erupt, lava coverage results include the following: SW rift zone flows threaten almost all of Mauna Loas SW flanks, even at effective effusion rates as low as 50 m³ s^–1 (the average effective effusion rate for SW rift zone eruptions since 1843 is close to 400 m³ s^–1). N flank eruptions, although rare in the recent geologic record, have the potential to threaten much of the coastline S of Keauhou with effective effusion rates of 50–100 m³ s^–1, and the coast near Anaehoomalu if effective effusion rates are 400–500 m³ s^–1 (the 1859 a flow reached this coast with an effective effusion rate of 400 m³ s^–1). If the NE rift zone continues to be active only at elevations >2,500 m, in order for a channel-fed flow to reach Hilo the effective effusion rate needs to be 400 m³ s^–1 (the 1984 flow by comparison, had an effective effusion rate of 200 m³ s^–1). Hilo could be threatened by NE rift zone channel-fed flows with lower effective effusion rates but only if they issue from vents at 2,000 m or lower. Populated areas on Mauna Loas SE flanks (e.g. Phala), could be threatened by SW rift zone eruptions with effective effusion rates of 100 m³ s^–1.Editorial responsibility: J Donnelly-Nolan     

EISCAT measurements of the solar wind

EISCAT observations of interplanetary scintillation have been used to measure the velocity of the solar wind at distances between 15 and 130R (solar radii) from the Sun. The results show that the solar wind consists of two distinct components, a fast stream with a velocity of 800 km s^–1 and a slow stream at 400 kms^–1. The fast stream appears to reach its final velocity much closer to the Sun than expected. The results presented here suggest that this is also true for the slow solar wind. Away from interaction regions the flow vector of the solar wind is purely radial to the Sun. Observations have been made of fast wind/slow wind interactions which show enhanced levels of scintillation in compression regions.     

Substorm associated radar auroral surges: a statistical study and possible generation model

Substorm-associated radar auroral surges (SARAS) are a short lived (15–90 minutes) and spatially localised (5° of latitude) perturbation of the plasma convection pattern observed within the auroral E-region. The understanding of such phenomena has important ramifications for the investigation of the larger scale plasma convection and ultimately the coupling of the solar wind, magnetosphere and ionosphere system. A statistical investigation is undertaken of SARAS, observed by the Sweden And Britain Radar Experiment (SABRE), in order to provide a more extensive examination of the local time occurrence and propagation characteristics of the events. The statistical analysis has determined a local time occurrence of observations between 1420 MLT and 2200 MLT with a maximum occurrence centred around 1700 MLT. The propagation velocity of the SARAS feature through the SABRE field of view was found to be predominately L-shell aligned with a velocity centred around 1750 ms^–1 and within the range 500 m s^–1 and 3500 m s^–1. This comprehensive examination of the SARAS provides the opportunity to discuss, qualitatively, a possible generation mechanism for SARAS based on a proposed model for the production of a similar phenomenon referred to as sub-auroral ion drifts (SAIDs). The results of the comparison suggests that SARAS may result from a similar geophysical mechanism to that which produces SAID events, but probably occurs at a different time in the evolution of the event.     

Two-dimensional kinematic and rheological modeling of the 1912 pyroclastic flow,Katmai, Alaska

Pyroclastic flow emplacement is strongly influenced by eruption column height. A surface along which kinetic energy is zero theoretically connects the loci of eruption column collapse with all coeval ignimbrite termini. This surface is reconstructed as a two-dimensional energy line for the 1912 Katmai pyroclastic flow in the Valley of Ten Thousand Smokes from mapped flow termini and the runup of the ignimbrite onto obstructions and through passes. Extrapolation of the energy line to the vicinity of the source vent at Novarupta suggests the eruption column which generated the ignimbrite eruption was approximately 425 m high. The 1912 pyroclastic flow travelled about 25 km downvalley. Empirical velocity data calculated from runup elevations and surveyed centrifugal superelevations indicate initial velocities near Novarupta were greater than 79–88 m s^–1. The flow progressively decelerated and was travelling only 2–8 m s^–1 when it crossed a moraine 16 km downvalley. The constant slope of the energy line away from Novarupta suggests the flow was systematically slowed by internal and basal friction. Using a simple physical model to calculate flow velocities and a constant kinetic friction coefficient (Heim coefficient) of 0.04 derived from the reconstructed energy line, the flow is estimated to have decelerated at an average rate of –0.16 m s^–2 and to have taken approximately 9.5 minutes to travel 25 km down the Valley of Ten Thousand Smokes. The shear strength of the flowing ignimbrite at the moraine was approximately 0.5 kPa, and its Bingham viscosity when it crossed the moraine was 3.5 × 10³ P. If the flow was Newtonian, its viscosity was 4.2 × 10³ P. Reynolds and Froude numbers at the moraine were only 41–62 and 0.84–1.04, respectively, indicating laminar, subcritical flow.     

Comparisons between Canadian prairie MF radars, FPI (green and OH lines) and UARS HRDI systems

Detailed comparisons have been completed between the MF radars (MFR) in the Canadian prairies and three other systems: two ground-based Fabry-Perot interferometers (FPI) and the UARS high resolution Doppler imager (HRDI) system. The radars were at Sylvan Lake (52°N, 114°W), Robsart (49°N, 109°W) and the main continuing facility is at Saskatoon (52°N, 107°W). Statistical comparisons of hourly mean winds (1988–1992) for the Saskatoon MFR and FPI (557.7 nm green line) using scatter plots, wind speed-ratios, and direction-difference histograms show excellent agreement for Saskatoon. No serious biases in speeds or directions occur at the height of best agreement, 98 km. If anything, the MFR speeds appear bigger. The same applies to the Sylvan Lake MFR and Calgary FPI, where the best height is 88 km. In both cases these are close to the preferred heights for the emission layers. Differences between measurements seen on individual days are likely related to the influence of gravity waves (GW) upon the optical and radar systems, each of which have inherent spatial averaging (350, 50 km respectively), as well as the spatial difference between the nominal measurement locations. For HRDI, similar statistical comparisons are made, using single-overpass satellite winds and hourly means (to improve data quality) from MFR. Heights of best agreement, based upon direction-difference histograms, are shown; there is a tendency, beginning near 87 km, for these MFR heights to be 2 or 3 km greater than the HRDI heights. Speeds at these heights are typically larger for the satellite (MFR/HRDI = 0.7-0.8). Reasons for the differences are investigated. It is shown that the estimated errors and short-term (90 min) differences are larger for HRDI than for the MFR, indicating more noise or GW contamination. This leads to modest but significant differences in median speed-ratio (MFR/HRDI < 1). Also, comparison of the two systems is made under conditions when they agree best and when they show large disagreement. For the latter cases both systems show higher relative errors, and the HRDI vectors are frequently small. It is suggested that spatial or temporal GW wind fluctuations are the likely cause of the larger HRDI-MFR disagreement when wind speeds are small. No satisfactory explanation exists for the overall discrepancy is speeds between the MFR and HRDI.     

The oceans as source of cloud-forming nuclei

Summary The number and size of salt particles produced by the bursting of air bubbles in sea water has been measured. Bubbles of diameters varying between 1/2 and 2 mm each produced about 300 nuclei the sizes of which, under the electron microscope, were mainly between 0.1 and 0.5 diameter. They appeared to consist mainly of sodium chloride, the smallest ones containing only 10^–15 g of salt.These results, together with measurements of the size distribution of salt nuclei collected over the oceans in areas of spray formation, indicate that the total concentrations of salt nuclei over the oceans in winds of up to 15 m sec^–1 probably do not exceed 100 cm^–3. The corresponding rate of production of salt nuclei at the sea surface is estimated to be 1000 cm^–2 sec^–1. It is therefore inferred that sea spray contributes perhaps only one-fifth of the nuclei involved in cloud formation, the majority being the products of combustion, either natural or man-made.     

Extensions of Leovy's linear mesospheric circulation model

Linear mesospheric circulation models can quite realistically simulate mid-latitude mesospheric circulation and structure. In doing so they provide estimates for dissipation (Rayleigh friction 2×10^–6s^–1) and vertical mixing (30 m²s^–1) appropriate near the mesopause.     

The lunar tides in the Antarctic mesosphere and lower thermosphere

Horizontal winds in the mesosphere and lower thermosphere over the Antarctic have been measured by a meteor radar at Rothera (67.5°S, 68.0°W) and MF radar at Davis (68.6°S, 78.0°E). Data from Rothera recorded over a 20-month interval in 2005–2006 and data from Davis recorded over the 13-year interval 1994–2006 are examined to investigate the monthly mean behaviour of the lunar semidiurnal tide. Both data sets show a clear signal of the 12.42-h lunar semidiurnal (M₂) tide. The amplitude reaches values as large as 8 m s⁻¹. The vertical wavelengths of the tide vary seasonally from 10 to 65 km. Comparisons of the phase of the tide measured over the two sites reveals that it does not purely consist of a migrating wavenumber 2 mode. This suggests that other, non-migrating, modes are likely to be present.     

Differences between Mean Sea Levels for the Pacific,Atlantic and Indian Oceans From Topex/Poseidon Altimetry

Geopotential values W of the mean equipotential surfaces representing the mean ocean topography were computed on the basis of four years (1993 - 1996) TOPEX/POSEIDON altimeter data: W = 62 636 854.10m ² s ^–2 for the Pacific (P), W = 62 636 858.20m ² s ^–2 for the Atlantic (A), W = 62 636 856.28m ²s^–2 for the Indian (I) Oceans. The corresponding mean separations between the ocean levels were obtained as follows: A – P = – 42 cm, I– P = – 22 cm, I – A = 20 cm, the rms errors came out at about 0.3 cm. No sea surface topography model was used in the solution.     

Effects of steep topography on the flow and stratification near Palmyra Atoll

Two interdisciplinary cruises aimed at relating the ecology of marine fish populations to oceanographic conditions were fielded during the late summer and late winter seasons near Palmyra Atoll (5.9°N, 162.1°W) in the Line Islands. Ocean current and hydrographic measurements revealed interaction of the flow with the steep topography. During the first cruise (August/September 1990) satellite-tracked surface drifters and acoustic Doppler current profiler (ADCP) measurements showed a strong eastward setting North Equatorial Counter Current (NECC) with maximum speeds exceeding 1 m s^–1 at 80 m depth approximately. This current turned southeastward on closer approach to Palmyra. The drifter paths exhibited excursions with zonal wavelength of approximately 250 km, meridional amplitude of 25 km and period of approximately 5 days. During the second cruise (February/March 1992), the ADCP-derived speeds of the NECC were weaker (maxima approximately 33 cm s^–1) while the relative geostrophic flow component was of magnitude similar to 1990 and the signal of zonal geostrophic currents reached much deeper to approximately 650 m depth (150 m in 1990). Doming isopycnals beneath the surface mixed layer as well as thick (10–25 m) internal mixed layers were found near Palmyra during both cruises, with slightly different positions relative to the island. The discontinuous vertical temperature profiles may have been a result of strong boundary mixing due to breaking internal waves on Palmyras steep slopes. In the immediate vicinity of the island variations in flow speed, stratification and mixing in both the alongshore and cross-isobath directions were observed. Overall, the current speeds were reduced during February/March 1992, the peak time of the 1991–1993 warm event in the tropical Pacific. While parameters of turbulent two-dimensional wake theory are suggestive of formation and shedding of eddies in the lee of the island, no direct observations of circular motions were made in either expedition.Responsible Editor: Hans Burchard     

Latitudinal variations of neutral wind structures in the lower thermosphere for the March equinox period

Neutral winds in the lower thermosphere (95–$\text{130}\text{km}$) measured during the March equinox period (1991–1992) by ground-based incoherent scatter radars at Arecibo (18°N), Millstone Hill (42.5°N), and Sondrestrom (67°N) and by the space-based wind imaging interferometer (WINDII) are compared and show overall good agreement but some differences. At 18°N, the wind field in the altitude region of 95–$\text{110}\text{km}$ displays prevailing upward propagating diurnal tides with wavelengths of about $\text{22}\text{km}$. The diurnal structure is affected by the semidiurnal tide resulting in regular minima separated by 11–$\text{12}\text{h}$. At altitudes above $\text{110}\text{km}$, the diurnal tide dominant wind structure changes to the semidiurnal tide dominant structure as illustrated clearly by WINDII data with $\text{24}\text{h}$ coverage. Winds at 42.5°N and 67°N show similar structures in which winds at 105–$\text{115}\text{km}$ are generally anti-sunward. Daytime ISR winds show prevailing upward propagating semidiurnal tides with wavelengths of 35–$\text{70}\text{km}$. Winds from WINDII reveal the existence of the in situ thermospheric diurnal tide with amplitudes comparable to those of the semidiurnal tide. The superimposition of the two tides result in a wind field stronger during daytime than during nighttime at mid- and high-latitudes. Geomagnetic influence on neutral winds is negligible at low- and mid-latitudes under solar quiet conditions, but is observed at high-latitudes, where wind vectors follow a clockwise one-cell pattern at altitudes above about $\text{118}\text{km}$ in geomagnetic coordinates. Most recent simulations for the three latitudes provided by the NCAR thermosphere/ionosphere/mesosphere electrodynamics general circulation model are compared to the observations. The results at low- and mid-latitudes agree well with the observed winds in both wind structures and magnitudes, and reveal details of wave transition. Simulations for high-latitudes are less satisfactory, and require further improvements.     

CUTLASS HF radar observations of high-latitude azimuthally propagating vortical currents in the nightside ionosphere during magnetospheric substorms

High-time resolution CUTLASS observations and ground-based magnetometers have been employed to study the occurrence of vortical flow structures propagating through the high-latitude ionosphere during magnetospheric substorms. Fast-moving flow vortices (800 m s^–1) associated with Hall currents flowing around upward directed field-aligned currents are frequently observed propagating at high speed (1 km s^–1) azimuthally away from the region of the ionosphere associated with the location of the substorm expansion phase onset. Furthermore, a statistical analysis drawn from over 1000 h of high-time resolution, nightside radar data has enabled the characterisation of the bulk properties of these vortical flow systems. Their occurrence with respect to substorm phase has been investigated and a possible generation mechanism has been suggested.