Atmospheric, oceanic and hydrological polar motion excitation mechanisms based on a combination of geometric and gravimetric space observations

The goal of our study is to determine accurate time series of geophysical Earth rotation excitations to learn more about global dynamic processes in the Earth system. For this purpose, we developed an adjustment model which allows to combine precise observations from space geodetic observation systems, such as Satellite Laser Ranging (SLR), Global Navigation Satellite Systems (GNSS), Very Long Baseline Interferometry (VLBI), Doppler Orbit determination and Radiopositioning Integrated on Satellite (DORIS), satellite altimetry and satellite gravimetry in order to separate geophysical excitation mechanisms of Earth rotation. Three polar motion time series are applied to derive the polar motion excitation functions (integral effect). Furthermore we use five time variable gravity field solutions from Gravity Recovery and Climate Experiment (GRACE) to determine not only the integral mass effect but also the oceanic and hydrological mass effects by applying suitable filter techniques and a land-ocean mask. For comparison the integral mass effect is also derived from degree 2 potential coefficients that are estimated from SLR observations. The oceanic mass effect is also determined from sea level anomalies observed by satellite altimetry by reducing the steric sea level anomalies derived from temperature and salinity fields of the oceans. Due to the combination of all geodetic estimated excitations the weaknesses of the individual processing strategies can be reduced and the technique-specific strengths can be accounted for. The formal errors of the adjusted geodetic solutions are smaller than the RMS differences of the geophysical model solutions. The improved excitation time series can be used to improve the geophysical modeling.

Análise de séries temporais de coordenadas estimadas com GPS: uma proposta metodológica para eliminação de efeitos sazonais

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Earth Rotation and Gravity Field Parameters from Satellite Laser Ranging

We present the results from a simultaneous estimation of the gravity field, Earth rotation parameters, and station coordinates from combined SLR solutions incorporating up to nine geodetic satellites: LAGEOS-1/2, Starlette, Stella, AJISAI, Beacon-C, Lares, Blits and LARES. These solutions cover all three pillars of satellite geodesy and ensure full consistency between the Earth rotation parameters, gravity field coefficients, and geometry-related parameters. We address benefits emerging from such an approach and discuss particular aspects and limitations of the gravity field recovery using SLR data. The current accuracy of SLR-derived polar motion, by the means of WRMS w.r.t. IERS-08-C04 series, is at a level of 118-149 μas, which corresponds to 4 to 5 mm on the Earth’s surface. The WRMS of SLR-derived Length-of-Day, when the gravity field parameters are simultaneously estimated, is 56 μs/day, corresponding to about 26 mm on the ground, and the mean bias of SLR-derived Length-of-Day is 6.3 μs/day, corresponding to 3 mm.

Impact of the arc length on GNSS analysis results

Homogeneously reprocessed combined GPS/GLONASS 1- and 3-day solutions from 1994 to 2013, generated by the Center for Orbit Determination in Europe (CODE) in the frame of the second reprocessing campaign REPRO-2 of the International GNSS Service, as well as GPS- and GLONASS-only 1- and 3-day solutions for the years 2009 to 2011 are analyzed to assess the impact of the arc length on the estimated Earth Orientation Parameters (EOP, namely polar motion and length of day), on the geocenter, and on the orbits. The conventional CODE 3-day solutions assume continuity of orbits, polar motion components, and of other parameters at the day boundaries. An experimental 3-day solution, which assumes continuity of the orbits, but independence from day to day for all other parameters, as well as a non-overlapping 3-day solution, is included into our analysis. The time series of EOPs, geocenter coordinates, and orbit misclosures, are analyzed. The long-arc solutions were found to be superior to the 1-day solutions: the RMS values of EOP and geocenter series are typically reduced between 10 and 40 %, except for the polar motion rates, where RMS reductions by factors of 2–3 with respect to the 1-day solutions are achieved for the overlapping and the non-overlapping 3-day solutions. In the low-frequency part of the spectrum, the reduction is even more important. The better performance of the orbits of 3-day solutions with respect to 1-day solutions is also confirmed by the validation with satellite laser ranging.

Earth’s Rotation: A Challenging Problem in Mathematics and Physics

A suitable knowledge of the orientation and motion of the Earth in space is a common need in various fields. That knowledge has been ever necessary to carry out astronomical observations, but with the advent of the space age, it became essential for making observations of satellites and predicting and determining their orbits, and for observing the Earth from space as well. Given the relevant role it plays in Space Geodesy, Earth rotation is considered as one of the three pillars of Geodesy, the other two being geometry and gravity. Besides, research on Earth rotation has fostered advances in many fields, such as Mathematics, Astronomy and Geophysics, for centuries. One remarkable feature of the problem is in the extreme requirements of accuracy that must be fulfilled in the near future, about a millimetre on the tangent plane to the planet surface, roughly speaking. That challenges all of the theories that have been devised and used to-date; the paper makes a short review of some of the most relevant methods, which can be envisaged as milestones in Earth rotation research, emphasizing the Hamiltonian approach developed by the authors. Some contemporary problems are presented, as well as the main lines of future research prospected by the International Astronomical Union/International Association of Geodesy Joint Working Group on Theory of Earth Rotation, created in 2013.

Motion of the dayside polar cap boundary during substorm cycles: II. Generation of poleward-moving events and polar cap patches by pulses in the magnetopause reconnection rate

Using data from the EISCAT (European Incoherent Scatter) VHF and CUTLASS (Co-operative UK Twin- Located Auroral Sounding System) HF radars, we study the formation of ionospheric polar cap patches and their relationship to the magnetopause reconnection pulses identified in the companion paper by Lockwood et al. (2005). It is shown that the poleward-moving, high-concentration plasma patches observed in the ionosphere by EISCAT on 23 November 1999, as reported by Davies et al. (2002), were often associated with corresponding reconnection rate pulses. However, not all such pulses generated a patch and only within a limited MLT range (11:00–12:00 MLT) did a patch result from a reconnection pulse. Three proposed mechanisms for the production of patches, and of the concentration minima that separate them, are analysed and evaluated: (1) concentration enhancement within the patches by cusp/cleft precipitation; (2) plasma depletion in the minima between the patches by fast plasma flows; and (3) intermittent injection of photoionisation-enhanced plasma into the polar cap. We devise a test to distinguish between the effects of these mechanisms. Some of the events repeat too frequently to apply the test. Others have sufficiently long repeat periods and mechanism (3) is shown to be the only explanation of three of the longer-lived patches seen on this day. However, effect (2) also appears to contribute to some events. We conclude that plasma concentration gradients on the edges of the larger patches arise mainly from local time variations in the subauroral plasma, via the mechanism proposed by Lockwood et al. (2000).

Motion of the dayside polar cap boundary during substorm cycles: I. Observations of pulses in the magnetopause reconnection rate

Using data from the EISCAT (European Incoherent Scatter) VHF radar and DMSP (Defense Meteorological Satellite Program) spacecraft passes, we study the motion of the dayside open-closed field line boundary during two substorm cycles. The satellite data show that the motions of ion and electron temperature boundaries in EISCAT data, as reported by Moen et al. (2004), are not localised around the radar; rather, they reflect motions of the open-closed field line boundary at all MLT throughout the dayside auroral ionosphere. The boundary is shown to erode equatorward when the IMF points southward, consistent with the effect of magnetopause reconnection. During the substorm expansion and recovery phases, the dayside boundary returns poleward, whether the IMF points northward or southward. However, the poleward retreat was much faster during the substorm for which the IMF had returned to northward than for the substorm for which the IMF remained southward – even though the former substorm is much the weaker of the two. These poleward retreats are consistent with the destruction of open flux at the tail current sheet. Application of a new analysis of the peak ion energies at the equatorward edge of the cleft/cusp/mantle dispersion seen by the DMSP satellites identifies the dayside reconnection merging gap to extend in MLT from about 9.5 to 15.5 h for most of the interval. Analysis of the boundary motion, and of the convection velocities seen near the boundary by EISCAT, allows calculation of the reconnection rate (mapped down to the ionosphere) from the flow component normal to the boundary in its own rest frame. This reconnection rate is not, in general, significantly different from zero before 06:45 UT (MLT<9.5 h) – indicating that the X line footprint expands over the EISCAT field-of-view to earlier MLT only occasionally and briefly. Between 06:45 UT and 12:45UT (9.5

Low-energy ion outflows from the ionosphere during a major polar cap expansion — evidence for equatorward motion of inverted-V structures

Data from the Dynamics Explorer 1 satellite and the EISCAT and Sondrestrom incoherent scatter radars, have allowed a study of low-energy ion outflows from the ionosphere into the magnetosphere during a rapid expansion of the polar cap. From the combined radar data, a 200kV increase in cross-cap potential is estimated. The upflowing ions show “X” signatures in the pitch angle-time spectrograms in the expanding midnight sector of the auroral oval. These signatures reveal low-energy (below about 60eV), light-ion beams sandwiched between two regions of ion conics and are associated with inverted-V electron precipitation. The lack of mass dispersion of the poleward edge of the event, despite great differences in the times of flight, reflects the equatorward expansion of the acceleration regions at velocities similar to those of the antisunward convection. In addition, a transient burst of upflow of 0+ is observed within the cap, possibly due to enhanced Joule heating during the event.

The cooling steam of the Polar Express : historical origins, properties and implications of performance capture

As the first academically rigorous interrogation of the generation of performance within the global frame of the motion capture volume, this research presents a historical contextualisation and develops and tests a set of first principles through an original series of theoretically informed, practical exercises to guide those working in the emergent space of performance capture. It contributes a new understanding of the framing of performance in The Omniscient Frame, and initiates and positions performance capture as a new and distinct interdisciplinary discourse in the fields of theatre, animation, performance studies and film.

Quantum ricochets: surface capture, release and energy loss of fast ions hitting a polar surface at grazing incidence

A diffraction mechanism is proposed for the capture, multiple bouncing and final escape of a fast ion (keV) impinging on the surface of a polarizable material at grazing incidence. Capture and escape are effected by elastic quantum diffraction consisting of the exchange of a parallel surface wave vector G= 2p/ a between the ion parallel momentum and the surface periodic potential of period a. Diffraction- assisted capture becomes possible for glancing angles F smaller than a critical value given by Fc 2- 2./ a-| Vim|/ E, where E is the kinetic energy of the ion,. = h/ Mv its de Broglie wavelength and Vim its average electronic image potential at the distance from the surface where diffraction takes place. For F< Fc, the ion can fall into a selected capture state in the quasi- continuous spectrum of its image potential and execute one or several ricochets before being released by the time reversed diffraction process. The capture, ricochet and escape are accompanied by a large, periodic energy loss of several tens of eV in the forward motion caused by the coherent emission of a giant number of quanta h. of Fuchs- Kliewer surface phonons characteristic of the polar material. An analytical calculation of the energy loss spectrum, based on the proposed diffraction process and using a model ion-phonon coupling developed earlier (Lucas et al 2013 J. Phys.: Condens. Matter 25 355009), is presented, which fully explains the experimental spectrum of Villette et al (2000 Phys. Rev. Lett. 85 3137) for Ne+ ions ricocheting on a LiF(001) surface.

Estimating Sensor Motion from Wide-Field Optical Flow on a Log-Dipolar Sensor

Log-polar image architectures, motivated by the structure of the human visual field, have long been investigated in computer vision for use in estimating motion parameters from an optical flow vector field. Practical problems with this approach have been: (i) dependence on assumed alignment of the visual and motion axes; (ii) sensitivity to occlusion form moving and stationary objects in the central visual field, where much of the numerical sensitivity is concentrated; and (iii) inaccuracy of the log-polar architecture (which is an approximation to the central 20°) for wide-field biological vision. In the present paper, we show that an algorithm based on generalization of the log-polar architecture; termed the log-dipolar sensor, provides a large improvement in performance relative to the usual log-polar sampling. Specifically, our algorithm: (i) is tolerant of large misalignmnet of the optical and motion axes; (ii) is insensitive to significant occlusion by objects of unknown motion; and (iii) represents a more correct analogy to the wide-field structure of human vision. Using the Helmholtz-Hodge decomposition to estimate the optical flow vector field on a log-dipolar sensor, we demonstrate these advantages, using synthetic optical flow maps as well as natural image sequences.

An objective climatology of the dynamical forcing of polar lows in the Nordic seas

A new objective climatology of polar lows in the Nordic (Norwegian and Barents) seas has been derived from a database of diagnostics of objectively identified cyclones spanning the period January 2000 to April 2004. There are two distinct parts to this study: the development of the objective climatology and a characterization of the dynamical forcing of the polar lows identified. Polar lows are an intense subset of polar mesocyclones. Polar mesocyclones are distinguished from other cyclones in the database as those that occur in cold air outbreaks over the open ocean. The difference between the wet-bulb potential temperature at 700 hPa and the sea surface temperature (SST) is found to be an effective discriminator between the atmospheric conditions associated with polar lows and other cyclones in the Nordic seas. A verification study shows that the objective identification method is reliable in the Nordic seas region. After demonstrating success at identifying polar lows using the above method, the dynamical forcing of the polar lows in the Nordic seas is characterized. Diagnostics of the ratio of mid-level vertical motion attributable to quasi-geostrophic forcing from upper and lower levels (U/L ratio) are used to determine the prevalence of a recently proposed category of extratropical cyclogenesis, type C, for which latent heat release is crucial to development. Thirty-one percent of the objectively identified polar low events (36 from 115) exceeded the U/L ratio of 4.0, previously identified as a threshold for type C cyclones. There is a contrast between polar lows to the north and south of the Nordic seas. In the southern Norwegian Sea, the population of polar low events is dominated by type C cyclones. These possess strong convection and weak low-level baroclinicity. Over the Barents and northern Norwegian seas, the well-known cyclogenesis types A and B dominate. These possess stronger low-level baroclinicity and weaker convection.

Direct observations of the evolution of polar cap ionization patches

Patches of ionization are common in the polar ionosphere where their motion and associated density gradients give variable disturbances to High Frequency (HF) radio communications, over-the-horizon radar location errors, and disruption and errors to satellite navigation and communication. Their formation and evolution are poorly understood, particularly under disturbed space weather conditions. We report direct observations of the full evolution of patches during a geomagnetic storm, including formation, polar cap entry, transpolar evolution, polar cap exit, and sunward return flow. Our observations show that modulation of nightside reconnection in the substorm cycle of the magnetosphere helps form the gaps between patches where steady convection would give a “tongue” of ionization (TOI).

Automated identification and tracking of polar-cap plasma patches at solar minimum

A method of automatically identifying and tracking polar-cap plasma patches, utilising data inversion and feature-tracking methods, is presented. A well-established and widely used 4-D ionospheric imaging algorithm, the Multi-Instrument Data Assimilation System (MIDAS), inverts slant total electron content (TEC) data from ground-based Global Navigation Satellite System (GNSS) receivers to produce images of the free electron distribution in the polar-cap ionosphere. These are integrated to form vertical TEC maps. A flexible feature-tracking algorithm, TRACK, previously used extensively in meteorological storm-tracking studies is used to identify and track maxima in the resulting 2-D data fields. Various criteria are used to discriminate between genuine patches and "false-positive" maxima such as the continuously moving day-side maximum, which results from the Earth's rotation rather than plasma motion. Results for a 12-month period at solar minimum, when extensive validation data are available, are presented. The method identifies 71 separate structures consistent with patch motion during this time. The limitations of solar minimum and the consequent small number of patches make climatological inferences difficult, but the feasibility of the method for patches larger than approximately 500 km in scale is demonstrated and a larger study incorporating other parts of the solar cycle is warranted. Possible further optimisation of discrimination criteria, particularly regarding the definition of a patch in terms of its plasma concentration enhancement over the surrounding background, may improve results.

The response of ionospheric convection in the polar cap to substorm activity

We report multi-instrument observations during an isolated substorm on 17 October 1989. The EISCAT radar operated in the SP-UK-POLI mode measuring ionospheric convection at latitudes 71°-78°. SAMNET and the EISCAT Magnetometer Cross provide information on the timing of substorm expansion phase onset and subsequent intensifications, as well as the location of the field aligned and ionospheric currents associated with the substorm current wedge. IMP-8 magnetic field data are also included. Evidence of a substorm growth phase is provided by the equatorward motion of a flow reversal boundary across the EISCAT radar field of view at 2130 MLT, following a southward turning of the interplanetary magnetic field (IMF). We infer that the polar cap expanded as a result of the addition of open magnetic flux to the tail lobes during this interval. The flow reversal boundary, which is a lower limit to the polar cap boundary, reached an invariant latitude equatorward of 71° by the time of the expansion phase onset. A westward electrojet, centred at 65.4°, occurred at the onset of the expansion phase. This electrojet subsequently moved poleward to a maximum of 68.1° at 2000 UT and also widened. During the expansion phase, there is evidence of bursts of plasma flow which are spatially localised at longitudes within the substorm current wedge and which occurred well poleward of the westward electrojet. We conclude that the substorm onset region in the ionosphere, defined by the westward electrojet, mapped to a part of the tail radially earthward of the boundary between open and closed magnetic flux, the “distant” neutral line. Thus the substorm was not initiated at the distant neutral line, although there is evidence that it remained active during the expansion phase. It is not obvious whether the electrojet mapped to a near-Earth neutral line, but at its most poleward, the expanded electrojet does not reach the estimated latitude of the polar cap boundary.