Using S-band dual polarized radar for convective/stratiform rain indexing and the correspondence with AMSR-E GSFC profiling algorithm

The separation of rain types in convective and stratiform regimes has long been a goal in microwave remote sensing of precipitation research. In this essence, a dual polarized radar based indexing scheme that provides information on convective and stratiform (C/S) rain regimes has been presented in correspondence with advanced microwave scanning radiometer – earth observing system (AMSR-E) GSFC profiling algorithm estimate of convective rain percentage. The dual polarized radar based C/S indexing scheme first retrieves the normalized gamma drop size distribution parameters, median volume drop diameter (D₀) and concentration parameter (N_w), from dual polarized radar measurements Z_H and Z_DR, representing reflectivity and differential reflectivity respectively, by means of the genetic programming approach. Next, the C/S rain index is calculated based on the formulation of an empirical relation in N_w–D₀ domain. The scheme has been inspected and applied on measurements from the S-band Chilbolton dual polarized radar. A considerable number of “coincident” cases from the radar and the AMSR-E observations are investigated. It has been revealed that the dual polarized radar based C/S rain indexing is in a similar pattern with the AMSR-E GSFC profiling algorithm estimate of convective rain percentage. Generally, as C/S rain index value increases, which signifies a stratiform to convective trend, the AMSR-E convective rain percentage also increases.     

0 °C isotherm height for satellite communication in Malaysia

The 0 °C isotherm height is an important parameter for prediction of rain attenuation of microwave and millimeter wave for Earth-space communication. The variations of 0 °C isotherm heights for two monsoon seasons have been studied based on an analysis of radiosonde over three stations. The exceedence probability statistics of rain height are compared between the two seasons. The results on the 0 °C isotherm height can be utilized for the estimation of attenuation of microwave and millimeter wave due to rain over Earth-space paths. Attenuations of radio wave due to rain at frequencies above 10 GHz and above have also been estimated using the 0 °C isotherm height so derived.     

Variation of rain drop size distribution with rain rate at a few coastal and high altitude stations in southern peninsular India

Rain drop size distribution (DSD) was measured at four places in Southern India {Thiruvananthapuram, Kochi, Munnar and Sriharikota (SHAR)} using a Joss–Waldvogel (JW) impact type disdrometer. The data for each minute were corrected for dead time errors and rain rate was computed from the corrected data. The data for a whole month were then sorted according to rain rate (R) into several classes ranging from 0.1 to >100 mm/h. The average DSD in each class was computed, and the lognormal distribution function was fitted to the average. In all the cases, the function fitted the data very well. The fit parameters were found to have dependence on rain rate. The total number of drops (N_T), the geometric mean diameter (D_g) and the standard geometric deviation (σ) were also computed from the fit parameters. The standard geometric deviation (σ) was found to be more or less constant with rain rate at all the sites and in all months. The other two parameters (N_T and D_g) were found to vary exponentially with rain rate except in Munnar, a high altitude station. At Thiruvananthapuram, in most of the months, N_T increased exponentially with rain rate up to some value of R, which was different in different months, and then remained more or less constant or decrease slightly. In all cases, the variation of N_T and D_g was such that N_TD_g³ increased linearly with rain rate.     

A multilayer perceptron and multiclass support vector machine based high accuracy technique for daily rainfall estimation from MSG SEVIRI data

The current paper introduces a new multilayer perceptron (MLP) and support vector machine (SVM) based approach to improve daily rainfall estimation from the Meteosat Second Generation (MSG) data. In this study, the precipitation is first detected and classified into convective and stratiform rain by two MLP models, and then four multi-class SVM algorithms were used for daily rainfall estimation. Relevant spectral and textural input features of the developed algorithms were derived from the spectral MSG SEVIRI radiometer channels. The models were trained using radar rainfall data set colected over north Algeria. Validation of the proposed daily rainfall estimation technique was performed by rain gauge network data set recorded over north Algeria. Thus, several statistical scores were calculated, such as correlation coefficient (r), root mean square error (RMSE), mean error (Bias), and mean absolute error (MAE). The findings given by: (r = 0.97, bias = 0.31 mm, RMSE = 2.20 mm and MAE = 1.07 mm), showed a quite satisfactory relationship between the estimation and the respective observed daily precipitation. Moreover, the comparison of the results with those of two advanced techniques based on random forests (RF) and weighted ‘k’ nearest neighbor (WkNN) showed higher accuracy obtained by the proposed model.     

An empirical model for the variation of rain drop size distribution with rain rate at a few locations in southern India

A Joss–Waldvogel impact type disdrometer was installed at four different locations in the Indian peninsula during various periods from 2001 till date. The data are analysed to study the nature of rain drop size distribution (DSD) in this region. Out of the three well known distributions that describe DSD, namely, the Marshall–Palmer, Gamma and Lognormal, it has been found that Lognormal distribution fits the DSD in this region better than the other ones. Lognormal distributions for different rain rates were then derived by fitting the lognormal function to the data using a curve fitting software. Then the variation of fit parameters with rain rate was evaluated. Incorporating these variations, into the Lognormal distribution, an empirical equation that describes the DSD in this region for different rain rates was derived. Then this equation was tested with sample data from each of these stations. The data used for validation were not used for fitting lognormal equation to derive the fit parameters. The correlation between the DSD measured and derived using the empirical model was found to be quite good (0.9) except in some cases where the coefficient dropped to 0.75. The empirical model can be updated when more data are available.     

Study of daytime vertical E × B drift velocities inferred from ground-based magnetometer observations of ΔH,at low latitudes under geomagnetically disturbed conditions

In this study, 30 storm sudden commencement (SSC) events during the period 2001–2007 for which daytime vertical E × B drift velocities from JULIA radar, Jicamarca (geographic latitude 11.91°S, geographic longitude 283.11°E, 0.81°N dip latitude), Peru and ΔH component of geomagnetic field measured as the difference between the magnitudes of the horizontal (H) components between two magnetometers deployed at two different locations Jicamarca (geographic latitude 11.91°S, geographic longitude 283.11°E, 0.81°N dip latitude) and Piura (geographic latitude 5.21°S, geographic longitude 279.41°E, 6.81°N dip latitude), in Peru, were considered. It is observed that a positive correlation exists between peak value of daytime vertical E × B drift velocity and peak value of ΔH for the three consecutive days of SSC. A qualitative analysis made after selecting the peak values of daytime vertical E × B drift velocity and ΔH showed that 57% of the events have daytime vertical E × B drift velocity peak in the magnitude range 20–30 m/s and 63% of the events have ΔH peak in the range 80–100 nT. The maximum probable (45%) range of time of occurrence of peak value for both vertical E × B drift velocity and ΔH during the daytime hours were found to be the same, i.e., 10:00–12:00 LT. A strong positive correlation was also found to exist between the daytime vertical E × B drift velocity and ΔH for all the three consecutive days of SSC, for all the events considered. To establish a quantitative relationship between day time vertical E × B drift velocity and ΔH, linear and polynomial (order 2 and 3) regression analysis (Least Square Method (LSM)) were carried out, considering the fully disturbed day after the commencement of the storm as ‘disturbed period’ for the SSC events selected for analysis. The formulae indicating the relationship between daytime vertical E × B drift velocity and ΔH, for the ‘disturbed periods’, obtained through the regression analysis were verified using the JULIA radar observed E × B drift velocity for 3 selected events. Root Mean Square (RMS) error analysis carried out for each case suggest that polynomial regression (order 3) analysis provides a better agreement with the observations from among the linear, polynomial (order 2 and 3) analysis.     

Impact of rain on wave retrieval from Sentinel-1 synthetic aperture radar images in tropical cyclones

The main objective of our work was to investigate the impact of rain on wave observations from C-band (~5.3 GHz) synthetic aperture radar (SAR) in tropical cyclones. In this study, 10 Sentinel-1 SAR images were available from the Satellite Hurricane Observation Campaign, which were taken under cyclonic conditions during the 2016 hurricane season. The third-generation wave model, known as Simulating WAves Nearshore (SWAN) (version 41.31), was used to simulate the wave fields corresponding to these Sentinel-1 SAR images. In addition, rainfall data from the Tropical Rainfall Measuring Mission satellite passing over the spatial coverage of the Sentinel-1 SAR images were collected. The simulated results were validated against significant wave heights (SWHs) from the Jason-2 altimeter and European Centre for Medium-Range Weather Forecasts data, revealing a root mean square error (RMSE) of ~0.5 m with a 0.25 scatter index. Winds retrieved from the VH-polarized Sentinel-1 SAR images using the Sentinel-1 Extra Wide-swath Mode Wind Speed Retrieval Model after Noise Removal were taken as prior information for wave retrieval. It was discovered that rain did indeed affect the SAR wave retrieval, as evidenced by the 3.21-m RMSE of SWHs between the SAR images and the SWAN model, which was obtained for the ~1000 match-ups with raindrops. The raindrops dampened the wave retrieval when the rain rate was < ~5 mm/hr; however, they enhanced wave retrieval for higher rain rates. It was also found that the portion of the rain-induced ring wave with a wave number > 0.05 rad/m (~125 m wavelength) was clearly observed in the SAR-derived wave spectra.     

Antarctic meteor observations using the Davis MST and meteor radars

This paper presents the meteor observations obtained using two radars installed at Davis (68.6°S, 78.0°E), Antarctica. The Davis MST radar was installed primarily for observation of polar mesosphere summer echoes, with additional transmit and receive antennas installed to allow all-sky interferometric meteor radar observations. The Davis meteor radar performs dedicated all-sky interferometric meteor radar observations. The annual count rate variation for both radars peaks in mid-summer and minimizes in early Spring. The height distribution shows significant annual variation, with minimum (maximum) peak heights and maximum (minimum) height widths in early Spring (mid-summer). Although the meteor radar count rate and height distribution variations are consistent with a similar frequency meteor radar operating at Andenes (69.3°N), the peak heights show a much larger variation than at Andenes, while the count rate maximum-to-minimum ratios show a much smaller variation. Investigation of the effects of the temporal sampling parameters suggests that these differences are consistent with the different temporal sampling strategies used by the Davis and Andenes meteor radars. The new radiant mapping procedure of [Jones, J., Jones, W., Meteor radiant activity mapping using single-station radar observations, Mon. Not. R. Astron. Soc., 367(3), 1050–1056, doi: 10.1111/j.1365-2966.2006.10025.x, 2006] is investigated. The technique is used to detect the Southern delta-Aquarid meteor shower, and a previously unknown weak shower. Meteoroid speeds obtained using the Fresnel transform are presented. The diurnal, annual, and height variation of meteoroid speeds are presented, with the results found to be consistent with those obtained using specular meteor radars. Meteoroid speed estimates for echoes identified as Southern delta-Aquarid and Sextantid meteor candidates show good agreement with the theoretical pre-atmospheric speeds of these showers (41 km s⁻¹ and 32 km s⁻¹, respectively). The meteoroid speeds estimated for these showers show decreasing speed with decreasing height, consistent with the effects of meteoroid deceleration. Finally, we illustrate how the new radiant mapping and meteoroid speed techniques can be combined for unambiguous meteor shower detection, and use these techniques to detect a previously unknown weak shower.     

Rainfall estimation over a Mediterranean region using a method based on various spectral parameters of SEVIRI-MSG

The ultimate objective of this paper is the estimation of rainfall over an area in Algeria using data from the SEVIRI radiometer (Spinning Enhanced Visible and Infrared Imager). To achieve this aim, we use a new Convective/Stratiform Rain Area Delineation Technique (CS-RADT). The satellite rainfall retrieval technique is based on various spectral parameters of SEVIRI that express microphysical and optical cloud properties. It uses a multispectral thresholding technique to distinguish between stratiform and convective clouds. This technique (CS-RADT) is applied to the complex situation of the Mediterranean climate of this region. The tests have been conducted during the rainy seasons of 2006/2007 and 2010/2011 where stratiform and convective precipitation is recorded. The developed scheme (CS-RADT) is calibrated by instantaneous meteorological radar data to determine thresholds, and then rain rates are assigned to each cloud type by using radar and rain gauge data. These calibration data are collocated with SEVIRI data in time and space.     

The gravitropic and phototropic responses of wheat grown in a space greenhouse prototype with hemispherical planting surface

The time course of gravicurvature of 3-day-old wheat (Triticum aestivum L., cv. Apogee) coleoptiles and 7-day-old wheat stems were studied in darkness and under red and red-blue light illumination after declination from the vertical at various angles. The experiments showed that the shortest gravitropic curvature corresponded to 30° initial angle of gravistimulation (IAG). The time course became longer as the IAG increased and with plant age. The effects of unilateral red (660 nm) and red-blue light (660 nm; 470 nm) at photosynthetic photon flux (PPF) of 30 μmol m⁻² s⁻¹ on the curvature of 3-day-old coleoptiles were evaluated. Red light did not produce phototropic bending of wheat coleoptiles in contrast with red-blue light. The analysis of experimental data showed that the curvature in response to a gravitropic stimulus or to combined gravity-light stimuli were not statistically different. Time course of gravitropic curvature were used to determine the acceptable crop rotation rate around the horizontal axis. Approximation of stem bending to a linear dynamic system described by a first-order aperiodic element with a lag allowed the determination of the dependence of the amplitude of apex oscillations on the rate of horizontal rotation under 1-g conditions. The calculated lowest minimal rotation rate (MRR) minimizing the gravitropic effects on wheat was about 1 revolution per hour (rph). Rotating the plant growth chamber (PGC) at a rate of more than MRR eliminated the effect of gravitropic curvature.     

Study of mid-latitude ionospheric convection during quiet and disturbed periods using the SuperDARN Hokkaido radar

Westward ionospheric convective flows around midnight are frequently observed at mid-latitudes. They can be generated by so-called disturbance dynamo mechanisms working mainly in the mid-latitudes. To understand the influence of disturbance dynamo effects in the mid-latitudes, we studied the latitudinal distribution of westward flows in association with several kinds of geomagnetic disturbances using the SuperDARN Hokkaido radar. This radar creates high temporal resolution (1 s to 2 min), two-dimensional observations measuring the line-of-sight velocities of ionospheric plasma irregularities, which can be regarded as line-of-sight velocities of ionospheric convection in the mid-latitude region from 40° to 50°. This region could not be monitored using preexisting SuperDARN radars. In this study, we used ionospheric echo data obtained by the SuperDARN Hokkaido radar over 5 years (December 2006 to November 2011). We identified westward flows around midnight at about 40° to 55° geomagnetic latitude. Additionally, the data showed that the westward flow around midnight intensified under high geomagnetic activity (high Kp). This suggests that the disturbance dynamo could affect the mid-latitude ionospheric convection. We performed Superposed Epoch Analysis (SEA) to study the influences from the geomagnetic disturbances on mid-latitude ionospheric convection. We found no obvious influence during major storms (minimum Dst below −60 nT). SEA was also used to study the temporal and latitudinal dependence on the influences from substorms. From analysis of 36 events of AL-defined substorms, we saw that the influence of substorms lasted from 5 to 20 h after the onset between 44° and 53° geomagnetic latitude. The westward flow at mid-latitude grew to a maximum at 12 h after the geomagnetic substorm onset. This is consistent with the results of past numerical simulation studies of the disturbance dynamo effects.     

Comparison of mesopause region meteor radar winds,medium frequency radar winds and low frequency drifts over Germany

During 2004 and 2005 measurements of mesospheric/lower thermospheric (80–100 km) winds have been carried out in Germany using three different ground-based systems, namely a meteor radar (36.2 MHz) at the Collm Observatory (51.3°N, 13°E), a MF radar (3.18 MHz) at Juliusruh (54.6°N, 13.4°E) and the LF D1 measurements using a transmitter (177 kHz) at Zehlendorf near Berlin and receivers at Collm with the reflection point at 52.1°N, 13.2°E. This provides the possibility of comparing the results of different radar systems in nearly the same measuring volume. Meteor radar winds are generally stronger than the winds observed by MF and especially by LF radars. This difference is small near 80 km but increases with height. The difference between meteor radar and medium frequency radar winds is larger during winter than during summer, which might indicate an indirect influence of gravity waves on spaced antenna measurements.     

Gravity wave propagation studies using the Indian MST radar observations

The Indian MST radar facility at Gadanki (13.5°N, 79.2°E) has been utilised to study the propagation of gravity waves from the troposphere/lower stratosphere to the mesosphere and their interaction with the radar backscattered signal variations. The main objective is to correlate vertically propagating gravity waves derived from the tropospheric velocity fields with the dynamics of mesospheric scattering centres. The tropospheric wind velocities and signal strengths over the entire height range have been subjected to power spectral and wavelet analysis to determine the predominant wave periods/amplitudes and the coupling between the lower atmosphere and mesosphere. Results show that (a) the gravity waves are clearly detectable near tropopause heights, (b) while relatively higher period gravity waves (20–50 min) interact with mesospheric scattering centres, the lower period waves (<20 min) are absorbed in the troposphere itself, (c) the mesospheric scattering layers are affected by gravity waves of complementary periods.     

E region electric fields at the dip equator and anomalous conductivity effects

Zonal and vertical electric fields were estimated at E region heights in the Brazilian sector. Zonal electric fields are obtained from the vertical electric fields based on their relation through the Hall-to-Pedersen ionospheric conductivities ratio. The technique for obtaining the vertical electric field is based on its proportionality to the Doppler velocities of type 2 irregularities as detected by coherent radars. The 50 MHz backscatter coherent (RESCO) radar was used to estimate the Doppler velocities of the type 2 irregularities embedded in the equatorial electrojet. A magnetic field-line integrated conductivity model was developed to provide the conductivities. It considers a multi-species ionosphere and a multi-species neutral atmosphere, and uses the IRI 2007, the MISIS 2000 and the IGRF 10 models as input parameters for ionosphere, neutral atmosphere and Earth’s magnetic field, respectively. The ion-neutral collision frequencies of all the species are combined through the momentum transfer collision frequency equation, and different percentages of electron-neutral collisions were artificially included for studying the implication of such increase in the zonal electric field, which resulted ranging from 0.13 to 0.49 mV/m between the 8 and 18 h (LT), under quiet magnetic conditions.     

Response of low latitude D-region ionosphere to the total solar eclipse of 22 July 2009 deduced from ELF/VLF analysis

Response of the D-region of the ionosphere to the total solar eclipse of 22 July 2009 at low latitude, Varanasi (Geog. lat., 25.27° N; Geog. long., 82.98° E; Geomag. lat. = 14° 55’ N) was investigated using ELF/VLF radio signal. Tweeks, a naturally occurring VLF signal and radio signals from various VLF navigational transmitters are first time used simultaneously to study the effect of total solar eclipse (TSE). Tweeks occurrence is a nighttime phenomena but the obscuration of solar disc during TSE in early morning leads to tweek occurrence. The changes in D-region ionospheric VLF reflection heights (h) and electron density (n_e: 22.6–24.6 cm⁻³) during eclipse have been estimated from tweek analysis. The reflection height increased from ∼89 km from the first occurrence of tweek to about ∼93 km at the totality and then decreased to ∼88 km at the end of the eclipse, suggesting significant increase in tweek reflection height of about 5.5 km during the eclipse. The reflection heights at the time of totality during TSE are found to be less by 2–3 km as compared to the usual nighttime tweek reflection heights. This is due to partial nighttime condition created by TSE. A significant increase of 3 dB in the strength of the amplitude of VLF signal of 22.2 kHz transmitted from JJI-Japan is observed around the time of the total solar eclipse (TSE) as compared to a normal day. The modeled electron density height profile of the lower ionosphere depicts linear variation in the electron density with respect to solar radiation as observed by tweek analysis also. These low latitude ionospheric perturbations on the eclipse day are discussed and compared with other normal days.     

Lunar ground penetrating radar: Minimizing potential data artifacts caused by signal interaction with a rover body

Ground-penetrating radar (GPR) is the leading geophysical candidate technology for future lunar missions aimed at mapping shallow stratigraphy (<5 m). The instrument’s exploration depth and resolution capabilities in lunar materials, as well as its small size and lightweight components, make it a very attractive option from both a scientific and engineering perspective. However, the interaction between a GPR signal and the rover body is poorly understood and must be investigated prior to a space mission. In doing so, engineering and survey design strategies should be developed to enhance GPR performance in the context of the scientific question being asked. This paper explores the effects of a rover (simulated with a vertical metal plate) on GPR results for a range of heights above the surface and antenna configurations at two sites: (i) a standard GPR testing site with targets of known position, size, and material properties, and; (ii) a frozen lake for surface reflectivity experiments. Our results demonstrate that the GPR antenna configuration is a key variable dictating instrument design, with the XX polarization considered optimal for minimizing data artifact generation. These findings could thus be used to help guide design requirements for an eventual flight instrument.     

Vertical wavenumber spectra of atmospheric ozone measured from ozonesonde observations

Vertical profiles of ozone have been measured at balloon altitudes. Our purpose is to examine the character of vertical wavenumber spectra of ozone fluctuations, to assess the possible roles of gravity wave field in ozone fluctuations, and to determine dominant vertical wavelengths of ozone spectra. Vertical wavenumber spectra of 12 ozone fluctuations obtained during June–August 2003 are presented. Results indicate that mean spectral slopes in the wavenumber range from 4.69 × 10⁻⁴ to 2.50 × 10⁻³ cyc/m are about −2.91 in the troposphere and −2.87 in the lower stratosphere, which is close to the slope of −3 predicted by current gravity wave saturation models. The consistency of the observed spectral slopes with the value of −3 predicted by current gravity wave saturation models suggests that the observed ozone fluctuations are due primarily to atmospheric gravity waves. At m = 1/(1000 m) the mean spectral amplitude is over 30 times larger in the lower stratosphere than in the troposphere. Mean vertical wavenumber spectra in area-preserving form reveal dominant vertical wavelengths of ∼2.6 km in the troposphere and ∼2.7 km in the lower stratosphere, which is consistent with the values varying between 1.5 and 3.0 km estimated from the velocity field and temperature field at these heights.     

An Artificial Neural Network based approach for estimation of rain intensity from spectral moments of a Doppler Weather Radar

By using a Doppler Weather Radar (DWR) at Shriharikota (13.66°N & 80.23°E), an Artificial Neural Network (ANN) based technique is proposed to improve the accuracy of rain intensity estimation. Three spectral moments of a Doppler spectra are utilized as an input data to an ANN. Rain intensity, as measured by the tipping bucket rain gauges around the DWR station, are considered as a target values for the given inputs. Rain intensity as estimated by the developed ANN model is validated by the rain gauges measurements. With the help of a developed technique, reasonable improvement in the estimation of rain intensity is observed. By using the developed technique, root mean square error and bias are reduced in the range of 34–18% and 17–3% respectively, compared to Z–R approach.     

Titan’s new pole: Implications for the Huygens entry and descent trajectory and landing coordinates

The European Space Agency’s Huygens probe separated from the NASA Cassini spacecraft on 25 December 2004, after having been attached for a 7-year interplanetary journey and three orbits around Saturn. The probe reached the predefined NASA/ESA interface point on 14 January 2005 at 09:05:52.523 (UTC). It performed a successful entry and descent sequence and softly landed on Titan’s surface on the same day at 11:38:10.77 (UTC) with a speed of about 4.54 m/s. Since the publication of the official project entry and descent trajectory reconstruction effort by the Descent Trajectory Working Group in 2007 (referred to as DTWG#4) various other efforts have been performed and published. This paper presents an overview of the most relevant reconstructions and compares their methodologies and results. Furthermore, the results of a new reconstruction effort (DTWG#5) are presented, which is based on the same methodology as DTWG#4 but takes into account new estimates of Titan’s pole coordinates which were derived from radar images of different Cassini Titan flybys. It can be shown that the primary effect can be observed in the meridional direction which is represented by a stark southward shift of the trajectory by about 0.3 deg. A much smaller effect is seen in the zonal direction (i.e., less than 0.01 deg in the west to east direction). The revised probe landing coordinates are 192.335 deg W and 10.573 deg S. A comparison of these coordinates with results of recent landing site investigations using visual and radar images of the Cassini VIMS instrument shows excellent agreement of the two independently derived landing coordinates, i.e., longitude and latitude residuals of respectively 0.035 deg and 0.007 deg.     

The characteristics of the semi-diurnal tides in mesosphere/low-thermosphere (MLT) during 2002 at Wuhan (30.6°N, 114.4°E) – using canonical correlation analysis technique

In this paper, we use canonical correlation analysis (CCA) method to investigate the semi-diurnal tidal winds in mesosphere and low thermosphere (MLT) region, observed by a newly installed meteor radar at Wuhan (30.6°N, 114.4°E), during the year 2002. In general, 4(3) effective semi-diurnal tidal pairs of patterns are obtained, which represent ∼2/3 total variances of the origin data set. These patterns are expected to be corresponding to the atmospheric oscillations within the semi-diurnal frequency band excited or modulated by different sources, i.e., the seasonal variations, the modulations by the planetary wave oscillations or the solar 27-day activity. Among all the patterns, the 1st pattern, which represents ∼1/3 of total variances, is the most notable. Its amplitudes show maximum values in spring and autumn, and the vertical wavelengths are longer in summer and shorter in winter, which is in line with the results obtained from traditional harmonic analysis. The vertical wavelengths of the higher order patterns (∼50 km) suggest the classic semi-tidal mode S(2, 4)/S(2, 5) is dominant.