Models of interacting pairs of thin,quasi-geostrophic vortices: steady-state solutions and nonlinear stability

We study pairwise interactions of elliptical quasi-geostrophic (QG) vortices as the limiting case of vanishingly thin uniform potential vorticity ellipsoids. In this limit, the product of the vertical extent of the ellipsoid and the potential vorticity within it is held fixed to a finite non-zero constant. Such elliptical “lenses” inherit the property that, in isolation, they steadily rotate without changing shape. Here, we use this property to extend both standard moment models and Hamiltonian ellipsoidal models to approximate the dynamical interaction of such elliptical lenses. By neglecting non-elliptical deformations, the simplified models reduce the dynamics to just four degrees of freedom per vortex. For simplicity, we focus on pairwise interactions between identical elliptical vortices initially separated in both the horizontal and vertical directions. The dynamics of the simplified models are compared with the full QG dynamics of the system, and show good agreement as expected for sufficiently distant lenses. The results reveal the existence of families of steadily rotating equilibria in the initial horizontal and vertical separation parameter space. For sufficiently large vertical separations, equilibria with varying shape exist for all horizontal separations. Below a critical vertical separation (stretched by the constant ratio of buoyancy to Coriolis frequencies N / f), comparable to the mean radius of either vortex, a gap opens in horizontal separation where no equilibria are possible. Solutions near the edge of this gap are unstable. In the full QG system, equilibria at the edge of the gap exhibit corners (infinite curvature) along their boundaries. Comparisons of the model results with the full nonlinear QG evolution show that the early stages of the instability are captured by the Hamiltonian elliptical model but not by the moment model that inaccurately estimates shorter-range interactions.     

A low temperature X-ray single-crystal diffraction and polarised infra-red study of epidote

The effects of low-temperature on the crystal structure of a natural epidote [Ca_1.925Fe_0.745Al_2.265Ti_0.004Si_3.037O₁₂(OH), a = 8.8924(7), b = 5.6214(3), c = 10.1547(6)? and β = 115.396(8)° at room conditions, Sp. Gr. P2₁ /m] have been investigated with a series of structure refinements down to 100 K on the basis of X-ray single-crystal diffraction data. The reflection conditions confirm that the space group is maintained within the T-range investigated. Structural refinements at all temperatures show the presence of Fe³⁺ at the octahedral M(3) site only [%Fe(M3) = 70.6(4)% at 295 K]. Only one independent proton site was located and two possible H-bonds occur, with O(10) as donor and O(4) and O(2) as acceptors. The H-bonding scheme is maintained down to 100 K and is supported by single crystal room-T polarised FTIR data. FTIR Spectra over the region 4,000–2,500 cm⁻¹ are dominated by the presence of a strongly pleochroic absorption feature which can be assigned to protonation of O(10)–O(4). Previously unobserved splitting of this absorption features is consistent with a NNN influence due to the presence of Al and Fe³⁺ on the nearby M(3) site. An additional relatively minor absorption feature in FTIR spectra can be tentatively assigned to protonation of O(10)–O(2). Low-T does not affect significantly the tetrahedral and octahedral bond distances and angles, even when distances are corrected for “rigid body motions”. A more significant effect is observed for the bond distances of the distorted Ca(1)- and Ca(2)-polyhedra, especially when corrected for “non-correlated motion”. The main low-T effect is observed on the vibrational regime of the atomic sites, and in particular for the two Ca-sites. A significant reduction of the magnitude of the thermal displacement ellipsoids, with a variation of U _eq (defined as one-third of the trace of the orthogonalised U _ij tensor) by ~40% is observed for the Ca-sites between 295 and 100 K. Within the same T-range, the U _eq of the octahedral and oxygen sites decrease similarly by ~35%, whereas those of the tetrahedral cations by ~22%.     

HTP21/c–C2/c phase transition and kinetics of Fe2+–Mg order–disorder of an Fe-poor pigeonite: implications for the cooling history of ureilites

A natural Ca-poor pigeonite (Wo₆En₇₆Fs₁₈) from the ureilite meteorite sample PCA82506-3, free of exsolved augite, was studied by in situ high-temperature single-crystal X-ray diffraction. The sample, monoclinic P2₁/c, was annealed up to 1,093°C to induce a phase transition from P2₁/c to C2/c symmetry. The variation with increasing temperature of the lattice parameters and of the intensity of the b-type reflections (h + k = 2n + 1, present only in the P2₁/c phase) showed a displacive phase transition P2₁/c to C2/c at a transition temperature T _Tr = 944°C, first order in character. The Fe–Mg exchange kinetics was studied by ex situ single-crystal X-ray diffraction in a range of temperatures between the closure temperature of the Fe–Mg exchange reaction and the transition temperature. Isothermal disordering annealing experiments, using the IW buffer, were performed on three crystals at 790, 840 and 865°C. Linear regression of ln k _D versus 1/T yielded the following equation: ln k_\textD = - 3717( ±416)/T(K) + 1.290( ±0.378);    (R² = 0.988) \ln \,k_{\text{D}} = - 3717( \pm 416)/T(K) + 1.290( \pm 0.378);\quad (R^{2} = 0.988) . The closure temperature (T _c) calculated using this equation was ∼740(±30)°C. Analysis of the kinetic data carried out taking into account the e.s.d.'s of the atomic fractions used to define the Fe–Mg degree of order, performed according to Mueller’s model, allowed us to retrieve the disordering rate constants C ₀ K _dis⁺ for all three temperatures yielding the following Arrhenius relation: ln( C₀ K_\textdis ⁺ ) = ln K₀ - Q/(RT) = 20.99( ±3.74) - 26406( ±4165)/T(K);    (R² = 0.988) \ln \left( {C_{0} K_{\text{dis}}^{ + } } \right) = \ln \,K_{0} - Q/(RT) = 20.99( \pm 3.74) - 26406( \pm 4165)/T(K);\quad (R^{2} = 0.988) . An activation energy of 52.5(±4) kcal/mol for the Fe–Mg exchange process was obtained. The above relation was used to calculate the following Arrhenius relation modified as a function of X _Fe (in the range of X _Fe = 0.20–0.50): ln( C₀ K_\textdis ⁺ ) = (21.185 - 1.47X_\textFe ) - \frac(27267 - 4170X_\textFe )T(K) \ln \left( {C_{0} K_{\text{dis}}^{ + } } \right) = (21.185 - 1.47X_{\text{Fe}} ) - {\frac{{(27267 - 4170X_{\text{Fe}} )}}{T(K)}} . The cooling time constant, η = 6 × 10⁻¹ K⁻¹ year⁻¹ calculated on the PCA82506-3 sample, provided a cooling rate of the order of 1°C/min consistent with the extremely fast late cooling history of the ureilite parent body after impact excavation.     

Magnetic flares and the optical variability of the X-ray transient XTE J1118+480

The simultaneous presence of a strong quasi-periodic oscillation, of period ∼10 s, in the optical and X-ray light curves of the X-ray transient XTE J1118+480 suggests that a significant fraction of the optical flux originates from the inner part of the accretion flow, where most of the X-rays are produced. We present a model of magnetic flares in an accretion disc corona where thermal cyclo-synchrotron emission contributes significantly to the optical emission, while the X-rays are produced by inverse Compton scattering of the soft photons produced by dissipation in the underlying disc and by the synchrotron process itself. Given the observational constraints, we estimate the values for the coronal temperature, optical depth and magnetic field intensity, as well as the accretion rate for the source. Within our model we predict a correlation between optical and hard X-ray variability and an anticorrelation between optical and soft X-rays. We also expect optical variability on flaring time-scales (∼tens of ms), with a power-density spectrum similar to that observed in the X-ray band. Finally, we use both the available optical/extreme-ultraviolet/X-ray spectral energy distribution and the low-frequency time variability to discuss limits on the inner radius of the optically thick disc.     

Runoff-generated debris flows: Observation of initiation conditions and erosion–deposition dynamics along the channel at Cancia (eastern Italian Alps)

In the Dolomitic region, abundant coarse hillslope sediment is commonly found at the toe of rocky cliffs. Ephemeral channels originate where lower permeability bedrock surfaces concentrate surface runoff. Debris flows initiate along such channels following intense rainfall and determine the progressive erosion and deepening of the channels. Sediment recharge mechanisms include rock fall, dry ravel processes and channel-bank failures. Here we document debris flow activity that took place in an active debris flow basin during the year 2015. The Cancia basin is located on the southwestern slope of Mount Antelao (3264 m a.s.l.) in the dolomitic region of the eastern Italian Alps. The 2.5 km² basin is incised in dolomitic limestone rocks. The data consist of repeated topographic surveys, distributed rainfall measurements, time-lapse (2 s) videos of two events and pore pressure measurements in the channel bed. During July and August 2015, two debris flow events occurred, following similarly intense rainstorms. We compared rainfall data to existing rainfall triggering thresholds and simulated the hydrological response of the headwater catchment with a distributed model in order to estimate the total and peak water discharge. Our data clearly illustrate how debris entrainment along the channel is the main contributor to the overall mobilized volume and that erosion is dominant when the channel slope exceeds 16°. Further downstream, sediment accumulation and depletion occurred alternately for the two successive events, indicating that sediment availability along the channel also influences the flow behaviour along the prevailing-transport reach. The comparison between monitoring data, topographical analysis and hydrological simulation allows the estimation of the average solid concentration of the two events and suggests that debris availability has a significant influence on the debris flow volume. © 2020 John Wiley & Sons, Ltd.     

Nonlinear analysis of masonry structures using fiber-section line elements

For seismic analysis of unreinforced masonry (URM) buildings characterized by a box-like behavior, a widely accepted model is based on the equivalent frame idealization of walls. The equivalent frame model uses 1D elements to represent the vertical piers and horizontal spandrels which are connected by rigid nodes. The mechanical characterization of the elements is one of the crucial aspects to predict reasonably the building seismic behavior. Through the comparison with pseudo-static and dynamic experimental tests performed on two-story full-scale buildings, this paper validates the frame modeling in the OpenSees framework, which includes a fiber-section force-based beam element for the axial-flexural behavior, coupled with a cyclic shear-deformation phenomenological law.     

Thermodynamics and phonon dispersion of pyrope and grossular silicate garnets from ab initio simulations

The phonon dispersion and thermodynamic properties of pyrope (\(\hbox {Mg}_3\hbox {Al}_2\hbox {Si}_3\hbox {O}_{12}\)) and grossular (\(\hbox {Ca}_3\hbox {Al}_2\hbox {Si}_3\hbox {O}_{12}\) ) have been computed by using an ab initio quantum mechanical approach, an all-electron variational Gaussian-type basis set and the B3LYP hybrid functional, as implemented in the Crystal program. Dispersion effects in the phonon bands have been simulated by using supercells of increasing size, containing 80, 160, 320, 640, 1280 and 2160 atoms, corresponding to 1, 2, 4, 8, 16 and 27 \(\mathbf {k}\) points in the first Brillouin zone. Phonon band structures, density of states and corresponding inelastic neutron scattering spectra are reported. Full convergence of the various thermodynamic properties, in particular entropy (S) and specific heat at constant volume (\(C_\mathrm{{V}}\)), with the number of \(\mathbf {k}\) points is achieved with 27 \(\mathbf {k}\) points. The very regular behavior of the S(T) and \(C_\mathrm{{V}}(T)\) curves as a function of the number of \(\mathbf {k}\) points, determined by high numerical stability of the code, permits extrapolation to an infinite number of \(\mathbf {k}\) points. The limiting value differs from the 27-\(\mathbf {k}\) case by only 0.40 % at 100 K for S (the difference decreasing to 0.11 % at 1000 K) and by 0.29 % (0.05 % at 1000 K) for \(C_\mathrm{{V}}\). The agreement with the experimental data is rather satisfactory. We also address the problem of the relative entropy of pyrope and grossular, a still debated question. Our lattice dynamical calculations correctly describe the larger entropy of pyrope than grossular by taking into account merely vibrational contributions and without invoking “static disorder” of the Mg ions in dodecahedral sites. However, as the computed entropy difference is found to be larger than the experimental one by a factor of 2–3, present calculations cannot exclude possible thermally induced structural changes, which could lead to further conformational contributions to the entropy.