ALBERT

Beschreibung: The numerical simulation on ionospheric perturbations in electric field before large earthquakes Annales Geophysicae, 32, 1487-1493, 2014 Author(s): S. F. Zhao, X. M. Zhang, Z. Y. Zhao, and X. H. Shen Many observational results have shown electromagnetic abnormality in the ionosphere before large earthquakes. The theoretical simulation can help us to understand the internal mechanism of these anomalous electromagnetic signals resulted from seismic regions. In this paper, the horizontal and vertical components of electric and magnetic field at the topside ionosphere are simulated by using the full wave method that is based on an improved transfer matrix method in the lossy anisotropic horizontally stratified ionosphere. Taken account into two earthquakes with electric field perturbations recorded by the DEMETER satellite, the numerical results reveal that the propagation and penetration of ULF (ultra-low-frequency) electromagnetic waves into the ionosphere is related to the spatial distribution of electron and ion densities at different time and locations, in which the ion density has less effect than electron density on the field intensity. Compared with different frequency signals, the minimum values of electric and magnetic field excited by earthquakes can be detected by satellite in current detection capability have also been calculated, and the lower frequency wave can be detected easier.

Beschreibung: A comparison between VEGA 1, 2 and Giotto flybys of comet 1P/Halley: implications for Rosetta Annales Geophysicae, 32, 1441-1453, 2014 Author(s): M. Volwerk, K.-H. Glassmeier, M. Delva, D. Schmid, C. Koenders, I. Richter, and K. Szegö Three flybys of comet 1P/Halley, by VEGA 1, 2 and Giotto, are investigated with respect to the occurrence of mirror mode waves in the cometosheath and field line draping in the magnetic pile-up region around the nucleus. The time interval covered by these flybys is approximately 8 days, which is also the approximate length of an orbit or flyby of Rosetta around comet 67P/Churyumov–Gerasimenko. Thus any significant changes observed around Halley are changes that might occur for Rosetta during one pass of 67P/CG. It is found that the occurrence of mirror mode waves in the cometosheath is strongly influenced by the dynamical pressure of the solar wind and the outgassing rate of the comet. Field line draping happens in the magnetic pile-up region. Changes in nested draping regions (i.e. regions with different B x directions) can occur within a few days, possibly influenced by changes in the outgassing rate of the comet and thereby the conductivity of the cometary ionosphere.

Beschreibung: Instability and change detection in exponential families and generalized linear models, with a study of Atlantic tropical storms Nonlinear Processes in Geophysics, 21, 1133-1143, 2014 Author(s): Y. Lu and S. Chatterjee Exponential family statistical distributions, including the well-known normal, binomial, Poisson, and exponential distributions, are overwhelmingly used in data analysis. In the presence of covariates, an exponential family distributional assumption for the response random variables results in a generalized linear model. However, it is rarely ensured that the parameters of the assumed distributions are stable through the entire duration of the data collection process. A failure of stability leads to nonsmoothness and nonlinearity in the physical processes that result in the data. In this paper, we propose testing for stability of parameters of exponential family distributions and generalized linear models. A rejection of the hypothesis of stable parameters leads to change detection. We derive the related likelihood ratio test statistic. We compare the performance of this test statistic to the popular normal distributional assumption dependent cumulative sum (Gaussian CUSUM) statistic in change detection problems. We study Atlantic tropical storms using the techniques developed here, so to understand whether the nature of these tropical storms has remained stable over the last few decades.

Beschreibung: Magnetic curvature driven Rayleigh-Taylor instability revisited Annales Geophysicae, 29, 411-413, 2011 Author(s): O. A. Pokhotelov and O. G. Onishchenko The problem of incomplete finite ion Larmor radius (FLR) stabilization of the magnetic curvature driven Rayleigh-Taylor instability (RTI) in low beta plasma with homogeneous ion temperature is investigated. For this purpose a model hydrodynamic description of nonlinear flute waves with arbitrary spatial scales compared to the ion Larmor radius is developed. It is shown that the RTI is not stabilized by FLR effects in a plasma with cold electrons when the ratio of characteristic spatial scale of the plasma inhomogeneity to local effective radius of curvature of the magnetic field lines is larger than 1/4. The crucial role in the absence of the complete FLR stabilization plays the contribution of the compressibility of the polarization part of the ion velocity.

Beschreibung: Forecasting relativistic electron flux using dynamic multiple regression models Annales Geophysicae, 29, 415-420, 2011 Author(s): H.-L. Wei, S. A. Billings, A. Surjalal Sharma, S. Wing, R. J. Boynton, and S. N. Walker The forecast of high energy electron fluxes in the radiation belts is important because the exposure of modern spacecraft to high energy particles can result in significant damage to onboard systems. A comprehensive physical model of processes related to electron energisation that can be used for such a forecast has not yet been developed. In the present paper a systems identification approach is exploited to deduce a dynamic multiple regression model that can be used to predict the daily maximum of high energy electron fluxes at geosynchronous orbit from data. It is shown that the model developed provides reliable predictions.

Beschreibung: Airglow observations over the equatorial ionization anomaly zone in Taiwan Annales Geophysicae, 29, 749-757, 2011 Author(s): J. Y. Liu, P. K. Rajesh, I. T. Lee, and T. C. Chow Airglow imaging at mid-latitude stations often show intensity modulations associated with medium scale travelling ionospheric disturbances (MSTID), while those carried out near the equatorial regions reveal depletions caused by equatorial plasma bubbles (EPB). Two all sky cameras are used to observe plasma depletions in the 630.0 nm emission over the equatorial ionization anomaly (EIA) region, Taiwan (23° N, 121° E; 13.5° N Magnetic) during 1998–2002 and 2006–2007. The results show EPB and MSTID depletions in different solar activity conditions. Several new features of the EPB depletions such as bifurcation, secondary structure on the walls, westward tilt, etc., are discussed in this paper. Evidence of tilted depletions with secondary structures developing on the eastern wall that later evolve to appear as bifurcations, are presented for the first time. Moreover, detail investigations are carried out using International Reference Ionosphere (IRI) model as well as the electron density from Ionosonde and Global Positioning System (GPS) Occultation Experiment (GOX) onboard FORMOSAT-3/COSMIC satellite, to understand the conditions that favor the propagation of MSTID to the latitude of Taiwan.

Beschreibung: 3-D mesoscale MHD simulations of a cusp-like magnetic configuration: method and first results Annales Geophysicae, 29, 759-770, 2011 Author(s): E. Adamson and K. Nykyri We present a local mesoscale model of the magnetospheric cusp region with high resolution (up to 300 km). We discuss the construction and implementation of the initial configuration and give a detailed description of the numerical simulation. An overview of simulation results for the case of strongly northward interplanetary magnetic field (IMF) is then presented and compared with data from Cluster 2 spacecraft from 14 February 2003. Results show a cusp diamagnetic cavity (CDC) with depth normal to the magnetospheric boundary on the order of 1–2 R E and a much larger extent of ~5–9 R E tangential to the boundary, bounded by a gradual inner boundary with the magnetospheric lobe and a more distinct exterior boundary with the magnetosheath. These results are qualitatively consistent with observational data.

Beschreibung: Polar cap ion beams during periods of northward IMF: Cluster statistical results Annales Geophysicae, 29, 771-787, 2011 Author(s): R. Maggiolo, M. Echim, J. De Keyser, D. Fontaine, C. Jacquey, and I. Dandouras Above the polar caps and during prolonged periods of northward IMF, the Cluster satellites detect upward accelerated ion beams with energies up to a few keV. They are associated with converging electric field structures indicating that the acceleration is caused by a quasi-static field-aligned electric field that can extend to altitudes higher than 7 R E (Maggiolo et al., 2006; Teste et al., 2007). Using the AMDA science analysis service provided by the Centre de Données de la Physique des Plasmas, we have been able to extract about 200 events of accelerated upgoing ion beams above the polar caps from the Cluster database. Most of these observations are taken at altitudes lower than 7 R E and in the Northern Hemisphere. We investigate the statistical properties of these ion beams. We analyze their geometry, the properties of the plasma populations and of the electric field inside and around the beams, as well as their dependence on solar wind and IMF conditions. We show that ~40 % of the ion beams are collocated with a relatively hot and isotropic plasma population. The density and temperature of the isotropic population are highly variable but suggest that this plasma originates from the plasma sheet. The ion beam properties do not change significantly when the isotropic, hot background population is present. Furthermore, during one single polar cap crossing by Cluster it is possible to detect upgoing ion beams both with and without an accompanying isotropic component. The analysis of the variation of the IMF B Z component prior to the detection of the beams indicates that the delay between a northward/southward turning of IMF and the appearance/disappearance of the beams is respectively ~2 h and 20 min. The observed electrodynamic characteristics of high altitude polar cap ion beams suggest that they are closely connected to polar cap auroral arcs. We discuss the implications of these Cluster observations above the polar cap on the magnetospheric dynamics and configuration during prolonged periods of northward IMF.

Beschreibung: Rank ordering multifractal analysis of the auroral electrojet index Nonlinear Processes in Geophysics, 18, 277-285, 2011 Author(s): G. Consolini and P. De Michelis In the second half of the 90s interest grew on the complex features of the magnetospheric dynamics in response to solar wind changes. An important series of papers were published on the occurrence of chaos, turbulence and complexity. Among them, particularly interesting was the study of the bursty and fractal/multifractal character of the high latitude energy release during geomagnetic storms, which was evidenced by analyzing the features of the Auroral Electrojet (AE) indices. Recently, the multifractal features of the small time-scale increments of AE-indices have been criticized in favor of a more simple fractal behavior. This is particularly true for the scaling features of the probability density functions (PDFs) of the AE index increments. Here, after a brief review of the nature of the fractal/multifractal features of the magnetospheric response to solar wind changes, we investigate the multifractal nature of the scaling features of the AE index increments PDFs using the Rank Ordering Multifractal Analysis (ROMA) technique. The ROMA results clearly demonstrate the existence of a hierarchy of scaling indices, depending on the increment amplitude, for the data collapsing of PDFs relative to increments at different time scales. Our results confirm the previous results by Consolini et al. (1996) and the more recent results by Rypdal and Rypdal (2010).

Beschreibung: Diurnal variation of short-period (20–120 min) gravity waves in the equatorial Mesosphere and Lower Thermosphere and its relation to deep tropical convection Annales Geophysicae, 29, 623-629, 2011 Author(s): N. Venkateswara Rao, Y. Shibagaki, and T. Tsuda We study short period gravity waves (20–120 min) in the equatorial Mesosphere and Lower Thermosphere (MLT) using a Medium Frequency (MF) radar at Pameungpeuk (7.4° S, 107.4° E), Indonesia. In particular, we study local time and seasonal variation of the gravity wave variance and its relation to tropical convection. The gravity wave variance at 88 km enhances between 20:00 LT and 07:00 LT, with a peak at 02:00–03:00 LT. The enhancement is mainly observed during February–April and September–October and shows inter-annual variability. Convective activity over the same location persists from 16:00–21:00 LT with a peak activity ~18:00 LT and enhances between November–April. Time delay between the peak of convection and that of gravity wave activity ranges 1–15 h, which is consistent with theoretical calculations and previous reports based on reverse ray tracing analysis.

Beschreibung: The features of electronics structure of the multichannel scintillation module for the EMMA experiment Astrophysics and Space Sciences Transactions, 7, 171-174, 2011 Author(s): V. Volchenko, G. Volchenko, E. Akhrameev, L. Bezrukov, I. Dzaparova, T. Enqvist, L. Inzhechik, A. Izmaylov, J. Joutsenvaara, M. Khabibullin, A. Khotjantsev, P. Kuusiniemi, B. Lubsandorzhiev, O. Mineev, V. Petkov, R. Poleshuk, B. Shaibonov, J. Sarkamo, A. Shaykhiev, W. Trzaska, A. Yanin, and N. Yershov A brief description of the developed structural electric diagrams of 16-channel scintillation module for the underground EMMA experiment, the basic characteristics and parameters of the electrical diagrams of this module are presented. Multi-pixel photodiodes operating in a limited Geiger mode are used for photoreadout of the scintillator detectors in 16-channel scintillation module. The method of the automatic tuning of the photosensors gain based on the stabilization of an average counting rate of the scintillation detectors from gamma rays of a natural radioactive background is described.

Beschreibung: Observations of AGW/TID propagation across the polar cap: a case study Annales Geophysicae, 29, 1355-1363, 2011 Author(s): H. T. Cai, F. Yin, S. Y. Ma, and I. W. McCrea In this paper, we present observational evidence for the trans-polar propagation of large-scale Traveling Ionospheric Disturbances (TIDs) from their nightside source region to the dayside. On 13 February 2001, the 32 m dish of EISCAT Svalbard Radar (ESR) was directing toward the geomagnetic pole at low elevation (30°) during the interval 06:00–12:00 UT (MLT ≈ UT + 3 h), providing an excellent opportunity to monitor the ionosphere F-region over the polar cap. The TIDs were first detected by the ESR over the dayside north polar cap, propagating equatorward, and were subsequently seen by the mainland UHF radar at auroral latitudes around geomagnetic local noon. The propagation properties of the observed ionization waves suggest the presence of a moderately large-scale TIDs, propagating across the northern polar cap from the night-time auroral source during substorm conditions. Our results agree with the theoretical simulations by Balthazor and Moffett (1999) in which poleward-propagating large-scale traveling atmospheric disturbances were found to be self-consistently driven by enhancements in auroral heating.

Beschreibung: Statistical evidence for O + energization and outflow caused by wave-particle interaction in the high altitude cusp and mantle Annales Geophysicae, 29, 945-954, 2011 Author(s): M. Waara, R. Slapak, H. Nilsson, G. Stenberg, M. André, and I. A. Barghouthi We present a statistical study of the low (

Beschreibung: O + heating associated with strong wave activity in the high altitude cusp and mantle Annales Geophysicae, 29, 931-944, 2011 Author(s): R. Slapak, H. Nilsson, M. Waara, M. André, G. Stenberg, and I. A. Barghouthi We use the Cluster spacecraft to study three events with intense waves and energetic oxygen ions (O + ) in the high altitude cusp and mantle. The ion energies considered are of the order 1000 eV and higher, observed above an altitude of 8 earth radii together with high wave power at the O + gyrofrequency. We show that heating by waves can explain the observed high perpendicular energy of O + ions, using a simple gyroresonance model and 25–45% of the observed wave spectral density at the gyrofrequency. This is in contrast to a recently published study where the wave intensity was too low to explain the observed high altitude ion energies. Long lasting cases (〉10 min) of high perpendicular-to-parallel temperature ratios are sometimes associated with low wave activity, suggesting that high perpendicular-to-parallel temperature ratio is not a good indicator of local heating. Using multiple spacecraft, we show that the regions of enhanced wave activity are at least one order of magnitude larger than the gyroradius of the heated ions.

Beschreibung: Simulation of aerosol optical properties over a tropical urban site in India using a global model and its comparison with ground measurements Annales Geophysicae, 29, 955-963, 2011 Author(s): D. Goto, K. V. S. Badarinath, T. Takemura, and T. Nakajima Aerosols have great impacts on atmospheric environment, human health, and earth's climate. Therefore, information on their spatial and temporal distribution is of paramount importance. Despite numerous studies have examined the variation and trends of BC and AOD over India, only very few have focused on their spatial distribution or even correlating the observations with model simulations. In the present study, a three-dimensional aerosol transport-radiation model coupled with a general circulation model. SPRINTARS, simulated atmospheric aerosol distributions including BC and aerosol optical properties, i.e., aerosol optical thickness (AOT), Ångström Exponent (AE), and single scattering albedo (SSA). The simulated results are compared with both BC measurements by aethalometer and aerosol optical properties measured by ground-based skyradiometer and by satellite sensor, MODIS/Terra over Hyderabad, which is a tropical urban area of India, for the year 2008. The simulated AOT and AE in Hyderabad are found to be comparable to ground-based measured ones. The simulated SSA tends to be higher than the ground-based measurements. Both these comparisons of aerosol optical properties between the simulations with different emission inventories and the measurements indicate that, firstly the model uncertainties derived from aerosol emission inventory cannot explain the gaps between the simulations and the measurements and secondly the vertical transport of BC and the treatment of BC-containing particles can be the main issue in the global model to solve the gap.

Beschreibung: Energy and flux variations across thin auroral arcs Annales Geophysicae, 29, 1699-1712, 2011 Author(s): H. Dahlgren, B. Gustavsson, B. S. Lanchester, N. Ivchenko, U. Brändström, D. K. Whiter, T. Sergienko, I. Sandahl, and G. Marklund Two discrete auroral arc filaments, with widths of less than 1 km, have been analysed using multi-station, multi-monochromatic optical observations from small and medium field-of-view imagers and the EISCAT radar. The energy and flux of the precipitating electrons, volume emission rates and local electric fields in the ionosphere have been determined at high temporal (up to 30 Hz) and spatial (down to tens of metres) resolution. A new time-dependent inversion model is used to derive energy spectra from EISCAT electron density profiles. The energy and flux are also derived independently from optical emissions combined with ion-chemistry modelling, and a good agreement is found. A robust method to obtain detailed 2-D maps of the average energy and number flux of small scale aurora is presented. The arcs are stretched in the north-south direction, and the lowest energies are found on the western, leading edges of the arcs. The large ionospheric electric fields (250 mV m −1 ) found from tristatic radar measurements are evidence of strong currents associated with the region close to the optical arcs. The different data sets indicate that the arcs appear on the boundaries between regions with different average energy of diffuse precipitation, caused by pitch-angle scattering. The two thin arcs on these boundaries are found to be related to an increase in number flux (and thus increased energy flux) without an increase in energy.

Beschreibung: Storm-time total electron content and its response to penetration electric fields over South America Annales Geophysicae, 29, 1765-1778, 2011 Author(s): P. M. de Siqueira, E. R. de Paula, M. T. A. H. Muella, L. F. C. Rezende, M. A. Abdu, and W. D. Gonzalez In this work the response of the ionosphere due to the severe magnetic storm of 7–10 November 2004 is investigated by analyzing GPS Total Electron Content (TEC) maps constructed for the South America sector. In order to verify the disturbed zonal electric fields in South America during the superstorm, ionospheric vertical drift data obtained from modeling results are used in the analysis. The vertical drifts were inferred from Δ H magnetometer data (Jicamarca-Piura) following the methodology presented by Anderson et al. (2004). Also used were vertical drifts measured by the Jicamarca ISR. Data from a digisonde located at São Luís, Brazil (2.33° S, 44.2° W, dip latitude 0.25°) are presented to complement the Jicamarca equatorial data. Penetration electric fields were observed by the comparison between the equatorial vertical drifts and the Interplanetary Electric Field (IEF). The TEC maps obtained from GPS data reflect the ionospheric response over the South America low-latitude and equatorial region. They reveal unexpected plasma distributions and TEC levels during the main phase of the superstorm on 7 November, which is coincident with the local post-sunset hours. At this time an increase in the pre-reversal enhancement was expected to develop the Equatorial Ionization Anomaly (EIA) but we observed the absence of EIA. The results also reveal well known characteristics of the plasma distributions on 8, 9, and 10 November. The emphasized features are the expansion and intensification of EIA due to prompt penetration electric fields on 9 November and the inhibition of EIA during post-sunset hours on 7, 8, and 10 November. One important result is that the TEC maps provided a bi-dimensional view of the ionospheric changes offering a spatial description of the electrodynamics involved, which is an advantage over TEC measured by isolated GPS receivers.

Beschreibung: Climatology of rapid geomagnetic variations at high latitudes over two solar cycles Annales Geophysicae, 29, 1783-1792, 2011 Author(s): A. Viljanen and E. Tanskanen We investigate the characteristics of rapid geomagnetic variations at high latitudes based on the occurrence of large time derivatives of the horizontal magnetic field ( d H / dt exceeding 1 nT s −1 ). Analysis of IMAGE magnetometer data from North Europe in 1983–2010, covering more than two solar cycles, confirms and specifies several previous findings. We show that d H / dt activity is high around the midnight and early morning hours, and nearly vanishes at noon and early afternoon. This happens during all seasons, although the midnight maximum is nearly invisible during summer. As indicated by modelled ionospheric equivalent currents, large d H / dt values occur predominantly during westward ionospheric electrojets. Before and around midnight, d H / dt tends to be north-south oriented, whereas in the morning hours, its direction is more west-east directed. d H / dt tends to be more strictly north-south oriented during winter than other seasons. The seasonal occurrence of large d H / dt values is similar to the variation of the maximum amplitude of westward equivalent currents. The yearly fraction of east-west directed large d H / dt vectors at the Kilpisjärvi station (MLAT 65.88) varies from 31 to 47 % without any clear correlation with the general geomagnetic activity nor with the yearly averages of solar wind parameters.

Beschreibung: Is thermospheric long-term cooling due to CO 2 or O 3 ? Annales Geophysicae, 29, 1779-1782, 2011 Author(s): P. L. Walsh and W. L. Oliver While greenhouse gases trap heat emanating from the Earth and thereby heat the surface atmosphere, they act as emitters in the high atmosphere and cool the air there. In 1989 Roble and Dickinson (1989) estimated the cooling that would occur in the thermosphere (250–500 km altitude) due to a doubling of greenhouse gas densities. Ever since, long-term data bases have been scoured for evidence of this thermospheric "global cooling." Here we show evidence that the thermosphere did indeed cool over the period 1966–1987, but the data suggest that the cooling accelerated at a "breakpoint year" around 1979 to a rate far larger than may be attributed to greenhouse cooling. This 1979 breakpoint year appears to coincide with a breakpoint year in ozone (O 3 ) column density. Further, the cooling was confined largely to the daytime thermosphere while the nighttime showed only a small trend. These results suggest, first, that the greenhouse cooling of the thermosphere may well not be detectable with current data sets and, second, that the long-term cooling that is clearly seen may be due largely to O 3 depletion.

Beschreibung: Characteristics of the terrestrial field-aligned current system Annales Geophysicae, 29, 1713-1729, 2011 Author(s): J. W. Gjerloev, S. Ohtani, T. Iijima, B. Anderson, J. Slavin, and G. Le We present the first ever comprehensive statistical study of the spatiotemporal characteristics of field-aligned currents in the terrestrial magnetosphere-ionosphere system using multi point measurements. We determine how the FAC density, variability and scale size are coupled. The three ST 5 satellites were in a pearls-on-a-string formation making measurements of the magnetic field with variable inter-spacecraft separations ranging from a few seconds to about 10 min. More than 4700 sets of satellite passes are analyzed using a robust correlation analysis aimed at determining the variability of the FAC system as a function of scale size and satellite spacing. We find significant differences between the FAC characteristics on the dayside and on the nightside in terms of dynamics of the current systems. On the dayside the FAC characteristics are found to be independent of IMF B z and geomagnetic activity while the nightside indicates increased variability during disturbed conditions. The boundary separating highly and poorly correlated FACs can be fitted by a linear line for satellite separations shorter than 60 s (dayside) and 160 s (nightside). We interpret this as the dayside and nightside magnetospheric reconfiguration times respectively. For times exceeding this the FAC characteristics are suggested to be controlled by the solar wind (dayside) and plasma sheet (nightside) dynamics. Finally, the characteristics of FAC system with scale sizes larger than ~200 km (at ionospheric altitude) appear to be stable and repeatable on time scales of the order of a minute (i.e. comparable to the low-altitude orbiting satellite's traverse time across the auroral belt). In this sense, our results effectively validate the Iijima and Potemra (1978) assumption that on average the large-scale currents with scale sizes of the Region1 and Region2 are quasi-persistently significant in the transport of energy and momentum between the magnetosphere and the ionosphere.

Beschreibung: Preface Large amplitude internal waves in the coastal ocean Nonlinear Processes in Geophysics, 18, 653-655, 2011 Author(s): R. Grimshaw, K. Helfrich, and A. Scotti No Abstract available.

Beschreibung: Whistler-Langmuir oscillitons and their relation to auroral hiss Annales Geophysicae, 29, 1739-1753, 2011 Author(s): K. Sauer and R. D. Sydora A new type of oscilliton (soliton with superimposed spatial oscillations) is described which arises in plasmas if the electron cyclotron frequency Ω e is larger than the electron plasma frequency ω e , which is a typical situation for auroral regions in planetary magnetospheres. Both high-frequency modes of concern, the Langmuir and the whistler wave, are completely decoupled if they propagate parallel to the magnetic field. However, for oblique propagation two mixed modes are created with longitudinal and transverse electric field components. The lower mode (in the literature commonly called the whistler mode, e.g. Gurnett et al., 1983) has whistler wave characteristics at small wave numbers and asymptotically transforms into the Langmuir mode. As a consequence of the coupling between these two modes, with different phase velocity dependence, a maximum in phase velocity appears at finite wave number. The occurrence of such a particular point where phase and group velocity coincide creates the condition for the existence of a new type of oscillating nonlinear stationary structure, which we call the whistler-Langmuir (WL) oscilliton. After determining, by means of stationary dispersion theory, the parameter regime in which WL oscillitons exist, their spatial profiles are calculated within the framework of cold (non-relativistic) fluid theory. Particle-in-cell (PIC) simulations are used to demonstrate the formation of WL oscillitons which seem to play an important role in understanding electron beam-excited plasma radiation that is observed as auroral hiss in planetary magnetospheres far away from the source region.

Beschreibung: Density and temperature of energetic electrons in the Earth's magnetotail derived from high-latitude GPS observations during the declining phase of the solar cycle Annales Geophysicae, 29, 1755-1763, 2011 Author(s): M. H. Denton and T. E. Cayton Single relativistic-Maxwellian fits are made to high-latitude GPS-satellite observations of energetic electrons for the period January 2006–November 2010; a constellation of 12 GPS space vehicles provides the observations. The derived fit parameters (for energies ~0.1–1.0 MeV), in combination with field-line mapping on the nightside of the magnetosphere, provide a survey of the energetic electron density and temperature distribution in the magnetotail between McIlwain L-values of L =6 and L =22. Analysis reveals the characteristics of the density-temperature distribution of energetic electrons and its variation as a function of solar wind speed and the Kp index. The density-temperature characteristics of the magnetotail energetic electrons are very similar to those found in the outer electron radiation belt as measured at geosynchronous orbit. The energetic electron density in the magnetotail is much greater during increased geomagnetic activity and during fast solar wind. The total electron density in the magnetotail is found to be strongly correlated with solar wind speed and is at least a factor of two greater for high-speed solar wind ( V SW =500–1000 km s −1 ) compared to low-speed solar wind ( V SW =100–400 km s −1 ). These results have important implications for understanding (a) how the solar wind may modulate entry into the magnetosphere during fast and slow solar wind, and (b) if the magnetotail is a source or a sink for the outer electron radiation belt.

Beschreibung: Solar activity and the associated ground level enhancements of solar cosmic rays during solar cycle 23 Astrophysics and Space Sciences Transactions, 7, 439-443, 2011 Author(s): M. Andriopoulou, H. Mavromichalaki, P. Preka-Papadema, C. Plainaki, A. Belov, and E. Eroshenko The solar cycle 23 seems to be of great interest for the researchers due to many peculiarities. A study of the parameters of the sixteen ground level enhancements recorded during the approximately 12-year period of it (1996–2008) together with the associated solar activity, including the main properties of the solar flares, the coronal mass ejections and the radio bursts has been realized, in an effort to understand the connection of these events. All studied cases seem to be connected with very intense flares of long duration, having a mean importance value of X 5.9 and a mean duration of 164.5 min, with either halo or partial halo coronal mass ejections with a mean linear velocity of 1876 km/sec, as well as with intense radio bursts. It is also noticed that the ground level enhancements of the 23rd solar cycle occurred after the onset time of the associated solar X -ray flares with a mean time delay of about 38 min, very useful result for their monitoring and prediction.

Beschreibung: Lower hybrid resonances stimulated by the four CLUSTER relaxation sounders deep inside the plasmasphere: observations and inferred plasma characteristics Annales Geophysicae, 29, 2003-2018, 2011 Author(s): S. Kougblénou, G. Lointier, P. M. E. Décréau, J.-G. Trotignon, J.-L. Rauch, X. Vallières, P. Canu, A. Masson, and J. Pickett The frequency range of the WHISPER relaxation sounder instrument on board CLUSTER, 4–80 kHz, has been chosen so as to encompass the electron gyro-frequency, F ce , and the electron plasma frequency, F p , in most regions to be explored. Measurement of those frequencies, which are triggered as resonances by the sounder, provides a direct estimation of in situ fundamental plasma characteristics: electron density and magnetic field intensity. In the late mission phase, CLUSTER penetrated regions deep inside the plasmasphere where F ce and F p are much higher than the upper frequency of the sounder's range. However, they are of the right order of magnitude as to place the lower hybrid frequency, F lh , in the 4–15 kHz band. This characteristic frequency, placed at a resonance of the medium, is triggered by the sounder's transmitter and shows up as an isolated peak in the received spectrum, not present in spectra of naturally occuring VLF waves. This paper illustrates, from analysis of case events, how measured F lh values give access to a plasma diagnostic novel of its kind. CLUSTER, travelling along its orbit, encounters favourable conditions where F ce is increasing and F p decreasing, such that F ce / F p increases from values below unity to values above unity. Measured F lh values thus give access, in turn, to the effective mass, M eff , indicative of plasma ion composition, and to the core plasmasphere electron density value, a parameter difficult to measure. The analysed case events indicate that the estimated quantities ( M eff in the 1.0–1.4 range, N e in the 5 × 10 2 –10 4 cm −3 range) are varying with external factors (altitude, L value, geomagnetic activity) in a plausible way. Although covering only a restricted region (mid-latitude, low altitude inner plasmasphere), these measurements are available, since late 2009, for all CLUSTER perigee passes not affected by eclipses (on average, roughly a third of a total of ~200 passes per year) and offer multipoint observations previously unavailable in this region.

Beschreibung: An explanation of the Pioneer anomaly involving accelerated atomic clocks Astrophysics and Space Sciences Transactions, 7, 487-494, 2011 Author(s): K. Wilhelm and B. N. Dwivedi The Pioneer anomaly stands for unexplained frequency shifts of the Doppler radio-tracking signals received at the ground stations from the Pioneer 10 and 11 spacecraft in disagreement with expectations based on model calculations. We consider here observations of Pioneer 10 at heliocentric distances between 40 ua and 70.5 ua over a time interval of 11.55 years from 1987 to 1998. The anomaly has been interpreted in the literature either as a Doppler shift caused by an apparent spacecraft deceleration not accounted for by known effects, or as an unexpected clock acceleration of the frequency standards at the ground stations. The reasons for the anomalous behaviour are not understood in both cases. Based on a gravitational impact model – requiring a secular mass increase of all massive bodies – a solution is proposed that implies a clock acceleration with a value close to that of the Hubble constant.

Beschreibung: Size distribution and structure of Barchan dune fields Nonlinear Processes in Geophysics, 18, 455-467, 2011 Author(s): O. Durán, V. Schwämmle, P. G. Lind, and H. J. Herrmann Barchans are isolated mobile dunes often organized in large dune fields. Dune fields seem to present a characteristic dune size and spacing, which suggests a cooperative behavior based on dune interaction. In Duran et al. (2009), we propose that the redistribution of sand by collisions between dunes is a key element for the stability and size selection of barchan dune fields. This approach was based on a mean-field model ignoring the spatial distribution of dune fields. Here, we present a simplified dune field model that includes the spatial evolution of individual dunes as well as their interaction through sand exchange and binary collisions. As a result, the dune field evolves towards a steady state that depends on the boundary conditions. Comparing our results with measurements of Moroccan dune fields, we find that the simulated fields have the same dune size distribution as in real fields but fail to reproduce their homogeneity along the wind direction.

Beschreibung: Observations of electron vorticity in the inner plasma sheet Annales Geophysicae, 29, 1517-1527, 2011 Author(s): C. Gurgiolo, M. L. Goldstein, A. F. Viñas, W. H. Matthaeus, and A. N. Fazakerley From a limited number of observations it appears that vorticity is a common feature in the inner plasma sheet. With the four Cluster spacecraft and the four PEACE instruments positioned in a tetrahedral configuration, for the first time it is possible to directly estimate the electron fluid vorticity in a space plasma. We show examples of electron fluid vorticity from multiple plasma sheet crossings. These include three time periods when Cluster passed through a reconnection ion diffusion region. Enhancements in vorticity are seen in association with each crossing of the ion diffusion region.

Beschreibung: Basic principles and aims of model order reduction in compliant mechanisms Mechanical Sciences, 2, 197-204, 2011 Author(s): M. Rösner and R. Lammering Model order reduction appears to be beneficial for the synthesis and simulation of compliant mechanisms due to computational costs. Model order reduction is an established method in many technical fields for the approximation of large-scale linear time-invariant dynamical systems described by ordinary differential equations. Based on system theory, underlying representations of the dynamical system are introduced from which the general reduced order model is derived by projection. During the last years, numerous new procedures were published and investigated appropriate to simulation, optimization and control. Singular value decomposition, condensation-based and Krylov subspace methods representing three order reduction methods are reviewed and their advantages and disadvantages are outlined in this paper. The convenience of applying model order reduction in compliant mechanisms is quoted. Moreover, the requested attributes for order reduction as a future research direction meeting the characteristics of compliant mechanisms are commented.

Beschreibung: Combining 2-m temperature nowcasting and short range ensemble forecasting Nonlinear Processes in Geophysics, 18, 903-910, 2011 Author(s): A. Kann, T. Haiden, and C. Wittmann During recent years, numerical ensemble prediction systems have become an important tool for estimating the uncertainties of dynamical and physical processes as represented in numerical weather models. The latest generation of limited area ensemble prediction systems (LAM-EPSs) allows for probabilistic forecasts at high resolution in both space and time. However, these systems still suffer from systematic deficiencies. Especially for nowcasting (0–6 h) applications the ensemble spread is smaller than the actual forecast error. This paper tries to generate probabilistic short range 2-m temperature forecasts by combining a state-of-the-art nowcasting method and a limited area ensemble system, and compares the results with statistical methods. The Integrated Nowcasting Through Comprehensive Analysis (INCA) system, which has been in operation at the Central Institute for Meteorology and Geodynamics (ZAMG) since 2006 (Haiden et al., 2011), provides short range deterministic forecasts at high temporal (15 min–60 min) and spatial (1 km) resolution. An INCA Ensemble (INCA-EPS) of 2-m temperature forecasts is constructed by applying a dynamical approach, a statistical approach, and a combined dynamic-statistical method. The dynamical method takes uncertainty information (i.e. ensemble variance) from the operational limited area ensemble system ALADIN-LAEF (Aire Limitée Adaptation Dynamique Développement InterNational Limited Area Ensemble Forecasting) which is running operationally at ZAMG (Wang et al., 2011). The purely statistical method assumes a well-calibrated spread-skill relation and applies ensemble spread according to the skill of the INCA forecast of the most recent past. The combined dynamic-statistical approach adapts the ensemble variance gained from ALADIN-LAEF with non-homogeneous Gaussian regression (NGR) which yields a statistical \mbox{correction} of the first and second moment (mean bias and dispersion) for Gaussian distributed continuous variables. Validation results indicate that all three methods produce sharp and reliable probabilistic 2-m temperature forecasts. However, the statistical and combined dynamic-statistical methods slightly outperform the pure dynamical approach, mainly due to the under-dispersive behavior of ALADIN-LAEF outside the nowcasting range. The training length does not have a pronounced impact on forecast skill, but a spread re-scaling improves the forecast skill substantially. Refinements of the statistical methods yield a slight further improvement.

Beschreibung: First detection of the O III 495.8911 and 500.6843 nm lines in the Earth's upper atmosphere Annales Geophysicae, 29, 2235-2238, 2011 Author(s): O. Witasse, T. Slanger, and R. Thissen We report the first detection of two emission lines of the atomic oxygen doubly-charged ion at 495.8911 and 500.6843 nm in the terrestrial upper atmosphere. They correspond to the transitions 1 D 2 - 3 P 1 and 1 D 2 - 3 P 2 of the O ++ ion, respectively. The measurements were performed on 30 October 2003 during the "Halloween" storms, with the Ultraviolet and Visual Echelle Spectrograph (UVES) mounted on the Very Large Telescope (VLT) in Chile. The intensities of these emissions are ~70 mRayleigh, and ~260 mRayleigh, respectively. These emissions constitute a new diagnostic of the state of the ionosphere.

Beschreibung: Solar wind and geomagnetism: toward a standard classification of geomagnetic activity from 1868 to 2009 Annales Geophysicae, 30, 421-426, 2012 Author(s): J. L. Zerbo, C. Amory Mazaudier, F. Ouattara, and J. D. Richardson We examined solar activity with a large series of geomagnetic data from 1868 to 2009. We have revisited the geomagnetic activity classification scheme of Legrand and Simon (1989) and improve their scheme by lowering the minimum Aa index value for shock and recurrent activity from 40 to 20 nT. This improved scheme allows us to clearly classify about 80% of the geomagnetic activity in this time period instead of only 60% for the previous Legrand and Simon classification.

Beschreibung: A statistical study of the performance of the Hakamada-Akasofu-Fry version 2 numerical model in predicting solar shock arrival times at Earth during different phases of solar cycle 23 Annales Geophysicae, 30, 405-419, 2012 Author(s): S. M. P. McKenna-Lawlor, C. D. Fry, M. Dryer, D. Heynderickx, K. Kecskemety, K. Kudela, and J. Balaz The performance of the Hakamada Akasofu-Fry, version 2 (HAFv.2) numerical model, which provides predictions of solar shock arrival times at Earth, was subjected to a statistical study to investigate those solar/interplanetary circumstances under which the model performed well/poorly during key phases (rise/maximum/decay) of solar cycle 23. In addition to analyzing elements of the overall data set (584 selected events) associated with particular cycle phases, subsets were formed such that those events making up a particular sub-set showed common characteristics. The statistical significance of the results obtained using the various sets/subsets was generally very low and these results were not significant as compared with the hit by chance rate (50%). This implies a low level of confidence in the predictions of the model with no compelling result encouraging its use. However, the data suggested that the success rates of HAFv.2 were higher when the background solar wind speed at the time of shock initiation was relatively fast. Thus, in scenarios where the background solar wind speed is elevated and the calculated success rate significantly exceeds the rate by chance, the forecasts could provide potential value to the customer. With the composite statistics available for solar cycle 23, the calculated success rate at high solar wind speed, although clearly above 50%, was indicative rather than conclusive. The RMS error estimated for shock arrival times for every cycle phase and for the composite sample was in each case significantly better than would be expected for a random data set. Also, the parameter "Probability of Detection, yes" (PODy) which presents the Proportion of Yes observations that were correctly forecast (i.e. the ratio between the shocks correctly predicted and all the shocks observed), yielded values for the rise/maximum/decay phases of the cycle and using the composite sample of 0.85, 0.64, 0.79 and 0.77, respectively. The statistical results obtained through detailed analysis of the available data provided insights into how changing circumstances on the Sun and in interplanetary space can affect the performance of the model. Since shock arrival predictions are widely utilized in making commercially significant decisions re. protecting space assets, the present detailed archival studies can be useful in future operational decision making during solar cycle 24. It would be of added value in this context to use Briggs-Rupert methodology to estimate the cost to an operator of acting on an incorrect forecast.

Beschreibung: Spatial patterns of linear and nonparametric long-term trends in Baltic sea-level variability Nonlinear Processes in Geophysics, 19, 95-111, 2012 Author(s): R. V. Donner, R. Ehrcke, S. M. Barbosa, J. Wagner, J. F. Donges, and J. Kurths The study of long-term trends in tide gauge data is important for understanding the present and future risk of changes in sea-level variability for coastal zones, particularly with respect to the ongoing debate on climate change impacts. Traditionally, most corresponding analyses have exclusively focused on trends in mean sea-level. However, such studies are not able to provide sufficient information about changes in the full probability distribution (especially in the more extreme quantiles). As an alternative, in this paper we apply quantile regression (QR) for studying changes in arbitrary quantiles of sea-level variability. For this purpose, we chose two different QR approaches and discuss the advantages and disadvantages of different settings. In particular, traditional linear QR poses very restrictive assumptions that are often not met in reality. For monthly data from 47 tide gauges from along the Baltic Sea coast, the spatial patterns of quantile trends obtained in linear and nonparametric (spline-based) frameworks display marked differences, which need to be understood in order to fully assess the impact of future changes in sea-level variability on coastal areas. In general, QR demonstrates that the general variability of Baltic sea-level has increased over the last decades. Linear quantile trends estimated for sliding windows in time reveal a wide-spread acceleration of trends in the median, but only localised changes in the rates of changes in the lower and upper quantiles.

Beschreibung: On closure parameter estimation in chaotic systems Nonlinear Processes in Geophysics, 19, 127-143, 2012 Author(s): J. Hakkarainen, A. Ilin, A. Solonen, M. Laine, H. Haario, J. Tamminen, E. Oja, and H. Järvinen Many dynamical models, such as numerical weather prediction and climate models, contain so called closure parameters. These parameters usually appear in physical parameterizations of sub-grid scale processes, and they act as "tuning handles" of the models. Currently, the values of these parameters are specified mostly manually, but the increasing complexity of the models calls for more algorithmic ways to perform the tuning. Traditionally, parameters of dynamical systems are estimated by directly comparing the model simulations to observed data using, for instance, a least squares approach. However, if the models are chaotic, the classical approach can be ineffective, since small errors in the initial conditions can lead to large, unpredictable deviations from the observations. In this paper, we study numerical methods available for estimating closure parameters in chaotic models. We discuss three techniques: off-line likelihood calculations using filtering methods, the state augmentation method, and the approach that utilizes summary statistics from long model simulations. The properties of the methods are studied using a modified version of the Lorenz 95 system, where the effect of fast variables are described using a simple parameterization.

Beschreibung: Circulation changes in the winter lower atmosphere and long-lasting solar/geomagnetic activity Annales Geophysicae, 30, 1719-1726, 2012 Author(s): J. Bochníček, H. Davídkovová, P. Hejda, and R. Huth The paper describes the association between high long-lasting solar/geomagnetic activity and geopotential height (GPH) changes in the winter lower atmosphere, based on their development in the Northern Hemisphere in the winter periods (December–March) of 1950–1969 and 1970–2002. Solar/geomagnetic activity is characterised by the 60-day mean of the sunspot number R/by the 60-day mean of the daily sum of the Kp index. The GPH distributions in the lower atmosphere are described by 60-day anomalies from their long-term daily average at 20 hPa/850 hPa. The data have been adopted from the NCEP/NCAR reanalysis. The 60-day mean values of solar/geomagnetic activity and GPH anomalies were calculated in five-day steps over the whole winter period. The analysis was carried out using composite maps which represent their distribution of the GPH anomalies during high solar activity ( R ≥ 100) and high geomagnetic activity (ΣKp ≥ 20). Analysis has shown that the distribution of GPH anomalies depends on solar activity, geomagnetic activity and the phase of winter period (early or late winter). The nature of this relationship then depends on the time interval involved, i.e. 1950–1969 or 1970–2002. Positive anomalies in the polar stratosphere (20 hPa) were detected during the whole winter periods of the years 1950–1969. Significant anomalies were detected in the lower troposphere (850 hPa) during the second half of the winter period. The distribution of GPH anomalies on the maps compiled with regard to solar activity was similar to the distribution on maps compiled with regard to geomagnetic activity. In the interval 1970–2002, significant negative GPH anomalies were detected in the stratosphere at high latitudes, and positive anomalies were detected in the region of low latitudes. The distribution of GPH anomalies in the lower troposphere was substantially affected by situations in which, together with high solar activity, also high geomagnetic activity occurred.

Beschreibung: Comprehensive analysis of tornado statistics in comparison to earthquakes: intensity and temporal behaviour Nonlinear Processes in Geophysics, 20, 47-57, 2013 Author(s): L. Schielicke and P. Névir Tornadoes and earthquakes are characterised by a high variability in their properties concerning intensity, geometric properties and temporal behaviour. Earthquakes are known for power-law behaviour in their intensity (Gutenberg–Richter law) and temporal statistics (e.g. Omori law and interevent waiting times). The observed similarity of high variability of these two phenomena motivated us to compare the statistical behaviour of tornadoes using seismological methods and quest for power-law behaviour. In general, the statistics of tornadoes show power-law behaviour partly coextensive with characteristic scales when the temporal resolution is high (10 to 60 min). These characteristic scales match with the typical diurnal behaviour of tornadoes, which is characterised by a maximum of tornado occurrences in the late afternoon hours. Furthermore, the distributions support the observation that tornadoes cluster in time. Finally, we shortly discuss a possible similar underlying structure composed of heterogeneous, coupled, interactive threshold oscillators that possibly explains the observed behaviour.

Beschreibung: Quantifying soil complexity using network models of soil porous structure Nonlinear Processes in Geophysics, 20, 41-45, 2013 Author(s): M. Samec, A. Santiago, J. P. Cárdenas, R. M. Benito, A. M. Tarquis, S. J. Mooney, and D. Korošak This paper describes an investigation into the properties of spatially embedded complex networks representing the porous architecture of soil systems. We suggest an approach to quantify the complexity of soil pore structure based on the node-node link correlation properties of the networks. We show that the complexity depends on the strength of spatial embedding of the network and that this is related to the transition from a non-compact to compact phase of the network.

Beschreibung: Analysis of asymmetries in propagating mode-2 waves Nonlinear Processes in Geophysics, 20, 59-69, 2013 Author(s): J. Olsthoorn, A. Baglaenko, and M. Stastna Using numerical simulations performed with a pseudo-spectral incompressible Navier–Stokes solver, we describe the asymmetries that arise in the recirculating core of mode-2 internal, solitary-like waves. The waves are generated in a manner consistent with many laboratory studies, namely via the collapse of a region of mixed fluid. Analysis of the simulations reveals that asymmetries across both the wave crest and the pycnocline centre develop in the spatial distribution of density, kinetic energy and a passive tracer transported by the mode-2 waves. The simulations are extended to three-dimensions to allow for the formation of spanwise instabilities. We find that three-dimensionalization modifies the structure and energetics of the core, but that the majority of the results obtained from two dimensional simulations remain valid. Taken together, our simulations demonstrate that the cores of solitary-like mode-2 waves are different then their counterparts for mode-1 waves and that their accurate characterization on both lab and field scales should account for the core asymmetry across the pycnocline centre.

Beschreibung: Earthquake dynamic response of large flexible multibody systems Mechanical Sciences, 4, 131-137, 2013 Author(s): E. V. Zahariev In the paper dynamics of large flexible structures imposed on earthquakes and high amplitude vibrations is regarded. Precise dynamic equations of flexible systems are the basis for reliable motion simulation and analysis of loading of the design scheme elements. Generalized Newton–Euler dynamic equations for rigid and flexible bodies are applied. The basement compulsory motion realized because of earthquake or wave propagation is presented in the dynamic equations as reonomic constraints. The dynamic equations, algebraic equations and reonomic constraints compile a system of differential algebraic equations which are transformed to a system of ordinary differential equations with respect to the generalized coordinates and the reactions due to the reonomic constraints. Examples of large flexible structures and wind power generator dynamic analysis are presented.

Beschreibung: Support vector machines for TEC seismo-ionospheric anomalies detection Annales Geophysicae, 31, 173-186, 2013 Author(s): M. Akhoondzadeh Using time series prediction methods, it is possible to pursue the behaviors of earthquake precursors in the future and to announce early warnings when the differences between the predicted value and the observed value exceed the predefined threshold value. Support Vector Machines (SVMs) are widely used due to their many advantages for classification and regression tasks. This study is concerned with investigating the Total Electron Content (TEC) time series by using a SVM to detect seismo-ionospheric anomalous variations induced by the three powerful earthquakes of Tohoku (11 March 2011), Haiti (12 January 2010) and Samoa (29 September 2009). The duration of TEC time series dataset is 49, 46 and 71 days, for Tohoku, Haiti and Samoa earthquakes, respectively, with each at time resolution of 2 h. In the case of Tohoku earthquake, the results show that the difference between the predicted value obtained from the SVM method and the observed value reaches the maximum value (i.e., 129.31 TECU) at earthquake time in a period of high geomagnetic activities. The SVM method detected a considerable number of anomalous occurrences 1 and 2 days prior to the Haiti earthquake and also 1 and 5 days before the Samoa earthquake in a period of low geomagnetic activities. In order to show that the method is acting sensibly with regard to the results extracted during nonevent and event TEC data, i.e., to perform some null-hypothesis tests in which the methods would also be calibrated, the same period of data from the previous year of the Samoa earthquake date has been taken into the account. Further to this, in this study, the detected TEC anomalies using the SVM method were compared to the previous results (Akhoondzadeh and Saradjian, 2011; Akhoondzadeh, 2012) obtained from the mean, median, wavelet and Kalman filter methods. The SVM detected anomalies are similar to those detected using the previous methods. It can be concluded that SVM can be a suitable learning method to detect the novelty changes of a nonlinear time series such as variations of earthquake precursors.

Beschreibung: Analysis of servo-constraint problems for underactuated multibody systems Mechanical Sciences, 4, 113-129, 2013 Author(s): R. Seifried and W. Blajer Underactuated multibody systems have fewer control inputs than degrees of freedom. In trajectory tracking control of such systems an accurate and efficient feedforward control is often necessary. For multibody systems feedforward control by model inversion can be designed using servo-constraints. So far servo-constraints have been mostly applied to differentially flat underactuated mechanical systems. Differentially flat systems can be inverted purely by algebraic manipulations and using a finite number of differentiations of the desired output trajectory. However, such algebraic solutions are often hard to find and therefore the servo-constraint approach provides an efficient and practical solution method. Recently first results on servo-constraint problems of non-flat underactuated multibody systems have been reported. Hereby additional dynamics arise, so-called internal dynamics, yielding a dynamical system as inverse model. In this paper the servo-constraint problem is analyzed for both, differentially flat and non-flat systems. Different arising important phenomena are demonstrated using two illustrative examples. Also strategies for the numerical solution of servo-constraint problems are discussed.

Beschreibung: Ion motion in the current sheet with sheared magnetic field – Part 1: Quasi-adiabatic theory Nonlinear Processes in Geophysics, 20, 163-178, 2013 Author(s): A. V. Artemyev, A. I. Neishtadt, and L. M. Zelenyi We present a theory of trapped ion motion in the magnetotail current sheet with a constant dawn–dusk component of the magnetic field. Particle trajectories are described analytically using the quasi-adiabatic invariant corresponding to averaging of fast oscillations around the tangential component of the magnetic field. We consider particle dynamics in the quasi-adiabatic approximation and demonstrate that the principal role is played by large (so called geometrical) jumps of the quasi-adiabatic invariant. These jumps appear due to the current sheet asymmetry related to the presence of the dawn–dusk magnetic field. The analytical description is compared with results of numerical integration. We show that there are four possible regimes of particle motion. Each regime is characterized by certain ranges of values of the dawn–dusk magnetic field and particle energy. We find the critical value of the dawn–dusk magnetic field, where jumps of the quasi-adiabatic invariant vanish.

Beschreibung: Magnetosheath dynamic pressure enhancements: occurrence and typical properties Annales Geophysicae, 31, 319-331, 2013 Author(s): M. O. Archer and T. S. Horbury The first comprehensive statistical study of large-amplitude (〉 100%) transient enhancements of the magnetosheath dynamic pressure reveals events of up to ~ 15 times the ambient dynamic pressure with durations up to 3 min and an average duration of around 30 s, predominantly downstream of the quasi-parallel shock. The dynamic pressure transients are most often dominated by velocity increases along with a small fractional increase in the density, though the velocity is generally only deflected by a few degrees. Superposed wavelet transforms of the magnetic field show that, whilst most enhancements exhibit changes in the magnetosheath magnetic field, the majority are not associated with changes in the Interplanetary Magnetic Field (IMF). However, there is a minority of enhancements that do appear to be associated with solar wind discontinuities which cannot be explained simply by random events. In general, it is found that during periods of magnetosheath dynamic pressure enhancements the IMF is steadier than usual. This suggests that a stable foreshock and hence foreshock structures or processes may be important in the generation of the majority of magnetosheath dynamic pressure enhancements.

Beschreibung: Conditioning model output statistics of regional climate model precipitation on circulation patterns Nonlinear Processes in Geophysics, 19, 623-633, 2012 Author(s): F. Wetterhall, F. Pappenberger, Y. He, J. Freer, and H. L. Cloke Dynamical downscaling of Global Climate Models (GCMs) through regional climate models (RCMs) potentially improves the usability of the output for hydrological impact studies. However, a further downscaling or interpolation of precipitation from RCMs is often needed to match the precipitation characteristics at the local scale. This study analysed three Model Output Statistics (MOS) techniques to adjust RCM precipitation; (1) a simple direct method (DM), (2) quantile-quantile mapping (QM) and (3) a distribution-based scaling (DBS) approach. The modelled precipitation was daily means from 16 RCMs driven by ERA40 reanalysis data over the 1961–2000 provided by the ENSEMBLES (ENSEMBLE-based Predictions of Climate Changes and their Impacts) project over a small catchment located in the Midlands, UK. All methods were conditioned on the entire time series, separate months and using an objective classification of Lamb's weather types. The performance of the MOS techniques were assessed regarding temporal and spatial characteristics of the precipitation fields, as well as modelled runoff using the HBV rainfall-runoff model. The results indicate that the DBS conditioned on classification patterns performed better than the other methods, however an ensemble approach in terms of both climate models and downscaling methods is recommended to account for uncertainties in the MOS methods.

Beschreibung: New type of ensemble of quasi-periodic, long-lasting VLF emissions at the auroral zone Annales Geophysicae, 30, 1655-1660, 2012 Author(s): J. Manninen, N. G. Kleimenova, and O. V. Kozyreva A new type of the series of quasi-periodic (QP) very low frequency (VLF) emissions in frequency range of 1–5 kHz, and not associated with geomagnetic pulsations, has been discovered at auroral latitudes ( L = 5.3) during the Finnish VLF campaign (held in December 2011). At least five unusually spectacular events, each with a duration of several hours, have been observed during the night under conditions of quiet geomagnetic activity (Kp = 0–1), although QPs usually occur during the daytime. Contrary to the QP emissions typically occurring during the day, the spectral structure of these QP events represented an extended, complicated sequence of repeated discrete rising VLF signals. Their duration was about 2–3 min each, with the repetition periods ranging from ~1 min to ~10 min. Two such nighttime non-typical events are reported in this paper. The fine structure of the separated QP elements may represent a mixture of the different frequency band signals, which seem to have independent origins. It was found that the periodic signals with lower frequency appear to trigger the strong dispersive upper frequency signals. The temporal dynamics of the spectral structure of the QPs studied were significantly controlled by some disturbances in the solar wind and interplanetary magnetic field (IMF). This finding is very important for future theoretical investigations because the generation mechanism of this new type of QP emissions is not yet understood.

Beschreibung: Spatial characteristics of wave-like structures in diffuse aurora obtained using optical observations Annales Geophysicae, 30, 1693-1701, 2012 Author(s): K. Axelsson, T. Sergienko, H. Nilsson, U. Brändström, Y. Ebihara, K. Asamura, and M. Hirahara We present the results of a statistical study using optical images from ALIS (Auroral Large Imaging System) to investigate the spatial and temporal variations of structures in diffuse aurora. Analysis of conjugate Reimei data shows that such fine structures are a result of modulation of high-energy precipitating electrons. Pitch angle diffusion into the loss cone due to interaction of whistler mode waves with plasma sheet electrons is the most feasible mechanism leading to high-energy electron precipitation. This suggests that the fine structure is an indication of modulations of the efficiency of the wave–particle interaction. The scale sizes and variations of these structures, mapped to the magnetosphere, can give us information about the characteristics of the modulating wave activity. We found the scale size of the auroral stripes and the spacing between them to be on average 13–14 km, which corresponds to 3–4 ion gyro radii for protons with an energy of 7 keV. The structures move southward with a speed close to zero in the plasma convection frame.

Beschreibung: A typical wave wake from high-speed vessels: its group structure and run-up Nonlinear Processes in Geophysics, 20, 179-188, 2013 Author(s): I. Didenkulova and A. Rodin High-amplitude water waves induced by high-speed vessels are regularly observed in Tallinn Bay, the Baltic Sea, causing intense beach erosion and disturbing marine habitants in the coastal zone. Such a strong impact on the coast may be a result of a certain group structure of the wave wake. In order to understand it, here we present an experimental study of the group structure of these wakes at Pikakari beach, Tallinn Bay. The most energetic vessel waves at this location (100 m from the coast at the water depth 2.7 m) have amplitudes of about 1 m and periods of 8–10 s and cause maximum run-up heights on a beach up to 1.4 m. These waves represent frequency modulated packets where the largest and longest waves propagate ahead of other smaller amplitude and period waves. Sometimes the groups of different heights and periods can be separated even within one wave wake event. The wave heights within a wake are well described by the Weibull distribution, which has different parameters for wakes from different vessels. Wave run-up heights can also be described by Weibull distribution and its parameters can be connected to the parameters of the distribution of wave heights 100 m from the coast. Finally, the run-up of individual waves within a packet is studied. It is shown that the specific structure of frequency modulated wave packets, induced by high-speed vessels, leads to a sequence of high wave run-ups at the coast, even when the original wave heights are rather moderate. This feature can be a key to understanding the significant impact on coasts caused by fast vessels.

Beschreibung: Response of polar mesosphere summer echoes to geomagnetic disturbances in the Southern and Northern Hemispheres: the importance of nitric oxide Annales Geophysicae, 31, 333-347, 2013 Author(s): S. Kirkwood, E. Belova, P. Dalin, M. Mihalikova, D. Mikhaylova, D. Murtagh, H. Nilsson, K. Satheesan, J. Urban, and I. Wolf The relationship between polar mesosphere summer echoes (PMSE) and geomagnetic disturbances (represented by magnetic K indices) is examined. Calibrated PMSE reflectivities for the period May 2006–February 2012 are used from two 52.0/54.5 MHz radars located in Arctic Sweden (68° N, geomagnetic latitude 65°) and at two different sites in Queen Maud Land, Antarctica (73°/72° S, geomagnetic latitudes 62°/63°). In both the Northern Hemisphere (NH) and the Southern Hemisphere (SH) there is a strong increase in mean PMSE reflectivity between quiet and disturbed geomagnetic conditions. Mean volume reflectivities are slightly lower at the SH locations compared to the NH, but the position of the peak in the lognormal distribution of PMSE reflectivities is close to the same at both NH and SH locations, and varies only slightly with magnetic disturbance level. Differences between the sites, and between geomagnetic disturbance levels, are primarily due to differences in the high-reflectivity tail of the distribution. PMSE occurrence rates are essentially the same at both NH and SH locations during most of the PMSE season when a sufficiently low detection threshold is used so that the peak in the lognormal distribution is included. When the local-time dependence of the PMSE response to geomagnetic disturbance level is considered, the response in the NH is found to be immediate at most local times, but delayed by several hours in the afternoon sector and absent in the early evening. At the SH sites, at lower magnetic latitude, there is a delayed response (by several hours) at almost all local times. At the NH (auroral zone) site, the dependence on magnetic disturbance is highest during evening-to-morning hours. At the SH (sub-auroral) sites the response to magnetic disturbance is weaker but persists throughout the day. While the immediate response to magnetic activity can be qualitatively explained by changes in electron density resulting from energetic particle precipitation, the delayed response can largely be explained by changes in nitric oxide concentrations. Observations of nitric oxide concentration at PMSE heights by the Odin satellite support this hypothesis. Sensitivity to geomagnetic disturbances, including nitric oxide produced during these disturbances, can explain previously reported differences between sites in the auroral zone and those at higher or lower magnetic latitudes. The several-day lifetime of nitric oxide can also explain earlier reported discrepancies between high correlations for average conditions (year-by-year PMSE reflectivities and K indices) and low correlations for minute-to-day timescales.

Beschreibung: Observations of poleward-propagating large-scale traveling ionospheric disturbances in southern China Annales Geophysicae, 31, 377-385, 2013 Author(s): F. Ding, W. Wan, B. Ning, B. Zhao, Q. Li, Y. Wang, L. Hu, R. Zhang, and B. Xiong We report here on two cases of poleward-propagating large-scale traveling ionospheric disturbances (LSTIDs) in China during a medium-scale storm between 27 May and 1 June 2011. The observations were conducted by making use of the Global Positioning System network and ionosondes in China and Southeast Asia. One northeastward-propagating LSTID occurred on the morning of 30 May, while the other was observed during the nighttime of 1 June. Both poleward-traveling LSTIDs occurred during the storm's recovery phase in southern China's low-latitude region (geomagnetic latitude ~ 7.3–24° N) and experienced severe dissipation during their propagation from south to north. Although the initial relative amplitude of the nighttime LSTID was ~ 60% larger than that of the morning event, the nighttime event dissipated more quickly than the morning event because of a strong nighttime enhancement in background total electronic content (TEC) during storm time, which led to strong ion-drag dissipation during the evening. The poleward-propagating LSTIDs exhibit a narrower latitudinal range, a smaller amplitude, and a slightly higher elevation compared with the equatorward-moving LSTIDs observed in the same region. Given these features, the poleward-propagating LSTIDs were likely excited by some local source near southern China. Excitation of secondary LSTIDs during the dissipation of some primary medium-scale disturbances from the lower atmosphere is a possible mechanism.

Beschreibung: A new method for solving the MHD equations in the magnetosheath Annales Geophysicae, 31, 419-437, 2013 Author(s): C. Nabert, K.-H. Glassmeier, and F. Plaschke We present a new analytical method to derive steady-state magnetohydrodynamic (MHD) solutions of the magnetosheath in different levels of approximation. With this method, we calculate the magnetosheath's density, velocity, and magnetic field distribution as well as its geometry. Thereby, the solution depends on the geomagnetic dipole moment and solar wind conditions only. To simplify the representation, we restrict our model to northward IMF with the solar wind flow along the stagnation streamline. The sheath's geometry, with its boundaries, bow shock and magnetopause, is determined self-consistently. Our model is stationary and time relaxation has not to be considered as in global MHD simulations. Our method uses series expansion to transfer the MHD equations into a new set of ordinary differential equations. The number of equations is related to the level of approximation considered including different physical processes. These equations can be solved numerically; however, an analytical approach for the lowest-order approximation is also presented. This yields explicit expressions, not only for the flow and field variations but also for the magnetosheath thickness, depending on the solar wind parameters. Results are compared to THEMIS data and offer a detailed explanation of, e.g., the pile-up process and the corresponding plasma depletion layer, the bow shock and magnetopause geometry, the magnetosheath thickness, and the flow deceleration.

Beschreibung: Investigation of a mesospheric bore event over northern China Annales Geophysicae, 31, 409-418, 2013 Author(s): Q. Li, J. Xu, J. Yue, X. Liu, W. Yuan, B. Ning, S. Guan, and J. P. Younger A mesospheric bore event was observed using an OH all-sky airglow imager (ASAI) at Xinglong (40.2° N, 117.4° E), in northern China, on the night of 8–9 January 2011. Simultaneous observations by a Doppler meteor radar, a broadband sodium lidar, and TIMED/SABER OH intensity and temperature measurements are used to investigate the characteristics and environment of the bore propagation and the possible relations with the Na density perturbations. The bore propagated from northeast to southwest and divided the sky into bright and dark halves. The calculations show that the bore has an average phase velocity of 68 m s −1 . The crests following the bore have a horizontal wavelength of ~ 22 km. These parameters are consistent with the hydraulic jump theory proposed by Dewan and Picard, as well as the previous bore reports. Simultaneous wind measurements from the Doppler meteor radar at Shisanling (40.3° N, 116.2° E) and temperature data from SABER on board the TIMED satellite are used to characterize the propagating environment of the bore. The result shows that a thermal-Doppler duct exists near the OH layer that supports the horizontal propagation of the bore. Simultaneous Na lidar observations at Yanqing (40.4° N, 116.0° E) suggest that there is a downward displacement of Na density during the passage of the mesospheric bore event.

Beschreibung: Near real-time estimation of water vapour in the troposphere using ground GNSS and the meteorological data Annales Geophysicae, 30, 1379-1391, 2012 Author(s): J. Bosy, J. Kaplon, W. Rohm, J. Sierny, and T. Hadas The near real-time (NRT) high resolution water vapour distribution models can be constructed based on GNSS observations delivered from Ground Base Augmentation Systems (GBAS) and ground meteorological data. Since 2008 in the territory of Poland, a GBAS system called ASG-EUPOS (Active Geodetic Network) has been operating. This paper addresses the problems concerning construction of the NRT model of water vapour distribution in the troposphere near Poland. The first section presents all available GNSS and ground meteorological stations in the area of Poland and neighbouring countries. In this section, data feeding scheme is discussed, together with timeline and time resolution. The high consistency between measured and interpolated temperature value is shown, whereas some discrepancy in the pressure is observed. In the second section, the NRT GNSS data processing strategy of ASG-EUPOS network is discussed. Preliminary results show fine alignment of the obtained Zenith Troposphere Delays (ZTDs) with reference data from European Permanent Network (EPN) processing center. The validation of NRT troposphere products against daily solution shows 15 mm standard deviation of obtained ZTD differences. The last section presents the first results of 2-D water vapour distribution above the GNSS network and application of the tomographic model to 3-D distribution of water vapour in the atmosphere. The GNSS tomography model, working on the simulated data from numerical forecast model, shows high consistency with the reference data (by means of standard deviation 4 mm km −1 or 4 ppm), however, noise analysis shows high solution sensitivity to errors in observations. The discrepancy for real data preliminary solution (measured as a mean standard deviation) between reference NWP data and tomography data was on the level of 9 mm km −1 (or 9 ppm) in terms of wet refractivity.

Beschreibung: Variations of the ionospheric TEC using simultaneous measurements from the China Crustal Movement Observation Network Annales Geophysicae, 30, 1423-1433, 2012 Author(s): Y. W. Wu, R. Y. Liu, B. C. Zhang, Z. S. Wu, J. S. Ping, J. M. Liu, and Z. J. Hu Variations of the ionospheric Total Electron Content (TEC) over China are investigated using the TEC data obtained from China Crustal Movement Observation Network in the year 2004. The results show a single-peak occurred in post-noon for the diurnal variation and two peaks exit around two equinox points, respectively, for the seasonal variation. Overall, the values of TEC increased from the north to the south of China. There were small but clear longitudinal differences in both sides of the longitudes with zero magnetic declination. The intensity of the day-to-day variation of TEC was not a monotonic change along the latitudes. It was usually weaker in the middle of China than that in the north or south. Comparing with the maximum F-layer electron density ( Nm F2) derived from the ionosonde stations in China, it is found that the day-to-day variation of TEC was less significant than that of Nm F2, and that the northern crest of the equatorial anomaly identified from the Nm F2 data can reach Guangzhou-region. While, the TEC crest was hardly observed in the same location. This is probably caused by the tilt of topside ionosphere near the northern anomaly crest region at lower latitudes.

Beschreibung: Collisionless reconnection: magnetic field line interaction Annales Geophysicae, 30, 1515-1528, 2012 Author(s): R. A. Treumann, W. Baumjohann, and W. D. Gonzalez Magnetic field lines are quantum objects carrying one quantum Φ 0 = 2π h/e of magnetic flux and have finite radius λ m . Here we argue that they possess a very specific dynamical interaction. Parallel field lines reject each other. When confined to a certain area they form two-dimensional lattices of hexagonal structure. We estimate the filling factor of such an area. Anti-parallel field lines, on the other hand, attract each other. We identify the physical mechanism as being due to the action of the gauge potential field, which we determine quantum mechanically for two parallel and two anti-parallel field lines. The distortion of the quantum electrodynamic vacuum causes a cloud of virtual pairs. We calculate the virtual pair production rate from quantum electrodynamics and estimate the virtual pair cloud density, pair current and Lorentz force density acting on the field lines via the pair cloud. These properties of field line dynamics become important in collisionless reconnection, consistently explaining why and how reconnection can spontaneously set on in the field-free centre of a current sheet below the electron-inertial scale.

Beschreibung: Observation of a tropopause fold by MARA VHF wind-profiler radar and ozonesonde at Wasa, Antarctica: comparison with ECMWF analysis and a WRF model simulation Annales Geophysicae, 30, 1411-1421, 2012 Author(s): M. Mihalikova, S. Kirkwood, J. Arnault, and D. Mikhaylova Tropopause folds are one of the mechanisms of stratosphere–troposphere exchange, which can bring ozone rich stratospheric air to low altitudes in the extra-tropical regions. They have been widely studied at northern mid- or high latitudes, but so far almost no studies have been made at mid- or high southern latitudes. The Moveable Atmospheric Radar for Antarctica (MARA), a 54.5 MHz wind-profiler radar, has operated at the Swedish summer station Wasa, Antarctica (73° S, 13.5° W) during austral summer seasons from 2007 to 2011 and has observed on several occasions signatures similar to those caused by tropopause folds at comparable Arctic latitudes. Here a case study is presented of one of these events when an ozonesonde successfully sampled the fold. Analysis from European Center for Medium Range Weather Forecasting (ECMWF) is used to study the circumstances surrounding the event, and as boundary conditions for a mesoscale simulation using the Weather Research and Forecasting (WRF) model. The fold is well resolved by the WRF simulation, and occurs on the poleward side of the polar jet stream. However, MARA resolves fine-scale layering associated with the fold better than the WRF simulation.

Beschreibung: Chorus wave-normal statistics in the Earth's radiation belts from ray tracing technique Annales Geophysicae, 30, 1223-1233, 2012 Author(s): H. Breuillard, Y. Zaliznyak, V. Krasnoselskikh, O. Agapitov, A. Artemyev, and G. Rolland Discrete ELF/VLF (Extremely Low Frequency/Very Low Frequency) chorus emissions are one of the most intense electromagnetic plasma waves observed in radiation belts and in the outer terrestrial magnetosphere. These waves play a crucial role in the dynamics of radiation belts, and are responsible for the loss and the acceleration of energetic electrons. The objective of our study is to reconstruct the realistic distribution of chorus wave-normals in radiation belts for all magnetic latitudes. To achieve this aim, the data from the electric and magnetic field measurements onboard Cluster satellite are used to determine the wave-vector distribution of the chorus signal around the equator region. Then the propagation of such a wave packet is modeled using three-dimensional ray tracing technique, which employs K. Rönnmark's WHAMP to solve hot plasma dispersion relation along the wave packet trajectory. The observed chorus wave distributions close to waves source are first fitted to form the initial conditions which then propagate numerically through the inner magnetosphere in the frame of the WKB approximation. Ray tracing technique allows one to reconstruct wave packet properties (electric and magnetic fields, width of the wave packet in k -space, etc.) along the propagation path. The calculations show the spatial spreading of the signal energy due to propagation in the inhomogeneous and anisotropic magnetized plasma. Comparison of wave-normal distribution obtained from ray tracing technique with Cluster observations up to 40° latitude demonstrates the reliability of our approach and applied numerical schemes.

Beschreibung: Fractal analysis of ULF electromagnetic emissions in possible association with earthquakes in China Nonlinear Processes in Geophysics, 19, 577-583, 2012 Author(s): Y. Ida, D. Yang, Q. Li, H. Sun, and M. Hayakawa The long-term data (during the period of 1 March 2003 through 31 December 2006) of ULF geomagnetic variations observed at Kashi station (geographic coordinates: 39.5° N, 76.0° E) in China have been used to investigate the long-term variation of fractal dimension of ULF emissions. We have studied the changes in fractal dimension in association with several earthquakes around the observation station. It is then found that a significant change (or decrease) in the fractal dimension of the Z component took place before the 1 September 2003 earthquake, which lends a further support to our previous finding based on our improved polarization analysis for the same earthquake. The results obtained are discussed in the contexts of a few aspects (detectability of seismogenic emissions, comparison with previous results by other analysis methods, the importance of fractal analysis in the nonlinear process of the lithosphere and earthquake prediction).

Beschreibung: High-precision measurement of satellite velocity using the EISCAT radar Annales Geophysicae, 30, 1555-1565, 2012 Author(s): T. Nygrén, J. Markkanen, A. Aikio, and M. Voiculescu This paper presents a method of measuring the velocity of a hard target using radar pulses reflected from the target flying through the radar beam. The method has two stages. First, the Doppler shifts of the echo pulses are calculated at a high accuracy with an algorithm which largely improves the accuracy given by the Fourier transform. The algorithm also calculates the standard deviations of the Doppler frequencies with Monte Carlo simulation. The second step is to fit the results from a sequence of radar pulses to a velocity model allowing linear variation of the second time derivative of target range. The achieved accuracies are demonstrated using radio pulses reflected by a satellite passing through the beam of the EISCAT UHF radar working at 930-MHz frequency. At high SNR levels, the standard deviations of the frequency from a single pulse reach typically down to 0.2 Hz. The best standard deviations of velocity fit are below 5 mm s −1 while those of the second time derivative of range are below 1 cm s −2 .

Beschreibung: Predictability of a low-order interactive ensemble Nonlinear Processes in Geophysics, 19, 273-282, 2012 Author(s): L. Siqueira and B. Kirtman In this paper, numerical and analytical studies were performed to uncover the mechanisms controlling the changes in ensemble spread of a low-order coupled model with multiple atmospheric realizations. An interactive ensemble approach was applied to a coupled dynamical system based on two versions of the Lorenz 63 model designed in order to imitate the behavior of a coupled system with different time scales. In the dynamic system used in this work the spread of ensemble members is highly dependent on the mean state corresponding to asymmetries in predictability. The slowness of the slow model and the intensity of the boundary forcing anomalies both contribute to the asymmetry and phase locking of both subsystems. The mechanisms controlling the fast model spread were uncovered revealing uncertainty dynamics depending on the location of ensemble members in the fast model phase space and implicitly on the slowness and magnitude of the slow model anomalies.

Beschreibung: Equatorial plasma bubbles and L-band scintillations in Africa during solar minimum Annales Geophysicae, 30, 675-682, 2012 Author(s): V. V. Paznukhov, C. S. Carrano, P. H. Doherty, K. M. Groves, R. G. Caton, C. E. Valladares, G. K. Seemala, C. T. Bridgwood, J. Adeniyi, L. L. N. Amaeshi, B. Damtie, F. D'Ujanga Mutonyi, J. O. H. Ndeda, P. Baki, O. K. Obrou, B. Okere, and G. M. Tsidu We report on the longitudinal, local time and seasonal occurrence of equatorial plasma bubbles (EPBs) and L band (GPS) scintillations over equatorial Africa. The measurements were made in 2010, as a first step toward establishing the climatology of ionospheric irregularities over Africa. The scintillation intensity is obtained by measuring the standard deviation of normalized GPS signal power. The EPBs are detected using an automated technique, where spectral analysis is used to extract and identify EPB events from the GPS TEC measurements. Overall, the observed seasonal climatology of the EPBs as well as GPS scintillations in equatorial Africa is adequately explained by geometric arguments, i.e., by the alignment of the solar terminator and local geomagnetic field, or STBA hypothesis (Tsunoda, 1985, 2010a). While plasma bubbles and scintillations are primarily observed during equinoctial periods, there are longitudinal differences in their seasonal occurrence statistics. The Atlantic sector has the most intense, longest lasting, and highest scintillation occurrence rate in-season. There is also a pronounced increase in the EPB occurrence rate during the June solstice moving west to east. In Africa, the seasonal occurrence shifts towards boreal summer solstice, with fewer occurrences and shorter durations in equinox seasons. Our results also suggest that the occurrence of plasma bubbles and GPS scintillations over Africa are well correlated, with scintillation intensity depending on depletion depth. A question remains about the possible physical mechanisms responsible for the difference in the occurrence phenomenology of EPBs and GPS scintillations between different regions in equatorial Africa.

Beschreibung: Understanding the drivers for the development of design rules for the synthesis of cylindrical flexures Mechanical Sciences, 3, 25-32, 2012 Author(s): M. J. Telleria and M. L. Culpepper Cylindrical flexures (CFs), defined as flexures with only one finite radius of curvature loaded normal to the plane of curvature, present an interesting research direction in compliant mechanisms. CFs are constructed out of a cylindrical stock which leads to geometry, manufacturability, and compatibility advantages. Synthesis rules must be developed to design these new systems effectively. Current knowledge in flexure design pertains to straight-beam flexures or curved flexures loaded along the plane of curvature. CFs present a challenge because their mechanics differ from those of straight beams, and although their modelling has been researched thoroughly it has yet to be distilled into element and system creation rules. This paper uses models and finite element analysis to demonstrate that current design rules for straight-beam flexures are insufficient and inadequate for the design of CF systems. The presented discussion will show that CFs differ both at the element and systems levels, and therefore future research will focus on developing the three components of the building block approach: (i) reworking of element mechanics models to reveal the parameters which cause the kinematics of the curved beam to differ from those of the straight beam, (ii) development of a visual stiffness representation, and (iii) formation of system creation rules.

Beschreibung: The ultra-fast Kelvin waves in the equatorial ionosphere: observations and modeling Annales Geophysicae, 31, 209-215, 2013 Author(s): A. N. Onohara, I. S. Batista, and H. Takahashi The main purpose of this study is to investigate the vertical coupling between the mesosphere and lower thermosphere (MLT) region and the ionosphere through ultra-fast Kelvin (UFK) waves in the equatorial atmosphere. The effect of UFK waves on the ionospheric parameters was estimated using an ionospheric model which calculates electrostatic potential in the E-region and solves coupled electrodynamics of the equatorial ionosphere in the E- and F-regions. The UFK wave was observed in the South American equatorial region during February–March 2005. The MLT wind data obtained by meteor radar at São João do Cariri (7.5° S, 37.5° W) and ionospheric F-layer bottom height ( h 'F) observed by ionosonde at Fortaleza (3.9° S; 38.4° W) were used in order to calculate the wave characteristics and amplitude of oscillation. The simulation results showed that the combined electrodynamical effect of tides and UFK waves in the MLT region could explain the oscillations observed in the ionospheric parameters.

Beschreibung: Payload charging events in the mesosphere and their impact on Langmuir type electric probes Annales Geophysicae, 31, 187-196, 2013 Author(s): T. A. Bekkeng, A. Barjatya, U.-P. Hoppe, A. Pedersen, J. I. Moen, M. Friedrich, and M. Rapp Three sounding rockets were launched from Andøya Rocket Range in the ECOMA campaign in December 2010. The aim was to study the evolution of meteoric smoke particles during a major meteor shower. Of the various instruments onboard the rocket payload, this paper presents the data from a multi-Needle Langmuir Probe (m-NLP) and a charged dust detector. The payload floating potential, as observed using the m-NLP instrument, shows charging events on two of the three flights. These charging events cannot be explained using a simple charging model, and have implications towards the use of fixed bias Langmuir probes on sounding rockets investigating mesospheric altitudes. We show that for a reliable use of a single fixed bias Langmuir probe as a high spatial resolution relative density measurement, each payload should also carry an additional instrument to measure payload floating potential, and an instrument that is immune to spacecraft charging and measures absolute plasma density.

Beschreibung: Sub-modeling approach for obtaining structural stress histories during dynamic analysis Mechanical Sciences, 4, 21-31, 2013 Author(s): T. T. Rantalainen, A. M. Mikkola, and T. J. Björk Modern machine structures are often fabricated by welding. From a fatigue point of view, the structural details and especially, the welded details are the most prone to fatigue damage and failure. Design against fatigue requires information on the fatigue resistance of a structure's critical details and the stress loads that act on each detail. Even though, dynamic simulation of flexible bodies is already current method for analyzing structures, obtaining the stress history of a structural detail during dynamic simulation is a challenging task; especially when the detail has a complex geometry. In particular, analyzing the stress history of every structural detail within a single finite element model can be overwhelming since the amount of nodal degrees of freedom needed in the model may require an impractical amount of computational effort. The purpose of computer simulation is to reduce amount of prototypes and speed up the product development process. Also, to take operator influence into account, real time models, i.e. simplified and computationally efficient models are required. This in turn, requires stress computation to be efficient if it will be performed during dynamic simulation. The research looks back at the theoretical background of multibody simulation and finite element method to find suitable parts to form a new approach for efficient stress calculation. This study proposes that, the problem of stress calculation during dynamic simulation can be greatly simplified by using a combination of Floating Frame of Reference Formulation with modal superposition and a sub-modeling approach. In practice, the proposed approach can be used to efficiently generate the relevant fatigue assessment stress history for a structural detail during or after dynamic simulation. Proposed approach is demonstrated in practice using one numerical example. Even though, examples are simplified the results show that approach is applicable and can be used as proposed.

Beschreibung: Flexible joints in structural and multibody dynamics Mechanical Sciences, 4, 65-77, 2013 Author(s): O. A. Bauchau and S. Han Flexible joints, sometimes called bushing elements or force elements, are found in all structural and multibody dynamics codes. In their simplest form, flexible joints simply consist of sets of three linear and three torsional springs placed between two nodes of the model. For infinitesimal deformations, the selection of the lumped spring constants is an easy task, which can be based on a numerical simulation of the joint or on experimental measurements. If the joint undergoes finite deformations, identification of its stiffness characteristics is not so simple, specially if the joint is itself a complex system. When finite deformations occur, the definition of deformation measures becomes a critical issue. This paper proposes a family of tensorial deformation measures suitable for elastic bodies of finite dimension. These families are generated by two parameters that can be used to modify the constitutive behavior of the joint, while maintaining the tensorial nature of the deformation measures. Numerical results demonstrate the objectivity of the deformations measures, a feature that is not shared by the deformations measures presently used in the literature. The impact of the choice of the two parameters on the constitutive behavior of the flexible joint is also investigated.

Beschreibung: Geometrically exact Cosserat rods with Kelvin–Voigt type viscous damping Mechanical Sciences, 4, 79-96, 2013 Author(s): J. Linn, H. Lang, and A. Tuganov We present the derivation of a simple viscous damping model of Kelvin–Voigt type for geometrically exact Cosserat rods from three-dimensional continuum theory. Assuming moderate curvature of the rod in its reference configuration, strains remaining small in its deformed configurations, strain rates that vary slowly compared to internal relaxation processes, and a homogeneous and isotropic material, we obtain explicit formulas for the damping parameters of the model in terms of the well known stiffness parameters of the rod and the retardation time constants defined as the ratios of bulk and shear viscosities to the respective elastic moduli. We briefly discuss the range of validity of the Kelvin–Voigt model and illustrate its behaviour for large bending deformations with a numerical example.

Beschreibung: Interaction of a monopole vortex with an isolated topographic feature in a three-layer geophysical flow Nonlinear Processes in Geophysics, 20, 107-119, 2013 Author(s): E. A. Ryzhov and K. V. Koshel In the frame of a three-layer, quasi-geostrophic analytical model of an f -plane geophysical flow, the Lagrangian advection induced by the interaction of a monopole vortex with an isolated topographic feature is addressed. Two different cases when the monopole is located either within the upper or the middle layer are of our interest. In the bottom layer, there is a delta-function topographic feature, which generates a closed recirculation region in its vicinity due to the background flow. This recirculation region extends to the middle and upper layers, and it plays the role of a topographic vortex. The interaction between the monopole and the topographic vortex causes a complex, including chaotic, advection of fluid particles. We show that the model's parameters, namely the monopole and topographic vortices' strengths and initial positions, and the layers' depths and densities, are responsible for the diverse advection patterns. While the patterns are rather complicated, one can single out two major processes, which mostly govern the fluid particle advection. The first one is the variation in time of the system's phase space structure, so that within the closed region of the topographic vortex, there appear periodically unclosed particle pathways by which the particles leave the topographic vortex. The second one is chaotic advection that arises from the nonstationarity of the monopole–topography interaction.

Beschreibung: Electric potential differences across auroral generator interfaces Annales Geophysicae, 31, 251-261, 2013 Author(s): J. De Keyser and M. Echim Strong localized high-altitude auroral electric fields, such as those observed by Cluster, are often associated with magnetospheric interfaces. The type of high-altitude electric field profile (monopolar, bipolar, or more complicated) depends on the properties of the plasmas on either side of the interface, as well as on the total electric potential difference across the structure. The present paper explores the role of this cross-field electric potential difference in the situation where the interface is a tangential discontinuity. A self-consistent Vlasov description is used to determine the equilibrium configuration for different values of the transverse potential difference. A major observation is that there exist limits to the potential difference, beyond which no equilibrium configuration of the interface can be sustained. It is further demonstrated how the plasma densities and temperatures affect the type of electric field profile in the transition, with monopolar electric fields appearing primarily when the temperature contrast is large. These findings strongly support the observed association of monopolar fields with the plasma sheet boundary. The role of shear flow tangent to the interface is also examined.

Beschreibung: Variability of ionospheric TEC during solar and geomagnetic minima (2008 and 2009): external high speed stream drivers Annales Geophysicae, 31, 263-276, 2013 Author(s): O. P. Verkhoglyadova, B. T. Tsurutani, A. J. Mannucci, M. G. Mlynczak, L. A. Hunt, and T. Runge We study solar wind–ionosphere coupling through the late declining phase/solar minimum and geomagnetic minimum phases during the last solar cycle (SC23) – 2008 and 2009. This interval was characterized by sequences of high-speed solar wind streams (HSSs). The concomitant geomagnetic response was moderate geomagnetic storms and high-intensity, long-duration continuous auroral activity (HILDCAA) events. The JPL Global Ionospheric Map (GIM) software and the GPS total electron content (TEC) database were used to calculate the vertical TEC (VTEC) and estimate daily averaged values in separate latitude and local time ranges. Our results show distinct low- and mid-latitude VTEC responses to HSSs during this interval, with the low-latitude daytime daily averaged values increasing by up to 33 TECU (annual average of ~20 TECU) near local noon (12:00 to 14:00 LT) in 2008. In 2009 during the minimum geomagnetic activity (MGA) interval, the response to HSSs was a maximum of ~30 TECU increases with a slightly lower average value than in 2008. There was a weak nighttime ionospheric response to the HSSs. A well-studied solar cycle declining phase interval, 10–22 October 2003, was analyzed for comparative purposes, with daytime low-latitude VTEC peak values of up to ~58 TECU (event average of ~55 TECU). The ionospheric VTEC changes during 2008–2009 were similar but ~60% less intense on average. There is an evidence of correlations of filtered daily averaged VTEC data with Ap index and solar wind speed. We use the infrared NO and CO 2 emission data obtained with SABER on TIMED as a proxy for the radiation balance of the thermosphere. It is shown that infrared emissions increase during HSS events possibly due to increased energy input into the auroral region associated with HILDCAAs. The 2008–2009 HSS intervals were ~85% less intense than the 2003 early declining phase event, with annual averages of daily infrared NO emission power of ~ 3.3 × 10 10 W and 2.7 × 10 10 W in 2008 and 2009, respectively. The roles of disturbance dynamos caused by high-latitude winds (due to particle precipitation and Joule heating in the auroral zones) and of prompt penetrating electric fields (PPEFs) in the solar wind–ionosphere coupling during these intervals are discussed. A correlation between geoeffective interplanetary electric field components and HSS intervals is shown. Both PPEF and disturbance dynamo mechanisms could play important roles in solar wind–ionosphere coupling during prolonged (up to days) external driving within HILDCAA intervals.

Beschreibung: CHRONO: a parallel multi-physics library for rigid-body, flexible-body, and fluid dynamics Mechanical Sciences, 4, 49-64, 2013 Author(s): H. Mazhar, T. Heyn, A. Pazouki, D. Melanz, A. Seidl, A. Bartholomew, A. Tasora, and D. Negrut The last decade witnessed a manifest shift in the microprocessor industry towards chip designs that promote parallel computing. Until recently the privilege of a select group of large research centers, Teraflop computing is becoming a commodity owing to inexpensive GPU cards and multi to many-core x86 processors. This paradigm shift towards large scale parallel computing has been leveraged in CHRONO, a freely available C++ multi-physics simulation package. CHRONO is made up of a collection of loosely coupled components that facilitate different aspects of multi-physics modeling, simulation, and visualization. This contribution provides an overview of CHRONO::Engine, CHRONO::Flex, CHRONO::Fluid, and CHRONO::Render, which are modules that can capitalize on the processing power of hundreds of parallel processors. Problems that can be tackled in CHRONO include but are not limited to granular material dynamics, tangled large flexible structures with self contact, particulate flows, and tracked vehicle mobility. The paper presents an overview of each of these modules and illustrates through several examples the potential of this multi-physics library.

Beschreibung: Algebraic analysis of kinematics of multibody systems Mechanical Sciences, 4, 33-47, 2013 Author(s): S. Piipponen and J. Tuomela The constructive commutative algebra is very useful in the kinematical analysis of the mechanisms because a large class of systems can be described using polynomial equations. We show that one can analyze quite complicated systems using a sort of divide and conquer strategy to decompose the system, and hence the configuration space, into simpler parts. The key observation is that it seems that typically systems indeed have a lot of distinct components, but usually only one of them is physically relevant. Hence if one finds the equations describing the component of interest the analysis of this system can be surprisingly simple compared to the original system. In particular typically the possible singularities of the original system disappear when one restricts the attention to the relevant component. On the technical side we show that some basic constraints used to define joints in 3 dimensional mechanisms can be decomposed to simpler parts. This has significant practical consequences because using these fundamental decompositions when writing the equations for complicated mechanisms decreases dramatically the complexity of the required computations.

Beschreibung: Case study of stratospheric gravity waves of convective origin over Arctic Scandinavia – VHF radar observations and numerical modelling Annales Geophysicae, 31, 239-250, 2013 Author(s): A. Réchou, J. Arnault, P. Dalin, and S. Kirkwood Orography is a well-known source of gravity and inertia-gravity waves in the atmosphere. Other sources, such as convection, are also known to be potentially important but the large amplitude of orographic waves over Scandinavia has generally precluded the possibility to study such other sources experimentally in this region. In order to better understand the origin of stratospheric gravity waves observed by the VHF radar ESRAD (Esrange MST radar) over Kiruna, in Arctic Sweden (67.88° N, 21.10° E), observations have been compared to simulations made using the Weather Research and Forecasting model (WRF) with and without the effects of orography and clouds. This case study concerns gravity waves observed from 00:00 UTC on 18 February to 12:00 UTC on 20 February 2007. We focus on the wave signatures in the static stability field and vertical wind deduced from the simulations and from the observations as these are the parameters which are provided by the observations with the best height coverage. As is common at this site, orographic gravity waves were produced over the Scandinavian mountains and observed by the radar. However, at the same time, southward propagation of fronts in the Barents Sea created short-period waves which propagated into the stratosphere and were transported, embedded in the cyclonic winds, over the radar site.

Beschreibung: Long-term trends in the ionospheric F2 region with different solar activity indices Annales Geophysicae, 31, 291-303, 2013 Author(s): J. Mielich and J. Bremer A new comprehensive data collection by Damboldt and Suessmann (2012a) with monthly fo F2 and M(3000)F2 median values is an excellent basis for the derivation of long-term trends in the ionospheric F2 region. Ionospheric trends have been derived only for stations with data series of at least 22 years (124 stations with fo F2 data and 113 stations with M(3000)F2 data) using a twofold regression analysis depending on solar and geomagnetic activity. Three main results have been derived: Firstly, it could be shown that the solar 10.7 cm radio flux F10.7 is a better index for the description of the solar activity than the relative solar sunspot number R as well as the solar EUV proxy E10.7. Secondly, the global mean fo F2 and hm F2 trends derived for the interval between 1948 and 2006 are in surprisingly good agreement with model calculations of an increasing atmospheric greenhouse effect (Rishbeth and Roble, 1992). Thirdly, during the years 2007 until 2009, the hm F2 values and to a smaller amount the fo F2 values strongly decrease. The reason for this effect is a reduction of the thermospheric density and ionization due to a markedly reduced solar EUV irradiation and extremely small geomagnetic activity during the solar cycle 23/24 minimum.

Beschreibung: Reaction Null Space of a multibody system with applications in robotics Mechanical Sciences, 4, 97-112, 2013 Author(s): D. N. Nenchev This paper provides an overview of implementation examples based on the Reaction Null Space formalism, developed initially to tackle the problem of satellite-base disturbance of a free-floating space robot, when the robot arm is activated. The method has been applied throughout the years to other unfixed-base systems, e.g. flexible-base and macro/mini robot systems, as well as to the balance control problem of humanoid robots. The paper also includes most recent results about complete dynamical decoupling of the end-link of a fixed-base robot, wherein the end-link is regarded as the unfixed-base. This interpretation is shown to be useful with regard to motion/force control scenarios. Respective implementation results are provided.

Beschreibung: Characterization of turbulence stability through the identification of multifractional Brownian motions Nonlinear Processes in Geophysics, 20, 97-106, 2013 Author(s): K. C. Lee Multifractional Brownian motions have become popular as flexible models in describing real-life signals of high-frequency features in geoscience, microeconomics, and turbulence, to name a few. The time-changing Hurst exponent, which describes regularity levels depending on time measurements, and variance, which relates to an energy level, are two parameters that characterize multifractional Brownian motions. This research suggests a combined method of estimating the time-changing Hurst exponent and variance using the local variation of sampled paths of signals. The method consists of two phases: initially estimating global variance and then accurately estimating the time-changing Hurst exponent. A simulation study shows its performance in estimation of the parameters. The proposed method is applied to characterization of atmospheric stability in which descriptive statistics from the estimated time-changing Hurst exponent and variance classify stable atmosphere flows from unstable ones.

Beschreibung: Simultaneous measurement of OI 557.7 nm, O 2 (0, 1) Atmospheric Band and OH (6, 2) Meinel Band nightglow at Kolhapur (17° N), India Annales Geophysicae, 31, 197-208, 2013 Author(s): N. Parihar, A. Taori, S. Gurubaran, and G. K. Mukherjee Near-simultaneous measurements of OI 557.7 nm, O 2 (0, 1) Atmospheric Band and OH (6, 2) Meinel Band nightglow were carried out at Kolhapur (17° N), India during February–March 2007. Atmospheric temperatures around 87 and 94 km were derived from the knowledge of intensity measurements of spectral features OH (6, 2) Meinel Band and O 2 Atmospheric Band, respectively. An account of the behaviour of derived temperatures has been presented. The nocturnal behaviour of OH and O 2 temperatures is governed by the waves of tidal origin, whereas the signatures of planetary wave-like oscillations is noted in the night-to-night variation of two temperatures. This is probably the first report of planetary waves observed in nightglow temperatures from the Indian subcontinent.

Beschreibung: Refilling process in the plasmasphere: a 3-D statistical characterization based on Cluster density observations Annales Geophysicae, 31, 217-237, 2013 Author(s): G. Lointier, F. Darrouzet, P. M. E. Décréau, X. Vallières, S. Kougblénou, J. G. Trotignon, and J.-L. Rauch The Cluster mission offers an excellent opportunity to investigate the evolution of the plasma population in a large part of the inner magnetosphere, explored near its orbit's perigee, over a complete solar cycle. The WHISPER sounder, on board each satellite of the mission, is particularly suitable to study the electron density in this region, between 0.2 and 80 cm −3 . Compiling WHISPER observations during 1339 perigee passes distributed over more than three years of the Cluster mission, we present first results of a statistical analysis dedicated to the study of the electron density morphology and dynamics along and across magnetic field lines between L = 2 and L = 10. In this study, we examine a specific topic: the refilling of the plasmasphere and trough regions during extended periods of quiet magnetic conditions. To do so, we survey the evolution of the ap index during the days preceding each perigee crossing and sort out electron density profiles along the orbit according to three classes, namely after respectively less than 2 days, between 2 and 4 days, and more than 4 days of quiet magnetic conditions (ap ≤ 15 nT) following an active episode (ap 〉 15 nT). This leads to three independent data subsets. Comparisons between density distributions in the 3-D plasmasphere and trough regions at the three stages of quiet magnetosphere provide novel views about the distribution of matter inside the inner magnetosphere during several days of low activity. Clear signatures of a refilling process inside an expended plasmasphere in formation are noted. A plasmapause-like boundary, at L ~ 6 for all MLT sectors, is formed after 3 to 4 days and expends somewhat further after that. In the outer part of the plasmasphere ( L ~ 8), latitudinal profiles of median density values vary essentially according to the MLT sector considered rather than according to the refilling duration. The shape of these density profiles indicates that magnetic flux tubes are not fully replenished after 6 days of quiet conditions. In addition, the outer plasmasphere in the night and dawn sectors (22:00 to 10:00 MLT range) maintains an overall clear deficit of ionospheric population, when compared to the situation in the noon and dusk sectors (10:00 to 22:00 MLT range).

Beschreibung: Spatiotemporal variability analysis of diffuse radiation in China during 1981–2010 Annales Geophysicae, 31, 277-289, 2013 Author(s): X. L. Ren, H. L. He, L. Zhang, L. Zhou, G. R. Yu, and J. W. Fan Solar radiation is the primary driver of terrestrial plant photosynthesis and the diffuse component can enhance canopy light use efficiency (LUE), which in turn influences the carbon balance of terrestrial ecosystems. In this study we calculated the spatial data of diffuse radiation in China from 1981 to 2010, using a radiation decomposition model and spatial interpolation method based on observational data. Furthermore, we explored the spatiotemporal characteristics of diffuse radiation using GIS and trend analysis techniques. The results show the following: (1) The spatial patterns of perennial average of annual diffuse radiation during 1981–2010 are complex and inhomogeneous in China, generally lower in the north and higher in the south and west. The perennial average ranges from 1730.20 to 3064.41 MJ m −2 yr −1 across the whole country. (2) There is an increasing trend of annual diffuse radiation in China from 1981 to 2010 on the whole, with mean increasing amplitude of 7.03 MJ m −2 yr −1 per decade. Whereas a significant downtrend was observed in the first 10 years, distinct anomalies in 1982, 1983, 1991 and 1992 occurred due to the eruptions of El Chinchon and Pinatubo. (3) The spatial distribution of the temporal variability of diffuse radiation showed significant regional heterogeneity in addition to the seasonal differences. Northwestern China has the most evident downtrend, with highest decreasing rate of 6% per decade, while the Tibetan Plateau has the most evident uptrend, with highest increasing rate of up to 9% per decade. Such quantitative spatiotemporal characteristics of diffuse radiation are essential in regional scale modeling of terrestrial carbon dynamics.

Beschreibung: A tri-stage cluster identification model for accurate analysis of seismic catalogs Nonlinear Processes in Geophysics, 20, 143-162, 2013 Author(s): S. J. Nanda, K. F. Tiampo, G. Panda, L. Mansinha, N. Cho, and A. Mignan In this paper we propose a tri-stage cluster identification model that is a combination of a simple single iteration distance algorithm and an iterative K-means algorithm. In this study of earthquake seismicity, the model considers event location, time and magnitude information from earthquake catalog data to efficiently classify events as either background or mainshock and aftershock sequences. Tests on a synthetic seismicity catalog demonstrate the efficiency of the proposed model in terms of accuracy percentage (94.81% for background and 89.46% for aftershocks). The close agreement between lambda and cumulative plots for the ideal synthetic catalog and that generated by the proposed model also supports the accuracy of the proposed technique. There is flexibility in the model design to allow for proper selection of location and magnitude ranges, depending upon the nature of the mainshocks present in the catalog. The effectiveness of the proposed model also is evaluated by the classification of events in three historic catalogs: California, Japan and Indonesia. As expected, for both synthetic and historic catalog analysis it is observed that the density of events classified as background is almost uniform throughout the region, whereas the density of aftershock events are higher near the mainshocks.

Beschreibung: Scaling for lobe and cleft patterns in particle-laden gravity currents Nonlinear Processes in Geophysics, 20, 121-130, 2013 Author(s): A. Jackson, B. Turnbull, and R. Munro Lobe and cleft patterns are frequently observed at the leading edge of gravity currents, including non-Boussinesq particle-laden currents such as powder snow avalanches. Despite the importance of the instability in driving air entrainment, little is known about its origin or the mechanisms behind its development. In this paper we seek to gain a better understanding of these mechanisms from a laboratory scale model of powder snow avalanches using lightweight granular material. The instability mechanisms in these flows appear to be a combination of those found in both homogeneous Boussinesq gravity currents and unsuspended granular flows, with the size of the granular particles playing a central role in determining the wavelength of the lobe and cleft pattern. When scaled by particle diameter a relationship between the Froude number and the wavelength of the lobe and cleft pattern is found, where the wavelength increases monotonically with the Froude number.

Beschreibung: Inversion of the western Pacific subtropical high dynamic model and analysis of dynamic characteristics for its abnormality Nonlinear Processes in Geophysics, 20, 131-142, 2013 Author(s): M. Hong, R. Zhang, J. X. Li, J. J. Ge, and K. F. Liu Based on time series data of 500 hPa potential field from NCEP/NCAR (National Center for Environmental Forecast of American/National Center for Atmospheric Research), a novel consideration of empirical orthogonal function (EOF) time–space separation and dynamic system reconstruction for time series is introduced. This method consists of two parts: first, the dynamical model inversion and model parameter optimization are carried out on the EOF time coefficient series using the genetic algorithm (GA), and, second, a nonlinear dynamic model representing the subtropical high (SH) activity and its abnormality is established. The SH activity and its abnormal mechanism is studied using the developed dynamical model. Results show that the configuration and diversification of the SH equilibriums have good correspondence with the actual short–medium term abnormal activity of the SH. Change of SH potential field brought by the combination of equilibriums is more complex than that by mutation, and their exhibition patterns are different. The mutation behavior from high-value to low-value equilibriums of the SH in summer corresponds with the southward drop of the SH in the observed weather process. The combination behavior of the two steady equilibriums corresponds with disappearance of the "double-ridge" phenomenon of the SH. Dynamical mechanisms of these phenomena are explained.

Beschreibung: Evolution of the DeNOC-based dynamic modelling for multibody systems Mechanical Sciences, 4, 1-20, 2013 Author(s): S. K. Saha, S. V. Shah, and P. V. Nandihal Dynamic modelling of a multibody system plays very essential role in its analyses. As a result, several methods for dynamic modelling have evolved over the years that allow one to analyse multibody systems in a very efficient manner. One such method of dynamic modelling is based on the concept of the Decoupled Natural Orthogonal Complement (DeNOC) matrices. The DeNOC-based methodology for dynamics modelling, since its introduction in 1995, has been applied to a variety of multibody systems such as serial, parallel, general closed-loop, flexible, legged, cam-follower, and space robots. The methodology has also proven useful for modelling of proteins and hyper-degree-of-freedom systems like ropes, chains, etc. This paper captures the evolution of the DeNOC-based dynamic modelling applied to different type of systems, and its benefits over other existing methodologies. It is shown that the DeNOC-based modelling provides deeper understanding of the dynamics of a multibody system. The power of the DeNOC-based modelling has been illustrated using several numerical examples.

Beschreibung: Wavefield decomposition and phase space dynamics of the seismic noise at Volcàn de Colima, Mexico: evidence of a two-state source process Nonlinear Processes in Geophysics, 20, 71-84, 2013 Author(s): M. Palo and P. Cusano We analyse the seismic noise recorded at the Colima Volcano (Mexico) in the period December 2005–May 2006 by four broadband three-component seismic stations. Specifically, we characterize the spectral content of the signal and follow its time evolution along all the data set. Moreover, we infer the properties of the attractor in the phase space by false nearest neighbours analysis and Grassberger–Procaccia algorithm, and adopt a time-domain decomposition method (independent component analysis) to find the basic constituents (independent components) of the system. Constraints on the seismic wavefield are inferred by the polarization analysis. We find two states of the background seismicity visible in different time-intervals that are Phase A and Phase B. Phase A has a spectrum with two peaks at 0.15 Hz and 0.3 Hz, with the latter dominating, an attractor of correlation dimension close to 3, three quasi-monochromatic independent components, and a relevant fraction of crater-pointing polarization solutions in the near-field. In Phase B, the spectrum is preserved but with the highest peak at 0.15 Hz, the attractor has a correlation dimension close to 2, two independent components are extracted, and the polarization solutions are dominated by Rayleigh waves incoming from the southwest direction. We depict two sources acting on the background seismicity that are the microseismic noise loading on the Pacific coastline and a low-energy volcanic tremor. A change in the amplitude of the microseismic noise can induce the switching from a state of the system to the other.

Beschreibung: Dynamics of a gravity car race with application to the Pinewood Derby Mechanical Sciences, 3, 73-84, 2012 Author(s): B. P. Mann, M. M. Gibbs, and S. M. Sah This paper investigates the underlying physics of a gravity car race. This work seeks to provide a sound theoretical basis to elucidate the design considerations that maximize performance while simultaneously dispelling false assertions that may arise from incomplete analyses. The governing equations are derived and solved analytically to predict race times; trend analyses are then performed along with a sensitivity analysis to ascertain the most important factors that influence performance. The inferences from a conservative energy balance are then compared with the predictions from the full set of differential equations, which include the dissipative terms associated with air resistance and friction.

Beschreibung: Comment on "Storming the Bastille: the effect of electric fields on the ionospheric F-layer" by Rishbeth et al. (2010) Annales Geophysicae, 31, 145-150, 2013 Author(s): B. T. Tsurutani, A. J. Mannuccci, O. P. Verkhoglyadova, and G. S. Lakhina No abstract available.

Beschreibung: Geosynchronous magnetic field responses to fast solar wind dynamic pressure enhancements: MHD field model Annales Geophysicae, 30, 1285-1295, 2012 Author(s): T. R. Sun, C. Wang, N. L. Borodkova, and G. N. Zastenker We performed global MHD simulations of the geosynchronous magnetic field in response to fast solar wind dynamic pressure ( P d ) enhancements. Taking three P d enhancement events in 2000 as examples, we found that the main features of the total field B and the dominant component B z can be efficiently predicted by the MHD model. The predicted B and B z varies with local time, with the highest level near noon and a slightly lower level around mid-night. However, it is more challenging to accurately predict the responses of the smaller component at the geosynchronous orbit (i.e., B x and B y ). In contrast, the limitations of T01 model in predicting responses to fast P d enhancements are presented.

Beschreibung: Low-frequency magnetic field fluctuations in Earth's plasma environment observed by THEMIS Annales Geophysicae, 30, 1271-1283, 2012 Author(s): L. Guicking, K.-H. Glassmeier, H.-U. Auster, Y. Narita, and G. Kleindienst Low-frequency magnetic wave activity in Earth's plasma environment was determined based on a statistical analysis of THEMIS magnetic field data. We observe that the spatial distribution of low-frequency magnetic field fluctuations reveals highest values in the magnetosheath, but the observations differ qualitatively from observations at Venus presented in a previous study since significant wave activity at Earth is also observed in the nightside magnetosheath. Outside the magnetosheath the low-frequency wave activity level is generally very low. By means of an analytical streamline model for the magnetosheath plasma flow, we are able to investigate the spatial and temporal evolution of wave intensity along particular streamlines in order to characterise possible wave generation mechanisms. We observe a decay of wave intensity along the streamlines, but contrary to the situation at Venus, we obtain good qualitative agreement with the theoretical concept of freely evolving/decaying turbulence. Differences between the dawn region and the dusk region can be observed only further away from the magnetopause. We conclude that wave generation mechanisms may be primarily attributed to processes at or in the vicinity of the bow shock. The difference with the observations of the Venusian magnetosheath we interpret to be the result of the different types of solar wind interaction processes since the Earth possesses a global magnetic field while Venus does not, and therefore the observed magnetic wave activities may be caused by diverse magnetic field controlled characteristics of wave generation processes.

Beschreibung: Variations of ionospheric profile parameters during solar maximum and comparison with IRI-2007 over Chung-Li, Taiwan Annales Geophysicae, 30, 1249-1257, 2012 Author(s): Y. J. Chuo This paper studies the seasonal changes in the diurnal variation of ionospheric bottomside slab thickness (B0), based on observations during high solar activities at the equatorial ionization anomaly (EIA) area station of Chung-Li (121.1° E, 24.9° N), Taiwan. The data examined in this investigation are derived from ionograms recorded at Chung-Li in 1999, and are compared with International Reference Ionosphere (IRI-2007) model values. In our data set B0 shows largest values and biggest changes during the daytime (06:00–12:00 LT) particularly in the summer. Moreover, the diurnal variation of B0 shows an abnormal peak during the pre-sunrise period, especially in the winter. The variation in the F-peak height ( hm F2) is related to a thermospheric wind traveling toward the equator, which also enhances B0 during the pre-sunrise period. The results of the comparison with the IRI model show that B0 is overestimated, in both the B0-table and the Gulyaeva option, after noon LT in the equinox (B0-table and Gulyaeva average values for the overestimation are 11 and 47 km, respectively) and summer (B0-table and Gulyaeva average values for the overestimation are 23 and 71 km, respectively) periods. Furthermore, the modeled values are underestimated at approximately 31 and 14 km for the table and Gulyaeva option during the daytime in the winter, respectively. The F2-layer maximum electron density ( Nm F2) data show reasonably favorable agreement with the model for a high correlation coefficient of approximately 0.97, with the major difference observed at approximately noon in the equinox and winter seasons. Regarding the hm F2 data, the model shows agreement with the observed values, and the largest discrepancy (average value is 39 km) was observed in the summer and the smallest (average value is 11 km) in the equinox season. This paper provides a comprehensive discussion on the relationship among B0, the Nm F2 and the hm F2 for geomagnetic storm events.

Beschreibung: Features of highly structured equatorial plasma irregularities deduced from CHAMP observations Annales Geophysicae, 30, 1259-1269, 2012 Author(s): C. Xiong, H. Lühr, S. Y. Ma, C. Stolle, and B. G. Fejer In this study five years of CHAMP (Challenging Mini-satellite Payload) fluxgate magnetometer (FGM) data is used to investigate the characteristics of Equatorial Plasma Bubbles (EPBs). We filtered the FGM data by using band-passes with four different cut-off periods to get the EPBs with different maximum spatial scale sizes in the meridional plane ranging from 76–608 km. Associated with the EPB observations at about 400 km, the typical altitude of CHAMP during the year 2000–2005, we also investigate the post-sunset equatorial vertical plasma drift data from ROCSAT-1 (Republic of China Satellite 1). Since the height of the F-layer is highly correlated with the vertical plasma drift and solar flux, we sorted the ROCSAT-1 data into different groups by F10.7. From the integrated vertical drift we have estimated the post-sunset uplift of the ionosphere. By comparing the properties of EPB occurrence for different scale sizes with the global distribution of plasma vertical uplift, we have found that EPBs reaching higher altitudes are more structured than those which are sampled by CHAMP near the top side of the depleted fluxtube. Such a result is in accord with 3-D model simulations (Aveiro and Hysell, 2010). Small-scale EPB structures are observed by CHAMP when the irregularities reach apex heights of 800 km and more. Such events are encountered primarily in the Brazilian sector during the months around November, when the post-sunset vertical plasma drift is high.

Beschreibung: Intensification of dayside diffuse auroral precipitation: contribution of dayside Whistler-mode chorus waves in realistic magnetic fields Annales Geophysicae, 30, 1297-1307, 2012 Author(s): R. Shi, D. Han, B. Ni, Z.-J. Hu, C. Zhou, and X. Gu Compared to the recently improved understanding of nightside diffuse aurora, the mechanism(s) responsible for dayside diffuse aurora remains poorly understood. While dayside chorus has been thought as a potential major contributor to dayside diffuse auroral precipitation, quantitative analyses of the role of chorus wave scattering have not been carefully performed. In this study we investigate a dayside diffuse auroral intensification event observed by the Chinese Arctic Yellow River Station (YRS) all-sky imagers (ASI) on 7 January 2005 and capture a substantial increase in diffuse auroral intensity at the 557.7 nm wavelength that occurred over almost the entire ASI field-of-view near 09:24 UT, i.e., ~12:24 MLT. Computation of bounce-averaged resonant scattering rates by dayside chorus emissions using realistic magnetic field models demonstrates that dayside chorus scattering can produce intense precipitation losses of plasma sheet electrons on timescales of hours (even approaching the strong diffusion limit) over a broad range of both energy and pitch angle, specifically, from ~1 keV to 50 keV with equatorial pitch angles from the loss cone to up to ~85° depending on electron energy. Subsequent estimate of loss cone filling index indicates that the loss cone can be substantially filled, due to dayside chorus driven pitch angle scattering, at a rate of ≥0.8 for electrons from ~500 eV to 50 keV that exactly covers the precipitating electrons for the excitation of green-line diffuse aurora. Estimate of electron precipitation flux at different energy levels, based on loss cone filling index profile and typical dayside electron distribution observed by THEMIS spacecraft under similar conditions, gives a total precipitation electron energy flux of the order of 0.1 erg cm −2 s −1 with ~1 keV characteristic energy (especially when using T01s), which can be very likely to cause intense green-line diffuse aurora activity on the dayside. Therefore, dayside chorus scattering in the realistic magnetic field can greatly contribute to the YRS ASI observed intensification of dayside green-line aurora. Besides wave induced scattering and changes in the ambient magnetic field, variations in associated electron flux can also contribute to enhanced diffuse aurora emissions, the possibility of which we cannot exactly rule out due to lack of simultaneous observations of magnetospheric particles. Since the geomagnetic activity level was rather low during the period of interest, it is reasonable to infer that changes in the associated electron flux in the magnetosphere should be small, and consequently its contribution to the observed enhanced diffuse auroral activity should be small as well. Our results support the scenario that dayside chorus could play a major role in the production of dayside diffuse aurora, and also demonstrate that changes in magnetospheric magnetic field should be considered to reasonably interpret observations of dayside diffuse aurora.

Beschreibung: Oxygen ion energization by waves in the high altitude cusp and mantle Annales Geophysicae, 30, 1309-1314, 2012 Author(s): M. Waara, H. Nilsson, R. Slapak, M. André, and G. Stenberg We present a comparative study of low frequency electric field spectral densities and temperatures observed by the Cluster spacecraft in the high altitude cusp/mantle region. We compare the relation between the O + temperature and wave intensity at the oxygen gyrofrequency at each measurement point and find a clear correlation. The trend of the correlation agrees with the predictions by both an asymptotic mean-particle theory and a test-particle approach. The perpendicular to parallel temperature ratio is also consistent with the predictions of the asymptotic mean-particle theory. At times the perpendicular temperature is significantly higher than predicted by the models. A simple study of the evolution of the particle distributions (conics) at these altitudes indicates that enhanced perpendicular temperatures would be observed over many R E after heating ceases. Therefore, sporadic intense heating is the likely explanation for cases with high temperature and comparably low wave activity. We observe waves of sufficient amplitude to explain the highest observed temperatures, while the theory in general overestimates the temperature associated with the highest observed wave activity, indicating that such high wave activity is very sporadic.

Beschreibung: The "step feature" of suprathermal ion distributions: a discriminator between acceleration processes? Annales Geophysicae, 30, 1315-1319, 2012 Author(s): H. J. Fahr and H. Fichtner The discussion of exactly which process is causing the preferred build-up of v −5 -power law tails of the velocity distribution of suprathermal particles in the solar wind is still ongoing. Criteria allowing one to discriminate between the various suggestions that have been made would be useful in order to clarify the physics behind these tails. With this study, we draw the attention to the so-called "step feature" of the velocity distributions and offer a criterion that allows one to distinguish between those scenarios that employ velocity diffusion, i.e. second-order Fermi processes, which are prime candidates in the present debate. With an analytical approximation to the self-consistently obtained velocity diffusion coefficient, we solve the transport equation for suprathermal particles. The numerical simulation reveals that this form of the diffusion coefficient naturally leads to the step feature of the velocity distributions. This finding favours – at least in regions of the appearance of the step feature (i.e. for heliocentric distances up to about 11 AU and at lower energies) – the standard velocity diffusion as a consequence of the particle's interactions with the plasma wave turbulence as opposed to that caused by velocity fluctuation-induced compressions and rarefactions.

Beschreibung: Plasma parameter analysis of the Langmuir decay process via Particle-in-Cell simulations Annales Geophysicae, 30, 1169-1183, 2012 Author(s): M. A. Diaz, M. Zettergren, J. L. Semeter, and M. Oppenheim The beam-plasma mechanism, based on the Langmuir decay process, has been proposed to explain naturally enhanced ion-acoustic lines (NEIALs), which are spectral distortions in incoherent scatter radar (ISR) data frequently observed in the vicinity of auroral arcs. In this work the effect of the Langmuir decay process on the ISR spectrum is studied and compared with an analytical model for different plasma parameters by using an electrostatic parallel particle-in-cell (EPPIC) code. Simulations show that the code is working in accordance with theory for a wide range of beam and plasma values and that the features of the spectrum are sensitive to changes of those values. These results suggest that the EPPIC code might be used to build a spectrum-plasma parameter model which will allow estimation of beam and plasma parameters from observed spectra. Simulations also confirm that background electron density ( n e ) plays an important role in determining the maximum detectable wavenumber of the enhancement. Specifically, results demonstrate that an increase in n e makes the enhancements of the ion acoustic more likely line at large wavenumbers, a finding consistent with statistical studies showing more frequent NEIAL occurrence near solar maximum. Finally, the simulations expose some inaccuracies of the current theoretical model in quantifying the energy passed from the beam to the Langmuir waves as well as with the range of enhanced wavenumbers. These differences may be attributable to the weak Langmuir turbulent regime assumption used in the theory.

Beschreibung: Review Article: "The Lagrangian description of aperiodic flows: a case study of the Kuroshio Current" Nonlinear Processes in Geophysics, 19, 449-472, 2012 Author(s): C. Mendoza and A. M. Mancho This article reviews several recently developed Lagrangian tools and shows how their combined use succeeds in obtaining a detailed description of purely advective transport events in general aperiodic flows. In particular, because of the climate impact of ocean transport processes, we illustrate a 2-D application on altimeter data sets over the area of the Kuroshio Current, although the proposed techniques are general and applicable to arbitrary time dependent aperiodic flows. The first challenge for describing transport in aperiodical time dependent flows is obtaining a representation of the phase portrait where the most relevant dynamical features may be identified. areas that are related to confinement regions. This representation is accomplished by using global Lagrangian descriptors that when applied for instance to the altimeter data sets retrieve over the ocean surface a phase portrait where the geometry of interconnected dynamical systems is visible. The phase portrait picture is essential because it evinces which transport routes are acting on the whole flow. Once these routes are roughly recognised, it is possible to complete a detailed description by the direct computation of the finite time stable and unstable manifolds of special hyperbolic trajectories that act as organising centres of the flow.

Beschreibung: Spatio-temporal structure of Alfvén waves excited by a sudden impulse localized on an L-shell Annales Geophysicae, 30, 1099-1106, 2012 Author(s): D. Yu. Klimushkin, P. N. Mager, and K.-H. Glassmeier This paper is concerned with the spatial structure and temporal evolution of the azimuthally small scale Alfvén wave generated by a sudden impulse concentrated on a given magnetic shell. At the outset, both poloidal and toroidal components are present in the wave's magnetic field. The oscillation in the poloidal component on a given magnetic shell is a superposition of two monochromatic oscillations, one with the local resonance frequency on this shell, and the other with the frequency corresponding to the resonance frequency on the source surface. The superposition of these two oscillations leads to beating. Due to phase mixing, the poloidal component of the oscillation decreases with time down to zero, transferring its energy to the toroidal component. Beating in the toroidal component is less pronounced. As time elapses, energy concentration near the source magnetic shell occurs with the frequency of the oscillation corresponding to the Alfvénic resonance frequency on this surface. Outside this thin region wave amplitudes become rather small at oscillation frequencies corresponding to the local resonance frequency of the respective magnetic shell.

Beschreibung: Predictive model of magnetosheath plasma flow and its validation against Cluster and THEMIS data Annales Geophysicae, 30, 973-982, 2012 Author(s): J. Soucek and C. P. Escoubet An analytical model of magnetosheath plasma flow is described and compared with a large dataset of magnetosheath ion flow velocity measurements from Cluster and THEMIS spacecraft. The model is based on previous works by Kobel and Flückiger (1994) and Génot et al. (2011) and has been modified to overcome the restrictions of these models on the shape of model magnetopause and bow shock. Our model is compatible with any parabolic bow shock model and arbitrary magnetopause model. The model is relatively simple to implement and computationally inexpensive, and its only inputs are upstream solar wind parameters. Comparison with observed data yields a good correspondence: median error in the direction of flow velocity is comparable with the instrumental error, and flow magnitude is predicted with a reasonable accuracy (relative error in flow speed was less than 25% for 86.5% of observations).

Beschreibung: Academician S. N. Vernov and cosmic ray variations research in the former USSR Astrophysics and Space Sciences Transactions, 8, 25-28, 2012 Author(s): L. Dorman During my long scientific career of about 60 years (from 1951), I was happy to meet and cooperate with Academician S. N. Vernov. During the first step (1950–1952), S. N. Vernov together with N. V. Pushkov organized the first Soviet network of cosmic ray (CR) stations equipped by large automatically operating ionization chambers of A. Compton type shielded by 10 cm of Pb. This was the beginning of the CR variations research in the former USSR. In the present paper I consider step by step the development of experimental and theoretical basis of CR variations research in the former USSR and the great role of S. N. Vernov in this process.

Beschreibung: Cosmic rays and stochastic magnetic reconnection in the heliotail Nonlinear Processes in Geophysics, 19, 351-364, 2012 Author(s): P. Desiati and A. Lazarian Galactic cosmic rays are believed to be generated by diffusive shock acceleration processes in Supernova Remnants, and the arrival direction is likely determined by the distribution of their sources throughout the Galaxy, in particular by the nearest and youngest ones. Transport to Earth through the interstellar medium is expected to affect the cosmic ray properties as well. However, the observed anisotropy of TeV cosmic rays and its energy dependence cannot be explained with diffusion models of particle propagation in the Galaxy. Within a distance of a few parsec, diffusion regime is not valid and particles with energy below about 100 TeV must be influenced by the heliosphere and its elongated tail. The observation of a highly significant localized excess region of cosmic rays from the apparent direction of the downstream interstellar flow at 1–10 TeV energies might provide the first experimental evidence that the heliotail can affect the transport of energetic particles. In particular, TeV cosmic rays propagating through the heliotail interact with the 100–300 AU wide magnetic field polarity domains generated by the 11 yr cycles. Since the strength of non-linear convective processes is expected to be larger than viscous damping, the plasma in the heliotail is turbulent. Where magnetic field domains converge on each other due to solar wind gradient, stochastic magnetic reconnection likely occurs. Such processes may be efficient enough to re-accelerate a fraction of TeV particles as long as scattering processes are not strong. Therefore, the fractional excess of TeV cosmic rays from the narrow region toward the heliotail direction traces sightlines with the lowest smearing scattering effects, that can also explain the observation of a harder than average energy spectrum.

Beschreibung: Calibration of a radiocarbon age Nonlinear Processes in Geophysics, 19, 345-350, 2012 Author(s): D. J. Keenan The calibration of a radiocarbon age to a calendar date is reviewed. It is shown that the commonly-used programs for calibration sometimes give results that are significantly in error.