Numerical modeling of inhomogeneous orographic wave influence on planetary waves in the middle atmosphere

Parameterization of dynamical and thermal effects of stationary orographic gravity waves (OGWs) generated by the Earth’s surface topography is incorporated into a numerical model of general circulation of the middle and upper atmosphere. Responses of atmospheric general circulation and characteristics of planetary waves at altitudes from the troposphere up to the thermosphere to the effects of OGWs propagating from the earth surface are studied. Changes in atmospheric circulation and amplitudes of planetary waves due to variations of OGW generation and propagation in different seasons are considered. It is shown that during solstices the main OGW dynamical and heat effects occur in the middle atmosphere of winter hemispheres, where changes in planetary wave amplitudes due to OGWs may reach up to 50%. During equinoxes OGW effects are distributed more homogeneously between northern and southern hemispheres.     

Atmosphere-ionosphere coupling from convectively generated gravity waves

Ionospheric variability impacts operational performances of a variety of technological systems, such as HF communication, Global Positioning System (GPS) navigation, and radar surveillance. The ionosphere is not only perturbed by geomagnetic inputs but is also influenced by atmospheric tides and other wave disturbances propagating from the troposphere to high altitudes. Atmospheric Gravity Waves (AGWs) excited by meteorological sources are one of the largest sources of mesoscale variability in the ionosphere. In this paper, Total Electron Content (TEC) data from networks of GPS receivers in the United States are analyzed to investigate AGWs in the ionosphere generated by convective thunderstorms. Two case studies of convectively generated gravity waves are presented. On April 4, 2014 two distinct large convective systems in Texas and Arkansas generated two sets of concentric AGWs that were observed in the ionosphere as Traveling Ionospheric Disturbances (TIDs). The period of the observed TIDs was 20.8 min, the horizontal wavelength was 182.4 km, and the horizontal phase speed was 146.4 m/s. The second case study shows TIDs generated from an extended squall line on December 23, 2015 stretching from the Gulf of Mexico to the Great Lakes in North America. Unlike the concentric wave features seen in the first case study, the extended squall line generated TIDs, which exhibited almost plane-parallel phase fronts. The TID period was 20.1 min, its horizontal wavelength was 209.6 km, and the horizontal phase speed was 180.1 m/s. The AGWs generated by both of these meteorological events have large vertical wavelength (>100 km), which are larger than the F2 layer thickness, thus allowing them to be discernible in the TEC dataset.     

The 16-day waves in the mesosphere and lower thermosphere over Wuhan (30.6°N, 114.5°E) and Adelaide (35°S, 138°E)

Winds from a meteor radar at Wuhan (30.6°N, 114.5°E) and a MF radar at Adelaide (35°S, 138°E) are used to study the 16-day waves in the mesosphere and lower thermosphere (MLT). The height range is 78–98 km at Wuhan and 70–98 km at Adelaide. By comparison, it is found that the zonal components at both sites are generally larger than the meridional ones, and eastward motion of the zonal background winds is favorable for the 16-day waves penetration to the MLT region. The zonal maximum amplitude appears in the autumn (September–October) around 86–98 km at Wuhan and in the winter months and early spring (July–October) around 72–82 km at Adelaide. Differences are found in wave amplitudes and time of appearance between the two years of 2002 and 2003. In 2003, the intensity of the wave amplitudes is relatively smaller than that for 2002 at both sites. The summer 16-day waves are comparatively weaker at Adelaide in both years, but stronger in 2002 at Wuhan near the mesopause and the lower thermosphere (86–98 km). The strong summer waves at Wuhan may come from the winter southern hemisphere.     

Development of a 2.8 μm film for scientific balloons

Development of a balloon to fly at higher altitudes is one of the most attractive challenges in scientific balloon technologies. After reaching the highest record setting balloon altitude of 53.0 km using the 3.4 μm film in 2002, we tried to make a thinner balloon film. In 2003, we developed a forming die and an air-ring and succeeded in forming a film with a thickness of 3.0 μm and a width of 220 cm. Using this film, we manufactured a balloon with a volume of 5000 m³ and succeeded in flying the balloon up to an altitude of 46.0 km. We then searched for a good combination of resins to make a thinner and wider film and obtained films with widths of 280 cm, and a thickness of 3.0 μm at first, and then 2.8 μm. In 2004, we performed balloon experiments making a 30,000 m³ balloon with the 3.0 μm film and a 5000 m³ balloon with the 2.8 μm film. Both balloons were well manufactured and reached the highest altitudes of 50.7 and 42.6 km, respectively.     

Seasonally steady planetary disturbances in the troposphere and stratosphere as seen in 30 years of NCEP reanalysis data

Zonal velocity and temperature daily global reanalysis data of 30 years are used to search seasonally steady planetary disturbances in the middle troposphere (400 hPa) and middle stratosphere (10 hPa). Significant wavenumber 1, 2 and 3 modes are found. Constant phase lines of zonal velocity 1 modes exhibit significant inclination angles with respect to the meridians. The winter hemisphere generally shows a more significant presence of structures. The Northern Hemisphere (NH) exhibits all over the year a larger amount of structures and more intense amplitudes than the Southern Hemisphere (SH). Middle latitudes exhibit the most significant cases and low latitudes the least significant ones. Longitudinally oriented land–sea transitions at ±$\pm$65° and −35° latitudes appear to play a significant role for the presence of steady planetary modes. The stratosphere exhibits a much simpler picture than the troposphere. Large scale structures with respectively NE–SW (NH) and NW–SE (SH) tilts in the observed temperature and zonal velocity constant phase lines recall the quasi-stationary Rossby wave trains that favor the poleward transport of angular momentum.     

Estimates of mesospheric gravity wave activity over convection from a global model

A new convective gravity wave source spectrum parameterization has been implemented in the Whole Atmosphere Community Climate Model version 2 (WACCM2). This parameterization specifies the momentum flux phase speed spectrum of gravity waves in the Tropics based on the properties of underlying convection; Hence, this parameterization provides realistic global estimates of gravity wave activity. In this paper, we show the estimated gravity wave phase speed spectra in the Tropics from a WACCM2 simulation, at the source level and at 85 km. Spatial distribution of gravity wave activity at 85 km is also presented. Subsequently, we discuss the factors that are primarily responsible for the estimated differences in gravity wave distribution across phase speeds with latitude and asymmetries in direction of gravity wave propagation in the mesosphere. We also examine which of the model assumptions can lead to uncertainties in our estimates of mesospheric gravity wave activity and we discuss how these assumptions provide challenges for comparison with observations of gravity waves in the mesosphere.     

Detection of long-term trends in the mesosphere/lower thermosphere from ground-based radio propagation measurements

Mesospheric temperature trends can be derived from LF phase-height observations in mid-latitudes supported by ionospheric absorption and ionosonde observations. Analysing the full observation period from 1959 until 2003, a mean yearly temperature trend has been derived with −0.25 K/yr for the height interval from 48 to 82 km. Subdividing the whole observation interval in two parts before and after 1979, the trend is markedly stronger in the second period with −0.38 K/yr compared with −0.20 K/yr in the first part before 1979. These differences can at least partly be explained by a steeper CO₂ increase and ozone decrease in the second interval. The differences in the mesospheric temperature trends are most evidently expressed during winter months and are markedly smaller during summer season. The reason of this seasonal difference is not quite clear; it may be related with detected ozone trends which are clearly stronger during winter months on both hemispheres.     

Multi-year behaviour of the midnight OH1 temperature according to observations at Zvenigorod over 2000–2016

Using spectral measurements of the hydroxyl airglow at the Zvenigorod station (56° N, 37° E), Moscow region, over 2000–2016, we obtained the long-term set of data comprising 1822 midnight values of the OH¹ temperature in the mesopause region. These data revealed a 17-year series of its mean annual values, as well as amplitudes and phases of the first two harmonics of its annual variation. The obtained parameters were analyzed to determine statistically relevant characteristics of their long-term variations. As a result, we found that the long-term behaviour of the mean annual OH¹ temperature features a small negative linear trend (?0.07 ± 0.03 K/year) over the addressed period. Besides, its dependence on solar activity is shown to be 4.1 ± 0.5 K/100 SFU. Regarding the long-term behaviour of the mean annual OH¹ temperature, we revealed the existence of two oscillations with 3-year (the amplitude being 1.3 ± 0.2 K) and 4.1-year (the amplitude being 0.6 ± 0.2 K) periods. We obtained empirical relations describing year-to-year variations in the amplitudes and phases of the annual and semi-annual harmonics.     

Analysis of electromagnetic and electrostatic effects of particle impacts on spacecraft

Particle impacts on spacecraft can cause considerable damage, even leading to complete failure. A theory for the resulting vehicle potential changes and the electromagnetic radiation from impact-induced plasma has been published by Close et al. (2010). Here we compare this theory to impacts registered by the Radio and Plasma Wave Science instrumentation on the Cassini spacecraft. We study both low-velocity (16 km/s) large particles (2.6 μm radius) detected in Saturn’s rings and high-velocity (450 km/s) small particles (1 nm radius) in the solar wind. The agreement with the theory is quite good. We also apply these results to earth orbit and conclude that both Electrostatic Discharge and Electromagnetic Pulse radiation are likely and could lead to spacecraft failure.     

On the 27-day variations of the galactic cosmic ray anisotropy and intensity for different periods of solar magnetic cycle

The average amplitude of the 27-day variation of the galactic cosmic ray anisotropy calculated based on the neutron monitors experimental data is larger in the qA > 0 period than in the qA < 0 period of solar magnetic cycle. The amplitudes of the 27-day variation of the galactic cosmic rays anisotropy do not depend on the tilt angles of the heliospheric neutral sheet for different the qA > 0 and the qA < 0 periods of solar magnetic cycle. A good correlation has been revealed between the changes of the amplitudes of the 27-day variations of the galactic cosmic ray anisotropy and intensity versus the qA > 0 and the qA < 0 periods of solar magnetic cycle.     

USV test flight by stratospheric balloon: Preliminary mission analysis

The Unmanned Space Vehicle test flights will use a 7 m 1300 kg aircraft. The first three launches will take place at the Italian Space Agency ASI base in Trapani–Milo, Sicily, through a stratospheric balloon that will drop the aircraft at a predefined height. After free fall acceleration to transonic velocities, the parachute deployment will allow a safe splash down in the central Mediterranean Sea. The goal of this article is to show the preliminary analysis results for the first USV flight.We carried out a statistical study for the year 2000–2003, evaluating the typical summer and winter launch windows of the Trapani–Milo base.First, in the center Mediterranean, we define safe recovery areas. They cannot be reached during the balloon ascending phase so, after a sufficiently long floating part able to catch the open sea, the balloon will go down to the release height (24 km). The simulation foresees a 400,000 m³ balloon and 3 valves for the altitude transfer.A safe splash down must occur far enough from the nearest coast: the minimum distance is considered around 25 km. The vehicle should be released at a distance, from the nearest coast, greater than this minimum amount plus the USV model maximum horizontal translation, during its own trajectory from balloon separation to splash down. In this way we define safe release areas for some possible translations.Winter stratospheric winds are less stable. The winter average flight duration is 7 h and it is probably too long for the diurnal recovery requirement and its scheduled procedures.Comparing past stratospheric balloons flights and trajectories computed using measured meteorological data (analysis), with their predictions made using forecast models and soundings, we obtain the standard deviation of the trajectory forecast uncertainty at the balloon–aircraft separation. Two cases are taken into account: predictions made 24 and 6 h before the launch.Assuming a Gaussian latitudinal uncertainty distribution for the prediction 6 h before the launch, we are able to identify the forecast trajectories that have a probability greater than 97% to reach the safe release areas.Simulating the summer windows trajectories for the years from 2000 to 2003 and for the favorable ground wind days, we obtain the number of trajectories with the desired forecast probabilities.     

Characterization of the semi-annual-oscillation in mesospheric temperatures at low-latitudes

Novel measurements of the seasonal variability in mesospheric temperature at low-latitudes have been obtained from Maui, Hawaii (20.8°N, 156.2°W) during a 25-month period from October 2001 to January 2004. Independent observations of the OH (6, 2) Meinel band (peak height ∼87 km) and the O₂ (0–1) atmospheric band emission (∼94 km) were made using the CEDAR Mesospheric Temperature Mapper. The data revealed a coherent oscillation in emission intensity and rotational temperature with a well-defined periodicity of 181 ± 7 days. The amplitude of this oscillation was determined to be ∼5–6 K in temperature and ∼8–9% in intensity for both the OH and O₂ data sets. In addition, a strong asymmetry in the shape of the oscillation was also observed with the spring maximum significantly larger than the fall peak. These data provide new evidence in support of a semi-annual-oscillation in mesospheric temperature (and airglow emission intensities) and help quantify its seasonal characteristics.     

Study of changes in the lineament structure,caused by earthquakes in South America by applying the lineament analysis to the Aster (Terra) satellite data

The region between Southern Peru and Northern Chile is one of the most seismically and volcanically active regions in South America. This is caused by a constant subduction of the South American Plate, converging with the Nazca Plate in the extreme North of Chile. We used the 15 and 30 m resolution satellite images, provided by the ASTER (VNIR and SWIR) instrument onboard the Terra satellite to study changes in the geological faults close to earthquake epicenters in southern Peru. Visible and infrared spectral bands were analysed using “The Lineament Extraction and Stripes Statistic Analysis” (LESSA) software package to examine changes in the lineament features and stripe density fields caused by seismic activity. We used the satellite images 128 and 48 days before and 73 days after a 5.2 Richter scale magnitude earthquake. The fact that the seasonal variations in the South of Peru and North of Chile are very small, and the vegetation is very limited, allowed us to establish substantial changes in the lineament and the stripe density field features. We develop a methodology that allows to evaluate the seismic risk in this region for the future.     

The lower ionosphere response to its disturbances by powerful radio waves

The paper presents data from some campaigns at Sura heating facility in 2011–1016. The experiments on probing of the artificial disturbed region of the lower ionosphere were carried out at two observation sites. One of them was located near Vasil’sursk 1 km from Sura facility (56.1°N; 46.1°E) and the other site was located at the Observatory (55.85°N; 48.8°E) of Kazan State University, 170 km to the East. Investigation of the features of the disturbed region of the lower ionosphere based on its diagnostics by the methods of the vertical sounding and oblique backscattering is the main goal of this paper. Ionosphere disturbance was fulfilled by the effect of the powerful radio wave of the ordinary or extraordinary polarization emitted by transmitters of the Sura facility with effective radiated power ERP = 50–120 MW at the frequency of 4.3, 4.7 and 5.6 MHz. Pumping waves were emitted with period from 30 s to 15 min. The disturbed region of the ionosphere in Vasil’sursk was probed by the vertical sounding technique using the partial reflexion radar at the frequency of 2.95 and 4.7 MHz. For the oblique sounding of the disturbed region the modified ionosonde Cyclon-M, operating at ten frequencies from 2.01 to 6.51 MHz was used at the Observatory site. On many heating sessions simultaneous variations of the probing partial reflection signals in Vasil’sursk and backscattered signals in Observatory were observed at the height at 40–100 km below the reflection height of the pumping wave. These observations were correlated with the pumping periods of the Sura facility. Possible mechanisms of the appearance of the disturbance in the lower ionosphere and its effect on the probing radio waves are discussed.     

Global MGS TES data and Mars-GRAM validation

Mars Global Reference Atmospheric Model (Mars-GRAM 2001) is an engineering-level Mars atmosphere model widely used for many Mars mission applications. From 0 to 80 km, it is based on NASA Ames Mars General Circulation Model (MGCM), while above 80 km it is based on University of Michigan Mars Thermospheric General Circulation Model. Mars-GRAM 2001 and MGCM use surface topography from Mars Global Surveyor Mars Orbiter Laser Altimeter (MOLA). Validation studies are described comparing Mars-GRAM with a global summary data set of Mars Global Surveyor Thermal Emission Spectrometer (TES) data. TES averages and standard deviations were assembled from binned TES data which covered surface to ∼40 km, over more than a full Mars year (February 1999–June 2001, just before start of a Mars global dust storm). TES data were binned in 10° × 10° latitude–longitude bins (36 longitude bins, centered at 5°–355°, by 18 latitude bins, centered at −85° to +85°), and 12 seasonal bins (based on 30° increments of L_s angle). Bin averages and standard deviations were assembled at 23 data levels (temperature at 21 pressure levels, plus surface temperature and surface pressure). Two time-of-day bins were used: local time near 2 or 14 h. Two dust optical depth bins were used: infrared optical depth, either less than or greater than 0.25 (which corresponds to visible optical depth less than or greater than about 0.5). For interests in aerocapture and precision entry and landing, comparisons focused on atmospheric density. TES densities versus height were computed from TES temperature versus pressure, using assumptions of perfect gas law and hydrostatics. Mars-GRAM validation studies used density ratio (TES/Mars-GRAM) evaluated at data bin center points in space and time. Observed average TES/Mars-GRAM density ratios were generally 1 ± 0.05, except at high altitudes (15–30 km, depending on season) and high latitudes (>45°N), or at most altitudes in the southern hemisphere at L_s ∼ 90° and 180°. Compared to TES averages for a given latitude and season, TES data had average density standard deviation about the mean of ∼2.5% for all data, or ∼1–4%, depending on time of day and dust optical depth. Average standard deviation of TES/Mars-GRAM density ratio was 8.9% for local time 2 h and 7.1% for local time 14 h. Thus standard deviation of observed TES/Mars-GRAM density ratio, evaluated at matching positions and times, is about three times the standard deviation of TES data about the TES mean value at a given position and season.     

On the production of highest energy solar protons at 20 January 2005

On January 20, 2005, 7:02–7:05 UT the Aragats Multidirectional Muon Monitor (AMMM) located at 3200 m a.s.l. registered enhancement of the high energy secondary muon flux (threshold ∼5 GeV). The enhancement, lasting for 3 min, has statistical significance of ∼4σ and is related to the X7.1 flare seen by the GOES, and very fast (>2500 km/s) CME seen by SOHO, and the Ground Level Enhancements (GLE) #69 detected by the world-wide network of neutron monitors and muon detectors. The energetic and temporal characteristics of the muon signal from the AMMM are compared with the characteristics of other monitors located at the Aragats Space-Environmental Center (ASEC) and with other neutron and muon detectors. Since secondary muons with energies >5 GeV are corresponding to solar proton primaries with energies 20–30 GeV we conclude that in the episode of the particle acceleration at 7:02–7:05 UT 20 January 2005 solar protons were accelerated up to energies in excess of 20 GeV.     

Multi precursors analysis associated with the powerful Ecuador (MW = 7.8) earthquake of 16 April 2016 using Swarm satellites data in conjunction with other multi-platform satellite and ground data

After DEMETER satellite mission (2004–2010), the launch of the Swarm satellites (Alpha (A), Bravo (B) and Charlie (C)) has created a new opportunity in the study of earthquake ionospheric precursors. Nowadays, there is no doubt that multi precursors analysis is a necessary phase to better understand the LAIC (Lithosphere Atmosphere Ionosphere Coupling) mechanism before large earthquakes. In this study, using absolute scalar magnetometer, vector field magnetometer and electric field instrument on board Swarm satellites, GPS (Global Positioning System) measurements, MODIS-Aqua satellite and ECMWF (European Centre for Medium-Range Weather Forecasts) data, the variations of the electron density and temperature, magnetic field, TEC (Total Electron Content), LST (Land Surface Temperature), AOD (Aerosol Optical Depth) and SKT (SKin Temperature) have been surveyed to find the potential seismic anomalies around the strong Ecuador (M_w = 7.8) earthquake of 16 April 2016. The four solar and geomagnetic indices: F10.7, D_st, K_p and a_p were investigated to distinguish whether the preliminary detected anomalies might be associated with the solar-geomagnetic activities instead of the seismo-ionospheric anomalies. The Swarm satellites (A, B and C) data analysis indicate the anomalies in time series of electron density variations on 7, 11 and 12 days before the event; the unusual variations in time series of electron temperature on 8 days preceding the earthquake; the analysis of the magnetic field scalar and vectors data show the considerable anomalies 52, 48, 23, 16, 11, 9 and 7 days before the main shock. A striking anomaly is detected in TEC variations on 1 day before earthquake at 9:00 UTC. The analysis of MODIS-Aqua night-time images shows that LST increase unusually on 11 days prior to main shock. In addition, the AOD variations obtained from MODIS measurements reach the maximum value on 10 days before the earthquake. The SKT around epicentral region presents anomalous higher value about 40 days before the earthquake. It should be noted that the different lead times of the observed anomalies could be acknowledged based on a reasonable LAIC earthquake mechanism. Our results emphasize that the Swarm satellites measurements play an undeniable role in progress the studies of the ionospheric precursors.     

OPTIS – A Satellite test of Special and General Relativity

OPTIS has been proposed as a small satellite platform in a high elliptical orbit (apogee 40,000 km, perigee 10,000 km) and is designed for high precision tests of foundations of Special and General Relativity. The experimental set-up consists of two ultrastable Nd:YAG lasers, three crossed optical resonators (monolithic cavities), an ensemble of atomic clocks, an optical comb generator, laser tracking devices and a drag-free control system. OPTIS enables improved tests of (1) the isotropy and (2) constancy of the speed of light, (3) special relativistic time dilation, (4) the universality of the gravitational redshift by comparison of various clocks, can measure (5) the absolute value of the gravitational redshift, (6) the Lense–Thirring effect and (7) the perigee advance and (8) can make a test of a hypothetical Yukawa part in the gravitational potential. To avoid any influence from atmospheric drag, solar radiation, or Earth albedo, the satellite needs drag-free control to depress the residual acceleration down to 10⁻¹⁴ m/s² in the frequency range between 10⁻² and 10⁻³ Hz. Precise thermal control must be used to stabilize the cavity temperature to within one part in 10⁷ at time scales of 100 s and to one part in 10⁵ on the orbit time scale.     

Study of >100 GeV electrons with BETS detector using a long duration balloon flight in Antarctica

We have observed cosmic-ray electrons from 10 to 1000 GeV by a long duration balloon flight using Polar Patrol Balloon (PPB) in Antarctica. The observation was carried out for 13 days at an altitude of 35 km in January 2004. The detector is an imaging calorimeter composed of scintillating-fiber belts and plastic scintillators inserted between lead plates. The geometrical factor of detector is about 600 cm²sr and the total thickness of lead absorber is 9 radiation lengths. The performance of the detector has been confirmed by the CERN-SPS beam test and also investigated by Monte-Carlo simulations. New telemetry system using a commercial satellite of iridium, power supply by solar batteries, and automatic level control using CPU have successfully been developed and operated during the flight. We have collected 5.7 × 10³ events over 100 GeV including nearly 100 candidates of primary electrons.     

On the quality of mutations in mammalian cells induced by high LET radiations

The deleterious effects of accelerated heavy ions as component of the space radiation environment on living cells are of increasing importance for long duration human space flight activities. The most important aspect of such densely ionizing particle radiation is attributed to the type and quality of biological damage induced by them. This issue is addressed by investigating cell inactivation and mutation induction at the Hprt locus (coding for hypoxanthine-guanine-phosphoribosyl-transferase) of cultured V79 Chinese hamster cells exposed to densely ionizing radiation (accelerated heavy ions with different LETs from oxygen to gold, specific energies ranging from 1.9 to 69.7 MeV/u, corresponding LET values range from 62 to 13,223 keV/μm) and to sparsely ionizing radiation (200 kV X-rays). 30 spontaneous, 40 X-ray induced and 196 heavy ion induced 6-thioguanine resistant Hprt mutant colonies were characterized by Southern technique using the restriction enzymes EcoRI, PstI and BglII and a full length Hprt cDNA probe isolated from the plasmid pHPT12. Restriction patterns of the spontaneous Hprt mutants were indistinguishable from the wild type pattern, as these mutants probably contain only small deletions or even point mutations in the Hprt locus. In contrast, the overall spectrum of heavy ion induced mutations revealed a majority of partial or total deletions of the Hprt gene. With constant particle fluence (3 × 10⁶ particles/cm²) the quality of heavy ion induced mutations in the Hprt locus depends on physical parameters of the beam (atomic number, specific energy, LET). This finding suggests a relationship between the type of DNA damage and track structure. The fraction of mutants with severe deletions in the Hprt locus after exposure to oxygen ions increases from 65% at 60 keV/μm up to a maximum (100%) at 300 keV/μm and declines with higher LET values to 75% at 750 keV/μm. With heavier ions (Ca- and Au-ions) and even higher LET-values this mutant fraction decreases to 58% at 13,200 keV/μm. Heavy ion induced DNA break points in the Hprt locus are not randomly distributed.