Wave-activity conservation laws and stability theorems for semi-geostrophic dynamics: part 1. pseudomomentum-based theory

There exists a well-developed body of theory based on quasi-geostrophic (QG) dynamics that is central to our present understanding of large-scale atmospheric and oceanic dynamics. An important question is the extent to which this body of theory may generalize to more accurate dynamical models. As a first step in this process, we here generalize a set of theoretical results, concerning the evolution of disturbances to prescribed basic states, to semi-geostrophic (SG) dynamics. SG dynamics, like QG dynamics, is a Hamiltonian balanced model whose evolution is described by the material conservation of potential vorticity, together with an invertibility principle relating the potential vorticity to the advecting fields. SG dynamics has features that make it a good prototype for balanced models that are more accurate than QG dynamics. In the first part of this two-part study, we derive a pseudomomentum invariant for the SG equations, and use it to obtain: (i) linear and nonlinear generalized Charney–Stern theorems for disturbances to parallel flows; (ii) a finite-amplitude local conservation law for the invariant, obeying the group-velocity property in the WKB limit; and (iii) a wave-mean-flow interaction theorem consisting of generalized Eliassen–Palm flux diagnostics, an elliptic equation for the stream-function tendency, and a non-acceleration theorem. All these results are analogous to their QG forms. The pseudomomentum invariant – a conserved second-order disturbance quantity that is associated with zonal symmetry – is constructed using a variational principle in a similar manner to the QG calculations. Such an approach is possible when the equations of motion under the geostrophic momentum approximation are transformed to isentropic and geostrophic coordinates, in which the ageostrophic advection terms are no longer explicit. Symmetry-related wave-activity invariants such as the pseudomomentum then arise naturally from the Hamiltonian structure of the SG equations. We avoid use of the so-called ‘massless layer’ approach to the modelling of isentropic gradients at the lower boundary, preferring instead to incorporate explicitly those boundary contributions into the wave-activity and stability results. This makes the analogy with QG dynamics most transparent. This paper treats the f-plane Boussinesq form of SG dynamics, and its recent extension to β-plane, compressible flow by Magnusdottir & Schubert. In the limit of small Rossby number, the results reduce to their respective QG forms. Novel features particular to SG dynamics include apparently unnoticed lateral boundary stability criteria in (i), and the necessity of including additional zonal-mean eddy correlation terms besides the zonal-mean potential vorticity fluxes in the wave-mean-flow balance in (iii). In the companion paper, wave-activity conservation laws and stability theorems based on the SG form of the pseudoenergy are presented.

Default perception of high-speed motion

When human observers are exposed to even slight motion signals followed by brief visual transients—stimuli containing no detectable coherent motion signals—they perceive large and salient illusory jumps. This novel effect, which we call “high phi”, challenges well-entrenched assumptions about the perception of motion, namely the minimal-motion principle and the breakdown of coherent motion perception with steps above an upper limit. Our experiments with transients such as texture randomization or contrast reversal show that the magnitude of the jump depends on spatial frequency and transient duration, but not on the speed of the inducing motion signals, and the direction of the jump depends on the duration of the inducer. Jump magnitude is robust across jump directions and different types of transient. In addition, when a texture is actually displaced by a large step beyond dmax, a breakdown of coherent motion perception is expected, but in the presence of an inducer observers again perceive coherent displacements at or just above dmax. In sum, across a large variety of stimuli, we find that when incoherent motion noise is preceded by a small bias, instead of perceiving little or no motion, as suggested by the minimal-motion principle, observers perceive jumps whose amplitude closely follows their own dmax limits.

Does grammatical aspect affect motion event cognition? A cross-linguistic comparison of English and Swedish speakers

In this article, we explore whether cross-linguistic differences in grammatical aspect encoding may give rise to differences in memory and cognition. We compared native speakers of two languages that encode aspect differently (English and Swedish) in four tasks that examined verbal descriptions of stimuli, online triads matching, and memory-based triads matching with and without verbal interference. Results showed between-group differences in verbal descriptions and in memory-based triads matching. However, no differences were found in online triads matching and in memory-based triads matching with verbal interference. These findings need to be interpreted in the context of the overall pattern of performance, which indicated that both groups based their similarity judgments on common perceptual characteristics of motion events. These results show for the first time a cross-linguistic difference in memory as a function of differences in grammatical aspect encoding, but they also contribute to the emerging view that language fine tunes rather than shapes perceptual processes that are likely to be universal and unchanging.

The role of the harmonic vector average in motion integration

The local speeds of object contours vary systematically with the cosine of the angle between the normal component of the local velocity and the global object motion direction. An array of Gabor elements whose speed changes with local spatial orientation in accordance with this pattern can appear to move as a single surface. The apparent direction of motion of plaids and Gabor arrays has variously been proposed to result from feature tracking, vector addition and vector averaging in addition to the geometrically correct global velocity as indicated by the intersection of constraints (IOC) solution. Here a new combination rule, the harmonic vector average (HVA), is introduced, as well as a new algorithm for computing the IOC solution. The vector sum can be discounted as an integration strategy as it increases with the number of elements. The vector average over local vectors that vary in direction always provides an underestimate of the true global speed. The HVA, however, provides the correct global speed and direction for an unbiased sample of local velocities with respect to the global motion direction, as is the case for a simple closed contour. The HVA over biased samples provides an aggregate velocity estimate that can still be combined through an IOC computation to give an accurate estimate of the global velocity, which is not true of the vector average. Psychophysical results for type II Gabor arrays show perceived direction and speed falls close to the IOC direction for Gabor arrays having a wide range of orientations but the IOC prediction fails as the mean orientation shifts away from the global motion direction and the orientation range narrows. In this case perceived velocity generally defaults to the HVA.

Dynamic portrait of the planetary 2/1 mean-motion resonance - II. Systems with a more massive inner planet

This paper presents the second part in our study of the global structure of the planar phase space of the planetary three-body problem, when both planets lie in the vicinity of a 2/1 mean-motion resonance. While Paper I was devoted to cases where the outer planet is the more massive body, the present work is devoted to the cases where the more massive body is the inner planet. As before, outside the well-known Apsidal Corotation Resonances (ACR), the phase space shows a complex picture marked by the presence of several distinct regimes of resonant and non-resonant motion, crossed by families of periodic orbits and separated by chaotic zones. When the chosen values of the integrals of motion lead to symmetric ACR, the global dynamics are generally similar to the structure presented in Paper I. However, for asymmetric ACR the resonant phase space is strikingly different and shows a galore of distinct dynamical states. This structure is shown with the help of dynamical maps constructed on two different representative planes, one centred on the unstable symmetric ACR and the other on the stable asymmetric equilibrium solution. Although the study described in the work may be applied to any mass ratio, we present a detailed analysis for mass values similar to the Jupiter-Saturn case. Results give a global view of the different dynamical states available to resonant planets with these characteristics. Some of these dynamical paths could have marked the evolution of the giant planets of our Solar system, assuming they suffered a temporary capture in the 2/1 resonance during the latest stages of the formation of our Solar system.

Bicomplexes and conservation laws in non-Abelian Toda models

A bicomplex structure is associated with the Leznov-Saveliev equation of integrable models. The linear problem associated with the zero-curvature condition is derived in terms of the bicomplex linear equation. The explicit example of a non-Abelian conformal affine Toda model is discussed in detail and its conservation laws are derived from the zero-curvature representation of its equation of motion.

Dynamics of two satellites in the 2/1 Mean-Motion resonance: application to the case of Enceladus and Dione

The dynamics of a pair of satellites similar to Enceladus-Dione is investigated with a two-degrees-of-freedom model written in the domain of the planar general three-body problem. Using surfaces of section and spectral analysis methods, we study the phase space of the system in terms of several parameters, including the most recent data. A detailed study of the main possible regimes of motion is presented, and in particular we show that, besides the two separated resonances, the phase space is replete of secondary resonances.

Basins of attraction changes by amplitude constraining of oscillators with limited power supply

We investigate the dynamics of a Duffing oscillator driven by a limited power supply, such that the source of forcing is considered to be another oscillator, coupled to the first one. The resulting dynamics come from the interaction between both systems. Moreover, the Duffing oscillator is subjected to collisions with a rigid wall (amplitude constraint). Newtonian laws of impact are combined with the equations of motion of the two coupled oscillators. Their solutions in phase space display periodic (and chaotic) attractors, whose amplitudes, especially when they are too large, can be controlled by choosing the wall position in suitable ways. Moreover, their basins of attraction are significantly modified, with effects on the final state system sensitivity. (c) 2005 Elsevier Ltd. All rights reserved.

Dynamics of Enceladus and Dione inside the 2:1 mean-motion resonance under tidal dissipation

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Intra-articular pressure profiles of the cadaveric equine fetlock joint in motion

The study of the influence of motion and initial intra-articular pressure (IAP) on intra-articular pressure profiles in equine cadaver metatarsophalangeal (MTP) joints was undertaken as a prelude to in vivo studies, Eleven equine cadaver MTP joints were submitted to 2 motion frequencies of 5 and 10 cycles/min of flexion and extension, simulating the condition of lower and higher (double) rates of passive motion. These frequencies were applied and pressure profiles generated with initial normal intra-articular pressure (-5 mmHg) and subsequently 30 mmHg intra-articular pressure obtained by injection of previously harvested synovial fluid.The 4 trials performed were 1) normal IAP; 5 cyles/min; 2) normal IAP; 10 cycles/min; 3) IAP at 30 mmHg; 5 cycles/min and 4) IAP at 30 mmHg; 10 cycles/min. The range of joint motion applied (mean +/- s.e.) was 67.6 +/- 1.61 degrees with an excursion from 12.2 +/- 1.2 degrees in extension to 56.2 +/- 2.6 degrees in flexion, Mean pressure recorded in mmHg for the first and last min of each trial, respectively, were 1) -5.7 +/- 0.9 and -6.3 +/- 1.1; 2) -5.3 +/- 1.1 and -6.2 +/- 1.1; 3) 58.8 +/- 8.0 and 42.3 +/- 7.2; 4) 56.6 +/- 3.7 and 40.3 +/- 4.6. Statistical analyses showed a trend for difference between the values for the first and last minute in trial 3 (0.05>P<0.1) with P = 0.1 and significant difference (P = 0.02) between the mean IAP of the first and last min in trial 4. The loss of intra-articular pressure associated with time and motion was 10.5, 16.9, 28.1 and 28.9% for trials 1-4, respectively. As initial intraarticular pressure and motion increased, the percent loss of intra-articular pressure increased.The angle of lowest pressure was 12.2 +/- 1.2 (mean +/- s.e.) in extension in trials 1 and 2, In trials 3 and 4, the lowest pressures were obtained in flexion with the joints at 18.5 +/- 2.0 degrees (mean +/- s.e.). This demonstrated that the joint angle of least pressure changed as the initial intra-articular pressure changed and there would not be a single angle of least pressure for a given joint.The volume of synovial fluid recovered from the MTP joints in trial 3 compared to 4 (trials in which fluid was injected to attain IAP of 30 mmHg) was not significantly different, supporting a soft tissue compliance change as a cause for the significant loss of intra-articular pressure during the 15 min of trial 4.The pressure profiles generated correlate well with in vivo values and demonstrated consistent pressure profiles. Our conclusions are summarised as follows:1. Clinically normal equine MTP joints which were frozen and then later thawed were found to have mostly negative baseline intra-articular pressures, as would be expected in living subjects,2. Alternate pressure profiles of the dorsal and plantar pouch at baseline intra-articular pressure document the presence of pressure forces that would support 'back and forth' fluid movement between joint compartments. This should result in movement of joint fluid during motion, assisting in lubrication and nutrition of articular cartilage,3. If joint pressure was initially greater than normal (30 mmHg), as occurs in diseased equine MTP joints, joint motion further increased joint capsule relaxation (compliance) and, therefore, reduced intra-articular pressure.4. Peak intra-articular pressures reached extremely high values (often >100 mmHg) in flexion when initial pressure was 30 mmHg. Joint effusion pressures recorded for clinical MCP joints are frequently 30 mmHg. These IAP values are expected to produce intermittent synovial ischaemia in clinical cases during joint flexion.5, Additional in vivo studies are necessary to confirm our conclusions from this study and to identify the contributions of fluid absorption and the presence of ischaemia in a vascularised joint.

Conservation laws for non-potential operators

A study was developed in order to build a function M invariant in time, by means of Hamiltonian's formulation, taking into account the equation associated to the problem, showing that starting from this function the equation of motion of the system with the contour conditions for non-conservative considered problems can be obtained. The Hamiltonian method is extended for these kind of systems in order to validate for non-potential operators through variational approach.

Coordination of spinal motion in the transverse and frontal planes during walking in people with and without recurrent low back pain

STUDY DESIGN. Observational cohort study. OBJECTIVE. To investigate spinal coordination during preferred and fast speed walking in pain-free subjects with and without a history of recurrent low back pain (LBP). SUMMARY OF BACKGROUND DATA. Dynamic motion of the spine during walking is compromised in the presence of back pain (LBP), but its analysis often presents some challenges. The coexistence of significant symptoms may change gait because of pain or adaptation of the musculoskeletal structures or both. A history of LBP without the overlay of a current symptomatic episode allows a better model in which to explore the impact on spinal coordination during walking. METHODS. Spinal and lower limb segmental motions were tracked using electromagnetic sensors. Analyses were conducted to explore the synchrony and spatial coordination of the segments and to compare the control and subjects with LBP. RESULTS. We found no apparent differences between the groups for either overall amplitude of motion or most indicators of coordination in the lumbar region; however, there were significant postural differences in the mid-stance phase and other indicators of less phase locking in controls compared with subjects with LBP. The lower thoracic spinal segment was more affected by the history of back pain than the lumbar segment. CONCLUSION. Although small, there were indicators that alterations in spinal movement and coordination in subjects with recurrent LBP were due to adaptive changes rather than the presence of pain. © 2013, Lippincott Williams & Wilkins.

Subject-specific musculoskeletal models of the lower limbs for the prediction of skeletal loads during motion

The determination of skeletal loading conditions in vivo and their relationship to the health of bone tissues, remain an open question. Computational modeling of the musculoskeletal system is the only practicable method providing a valuable approach to muscle and joint loading analyses, although crucial shortcomings limit the translation process of computational methods into the orthopedic and neurological practice. A growing attention focused on subject-specific modeling, particularly when pathological musculoskeletal conditions need to be studied. Nevertheless, subject-specific data cannot be always collected in the research and clinical practice, and there is a lack of efficient methods and frameworks for building models and incorporating them in simulations of motion. The overall aim of the present PhD thesis was to introduce improvements to the state-of-the-art musculoskeletal modeling for the prediction of physiological muscle and joint loads during motion. A threefold goal was articulated as follows: (i) develop state-of-the art subject-specific models and analyze skeletal load predictions; (ii) analyze the sensitivity of model predictions to relevant musculotendon model parameters and kinematic uncertainties; (iii) design an efficient software framework simplifying the effort-intensive phases of subject-specific modeling pre-processing. The first goal underlined the relevance of subject-specific musculoskeletal modeling to determine physiological skeletal loads during gait, corroborating the choice of full subject-specific modeling for the analyses of pathological conditions. The second goal characterized the sensitivity of skeletal load predictions to major musculotendon parameters and kinematic uncertainties, and robust probabilistic methods were applied for methodological and clinical purposes. The last goal created an efficient software framework for subject-specific modeling and simulation, which is practical, user friendly and effort effective. Future research development aims at the implementation of more accurate models describing lower-limb joint mechanics and musculotendon paths, and the assessment of an overall scenario of the crucial model parameters affecting the skeletal load predictions through probabilistic modeling.

Effects of anatomic conformation on three-dimensional motion of the caudal lumbar and lumbosacral portions of the vertebral column of dogs

OBJECTIVE: To determine the association between the 3-dimensional (3-D) motion pattern of the caudal lumbar and lumbosacral portions of the canine vertebral column and the morphology of vertebrae, facet joints, and intervertebral disks. SAMPLE POPULATION: Vertebral columns of 9 German Shepherd Dogs and 16 dogs of other breeds with similar body weights and body conditions. PROCEDURE: Different morphometric parameters of the vertebral column were assessed by computed tomography (CT) and magnetic resonance imaging. Anatomic conformation and the 3-D motion pattern were compared, and correlation coefficients were calculated. RESULTS: Total range of motion for flexion and extension was mainly associated with the facet joint angle, the facet joint angle difference between levels of the vertebral column in the transverse plane on CT images, disk height, and lever arm length. CONCLUSIONS AND CLINICAL RELEVANCE: Motion is a complex process that is influenced by the entire 3-D conformation of the lumbar portion of the vertebral column. In vivo dynamic measurements of the 3-D motion pattern of the lumbar and lumbosacral portions of the vertebral column will be necessary to further assess biomechanics that could lead to disk degeneration in dogs.

Application and assessment of a robust elastic motion correction algorithm to dynamic MRI

The purpose of this study was to assess the performance of a new motion correction algorithm. Twenty-five dynamic MR mammography (MRM) data sets and 25 contrast-enhanced three-dimensional peripheral MR angiographic (MRA) data sets which were affected by patient motion of varying severeness were selected retrospectively from routine examinations. Anonymized data were registered by a new experimental elastic motion correction algorithm. The algorithm works by computing a similarity measure for the two volumes that takes into account expected signal changes due to the presence of a contrast agent while penalizing other signal changes caused by patient motion. A conjugate gradient method is used to find the best possible set of motion parameters that maximizes the similarity measures across the entire volume. Images before and after correction were visually evaluated and scored by experienced radiologists with respect to reduction of motion, improvement of image quality, disappearance of existing lesions or creation of artifactual lesions. It was found that the correction improves image quality (76% for MRM and 96% for MRA) and diagnosability (60% for MRM and 96% for MRA).