Ring Currents and Internal Plasma Sources

The discovery of terrestrial O⁺ and other heavy ions in magnetospheric hot plasmas, combined with the association of energetic ionospheric outflows with geomagnetic activity, led to the conclusion that increasing geomagnetic activity is responsible for filling the magnetosphere with ionospheric plasma. Recently it has been discovered that a major source of ionospheric heavy ion plasma outflow is responsive to the earliest impact of coronal mass ejecta upon the dayside ionosphere. Thus a large increase in ionospheric outflows begins promptly during the initial phase of geomagnetic storms, and is already present during the main phase development of such storms. We hypothesize that enhancement of the internal source of plasma actually supports the transition from substorm enhancements of aurora to storm-time ring current development in the inner magnetosphere. Other planets known to have ring current-like plasmas also have substantial internal sources of plasma, notably Jupiter and Saturn. One planet having a small magnetosphere, but very little internal source of plasma, is Mercury. Observations suggest that Mercury has substorms, but are ambiguous with regard to the possibility of magnetic storms of the planet. The Messenger mission to Mercury should provide an interesting test of our hypothesis. Mercury should support at most a modest ring current if its internal plasma source is as small as is currently believed. If substantiated, this hypothesis would support a general conclusion that the magnetospheric inflationary response is a characteristic of magnetospheres with substantial internal plasma sources. We quantitatively define this hypothesis and pose it as a problem in comparative magnetospheres.     

Scientific Planning and Commanding of the Rosetta Payload

ESA’s Rosetta mission was launched in March 2004 and is on its way to comet 67P/Churyumov-Gerasimenko, where it is scheduled to arrive in summer 2014. It comprises a payload of 12 scientific instruments and a Lander. All instruments are provided by Principal Investigators, which are responsible for their operations. As for most ESA science missions, the ground segment of the mission consists of a Mission Operations Centre (MOC) and a Science Operations Centre (SOC). While the MOC is responsible for all spacecraft-related aspects and the final uplink of all command timelines to the spacecraft, the scientific operations of the instruments and the collection of the data and ingestion into the Planetary Science Archive are coordinated by the SOC. This paper focuses on the tasks of the SOC and in particular on the methodology and constraints to convert the scientific goals of the Rosetta mission to operational timelines.     

Integral-projection method and singularities of its application in mating shells with inconsistent parameters of approximation

We describe a numerical integral-projection method used by the authors for the approximate solution of systems of interrelated two-dimensional linear boundary-value problems in mechanics of composite shell systems. The method is based on discretization in each shell substructure of a two-dimensional problem along one of coordinates using a projection-grid variant of the Galerkin-Petrov method and its subsequent transformation to a system of ordinary differential equations; by integration and introduction of sought functions as unknown derivatives, the system is reduced to a system of integral equations being solved by the method of mechanical quadratures. The method is characterized by the fact that its application requires no additional conditions of conformity with discretization parameters of substructures being mated.     

How do interplanetary shock impact angles control the size of the geoeffective magnetosphere?

In this paper, we investigate temporal and spatial magnetosphere response to the impact of interplanetary (IP) shocks with different inclinations and speeds on the Earth’s magnetosphere. A data set with more than 500 IP shocks is used to identify positive sudden impulse (SI⁺) events as expressed by the SuperMAG partial ring current index. The SI⁺ rise time (RT), defined as the time interval between compression onset and maximum SI⁺ signature, is obtained for each event. We use RT and a model suggested by Takeuchi et al. (2002) to calculate the geoeffective magnetospheric distance (GMD) in the shock propagation direction as a function of shock impact angle and speed for each event. GMD is a generalization of the geoeffective magnetosphere length (GML) suggested by Takeuchi et al. (2002), defined from the subsolar point along the X line toward the tail. We estimate statistical GMD and GML values which are then reported for the first time. We also show that, similarly to well-known results for RT, the highest correlation coefficient for the GMD and impact angle is found for shocks with high speeds and small impact angles, and the faster and more frontal the shock, the smaller the GMD. This result indicates that the magnetospheric response depends heavily on shock impact angle. With these results, we argue that the prediction and forecasting of space weather events, such as those caused by coronal mass ejections, will not be accurately accomplished if the disturbances’ angles of impact are not considered as an important parameter within model and observation scheme capabilities.     

Diurnal and seasonal variations in long-term changes in the E-layer critical frequency

Long-term changes in the E-layer critical frequency, foE, at three stations of the European region (Juliusruh, Slough and Rome) and also at Moscow and Wakkanai stations are analyzed by the method developed by the authors and described in detail in the previous papers. It is found that Juliusruh and Slough stations demonstrate a well-pronounced change in foE (a trend) during two previous decades. At the same time, the same features of the behavior of the aforementioned trend k(foE) are obtained. The trend is positive and negative in the morning and evening hours, respectively. Similar diurnal behavior of k(foE) is found also for Moscow station but with lower absolute values of the trends. A well-pronounced seasonal behavior of k(foE) is detected at Juliusruh and Slough: the trend is minimal and maximal in the summer period and at the end of fall—beginning of winter, respectively. The maximal amplitude in the morning hours reaches +0.04?MHz per year, whereas the minimal amplitude in evening hours is ?0.06?MHz per year. No systematic changes exceeding by the magnitude 0.01?MHz per year are found for Rome and Wakkanai stations. It is assumed that the observed trends are related to changes (trends) in the meridional wind bringing NO molecules from the auroral oval to lower latitudes.     

Longitudinal variation of the equatorial ionosphere: Modeling and experimental results

We describe a new version of the Parameterized Regional Ionospheric Model (PARIM) which has been modified to include the longitudinal dependences. This model has been reconstructed using multidimensional Fourier series. To validate PARIM results, the South America maps of critical frequencies for the E (foE) and F (foF2) regions were compared with the values calculated by Sheffield Plasmasphere-Ionosphere Model (SUPIM) and IRI representations. PARIM presents very good results, the general characteristics of both regions, mainly the presence of the equatorial ionization anomaly, were well reproduced for equinoctial conditions of solar minimum and maximum. The values of foF2 and hmF2 recorded over Jicamarca (12°S; 77°W; dip lat. 1°N; mag. declination 0.3°) and sites of the conjugate point equatorial experiment (COPEX) campaign Boa Vista (2.8°N; 60.7°W; dip lat. 11.4°; mag. declination −13.1°), Cachimbo (9.5°S; 54.8°W; dip lat. −1.8°; mag. declination −15.5°), and Campo Grande (20.4°S; 54.6°W; dip lat. −11.1°; mag. declination −14.0°) have been used in this work. foF2 calculated by PARIM show good agreement with the observations, except during morning over Boa Vista and midnight-morning over Campo Grande. Some discrepancies were also found for the F-region peak height (hmF2) near the geomagnetic equator during times of F₃ layer occurrences. IRI has underestimated both foF2 and hmF2 over equatorial and low latitude sectors during evening-nighttimes, except for Jicamarca where foF2 values were overestimated.     

Vertical motions in Thailand after the 2004 Sumatra–Andaman Earthquake from GPS observations and its geophysical modelling

Following previous findings from ongoing GPS research in Thailand since 2004 we continue to exploit the GPS technique to monitor and model land motions induced by the Sumatra–Andaman Earthquake. Our latest results show that up to the end of 2010, Thailand has been co-seismically displaced and is subsequently undergoing a post-seismic horizontal deformation with total displacements (co-seismic plus post-seismic) ranging from 10.5 to 74.7 cm. We observed the largest horizontal displacements in the southern part of Thailand and moderate and small displacements in the central and northern parts. In addition to horizontal displacements throughout Thailand, continuous GPS measurements show that large parts of Thailand are subsiding at rates up to 1 cm/yr. It is the first time that such vertical post-seismic deformations at large distances (650–1500 km away from the Earthquake’s epicentre) have been recorded. We have investigated the physical processes leading to the observed subsidence. While after-slip on the subduction interface induces negligible or even slightly positive vertical motions, relaxation in the asthenosphere is associated with a sizable subsidence. Predictions from a 3D finite element model feature an asthenosphere with an effective viscosity of the order of 3 * 10¹⁸ Pas, fit the horizontal post-seismic data and the observed subsidence well. This model is then used to predict the subsidence over the whole seismic cycle. The subsidence should go on with a diminishing rate through the next two decades and its final magnitude should not exceed 10 cm in the Bangkok area.     

The effect of geocenter motion on Jason-2 orbits and the mean sea level

We compute a series of Jason-2 GPS and SLR/DORIS-based orbits using ITRF2005 and the std0905 standards ( Lemoine et al., 2010). Our GPS and SLR/DORIS orbit data sets span a period of 2 years from cycle 3 (July 2008) to cycle 74 (July 2010). We extract the Jason-2 orbit frame translational parameters per cycle by the means of a Helmert transformation between a set of reference orbits and a set of test orbits. We compare the annual terms of these time-series to the annual terms of two different geocenter motion models where biases and trends have been removed. Subsequently, we include the annual terms of the modeled geocenter motion as a degree-1 loading displacement correction to the GPS and SLR/DORIS tracking network of the POD process. Although the annual geocenter motion correction would reflect a stationary signal in time, under ideal conditions, the whole geocenter motion is a non-stationary process that includes secular trends. Our results suggest that our GSFC Jason-2 GPS-based orbits are closely tied to the center of mass (CM) of the Earth consistent with our current force modeling, whereas GSFC’s SLR/DORIS-based orbits are tied to the origin of ITRF2005, which is the center of figure (CF) for sub-secular scales. We quantify the GPS and SLR/DORIS orbit centering and how this impacts the orbit radial error over the globe, which is assimilated into mean sea level (MSL) error, from the omission of the annual term of the geocenter correction. We find that for the SLR/DORIS std0905 orbits, currently used by the oceanographic community, only the negligence of the annual term of the geocenter motion correction results in a – 4.67 ± 3.40 mm error in the Z-component of the orbit frame which creates 1.06 ± 2.66 mm of systematic error in the MSL estimates, mainly due to the uneven distribution of the oceans between the North and South hemisphere.     

Local-time variations of geomagnetic disturbances during intense geomagnetic storms and possible association with their interplanetary causes

In the present paper the local-time variations in the disturbance of the geomagnetic-field horizontal component (H) for eight intense geomagnetic storms that occurred during the descending phase of solar cycle 23 have been analyzed. The study was based on the plot of contour lines of the H-depletion intensity in the plane local time versus universal time (LT–UT maps) with the objective of observing how the morphology and evolution of the ring current is mapped into the surface of the Earth in presence of intense geomagnetic storms.