Search

While knowledge of the energy inputs from the Sun (as it is the primary energy source) is important for understanding the solar-terrestrial system, of equal importance is the manner in which the terrestrial part of the system organizes itself in a quasi-equilibrium state to accommodate and re-emit this energy. The ROSMIC project (2014–2018 inclusive) was the component of SCOSTEP’s Variability of the Sun and Its Terrestrial Impact (VarSITI) program which supported research into the terrestrial component of this system. The four themes supported under ROSMIC are solar influence on climate, coupling by dynamics, trends in the mesosphere lower thermosphere, and trends and solar influence in the thermosphere. Over the course of the VarSITI program, scientific advances were made in all four themes. This included improvements in understanding (1) the transport of photochemically produced species from the thermosphere into the lower atmosphere; (2) the manner in which waves produced in the lower atmosphere propagate upward and influence the winds, dynamical variability, and transport of constituents in the mesosphere, ionosphere, and thermosphere; (3) the character of the long-term trends in the mesosphere and lower thermosphere; and (4) the trends and structural changes taking place in the thermosphere. This paper reviews the progress made in these four areas over the past 5 years and summarizes the anticipated research directions in these areas in the future. It also provides a physical context of the elements which maintain the structure of the terrestrial component of this system. The effects that changes to the atmosphere (such as those currently occurring as a result of anthropogenic influences) as well as plausible variations in solar activity may have on the solar terrestrial system need to be understood to support and guide future human activities on Earth. [Figure not available: see fulltext.][Figure not available: see fulltext.][Figure not available: see fulltext.] © 2021, The Author(s)., WW and AS gratefully acknowledge funding from the ISEE (Institute for Space-Earth Environmental Research), Nagoya University, for support to attend the VarSITI summarizing workshop where the results discussed in this review paper were presented. With funding from the Spanish government through the Severo Ochoa Centre of Excellence accreditation SEV-2017-0709.

Accepted: 2019-04-28, 資料番号: SA1190017000

We have performed a statistical data analysis on relationship between simultaneous Na density data and cosmic noise absorption (CNA) data, which is an indicator for energetic particle precipitation, obtained at Syowa, Antarctic in 2000–2002. It is found that the Na densities around the topside of Na layers (above $$\sim 95 \,\hbox {km}$$ ∼ 95 km height) tended to be smaller when the CNA was larger. The amounts of Na density responses, i.e., Na density decrease or Na loss, were increasing with magnetic local time (MLT) from dusk hours to dawn hours, and those of CNA responses, i.e., CNA increase, were also increasing with MLT. Thus, there were clear negative correlations between the Na density and CNA variations. These results indicate that the Na loss around the topside of Na layer would be induced by the energetic particle precipitation, and its effect would be more severe in dawn hours.

A mountain wave with a significant brightness temperature amplitude and ~500 km horizontal wavelength was observed over the Andes on 24–25 July 2017 in Atmospheric Infrared Sounder (AIRS)/Aqua satellite data. In the Modern‐Era Retrospective Analysis for Research and Applications, version 2 (MERRA‐2), reanalysis data, the intense eastward wind flowed over the Andes. Visible/Infrared Imaging Radiometer Suite (VIIRS)/Suomi‐NPP (National Polar‐orbiting Partnership) did not detect the mountain waves; however, it observed concentric ring‐like waves in the nightglow emissions at ~87 km with ~100 km wavelengths on the same night over and leeward of the Southern Andes. A ray tracing analysis showed that the mountain waves propagated to the east of the Andes, where concentric ring‐like waves appeared above a region of mountain wave breaking. Therefore, the concentric ring‐like waves were likely secondary waves generated by momentum deposition that accompanied mountain wave breaking. These results provide the first direct evidence for secondary gravity waves generated by momentum deposition.

We analyzed OH airglow images observed from an all-sky camera at King Sejong Station, Antarctica for the period of 2012–2016. Using M-transform method, 2D-power spectra of short period waves ...

We have studied atmospheric gravity waves (AGWs) and nighttime medium‐scale traveling ionospheric disturbances (MSTIDs) by applying three‐dimensional spectral analysis technique to 557.7‐ and 630.0‐nm airglow images at Shigaraki (SGK) (35°N, 136°E, 1999–2017) and Rikubetsu (RIK) (44°N, 144°E, 1999–2017), Japan, Athabasca (ATH), Canada, (55°N, 247°E, 2005–2017), and Magadan (MGD), Russia (60°N, 151°E, 2008–2017), focusing on their horizontal wavenumber spectra. For the AGWs in 557.7‐nm images, the power spectra in summer are stronger than in other seasons, probably due to stronger tropospheric convection. The highest energy content of the waves are mostly at wavelengths between 20 and 300 km at MGD, ATH, and RIK, while it is above 200 km at SGK. The largest power spectral density is obtained at RIK at wavelengths of 30–100 km and then ATH. The slopes of the horizontal wavenumber spectra varies from −2.77 to −3.22. From the MSTIDs in 630.0‐nm images, the power spectra in summer at RIK and SGK are stronger than those in other seasons regardless of solar activity. The power spectra in solar quiet time are stronger than those in solar active time at all four stations. These features can be explained by the Perkins instability with coupling between sporadic E and F layers. The spectral slope decreases with increasing latitudes. Weak positive correlations were obtained between the daily wave power of AGWs in 557.7‐nm images and MSTIDs in 630.0‐nm images, suggesting that the MSTIDs in the thermosphere may be partially generated by the AGWs from the mesopause region.

We developed user-friendly software based on Matsuda et al.'s (2014) 3D-FFT method (Matsuda-transform, M-transform) for airglow imaging data analysis as a function of Interactive Data Language (IDL). Users can customize the range of wave parameters to process when executing the program. The input for this function is a 3-D array of a time series of a 2-D airglow image in geographical coordinates. We applied this new function to mesospheric airglow imaging data with slightly different observation parameters obtained for the period of April–May at three different latitudes: Syowa Station, the Antarctic (69∘ S, 40∘ E); Shigaraki, Japan (35∘ N, 136∘ E); and Tomohon, Indonesia (1∘ N, 122∘ E). The day-to-day variation of the phase velocity spectrum at the Syowa Station is smaller and the propagation direction is mainly westward. In Shigaraki, the day-to-day variation of the horizontal propagation direction is larger than that at the Syowa Station; the variation in Tomohon is even larger. In Tomohon, the variation of the nightly power spectrum magnitude is remarkable, which indicates the intermittency of atmospheric gravity waves (AGWs). The average nightly spectrum obtained from April–May shows that the dominant propagation is westward with a phase speed  m s−1 at the Syowa Station and east-southeastward with a phase speed of up to ∼80 m s−1 in Shigaraki. The day-to-day variation in Tomohon is too strong to discuss average characteristics; however, a phase speed of up to ∼100 m s−1 and faster is observed. The corresponding background wind profiles derived from MERRA-2 indicate that wind filtering plays a significant role in filtering out waves that propagate eastward at the Syowa Station. On the other hand, the background wind is not strong enough to filter out relatively high-speed AGWs in Shigaraki and Tomohon and the dominant propagation direction is likely related to the distribution and characteristics of the source region, at least in April and May.

Simultaneous observation of gravity waves (GWs) in the polar summer mesosphere over Syowa (69°S, 39°E) by a ground-based radar and satellite instrument are presented. On 21 January 2016, at 2.3 UT, the CIPS instrument on the AIM satellite observed Polar Mesospheric Clouds (PMCs) with GW structures over Syowa. The orientation of the wave crests suggests north-west propagation direction. A periodogram analysis indicates GWs with horizontal wavelengths of 105 ± 5 and 60 ± 5 km. Cross-spectral analysis of the PANSY radar horizontal wind components at the same time, location, and altitude indicate two dominant waves with apparent periods of 15 and 9 h, with respective horizontal wavelengths of 100 ± 5 km and 302 ± 130 km. At the PMC altitude of ∼83 km, the mean zonal component of the vertical momentum flux spectra of the observed GWs is westward and the meridional component is northward. The comparable horizontal wavelengths of the 15 h wave observed by PANSY radar and the CIPS instrument, along with the agreement between the direction of vertical flux of horizontal momentum from PANSY, the propagation direction estimated from the cross-spectral analysis, and the orientation of the wave crests from CIPS i.e. north-west propagation direction suggests that the same wave is observed by both instruments. This simultaneous observation of the same wave by CIPS and PANSY radar provide a distinctive view of both the horizontal and vertical extent of a GW in the Antarctic summer mesosphere.

We have obtained horizontal phase velocity distributions of the gravity waves around 90 km from four Antarctic airglow imagers, which belong to an international airglow imager/instrument network known as ANGWIN (Antarctic Gravity Wave Instrument Network). Results from the airglow imagers at Syowa (69°S, 40°E), Halley (76°S, 27°W), Davis (69°S, 78°E), and McMurdo (78°S, 167°E) were compared, using a new statistical analysis method based on 3-D Fourier transform (Matsuda et al., 2014) for the observation period between 7 April and 21 May 2013. Significant day-to-day and site-to-site differences were found. The averaged phase velocity spectrum during the observation period showed preferential westward direction at Syowa, McMurdo, and Halley, but no preferential direction at Davis. Gravity wave energy estimated by I'/I was ~5 times larger at Davis and Syowa than at McMurdo and Halley. We also compared the phase velocity spectrum at Syowa and Davis with the background wind field and found that the directionality only over Syowa could be explained by critical level filtering of the waves. This suggests that the eastward propagating gravity waves over Davis could have been generated above the polar night jet. Comparison of nighttime variations of the phase velocity spectra with background wind measurements suggested that the effect of critical level filtering could not explain the temporal variation of gravity wave directionality well, and other reasons such as variation of wave sources should be taken into account. Directionality was determined to be dependent on the gravity wave periods. ©2017. The Authors.

The potential energy of gravity waves (GWs) per unit mass (Ep), at altitudes of 15–70 km, has been examined from temperature profiles obtained by a Rayleigh/Raman (RR) lidar at Syowa Station (69°S, 40°E) from May 2011 to October 2013, with the exception of the summer months. The GWs with ground-based wave periods longer than 2 h and vertical wavelengths between 1.8 and 16 km were extracted from the temperature profiles. Ep was larger in winter than in spring and fall, although in 2012, at altitudes below 30 km, Ep was larger in spring than in winter and fall. Ep increased with a mean scale height of 11.3 km. Ep profiles showed a local maximum at an altitude of 20 km and a minimum at 25 km in almost every month, which has not been reported by previous studies observed by radiosondes. The values of Ep in October of 2012 were smaller at 35–60 km and larger at 20–35 km than those in October of 2011 and 2013. This difference in the Ep profile is most probably caused by different seasonal variations of zonal winds. The larger and smaller Ep values seem to be observed both below and above the altitude at which the zonal wind speed reached 0 m s-1. This result suggests that wind filtering of gravity waves with small phase speeds is significantly important in early spring.

We analyzed the horizontal phase velocity of gravity waves and medium-scale traveling ionospheric disturbances (MSTIDs) by using the three-dimensional fast Fourier transform method developed by Matsuda et al. [2014] for 557.7-nm (altitude: 90–100 km) and 630.0-nm (altitude: 200–300 km) airglow images obtained at Shigaraki MU Observatory (34.8∘N, 136.1∘E, dip angle: 49∘) over ∼16 years from 16 March 1999 to 20 February 2015. The analysis of 557.7-nm airglow images shows clear seasonal variation of the propagation direction of gravity waves in the mesopause region. In spring, summer, fall, and winter, the peak directions are northeastward, northeastward, northwestward, and southwestward, respectively. The difference in east-west propagation direction between summer and winter is probably caused by the wind filtering effect due to the zonal mesospheric jet. Comparison with tropospheric re-analysis data shows the difference in north-south propagation direction between summer and winter is caused by differences in the latitudinal location of wave sources due to convective activity in the troposphere relative to Shigaraki. The analysis of 630.0-nm airglow images shows that the propagation direction of MSTIDs is mainly southwestward with a minor northeastward component throughout the 16 years. A clear negative correlation is seen between the yearly power spectral density of MSTIDs and F10.7 solar flux. This negative correlation with solar activity may be explained by the linear growth rate of the Perkins instability and secondary wave generation of gravity waves in the thermosphere.

Small-scale ripple structures observed in OH airglow images are most likely induced by either dynamic instability due to large wind shear or convective instability due to superadiabatic lapse rate. Using the data set taken in the mesopause region with an OH all-sky imager at Yucca Ridge Field Station, Colorado (40.7°N, 104.9°W), from September 2003 to December 2005, we study the characteristics and seasonal variations of ripple structures. By analyzing the simultaneous background wind and temperature observed by the nearby sodium temperature/wind lidar at Fort Collins, Colorado (40.6°N, 105°W), and a nearby medium-frequency radar at Platteville, Colorado (40.2°N, 105.8°W), we are able to statistically study the possible relation between ripples and the background atmosphere conditions. Characteristics and seasonal variations of ripples are presented in detail in this study. The occurrence frequency of ripples exhibits clear seasonal variability, with peak in autumn. The occurrence of ripples shows a local time dependence, which is most likely associated with the solar tides. The lifetime and spatial scale of these ripples are typically 5–20 min and 5–10 km, respectively, and most of the ripples move preferentially either southward or northward. However, more than half of the observed ripples do not advect with background flow; they have higher Richardson numbers than those ripples that advect with background flow. It is possible that they are not instability features but wave structures that are hard to be distinguished from the real instability features.

A joint special issue explores the potential of collaboration to help understand atmospheric gravity waves in the Polar Regions and their effect on global circulation.

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 t-tests 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 6ms–1 for the Shigaraki meteor radar, 15ms–1 for the Jakarta meteor radar, and 20ms–1 for HRDI at 90-km altitude. These values can be used to confirm the statistical significance of the wind field obtained by averaging the observed winds.

A Rayleigh–Raman lidar system was installed in January 2011 at Syowa Station, Antarctica (69.0° S, 39.6° E). Polar mesospheric clouds (PMCs) were detected by lidar at around 22:30 UTC (LT −3 h) on 4 February 2011, which was the first day of observation. This was the first detection of PMCs over Syowa Station by lidar. On the same day, a Super Dual Auroral Radar Network (SuperDARN) HF radar with oblique-incidence beams detected polar mesospheric summer echoes (PMSE) between 21:30 and 23:00 UTC. This event is regarded as the last PMC activity around Syowa Station during the austral summer season (2010–2011), since no other PMC signals were detected by lidar in February 2011. This is consistent with results of PMC and mesopause temperature observations by satellite-born instruments of AIM (Aeronomy of Ice in the Mesosphere)/CIPS (Cloud Imaging and Particle Size) and AURA/MLS (Microwave Limb Sounder) and horizontal wind measurements taken by a separate MF radar. Doppler velocity of PMSE observed by the HF radar showed motion toward Syowa Station (westward). This westward motion is consistent with the wind velocities obtained by the MF radar. However, the PMSE region showed horizontal motion from a north-to-south direction during the PMC event. This event indicates that the apparent horizontal motion of the PMSE region can deviate from neutral wind directions and observed Doppler velocities.

Mesospheric data were analyzed by a composite method combining phase and group velocity tracing technique and the spectra method of Stokes parameter analysis to obtain the propagation parameters of atmospheric gravity waves (AGW) in the height ranges between 63.6 and 99.3 km, observed using the MU radar at Shigaraki in Japan in the months of November and July in the years 1986, 1988 and 1989. The data of waves with downward phase velocity and the data of waves with upward phase velocity were independently treated. First, the vertical phase velocity and vertical group velocity as well as the characteristic wave period for each wave packet were obtained by phase and group velocity tracing technique. Then its horizontal wavelength, intrinsic wave period and horizontal group velocity were obtained by the dispersion relation. The intrinsic frequency and azimuth of wave vector of each wave packet were checked by Stokes parameters analysis. The results showed that the waves with intrinsic periods in the range 30 min–4.5 h had horizontal wavelength ranging from 25 to 240 km, vertical wavelength from 2.5 to 12 km, and horizontal group velocities from 15 to 60 m s−1. Both upward moving wave packets and downward moving wave packets had horizontal group velocities mostly directed in the sector between directions NNE (north-north-east) and SEE in the month of November, and mostly in the sector between directions NW and SWS in the month of July. Comparing with mean wind directions, the gravity waves appeared to be more likely to propagate along with mean wind than against it. This apparent prevalence for downstream wave packets was found to be caused by a systematic filtering effect existing in the process of phase and group velocity tracing analysis: A significant portion of upstream wave packets might have been Doppler shifted out of the vertical range in phase and group velocity tracing analysis.

Previous investigations of ionospheric storm-enhanced density (SED) and tongue of ionization (TOI) focused mostly on the behavior of TOI during intense geomagnetic storms. Little attention has been paid to the spatial and temporal variations of TOI during weak to moderate geomagnetic disturbed conditions. In this paper we investigate the source and development of TOI during a moderate geomagnetic storm on 14 October 2012. Multi-instrumental observations including GPS total electron content (TEC), Defense Meteorological Satellite Program (DMSP) in situ measured total ion concentration and ion drift velocity, SuperDARN measured polar ion convection patterns, and electron density profiles from the Poker Flat Incoherent Scatter Radar (PFISR) have been utilized in the current analysis. GPS TEC maps show salient TOI structures persisting for about 5 h over high latitudes of North America on 14 October 2012 in the later recovery phase of the storm when the magnitudes of IMF By and Bz were less than 5 nT. The PFISR electron density profiles indicate that the extra ionization for TEC enhancements mainly occurred in the topside ionosphere with no obvious changes in the bottomside ionosphere and vertical plasma drifts. Additionally, there were no signatures of penetration electric fields in the equatorial electrojet data and upward ion drifts at high latitudes. At the same time, strong subauroral polarization streams with ion drift speeds exceeding 2.5 km/s carried sunward fluxes and migrated toward lower latitudes for about 5° based on the DMSP cross-track drift measurements. Based on those measurements, we postulate that the combined effects of initial build-up of ionization at midlatitudes through daytime production of ionization and equatorward (or less poleward than normal daytime) neutral wind reducing downward diffusion along the inclined filed lines, and an expanded polar ion convection pattern and its associated horizontal plasma transport are important in the formation of the TOI.

Multiple instrumental observations including GPS total electron content (TEC), foF2 and hmF2 from ionosondes, vertical ion drift measurements from Communication/Navigation Outage Forecasting System, magnetometer data, and far ultraviolet airglow measured by Thermosphere, Ionosphere, Mesosphere Energetics and Dynamics/Global Ultraviolet Imager (TIMED/GUVI) are used to investigate the profound ionospheric disturbances at midlatitude and low latitude during the 14–17 July 2012 geomagnetic storm event, which was featured by prolonged southward interplanetary geomagnetic field component for about 30 h below −10 nT. In the East Asian/Australian sector, latitudinal profile of TEC variations in the main phase were characterized by three bands of increments and separated by weak depressions in the equatorial ionospheric anomaly (EIA) crest regions, which were caused by the combined effects of disturbance dynamo electric fields (DDEF) and equatorward neutral winds. In the recovery phase, strong inhibition of EIA occurred and the summer crest of EIA disappeared on 16 July due to the combined effects of intrusion of neutral composition disturbance zone as shown by the TIMED/GUVI O/N2 measurements and long-lasting daytime westward DDEF inferred from the equatorial electrojet observations. The transit time of DDEF over the dip equator from westward to eastward is around 2200 LT. In the American longitude, the salient ionospheric disturbances in the summer hemisphere were characterized by daytime periodical intrusion of negative phase for three consecutive days in the recovery phase, preceded by storm-enhanced density plume in the initial phase. In addition, multiple short-lived prompt penetration electric fields appeared during stable southward interplanetary magnetic field (IMF) Bz in the recovery phase and were responsible for enhanced the EIA and equatorial ionospheric uplift around sunset.

We have developed a new analysis method for obtaining the power spectrum in the horizontal phase velocity domain from airglow intensity image data to study atmospheric gravity waves. This method can deal with extensive amounts of imaging data obtained on different years and at various observation sites without bias caused by different event extraction criteria for the person processing the data. The new method was applied to sodium airglow data obtained in 2011 at Syowa Station (69°S, 40°E), Antarctica. The results were compared with those obtained from a conventional event analysis in which the phase fronts were traced manually in order to estimate horizontal characteristics, such as wavelengths, phase velocities, and wave periods. The horizontal phase velocity of each wave event in the airglow images corresponded closely to a peak in the spectrum. The statistical results of spectral analysis showed an eastward offset of the horizontal phase velocity distribution. This could be interpreted as the existence of wave sources around the stratospheric eastward jet. Similar zonal anisotropy was also seen in the horizontal phase velocity distribution of the gravity waves by the event analysis. Both methods produce similar statistical results about directionality of atmospheric gravity waves. Galactic contamination of the spectrum was examined by calculating the apparent velocity of the stars and found to be limited for phase speeds lower than 30 m/s. In conclusion, our new method is suitable for deriving the horizontal phase velocity characteristics of atmospheric gravity waves from an extensive amount of imaging data.

We report extreme enhancements of the nitric oxide (NO) column density observed with the ground-based millimeter-wave spectroscopic radiometer installed at Syowa Station, Antarctica, during a large geomagnetic storm in April 2012. From the NO spectrum line shape and NO column density relationship with solar radiation, we concluded that the NO was emitted in the altitude range between 75 km and 100 km. The column density of NO gradually increased during the recovery phase. In addition to variations on a time frame of several days, we found diurnal variations. The increase of NO was related to precipitated electrons in the energy range of 30–300 keV observed by Polar-orbiting Operational Environmental Satellite (POES)/The Meteorological Operational (METOP). We found a rapid response (within 1 h) and a one-to-one correspondence between them. For the first time, we show that a remarkable increase of the column density of NO is caused by dawn-dusk asymmetry of the plasma sheet electrons.

Using 5 years of all-sky OH airglow imager data over Yucca Ridge Field Station, CO (40.7°N, 104.9°W), from September 2003 to September 2008, we extract and deduce quasi-monochromatic gravity wave (GW) characteristics in the mesopause region. The intrinsic periods are clustered between approximately 4 and 10 min, and many of them are unstable and evanescent. GW occurrence frequency exhibits a clear semiannual variation with equinoctial minima, which is likely related to the seasonal variation of background wind. The anomalous propagation direction in January 2006, with strong southward before major warming starting in 21 January and weak southward propagation afterward, was most likely affected by stratospheric sudden warming. The momentum fluxes show strongly anticorrelated with the tides, with ~180° out of phase in the zonal component. While in the meridional component, the easterly maximum occurred approximately 2–6 h after maximum easterly tidal wind. However, the anticorrelations are both weakest during the summer. The dissipating and breaking of small-scale and high-frequency GW's components could have a potential impact on the general circulation in the mesopause region.

We report Geminid meteor head echo observations with the high-power large-aperture (HPLA) Shigaraki middle and upper atmosphere (MU) radar in Japan (34.85° N, 136.10° E). The MU radar observation campaign was conducted from 13 December 2010, 08:00 UTC to 15 December, 20:00 UTC and resulted in 48 h of radar data. A total of ~ 270 Geminids were observed among ~ 8800 meteor head echoes with precisely determined orbits. The Geminid head echo activity is consistent with an earlier peak than the visual Geminid activity determined by the International Meteor Organization (IMO). The observed flux of Geminids is a factor of ~ 3 lower than the previously reported flux of the 2009 Orionids measured with an identical MU~radar setup. We use the observed flux ratio to discuss the relation between the head echo mass–velocity selection effect, the mass distribution indices of meteor showers and the mass threshold of the MU radar.

Accepted: 2016-11-10, 資料番号: SA1160279000

We have developed an automated analysis scheme for meteor head echo observations by the 46.5 MHz Middle and Upper atmosphere (MU) radar near Shigaraki, Japan (34.85° N, 136.10° E). The analysis procedure computes meteoroid range, velocity and deceleration as functions of time with unprecedented accuracy and precision. This is crucial for estimations of meteoroid mass and orbital parameters as well as investigations of the meteoroid-atmosphere interaction processes. In this paper we present this analysis procedure in detail. The algorithms use a combination of single-pulse-Doppler, time-of-flight and pulse-to-pulse phase correlation measurements to determine the radial velocity to within a few tens of metres per second with 3.12 ms time resolution. Equivalently, the precision improvement is at least a factor of 20 compared to previous single-pulse measurements. Such a precision reveals that the deceleration increases significantly during the intense part of a meteoroid's ablation process in the atmosphere. From each received pulse, the target range is determined to within a few tens of meters, or the order of a few hundredths of the 900 m long range gates. This is achieved by transmitting a 13-bit Barker code oversampled by a factor of two at reception and using a novel range interpolation technique. The meteoroid velocity vector is determined from the estimated radial velocity by carefully taking the location of the meteor target and the angle from its trajectory to the radar beam into account. The latter is determined from target range and bore axis offset. We have identified and solved the signal processing issue giving rise to the peculiar signature in signal to noise ratio plots reported by Galindo et al. (2011), and show how to use the range interpolation technique to differentiate the effect of signal processing from physical processes.

Ripples as seen in airglow imagers are small wavy structures with short horizontal wavelengths (

Based on 5 years of OH imager data between September 2003 and September 2008 over Yucca Ridge Field Station, CO (40.7oN, 104.9oW), we presented the variation of gravity wave (GW) occurrence frequency and propagation direction in the upper mesosphere. In summer the GW occurrence frequency was extremely high at above 95% compared to other seasons (around 85%). The GW propagation direction showed a strong northward (poleward) preference in summer and a southward (equatorward) preference in winter. This could be possibly due to ducting of waves in the mesopause thermal structure and wave generation by the strong deep convection located at south side in summer and possible storms located at north side in winter. Westward traveling waves were rare, but eastward were frequent. In addition to seasonal variability, significant interannual variability was also observed.

OH rotational temperatures have been observed at the Syowa Station, Antarctica (69° S, 39° E), which is located in the middle of the auroral zone and has a high-sensitivity spectrometer for the spectral region of the OH 8-4 band. A dataset of 153 nights was acquired during the 2008 austral winter season. Of the 153 nights, the weather and aurora conditions were only suitable on 6 nights to study the relationship between auroral activity and OH airglow variation. Of these 6 nights, a significant increase in the rotational temperature and a decrease in the intensity related to an aurora activity were identified on the night of 27/28 March 2008, but no such variations were seen during the other nights. The horizontal magnetic field disturbance on the night of 27/28 March was the largest of that winter, while the cosmic radio noise absorption was also very strong. These facts indicate that, compared with the other nights, a large flux of high-energy auroral particles precipitated during the night. It is suggested that the observed variations in the OH rotational temperature and airglow intensity were caused by a lowering of the average airglow height as a result of OH depletion in the upper part of the layer where high-energy auroral particles can reach.

Turbulence generation mechanisms prevalent in the atmosphere are mainly shear instabilities, breaking of internal buoyancy waves, and convective instabilities such as thermal convection due to heating of the ground. In the present work, clear-air turbulence underneath a cirrus cloud base is described owing to coincident observations from the VHF (46.5 MHz) middle and upper atmosphere (MU) radar, a Rayleigh–Mie–Raman (RMR) lidar, and a balloon radiosonde on 7–8 June 2006 (at Shigaraki, Japan; 34.85°N, 136.10°E). Time–height cross section of lidar backscatter ratio obtained at 2206 LT 7 June 2006 showed the presence of a cirrus layer between 8.0 and 12.5 km MSL. Downward-penetrating structures of ice crystals with horizontal and vertical extents of 1.0–4.0 km and 200–800 m, respectively, have been detected at the cirrus cloud base for about 35 min. At the same time, the MU radar data revealed clear-air turbulence layers developing downward from the cloud base in the environment of the protuberances detected by the RMR lidar. Their maximum depth was about 2.0 km for about 1.5 h. They were associated with oscillatory vertical wind perturbations of up to ±1.5 m s−1 and variances of Doppler spectrum of 0.2–1.5 m−2 s−2. Analysis of the data suggests that the turbulence and the downward penetration of cloudy air were possibly the consequence of a convective instability (rather than a dynamical shear instability) that was likely due to sublimation of ice crystals in the subcloud region. Downward clear-air motions measured by the MU radar were associated with the descending protuberances, and updrafts were observed between them. These observations suggest that the cloudy air might have been pushed down by the downdrafts of the convective instability and pushed up by the updrafts to form the observed protuberances at the cloud base. These structures may be virga or perhaps more likely mamma as reported by recent observations of cirrus mamma with similar instruments and by numerical simulations.

On 9 October 2007, long-horizontal-wavelength gravity waves were observed for the first time to steepen and form mesospheric bores at the altitude of ∼87 km, by an all-sky OH imager located at Fort Collins (41°N, 105°W), Colorado. The collocated sodium lidar simultaneously observed the presence of a temperature inversion layer as the ducting region. One mesospheric bore uniquely later evolved into a large-amplitude soliton-like perturbation. When the gravity wave and the associated soliton-like perturbation passed through the lidar beams, the lidar detected strong vertical disturbance at 90 km, indicating convective instability. A large cold front system recorded several hours before in the troposphere was aligned to phase fronts of these large gravity waves. For all of the 7 mesospheric bores observed over a 5 year period, we found a similar alignment with a cold front 1000–1500 km away as the likely source of these large-scale gravity waves.

Significant attention is being paid in recent times by several observational and modeling studies to quantify the spatial and temporal variabilities of diurnal tide in the mesosphere and lower thermosphere (MLT) region. These variabilities are ascribed to spatial and temporal variations in the tidal forcing or interactions between the propagating tides and background wind, planetary waves or gravity waves. The present work makes use of simultaneous ground-based radar wind observations of different durations from five equatorial/low latitude sites in the Indian, Indonesian and Pacific sectors: Tirunelveli (8.7°N, 77.8°E), Jakarta (6.4°S, 106.7°E), Pontianak (0.03°N, 109°E), Kauai (22°N, 160°W) and Christmas Island (2°N, 157°W). This study delineates the longitudinal differences in the tidal characteristics in (i) interannual time scales over Tirunelveli and Kauai during 1993–2002, (ii) seasonal time scales over Christmas Island, Jakarta and Tirunelveli for the years 1993–1997 and (iii) shorter than seasonal time scales over Christmas Island, Pontianak and Tirunelveli during 1996–1997. An important observational feature noticed in this work is the differing behavior of the long-term tidal fields over Tirunelveli and Kauai. The monthly tidal amplitudes over Tirunelveli reveal a strong QBO signature whereas a similar, strong QBO signal could not be traced in the long-term observations from Kauai.

We present the characteristics of small-scale (

Using Global Positioning System (GPS) satellites, we have been conducting equatorial ionospheric scintillation observations at Kototabang, Indonesia since January 2003. Scintillations caused by equatorial plasma bubbles appear between 2000 and 0100 LT in equinoctial months with a seasonal asymmetry, and their activity decreases with decreasing solar activity. A comparison between scintillation index (S 4) and Earth’s brightness temperature (T bb) variations suggests that the scintillation activity can be related to tropospheric disturbances over the Indian Ocean to the west of Kototabang. To understand better the reasons of day-to-day variability of S 4, we analyze S 4, T bband lower thermospheric neutral wind ("Equation missing") data. The results show that S 4 fluctuates with periods of about 2.5, 5, 8, 14 and 25 days, possibly due to atmospheric waves from below and that similar periods are also found in the T bband "Equation missing" variations. Using a general circulation model, we made numerical simulations to determine the behavior of neutral wind in the equatorial thermosphere. The results indicate the following: (1) 2- to 20-day waves dissipate rapidly above about an altitude of 125 km, and 0.5- to 3-hour waves become predominant above 100 km, (2) zonal winds above 200 km altitude are, on the whole, eastward during sunset-sunrise, (3) zonal wind patterns due to short-period (1–4 h) atmospheric gravity waves (AGWs) above 120 km altitude change day by day, exhibit wavy structures with scale lengths of about 30–1000 km and, as a whole, move eastward at about 100−1 while changing patterns over time. These simulations suggest that the Rayleigh-Taylor instability responsible for plasma bubble generation can be seeded by AGWs with short periods of about 0.5–3 h, and that background conditions necessary for this instability are modulated by planetary-scale atmospheric waves propagating up to an altitude of about 120 km from below.

Gravity wave signatures were extracted from OH airglow observations using all-sky CCD imagers at four different stations: Cachoeira Paulista (CP) (22.7° S, 45° W) and São João do Cariri (7.4° S, 36.5° W), Brazil; Tanjungsari (TJS) (6.9° S, 107.9° E), Indonesia and Shigaraki (34.9° N, 136° E), Japan. The gravity wave parameters are used as an input in a reverse ray tracing model to study the gravity wave vertical propagation trajectory and to estimate the wave source region. Gravity waves observed near the equator showed a shorter period and a larger phase velocity than those waves observed at low-middle latitudes. The waves ray traced down into the troposphere showed the largest horizontal wavelength and phase speed. The ray tracing results also showed that at CP, Cariri and Shigaraki the majority of the ray paths stopped in the mesosphere due to the condition of m2m2m|→∞), which suggests the presence of ducting waves and/or waves generated in-situ. In the troposphere, the possible gravity wave sources are related to meteorological front activities and cloud convections at CP, while at Cariri and TJS tropical cloud convections near the equator are the most probable gravity wave sources. The tropospheric jet stream and the orography are thought to be the major responsible sources for the waves observed at Shigaraki.

Radiosonde observations of winds and temperature over several sites in Southeast Asia during the CPEA radiosonde observation period (10 Apri1-9 May, 2004) reveal the presence of a 7-day wave in the upper troposphere and lower stratosphere (UTLS), especially during the first half of the radiosonde observation period. Many of the characteristics of the wave resemble those of Kelvin waves. The wave amplitude peak is observed at altitudes of 20-21 km. The vertical phase structure of the wave is consistent with altitude over the radiosonde sites considered. The vertical wavelength of the wave is found to be in the range of 5.5-6.5 km. The correlative analysis among the different radiosonde sites and analysis of TIMED-SABER data sets reveals that the wave is more active in the longitude band of 0-180°E and has a zonal wave number of 3. A wave of similar periodicity and zonal structure observed in OLR suggests that tropical convection could be a source for these waves. A slight shift is observed in the peak of wave amplitude toward the northern hemisphere. Moreover, the latitudinal width of the wave is narrower than that predicted by theory. This is discussed based on the existing model results, which indicate that the latitudinal structure of the wave may be affected by strong vertical shear in zonal wind, such as that reported in the present work.

The MF and Meteor radar observations of horizontal winds over Koto Tabang (0.2°S, 100.3°E), Pontianak (0°, 109.3°E), and Pameungpeuk (7.5°S, 107.5°E) are utilized to study the characteristics of 5-8-day waves in the mesosphere and lower thermosphere (MLT) region over Indonesia. The wave activity is larger in zonal wind than in meridional wind and in the year 2003, it maximizes before and after northern hemispheric spring equinox. In the year 2004, the wave activity in zonal wind is larger during June-August 2004. During this period, the phase of the wave over Pontianak leads that over Koto Tabang, indicating that the wave is westward propagating with zonal wave number 1. The wave period ranges from 6.3 to 7.0 days over the equatorial sites Koto Tabang and Pontianak, whereas over the southern hemispheric site, Pameungpeuk, it is in the range 6.2-6.4 days. This shows that the dominant period of the wave changes with latitude. The amplitude and phase structures of the wave with altitude are similar among the three sites considered. The wave dominates the MLT dynamics over Pameungpeuk even during second half of CPEA (Coupling Processes in the Equatorial Atmosphere) campaign period (10 April-9 May 2004) with similar dominant wave period near 6.4-days. The wave amplitude is larger during last quarter of the campaign period.

The mechanism to explain the variations of tropopause and tropopause inversion layer (TIL) in the Arctic region during a sudden stratospheric warming (SSW) in 2009 was studied with the Modern-Era Retrospective analysis for Research and Applications reanalysis data and GPS/Constellation Observing system for Meteorology, Ionosphere, and Climate (COSMIC) temperature data. During the prominent SSW in 2009, the cyclonic system changed to the anticyclonic system due to the planetary wave with wave number 2 (wave2). The GPS/COSMIC temperature data showed that during the SSW in 2009, the tropopause height in the Arctic decreased accompanied with the tropopause temperature increase and the TIL enhancement. The variations of the tropopause and TIL were larger in higher latitudes. A static stability analysis showed that the variations of the tropopause and TIL were associated with the variations of the residual circulation and the static stability due to the SSW. Larger static stability appeared in the upper stratosphere and moved downward to the narrow region just above the tropopause. The descent of strong downward flow was faster in higher latitudes. The static stability tendency analysis showed that the strong downward residual flow induced the static stability change in the stratosphere and around the tropopause. The strong downwelling in the stratosphere was mainly induced by wave2, which led to the tropopause height and temperature changes due to the adiabatic heating. Around the tropopause, a pair of downwelling above the tropopause and upwelling below the tropopause due to wave2 contributed to the enhancement of static stability in the TIL immediately after the SSW.

The horizontal structure of the horizontal wind velocity is derived by the meteor echo observation with the MU radar. By dividing the meteor echoes into four azimuth sectors, horizontal wind velocities at four locations with about 100 km horizontal separation are determined with a time and height resolutions of 1 h and 2 km . The typical wind difference was about 20 m / s for 100 km separation with a significant temporal variation. Examples of the discrimination of tides and gravity waves with a same wave period was demonstrated using horizontal wind structure analysis. A simultaneous observation with the OH airglow imaging showed that the uniformity/non-uniformity of the wind field well corresponds to the uniformity/non-uniformity in the gravity wave in the airglow layer.

We have analyzed long-term variations of the wind velocity variance of ultra-fast Kelvin (UFK) waves (3.0– 3.8 d ) and gravity waves in two wave period ranges; 30– 35 h (0.7– 0.8 c/d ) and 8– 11 h (2.1–2.9 c/d ), in the equatorial mesosphere and lower thermosphere region using meteor wind radar (MWR) observations in Jakarta (JKT) ( 6° S 107° E ) from November 1992 to June 1998, as well as two medium frequency (MF) radars at Christmas Island ( 2° N 157° W ) and Pontianak (0.03° N 109° E ) from January 1996 to July 1997. The long-term variations of the wind variance of these waves at JKT showed enhancements twice a year in 1993–1995, mostly coinciding with the westward wind phase of MSAO or the transition period of MSAO from eastward to westward wind. The activity of both UFK and gravity waves became significantly weak in 1996–1997. In particular, the UFK wave variance became simultaneously smaller at JKT and the two MF radar sites, indicating that this peculiar event had a global scale.

Simultaneous observations of short-period gravity waves were carried out using three all-sky cooled-CCD imagers of OH airglow at Moshiri (northern edge of Japan) and at Shigaraki and Bisei (middle part of Japan) during the FRONT campaign of May 1998. These stations were separated by horizontal distances of 250–1300 km, and they offered a unique opportunity to investigate the spatial extent of gravity waves in the mesopause region. The images from Shigaraki and Bisei showed short-period gravity waves (horizontal wavelength: 20–40 km, phase front: ∼east-west aligned) moving northward (phase velocity: 20–40 m/s) on May 19, 21 and 22. At Moshiri, similar east-west waves moving northward were observed on May 19 and 21, indicating a broad horizontal extent of the gravity waves. Waves at Moshiri were not evident on May 22, however. The horizontal wind measured by the MF radar at Wakkanai (near Moshiri) showed height profiles that were similar to those at Yamagawa (southern edge of Japan) on May 19 and 21 but quite dissimilar on May 22. On the basis of these data, we suggest that short-period gravity waves in the mesopause region can have a broad spatial extent of more than 1000 km. The generation and propagation of internal gravity waves are discussed as a cause of the spatial extent of the waves.

This paper reports the development of coordinated observations carried out with the use of the MU radar (middle and upper atmosphere radar) and optical measurement techniques in Shigaraki. Since the start of the MU radar operation, mesosphere wind and turbulence were studied using mesospheric turbulence echo in the daytime. More recently, meteor echoes have been used to observe wind and temperature around the mesopause, during day and night. Comparative observations with the MU radar meteor mode and sodium lidars were used to study gravity waves and variation of sodium layers. Airglow observations were also carried out as a campaign basis, with which small scale and large scale gravity waves have been studied. We have also developed two CCD imagers to observe airglow images at the MU radar site for long term observations with the radar. Now Shigaraki has become an integrated optical/radar site ready to contribute to the PSMOS (Planetary Scale Mesopause Observing System) project being carried out between 1998–2002.

In order to clarify the horizontal structure of the wavelike oscillation frequently observed in the night time sodium density profile with a small Gaussian half-width in the middle of the night and a broad distribution at dusk, simultaneous observations with two lidars at Shigaraki (34.9°N, 136.1°E) and Hachioji (35.6°N, 139.4°E) and the MU radar at Shigaraki has been carried out. In the campaign of 35 nights, simultaneous observation was successful in four nights. On December 27–28, 1995, a large scale wave motion was observed by two lidars and the MU radar with meteor observation mode. Phase velocities were almost the same at the two sites and there were only slight differences in phase between the two sites. The wave motion was inconsistent with the component of the tidal wave and similar with the results of hodograph analyses. It is possible that the wave observed by lidars on December 27–28, 1995 was a gravity wave. The results of analyses suggest the possibility of gravity waves which were observed at fixed local time.

Simultaneous observations of the night airglow OH (6, 2) band emission intensity and rotational temperature, by a sky scanning airglow spectrophotometer, and meteor winds, by a middle and upper atmosphere radar (MU radar), were carried out at Shigaraki (34.9°N, 136.1°E), Japan, from October 29 to November 11, 1994, as the first phase of a campaign, and from July 25 to July 31, 1995 as the second phase. Horizontal structures in the OH emission intensity and rotational temperature were monitored optically, together with the background wind and its wave induced fluctuations, measured by MU radar. Since the MU radar makes a direct measurement of the vertical wavelength, and the OH spectrophotometer makes a direct measurement of the horizontal wavelength, the two techniques are mutually complementary to determine intrinsic wave parameters. Gravity waves with intrinsic periods of 2 to 9 hours, horizontal wavelengths of 500 to 3000 km and vertical wavelengths of 12 to 75 km were identified. Between the two different observation techniques, there is a reasonable agreement in the inferred wave characteristics.

We conducted meteor echo observations using the Buckland Park MF radar (35°S, 138°E) at 00:40–05:45 LT on October 22, 1997. In addition to the usual full correlation analysis (FCA) technique to measure horizontal wind velocities from 60 to 100 km MF radars have a potential to detect meteor echoes and infer winds through their Doppler frequency shifts. Because of the relatively low radio frequency employed MF radars have a great advantage of providing meteor wind well above 100 km altitude, where very few techniques can measure wind velocities. There is a limitation which should be noted as well. The observations are possible only during night time when the electron density of E-region is low enough for the radio wave to penetrate into the upper region. We detected 233 underdense meteor echoes from 80 km to 120 km with a mean height of 104.4 km. Although the transmitting antenna beams were steered toward off-zenith angles of 25°, almost all the echoes were received outside of the main lobe, indicating that conventional MF radar systems with a broad transmitting beam can work well for meteor observations. Bi-hourly wind profiles were obtained from 94 to 114 km altitudes. The profiles revealed a clear wave structure with a downward phase progression with time. FCA winds from 80 to 100 km were also estimated, and a continuous wind structure was obtained from FCA to meteor heights. Note that the present observations happened to be conducted during a major meteor shower activity. However, a majority of the underdense echoes were from non-shower meteors, and observations during non-shower periods will also yield enough information.