Showing papers in "Radio Science in 2012"

Relative detection efficiency of the World Wide Lightning Location Network

Abstract: [1] Using the detected energy per strokes of the World Wide Lightning Location Network (WWLLN) we calculate the relative detection efficiency for the network as if it had a uniform detection efficiency. The model uses the energy statistics of located strokes to determine which stations are sensitive to what stroke energies. We are then able to estimate the number of strokes that may be missing from any given regions as compared to the best, most sensitive regions of the WWLLN network. Stroke density maps can be corrected with the knowledge of how sensitive various regions of the network are operating. This new model for the relative WWLLN detection efficiency compensates for the uneven global coverage of the network sensors as well as variations in very low frequency (VLF) propagation. The model gives a way to represent the global distribution of strokes as if observed by a globally uniform network. The model results are analyzed in spatial and temporal regimes, and the effects of a single VLF detector going offline are investigated in areas of sparse and dense detector coverage. The results are also used to show spatial, temporal and energy distributions as seen by the detection efficiency corrected WWLLN.

Ionospheric Signatures of Tohoku-Oki Tsunami of March 11, 2011: Model Comparisons Near the Epicenter

Abstract: [1] We observe ionospheric perturbations caused by the Tohoku earthquake and tsunami of March 11, 2011. Perturbations near the epicenter were found in measurements of ionospheric total electron content (TEC) from 1198 GPS receivers in the Japanese GEONET network. For the first time for this event, we compare these observations with the estimated magnitude and speed of a tsunami-driven atmospheric gravity wave, using an atmosphere-ionosphere-coupling model and a tsunami model of sea-surface height, respectively. Traveling ionospheric disturbances (TIDs) were observed moving away from the epicenter at approximate speeds of 3400 m/s, 1000 m/s and 200–300 m/s, consistent with Rayleigh waves, acoustic waves, and gravity waves, respectively. We focus our analysis on gravity waves moving south and east of the epicenter, since tsunamis propagating in the deep ocean have been shown to produce gravity waves detectable in ionospheric TEC in the past. Observed southeastward gravity wave perturbations, seen ∼60 min after the earthquake, are mostly between 0.5 to 1.5 TECU, representing up to ∼5% of the background vertical TEC (VTEC). Comparisons of observed TID gravity waves with the modeled tsunami speed in the ocean and the predicted VTEC perturbation amplitudes from an atmosphere-ionosphere-coupling model show the measurements and models to be in close agreement. Due to the dense GPS network and high earthquake magnitude, these are the clearest observations to date of the effect of a major earthquake and tsunami on the ionosphere near the epicenter. Such observations from a future real-time GPS receiver network could be used to validate tsunami models, confirm the existence of a tsunami, or track its motion where in situ buoy data is not available.

Assimilation of GIRO data into a real‐time IRI

Abstract: Increasingly accurate and detailed global 3-D specification of the Earth's space plasma environment is required to further understand its intricate organization and behavior. For a long time space physics and aeronomy research has been data starved due to the great variety of natural time scales involved in the plasma phenomenology. We have started developing a new approach to the global ionospheric specification called Real-Time Assimilative Mapping (RTAM). The IRI-RTAM will use data from the Global Ionospheric Radio Observatory (GIRO) to smoothly transform International Reference Ionosphere's (IRI's) background empirical maps of the ionospheric characteristics to match the observations. Such empirical assimilative modeling will provide a high-resolution, global picture of the ionospheric response to various short-term events observed during periods of storm activity or the impact of gravity waves coupling the ionosphere to the lower atmosphere, including timelines of the vertical restructuring of the plasma distribution. It will also contribute to the challenging task of providing a rapid insight into the temporal and spatial space weather development using the real-time GIRO data streams. The new assimilation technique "updates" the IRI electron density distribution while preserving the overall integrity of IRI s typical ionospheric feature representations. The technique adjusts the coefficients of the spherical/diurnal expansions used by the CCIR and URSI-88 model to obtain the global sub-peak electron density distribution. The set of global corrected coefficients can be generated as frequently as every 15 min and easily disseminated using a single real-time RTAM server operated by GIRO.

Sharpening our thinking about polar cap ionospheric patch morphology, research, and mitigation techniques

Abstract: [1] Since polar cap patches were discovered, their nature, physics, and impact on navigation and communication signals has been repeatedly addressed. Both terminology and inference of physical processes from diverse instruments have introduced confusion. Poleward moving auroral form is a morphological descriptor but cannot be equated to an island of high-density plasma. Particle precipitation produces low-density patches, but high-density patches derive from solar produced plasma. The challenge of patches is finding the dominant mechanism for chopping entering ionization into islands (∼100–1000 km in size). Velocity-dependent recombination physics is valid in principle but not relevant to patches formed in at least the European sector. Most patches are formed by transient magnetic reconnection events. While the plasma is at too high an altitude for the strong velocity shears to erode plasma densities, the shears become the dominant plasma structuring mechanism until the initial magnetic tension force is relaxed. Initial patch structuring is not by gradient drift as believed for decades but rather by the shear driven instability, impacting mitigation techniques. Large-scale shears, driven to 2–3 km/s, impact satellite drag through thermospheric heating. The study here is intended to sharpen understanding of patches for future research and development of techniques for mitigation of their effects on navigation and other systems dependent on receiving radio frequency signals from satellites and understanding reconnection driven impact on thermospheric density and satellite drag.

Simulating the impacts of ionospheric scintillation on L band SAR image formation

Abstract: [1] We develop a phase screen model called the Synthetic Aperture Radar (SAR) Scintillation Simulator (SAR-SS) for predicting the impacts of ionospheric scintillation on SAR image formation. SAR-SS consists of a phase screen generator and a propagator. The screen generator creates a 2-D random realization of spatial phase fluctuations resulting from the traversal of small-scale field-aligned irregularities in the ionosphere. It accounts for the motion of the radar platform, the drift of the ionospheric irregularities, and the oblique angle of propagation, all of which determine the scale sizes of the irregularities sampled by the radar beam. The propagator solves the 3-D parabolic wave equation using the split step technique to compute the ionospheric transfer function for two-way propagation. This ionospheric transfer function is used to modulate the SAR signal due to terrestrial features in order to assess the ionospheric impact on SAR image formation in the small target approximation. We compare simulated and observed PALSAR imagery over Brazil during disturbed ionospheric conditions. We demonstrate that SAR-SS can reproduce the field-aligned streaks in PALSAR imagery caused by irregularities in the equatorial ionosphere that have been observed by previous authors. The field-aligned streaks exhibited a dominant wavelength larger than the Fresnel break scale, which suggests that refractive scatter was dominant over diffraction as the physical mechanism responsible for the scintillation of the radar signal in this case. The spectral index of phase fluctuations in the screen was quite large (9.0), suggesting that these irregularities were possibly associated with bottomside sinudoidal irregularities rather than equatorial plasma bubbles.

MAARSY: The new MST radar on Andøya—System description and first results

Abstract: [1] The Middle Atmosphere Alomar Radar System (MAARSY) on the North-Norwegian island Andoya is a 535 MHz monostatic radar with an active phased array antenna consisting of 433 Yagi antennas The 3-element Yagi antennas are arranged in an equilateral triangle grid forming a circular aperture of approximately 6300 m2 Each individual antenna is connected to its own transceiver with independent phase control and a scalable power output up to 2 kW This arrangement provides a very high flexibility of beam forming and beam steering with a symmetric radar beam of a minimum beam width of 36° allowing classical beam swinging operation as well as experiments with simultaneous multiple beams and the use of interferometric applications for improved studies of the Arctic atmosphere from the troposphere up to the lower thermosphere with high spatio-temporal resolution The installation of the antenna array was completed in August 2009 The radar control and data acquisition hardware as well as an initial expansion stage of 196 transceiver modules was installed in spring 2010 and upgraded to 343 transceiver modules in November 2010 The final extension to 433 transceiver modules has recently been completed in May 2011 Beside standard observations of tropospheric winds and Polar Mesosphere Summer Echoes, the first multi-beam experiments using up to 97 quasi-simultaneous beams in the mesosphere have been carried out in 2010 and 2011 These results provide a first insight into the horizontal variability of polar mesosphere summer and winter echoes with time resolutions between 3 and 9 minutes In addition, first meteor head echo observations were conducted during the Geminid meteor shower in December 2010

A fast algorithm for multiscale electromagnetic problems using interpolative decomposition and multilevel fast multipole algorithm

Abstract: [1] The interpolative decomposition (ID) is combined with the multilevel fast multipole algorithm (MLFMA), denoted by ID-MLFMA, to handle multiscale problems. The ID-MLFMA first generates ID levels by recursively dividing the boxes at the finest MLFMA level into smaller boxes. It is specifically shown that near-field interactions with respect to the MLFMA, in the form of the matrix vector multiplication (MVM), are efficiently approximated at the ID levels. Meanwhile, computations on far-field interactions at the MLFMA levels remain unchanged. Only a small portion of matrix entries are required to approximate coupling among well-separated boxes at the ID levels, and these submatrices can be filled without computing the complete original coupling matrix. It follows that the matrix filling in the ID-MLFMA becomes much less expensive. The memory consumed is thus greatly reduced and the MVM is accelerated as well. Several factors that may influence the accuracy, efficiency and reliability of the proposed ID-MLFMA are investigated by numerical experiments. Complex targets are calculated to demonstrate the capability of the ID-MLFMA algorithm.

Propagation of medium scale traveling ionospheric disturbances at different latitudes and solar cycle conditions

Abstract: [1] In this work, an extension in latitude range and time span with respect previous studies on Medium Scale Traveling Ionospheric Disturbances (MSTID) propagation, is presented. So far they have been basically studied at mid latitude and for limited periods (less than few years) at solar maximum conditions. This extension has been possible due to the availability of local Global Positioning System (GPS) networks at mid-north hemisphere (California), mid-south hemisphere (New Zealand), high and low latitudes (Alaska and Hawaii), for the last 13, 11, 7 and 4 years respectively. Optimal algorithms specially suitable for mass data processing have been used, such as the Single Receiver Medium Scale Traveling Ionospheric activity index (SRMTID) and the phase difference method for MSTID propagation estimation. The results reveal that several of the main MSTID climatological trends at mid latitude are also shared at low and high latitude, also modulated in intensity also by the Solar Cycle. This is the case for local fall/winter day-time equatorward propagated MSTIDs with typical velocities and wavelengths of 150–250 m/s and 100–300 km respectively. Moreover the comparison of MSTID propagation estimation using different techniques, and their implications in terms of potential origins of MSTIDs, are also discussed.

The Low‐Latitude Ionosphere Sensor Network: Initial results

Abstract: [1] The Low-Latitude Ionospheric Sensor Network (LISN) is a distributed observatory designed to nowcast the state and dynamics of the low-latitude ionosphere and to develop forecasts of the electric fields, densities, and equatorial spread F over the South American continent. The LISN observatory consists of three different types of instruments: GPS receivers, fluxgate magnetometers, and vertical incidence pulsed ionospheric radar (VIPIR) ionosondes. This report provides a succinct summary of recent observations obtained using the LISN GPS receivers and complemented with measurements from other instruments and GPS receivers that operate in South America. More specifically, the following are shown here: (1) observations of total electron content (TEC) enhancements that occur near local midnight, (2) maps of TEC perturbations associated with the passage of traveling ionospheric disturbances over South America, and (3) statistics of TEC depletions for 2 years of low solar activity. Near-midnight TEC enhancements consist of sudden increases in TEC that occur after sunset at low latitudes on 30% of the days. These TEC enhancements last for several hours and can have amplitudes between 1 and 50 TEC units. On 11–12 March 2011 the largest TEC enhancement was observed in South America at times when the Jicamarca incoherent scatter radar operated and observed peak densities above 106 el/cc at 300 km altitude. It is suggested that a combination of zonal electric fields and meridional neutral winds are able to redistribute the plasma along the field lines and create regions of enhanced TEC. Maps of TEC perturbations associated with the passage of gravity waves (GWs) over South America have been used to measure the phase velocity and direction of propagation of GWs. The large number of GPS receivers over South America has allowed us to record bubble events for every day during 2008 and 2009. It was found that the number of TEC depletion detections varies with a periodicity of 28 days. It is mentioned how these new observations and the installation of the last four VIPIR ionosondes will lead to new discoveries in the near future.

Contrasting the responses of three different ground‐based instruments to energetic electron precipitation

Abstract: [1] In order to make best use of the opportunities provided by space missions such as the Radiation Belt Storm Probes, we determine the response of complementary subionospheric radiowave propagation measurements (VLF), riometer absorption measurements, cosmic noise absorption, and GPS-produced total electron content (vTEC) to different energetic electron precipitation (EEP). We model the relative sensitivity and responses of these instruments to idealized monoenergetic beams of precipitating electrons, and more realistic EEP spectra chosen to represent radiation belts and substorm precipitation. In the monoenergetic beam case, we find riometers are more sensitive to the same EEP event occurring during the day than during the night, while subionospheric VLF shows the opposite relationship, and the change in vTEC is independent. In general, the subionospheric VLF measurements are much more sensitive than the other two techniques for EEP over 200 keV, responding to flux magnitudes two-three orders of magnitude smaller than detectable by a riometer. Detectable TEC changes only occur for extreme monoenergetic fluxes. For the radiation belt EEP case, clearly detectable subionospheric VLF responses are produced by daytime fluxes that are ∼10 times lower than required for riometers, while nighttime fluxes can be 10,000 times lower. Riometers are likely to respond only to radiation belt fluxes during the largest EEP events and vTEC is unlikely to be significantly disturbed by radiation belt EEP. For the substorm EEP case both the riometer absorption and the subionospheric VLF technique respond significantly, as does the change in vTEC, which is likely to be detectable at ∼3–4 total electron content units.

Resolution of the sprite polarity paradox: The role of halos

Abstract: [1] This study revisits the sprite polarity paradox, first manifest by observations that exceptional cloud-to-ground flashes with negative polarity generally did not produce detectable sprites. The paradox is here resolved by the Transient Luminous Event (TLE) known as the halo, which on account of its inferior brightness (0.3 MR versus 1.5 MR) and substantially shorter duration (1 ms versus 10–100 ms) in comparison with the sprite, is not readily detectable in ground-based video cameras with standard field duration (16.7–20 ms). Observations with improved temporal resolution (ISUAL (Imager of Sprites and Upper Atmospheric Lightnings) from space and PIPER (Photometric Imager of Precipitated Electron Radiation) observations from the ground) provide evidence that flashes with negative polarity dominate the global halo population, and that the halo numbers are more than sufficient to account for the previously missing TLEs. The evidence for lightning polarity-dependent TLEs (sprites, positive and halos, negative) is attributable to the well established but incompletely understood contrast in the behavior of negative and positive lightning flashes to ground.

VLSS redux: Software improvements applied to the Very Large Array Low‐Frequency Sky Survey

Abstract: [1] We present details of improvements to data processing and analysis which were recently used for a re-reduction of the Very Large Array (VLA) Low-Frequency Sky Survey (VLSS) data. Algorithms described are implemented in the data-reduction package Obit, and include smart-windowing to reduce clean bias, improved automatic radio frequency interference removal, improved bright-source peeling, and higher-order Zernike fits to model the ionospheric phase contributions. An additional, but less technical improvement was using the original VLSS catalog as a same-frequency/same-resolution reference for calculating ionospheric corrections, allowing more accuracy and a higher percentage of data for which solutions are found. We also discuss new algorithms for extracting a source catalog and analyzing ionospheric fluctuations present in the data. The improved reduction techniques led to substantial improvements including images of six previously unpublished fields (1% of the survey area) and reducing the clean bias by 50%. The largest angular size imaged has been roughly doubled, and the number of cataloged sources is increased by 35% to 95,000.

Long-term ionospheric anomaly monitoring for ground based augmentation systems

Abstract: [1] Extreme ionospheric anomalies can pose a potential integrity threat to ground-based augmentation of the Global Positioning System (GPS), and thus the development of ionospheric anomaly threat models for each region of operation is essential for system design and operation. This paper presents a methodology for automated long-term ionospheric anomaly monitoring, which will be used to build an ionospheric anomaly threat model, evaluate its validity over the life cycle of the system, continuously monitor ionospheric anomalies, and update the threat model if necessary. This procedure automatically processes GPS data collected from external networks and estimates ionospheric gradients at regular intervals. If ionospheric gradients large enough to be potentially hazardous to users are identified, manual data examination is triggered. This paper also develops a simplified truth processing method to create precise ionospheric delay estimates in near real-time, which is the key to automating the ionospheric monitoring procedure. The performance of the method is examined using data from the 20 November 2003 and 9 November 2004 ionospheric storms. These results demonstrate the effectiveness of simplified truth processing within long-term ionosphere monitoring. From the case studies, the automated procedure successfully identified extreme ionospheric anomalies, including the two worst ionospheric gradients observed and validated previously based on manual analysis. The automation of data processing enables us to analyze ionospheric data continuously going forward and to more accurately categorize ionospheric behavior under both nominal and anomalous conditions.

Introducing a disturbance ionosphere index

Abstract: [1] Although ionospheric perturbations such as traveling ionospheric disturbances have a strong impact on Global Navigation Satellite Systems (GNSS) and other space-based radio systems, the description of individual perturbations is difficult. To overcome this problem, it is suggested to use a disturbance ionosphere index (DIX) that describes the perturbation degree of the ionosphere in a less specific form as a proxy. Although such an index does not describe the exact propagation conditions at the measurement site, the estimated index number indicates the probability of a potential impact on radio systems used in communication, navigation, and remote sensing. The definition of such a DIX must take into account the following major requirements: relevance to practical needs, objective measure of ionospheric conditions, easy and reproducible computation, and availability of a reliable database. Since the total electron content has been shown in many publications to act as an outstanding parameter for quantifying the range error and also the strength of ionospheric perturbations, we propose a DIX that is based on GNSS measurements. To illustrate the use of the index, recent storms monitored in 2011 and the Halloween storm are discussed. The proposed index is a robust and objective measure of the ionospheric state, applicable to radio systems which are impacted by a highly variable perturbed ionosphere.

Absolute calibration of a wideband antenna and spectrometer for accurate sky noise temperature measurements

Abstract: [1] A new method of absolute calibration of sky noise temperature using a three-position switched spectrometer, measurements of antenna and low noise amplifier impedance with a vector network analyzer, and ancillary measurements of the amplifier noise waves is described and compared with other methods of calibration. An initial test of the method has been made using brief observations from 55 to 110 MHz at West Forks, Maine, to estimate the sky noise spectral index. Estimates are made of the accuracy that might ultimately be achieved with observations to detect or set limits on the red-shifted 21 cm line at a radio quiet site. It is concluded that an antenna reflection coefficient better than about −20 dB is required to avoid being limited by the accuracy of the antenna reflection coefficient measurements using a vector network analyzer.

Characteristics of traveling ionospheric disturbances observed by the TIDDBIT sounder

Abstract: [1] HF Doppler sounders represent a low-cost and low-maintenance solution for monitoring wave activity in theF region ionosphere. HF Doppler sounders together with modern data analysis techniques can provide comprehensive traveling ionospheric disturbance (TID) characteristics, including both horizontal and vertical TID velocities and wavelengths across the entire spectrum from periods of 1 min to over an hour. Atmospheric and Space Technology Research Associates LLC has developed a new system called “TIDDBIT” (TID Detector Built in Texas), and data will be presented from a TIDDBIT system deployed in Virginia. These results reinforce the relationship between atmospheric gravity waves (AGWs) and TIDs. The TID propagation azimuths rotate through 360° in 24 h, mimicking the rotation of the thermospheric winds but with approximately a 90° offset. The rotation of TID azimuths and thermospheric winds in Virginia is similar to that observed previously by other Northern Hemisphere systems and opposite from the direction observed in Antarctica. These results illustrate the filtering effects that thermospheric neutral winds can have on the propagation of AGW. The completeness of the wave information obtained from the TIDDBIT system makes it possible to reconstruct the vertical displacement of isoionic contours over the ∼200 km horizontal dimension of the sounder array. Such information will be relevant for understanding the seeding of irregularities, as well as for several operational needs involving navigation, communication, and surveillance systems.

Wide Area RTK: A satellite navigation system based on precise real‐time ionospheric modelling

Abstract: [1] The Wide Area Real Time Kinematic (WARTK) is an augmentation system concept for multi-frequency users based on precise real-time ionospheric modeling It is able to provide a high accuracy and integrity GNSS positioning service over continental areas using the infrastructure of a network of permanent ground monitor stations, such as the European Geostationary Navigation Overlay Service (EGNOS) network of Ranging and Integrity Monitoring Stations (RIMS) in Europe In this way, it allows an additional benefit to be obtained from these reference stations, that is, the network has the potential to support two independent systems: a satellite-based augmentation system, such as EGNOS, and a high-precision positioning service, based on WARTK Indeed, thanks to the accuracy of the ionospheric corrections provided, WARTK users have available in real-time an extra constraint per satellite between the carrier phase ambiguities, which helps solve them quickly Once such ambiguities have been solved, the GNSS user obtains navigation accurate to within 20 cm at the 95th percentile (about 10 cm RMS) Moreover, this precise positioning is achieved in a few minutes (with two frequency signals) or in a single epoch, after initial convergence of the tropospheric delay (with three frequency signals), even up to hundreds of kilometers away from the nearest reference station While previous WARTK research has been devoted to implementing the concept and assessing its feasibility, considering in particular the accuracy achievable, the work reported in this paper focused on consolidating the results by analyzing a large and representative data set, and on deeper analysis of the integrity issue It was carried out in the context of the Multi-constellation Regional System (MRS) project, within the European Space Agency GNSS Evolution Programme, with the aim of designing a high accuracy service for GPS and/or Galileo Three months of actual data, from more than 25 permanent GPS stations in Europe, have been processed (some of them as a roving user), for high-, mid- and low-solar cycle conditions (in 2002, 2004 and 2006 respectively) In addition, several ionospheric storms occurred during the selected periods, with Dst values reaching up to � 150 nT Results based on these data show that user domain integrity was maintained for baselines of up to 400 km At the 95th percentile, the daily horizontal and vertical position errors were 20 and 30 cm, respectively, and the corresponding protection levels were about 1 and 2 m The convergence time was around 5 minutes with actual GPS constellation data The benefits of using a multi-constellation system were also studied, with simulated GPS and three-frequency Galileo data, showing that it is possible to reduce the convergence time to a few seconds

A comprehensive analysis of free-space and guided-wave techniques for extracting the permeability and permittivity of materials using reflection-only measurements

Abstract: [1] The electromagnetic characterization of layered materials is often performed using reflection-only measurements, and through the years several different techniques have been developed to accomplish this both in the laboratory and in situ. It is shown in this paper that the majority of these methods are variations on a single approach, and that closed-form expressions are available to obtain the permittivity and permeability from the measurements. The approach is applicable to free-space (plane wave) measurements, waveguide (non-TEM) measurements, and transmission line (TEM) measurements, and is based on varying the underlying structure of the material stack to obtain independent measurements. Because the formulation is based on specifying the impedances of adjacent regions, the same formula for permittivity and permeability is valid regardless of the complexity of the stack containing the material under test. Examples of several techniques are given using free-space, waveguide, and transmission-line measurements of both dielectric materials and magnetic absorbing materials.

A fast direct matrix solver for surface integral equation methods for electromagnetic wave scattering from non-penetrable targets

Abstract: [1] The implementation details of a fast direct solver is described herein for solving dense matrix equations from the application of surface integral equation methods for electromagnetic field scatterings from non-penetrable targets. The proposed algorithm exploits the smoothness of the far field and computes a low rank decomposition of the off-diagonal coupling blocks of the matrices through a set of skeletonization processes. Moreover, an artificial surface (the Huygens' surface) is introduced for each clustering group to efficiently account for the couplings between well-separated groups. Furthermore, a recursive multilevel version of the algorithm is presented. Although asymptotically the algorithm would not alter the bleak outlook of the complexity of the worst case scenario,O(N3) for required CPU time where N denotes the number of unknowns, for electrically large electromagnetic (EM) problems; through numerical examples, we found that the proposed multilevel direct solver can scale as good as O(N1.3) in memory consumption and O(N1.8) in CPU time for moderate-sized EM problems. Note that our conclusions are drawn based on a few sample examples that we have conducted and should not be taken as a true complexity analysis for general electrodynamic applications. However, for the fixed frequency (h-refinement) scenario, where the discretization size decreases, the computational complexities observed agree well with the theoretical predictions. Namely, the algorithm exhibits O(N) and O(N1.5) complexities for memory consumption and CPU time, respectively.

Impacts of ionospheric scintillations on GPS receivers intended for equatorial aviation applications

Abstract: [1] This study examines the impacts of ionospheric scintillations on GPS receivers that are intended for equatorial or transequatorial aviation applications. We analyzed GPS data that were acquired at Ascension Island during the Air Force Research Laboratory (AFRL) campaign of the solar maximum year of 2002. Strong scintillations impacted the receiver-satellite geometry, leading to poor dilution of precisions and positioning accuracy. In addition, deep signal fades (>20 dB-Hz), leading to navigation outages were observed during most of the nights of the campaign. Under quiescent conditions, the C/No of satellites fluctuated slowly between 50 dB-Hz and 35 dB-Hz baselines for both L1 (1.5754 GHz) and L2 (1.2276 GHz) signals, depending on the satellite's elevation angle. The satellite's elevation angle and the effective scan velocity of the satellite's ionospheric penetration point (IPP) with respect to the magnetic field and plasma drift influenced the rate of fading of satellite signals.

A methodology to compute GPS slant total delays in a numerical weather model

Abstract: [1] A numerical algorithm based on Fermat's Principle was developed to simulate the propagation of Global Positioning System (GPS) radio signals in the refractivity field of a numerical weather model. The unique in the proposed algorithm is that the ray-trajectory automatically involves the location of the ground-based receiver and the satellite, i.e. the posed two-point boundary value problem is solved by an implicit finite difference scheme. This feature of the algorithm allows the fast and accurate computation of the signal travel-time delay, referred to as Slant Total Delay (STD), between a satellite and a ground-based receiver. We provide a technical description of the algorithm and estimate the uncertainty of STDs due to simplifying assumptions in the algorithm and due to the uncertainty of the refractivity field. In a first application, we compare STDs retrieved from GPS phase-observations at the German Research Centre for Geosciences Potsdam (GFZ STDs) with STDs derived from the European Center for Medium-Range Weather Forecasts analyses (ECMWF STDs). The statistical comparison for one month (August 2007) for a large and continuously operating network of ground-based receivers in Germany indicates good agreement between GFZ STDs and ECMWF STDs; the standard deviation is 0.5% and the mean deviation is 0.1%.

Ewald method for 3D periodic dyadic Green's functions and complex modes in composite materials made of spherical particles under the dual dipole approximation

Abstract: [1] We derive the Ewald representation for the dyadic periodic Green's functions to represent the electromagnetic field in a three dimensional (3D) periodic array of electric and magnetic dipoles. Then we use the developed theory to analyze the modes with real and complex wave number in a 3D periodic lattice of lead telluride (PbTe) microspheres at infrared frequencies and in a 3D periodic lattice of titanium dioxide (TiO2) microspheres at millimeter waves. Each microsphere is equivalently modeled with both an electric and a magnetic dipole, via a method here called the dual dipole approximation (DDA). The 3D lattices exhibit first a magnetic-induced then an electric-induced feature determined by microsphere magnetic and electric resonances. The DDA wave number results are compared to the ones computed with single electric or single magnetic dipole approximation and to the ones retrieved by using the Nicolson-Ross-Weir (NRW) retrieval method from reflection and transmission of finite thickness slabs computed by a full-wave simulation. It is shown that the DDA method is in very good agreement with NRW, in contrast to the previously reported single dipole approximation methods that fail to predict one of the two features (either electric or magnetic). A mode with transverse polarization is found to be dominant and able to propagate inside the lattice, and therefore the composite material can be treated as a homogeneous one with effective refractive index. This is obtained by adopting five different retrieval procedures for each lattice, and their agreement or disagreement is discussed.

Variability of total electron content over an equatorial West African station during low solar activity

Abstract: [1] Slant total electron content (STEC) data measured by the Global Positioning System receiver at Ilorin, Nigeria, with geographical coordinates 8.47°N, 4.68°E for the year 2009 (a low-activity year) was used to study the diurnal, monthly standard deviation and monthly median value of total electron content (TEC). The vertical total electron content (VTEC) values are estimated from the STEC data. The thin shell approximation with an ionospheric shell height of 350 km was used for the analysis. The diurnal variation of VTEC (DTEC) and its corresponding monthly median variation (MTEC) shows a minimum at presunrise between the hours of 05:00 and 06:00 LT. The DTEC values show a maximum variation range from ∼24 to ∼34 total electron content unit (TECU). The daytime maximum TEC values observed in all the months were broad with a slight daytime depression in May, June, July, and November. The maximum variation of MTEC after slight daytime depression is greater than its variation before the slight daytime depression in the months affected with the month of July as exception. The slight daytime depression was lowest in the month of May and has a value of 0.99 TECU. A postsunset decrease at 20:00 LT with corresponding enhancement 2 h later was observed in the month of March. This post sunset decrease and enhancement in the month of March could be a strong indicator of the abrupt onset of scintillations, plasma bubbles, and spread F phenomenon. The monthly standard deviation depicts summary behavior of all the diurnal variations in each month. Annual and seasonal variations were also investigated.

New Vary‐Chap profile of the topside ionosphere electron density distribution for use with the IRI model and the GIRO real time data

Abstract: A new Vary-Chap function is introduced for the empirical modeling of the electron density N(h) profile in the topside ionosphere that uses a shape function S(h) in the generalized Chapman function. The Vary-Chap profile extends the bottomside profile that is specified by the IRI model or measured by the Global Ionospheric Radio Observatory (GIRO) to the altitude of the ISIS-2 satellite. Some 80,000 topside profiles, measured by the topside sounder on the ISIS-2 satellite were analyzed, and the shape function S(h) was calculated for each profile. A parameterized function S*(h), composed of two sub-functions S1(h) and S2(h), is fitted to the measured S(h) profile using three free parameters. At altitudes just above the F2 layer peak height hmF2, the shape function S1 controls S(h), and at greater altitudes S2 controls S(h). The height of the intersection of S1 and S2 is defined as the transition height h(sub T) indicating the transition from an O(+) to an H(+)-dominated profile shape. The observed transition heights range from approx.500 km to 800 km.

Global observations of L band scintillation at solar minimum made by COSMIC

Abstract: [1] We report observations of the L band scintillation climatology made during the minimum of solar cycle 23/24 using the Constellation Observing System for Meteorology, Ionosphere, and Climate (COSMIC, also known as FORMOSAT-3). The measurements of the S4 scintillation index were made using the GPS Occultation Experiment instrument, which is a dual-frequency GPS receiver used for routine ionospheric measurements via occultation of satellites in the GPS constellation. A method for geolocating the S4 measurements is presented and discussed. The geolocated COSMIC measurements are used to construct monthly climatologies of the S4 scintillation index, which are in good agreement with previous climatologies based on ground-based measurements. COSMIC measurements of the electron density and S4 index permit an assessment of the correlation of the S4 index with the electron density. The scintillation climatology for March and April 2007 is compared to the coincident monthly climatology of the electron density, and a correlation is observed, but this correlation does not appear to hold globally. The COSMIC measurements are shown to be a powerful data set for studying L band scintillation and its correlation with the underlying ionospheric morphology.

Sounding of HF heating-induced artificial ionospheric disturbances by navigational satellite radio transmissions

Abstract: [1] During experiments carried out in 2009–2011 the midlatitude ionosphere was modified by powerful HF pulses from the Sura heating facility located near Nizhny Novgorod (Russia) and operated by the Radio Physical Research Institute. GPS/GLONASS and Parus/Tsikada satellite radio transmissions responding to the heating-induced disturbances in electron density were analyzed. The variations in the total electron content (TEC), which are proportional to the reduced phase of navigational signals, were studied for various schemes of radiation of the heating wave. The variations in TEC (their amplitudes and temporal behavior) caused by HF heating are identified in several examples. The TEC spectra contain frequency components corresponding to the modulation periods of the heating wave. For the first time, the spatial structure of the wave disturbances generated in the ionosphere by high-power radio waves radiated by the Sura heating facility with a square wave modulation of the effective radiated power at a frequency lower than or of the order of the Brunt-Vaisala frequency of the neutral atmosphere is imaged using the method of low-orbital radio tomography and GPS/GLONASS data.

Ionospheric and auroral clutter models for HF surface wave and over‐the‐horizon radar systems

Abstract: [1] The detection performance of high frequency surface wave radar (HFSWR) and high frequency over-the-horizon radar (OTHR) systems is heavily influenced by the presence of radar clutter. In HFSWR systems, the clutter has its origins in vertical-incidence ionospheric reflections, whereas in OTHR systems, the origin is Bragg backscatter from plasma structures in the auroral zone. This paper models the spreading of the radar clutter signal in the Doppler and angle-of-arrival domains that arises from forward-scattering effects as the radar pulse propagates through regions of ionospheric plasma irregularities. The models use a geometric optics approach to determine the power spectrum of the radar signal phase. This power spectrum is then used to simulate three-dimensional space-time-range radar data cubes. The accuracy of the models is tested by comparing the simulated data to measured data cubes. As an application, the data are then used to evaluate the performance of the newly developed fast fully adaptive (FFA) space-time adaptive processing (STAP) scheme to improve the extraction of target echoes from a clutter background.

Ray-traced tropospheric delays in VLBI analysis

Abstract: [1] We develop a ray-tracing package for the calculation of path delays of microwave signals in the troposphere based on numerical weather models which we use for the determination of the delays of geodetic Very Long Baseline Interferometry (VLBI) observations. We show results for a two-week campaign of continuous VLBI sessions in 2008 (CONT08), where we apply those ray-traced delays and analyze the repeatability of baseline lengths in comparison to a standard approach with zenith delays and mapping functions. We find improvement in baseline length repeatabilities when no tropospheric gradients are estimated in the analysis. Furthermore, ray-traced delays are applied for Intensive sessions containing the stations Tsukuba (Japan) and Wettzell (Germany) for the determination of Universal Time (UT1). We perform an external validation using GPS-derived length-of-day values and find an improvement for UT1 with ray-traced delays by up to 4.5%.

Time‐reversal techniques for MISO and MIMO wireless communication systems

Abstract: [1] We consider the application of different time-reversal (TR) signal processing and beamforming techniques to multiple-input single-output (MISO) and multiple-input multiple-output (MIMO) wireless communication systems. Conventional TR beamforming provides spatial focusing at the intended receiver; however, it does not yield perfect channel equalization. Time-reversed pilot can be normalized to provide perfect equalization at the expense of power level. This equalization is particularly important for high data rates where the bit error rate performance is dominated by internal noise due to intersymbol interference. To increase physical layer covertness, TR beamforming is combined with the multiple-signal-classification (MUSIC) technique to produce null fields at eavesdroppers. This technique is also applied to MIMO setups to eliminate interuser interference and hence increase system capacity. Differential TR is used to obtain and update pilot signals for passive moving receivers, i.e., those that cannot (or do not) transmit pilot signals. Time-reversed differential backscattered signal is able to provide satisfactory spatial and temporal focusing at the moving receiver.

Second-kind integral solvers for TE and TM problems of diffraction by open arcs

Abstract: We present a novel approach for the numerical solution of problems of diffraction by open arcs in two dimensional space. Our methodology relies on composition of weighted versions of the classical integral operators associated with the Dirichlet and Neumann problems (TE and TM polarizations, respectively) together with a generalization to the open-arc case of the well known closed-surface Calderon formulae. When used in conjunction with spectrally accurate discretization rules and Krylov-subspace linear algebra solvers such as GMRES, the new second-kind TE and TM formulations for open arcs produce results of high accuracy in small numbers of iterations—for low and high frequencies alike.