Analysis of Modulated Multivariate Oscillations

The concept of a common modulated oscillation spanning multiple time series is formalized, a method for the recovery of such a signal from potentially noisy observations is proposed, and the time-varying bias properties of the recovery method are derived. The method, an extension of wavelet ridge analysis to the multivariate case, identifies the common oscillation by seeking, at each point in time, a frequency for which a bandpassed version of the signal obtains a local maximum in power. The lowest-order bias is shown to involve a quantity, termed the instantaneous curvature, which measures the strength of local quadratic modulation of the signal after demodulation by the common oscillation frequency. The bias can be made to be small if the analysis filter, or wavelet, can be chosen such that the signal's instantaneous curvature changes little over the filter time scale. An application is presented to the detection of vortex motions in a set of freely-drifting oceanographic instruments tracking the ocean currents

Modulated Oscillations in Three Dimensions

The analysis of the fully three-dimensional and time-varying polarization characteristics of a modulated trivariate, or three-component, oscillation is addressed. The use of the analytic operator enables the instantaneous three-dimensional polarization state of any square-integrable trivariate signal to be uniquely defined. Straightforward expressions are given which permit the ellipse parameters to be recovered from data. The notions of instantaneous frequency and instantaneous bandwidth, generalized to the trivariate case, are related to variations in the ellipse properties. Rates of change of the ellipse parameters are found to be intimately linked to the first few moments of the signal's spectrum, averaged over the three signal components. In particular, the trivariate instantaneous bandwidth---a measure of the instantaneous departure of the signal from a single pure sinusoidal oscillation---is found to contain five contributions: three essentially two-dimensional effects due to the motion of the ellipse within a fixed plane, and two effects due to the motion of the plane containing the ellipse. The resulting analysis method is an informative means of describing nonstationary trivariate signals, as is illustrated with an application to a seismic record.Comment: IEEE Transactions on Signal Processing, 201

Multi-scaled analysis of the damped dynamics of an elastic rod with an essentially nonlinear end attachment

We study multi-frequency transitions in the transient dynamics of a viscously damped dispersive finite rod with an essentially nonlinear end attachment. The attachment consists of a small mass connected to the rod by means of an essentially nonlinear stiffness in parallel to a viscous damper. First, the periodic orbits of the underlying hamiltonian system with no damping are computed, and depicted in a frequency–energy plot (FEP). This representation enables one to clearly distinguish between the different types of periodic motions, forming back bone curves and subharmonic tongues. Then the damped dynamics of the system is computed; the rod and attachment responses are initially analyzed by the numerical Morlet wavelet transform (WT), and then by the empirical mode decomposition (EMD) or Hilbert–Huang transform (HTT), whereby, the time series are decomposed in terms of intrinsic mode functions (IMFs) at different characteristic time scales (or, equivalently, frequency scales). Comparisons of the evolutions of the instantaneous frequencies of the IMFs to the WT spectra of the time series enables one to identify the dominant IMFs of the signals, as well as, the time scales at which the dominant dynamics evolve at different time windows of the responses; hence, it is possible to reconstruct complex transient responses as superposition of the dominant IMFs involving different time scales of the dynamical response. Moreover, by superimposing the WT spectra and the instantaneous frequencies of the IMFs to the FEPs of the underlying hamiltonian system, one is able to clearly identify the multi-scaled transitions that occur in the transient damped dynamics, and to interpret them as ‘jumps’ between different branches of periodic orbits of the underlying hamiltonian system. As a result, this work develops a physics-based, multi-scaled framework and provides the necessary computational tools for multi-scaled analysis of complex multi-frequency transitions of essentially nonlinear dynamical systems

Horizontal rotation signals detected by "G-Pisa" ring laser for the Mw=9.0, March 2011, Japan earthquake

We report the observation of the ground rotation induced by the Mw=9.0, 11th of March 2011, Japan earthquake. The rotation measurements have been conducted with a ring laser gyroscope operating in a vertical plane, thus detecting rotations around the horizontal axis. Comparison of ground rotations with vertical accelerations from a co-located force-balance accelerometer shows excellent ring laser coupling at periods longer than 100s. Under the plane wave assumption, we derive a theoretical relationship between horizontal rotation and vertical acceleration for Rayleigh waves. Due to the oblique mounting of the gyroscope with respect to the wave direction-of-arrival, apparent velocities derived from the acceleration / rotation rate ratio are expected to be always larger than, or equal to the true wave propagation velocity. This hypothesis is confirmed through comparison with fundamental-mode, Rayleigh wave phase velocities predicted for a standard Earth model.Comment: Accepted for publication in Journal of Seismolog

A New Strategy for Deep Wide-Field High Resolution Optical Imaging

We propose a new strategy for obtaining enhanced resolution (FWHM = 0.12 arcsec) deep optical images over a wide field of view. As is well known, this type of image quality can be obtained in principle simply by fast guiding on a small (D = 1.5m) telescope at a good site, but only for target objects which lie within a limited angular distance of a suitably bright guide star. For high altitude turbulence this 'isokinetic angle' is approximately 1 arcminute. With a 1 degree field say one would need to track and correct the motions of thousands of isokinetic patches, yet there are typically too few sufficiently bright guide stars to provide the necessary guiding information. Our proposed solution to these problems has two novel features. The first is to use orthogonal transfer charge-coupled device (OTCCD) technology to effectively implement a wide field 'rubber focal plane' detector composed of an array of cells which can be guided independently. The second is to combine measured motions of a set of guide stars made with an array of telescopes to provide the extra information needed to fully determine the deflection field. We discuss the performance, feasibility and design constraints on a system which would provide the collecting area equivalent to a single 9m telescope, a 1 degree square field and 0.12 arcsec FWHM image quality.Comment: 46 pages, 22 figures, submitted to PASP, a version with higher resolution images and other supplementary material can be found at http://www.ifa.hawaii.edu/~kaiser/wfhr

Kinematics nomenclature for physiological accelerations with special reference to vestibular applications

Kinematics nomenclature for physiological accelerations and special reference to vestibular apparatu