On spectral scaling laws for incompressible anisotropic MHD turbulence

A heuristic model is given for anisotropic magnetohydrodynamics (MHD) turbulence in the presence of a uniform external magnetic field B_0 {\bf {\hat e}_{\pa}}. The model is valid for both moderate and strong $B_0$ and is able to describe both the strong and weak wave turbulence regimes as well as the transition between them. The main ingredient of the model is the assumption of constant ratio at all scales between \add{the} linear wave period and \add{the} nonlinear turnover timescale. Contrary to the model of critical balance introduced by Goldreich and Sridhar [P. Goldreich and S. Sridhar, ApJ {\bf 438}, 763 (1995)], it is not assumed in addition that this ratio be equal to unity at all scales which allows us to use the Iroshnikov-Kraichnan phenomenology. It is then possible to recover the widely observed anisotropic scaling law \kpa \propto \kpe^{2/3} between parallel and perpendicular wavenumbers (with reference to B_0 {\bf {\hat e}_{\pa}}) and to obtain the universal prediction, $3\alpha + 2\beta = 7$, for the total energy spectrum E(\kpe,\kpa) \sim \kpe^{-\alpha} \kpa^{-\beta}. In particular, with such a prediction the weak Alfv\'en wave turbulence constant-flux solution is recovered and, for the first time, a possible explanation to its precursor found numerically by Galtier et al [S. Galtier et al., J. Plasma Phys. {\bf 63}, 447 (2000)] is given

Intermittency in MHD turbulence and coronal nanoflares modelling

High resolution numerical simulations, solar wind data analysis, and measurements at the edges of laboratory plasma devices have allowed for a huge progress in our understanding of MHD turbulence. The high resolution of solar wind measurements has allowed to characterize the intermittency observed at small scales. We are now able to set up a consistent and convincing view of the main properties of MHD turbulence, which in turn constitutes an extremely efficient tool in understanding the behaviour of turbulent plasmas, like those in solar corona, where in situ observations are not available. Using this knowledge a model to describe injection, due to foot-point motions, storage and dissipation of MHD turbulence in coronal loops, is built where we assume strong longitudinal magnetic field, low beta and high aspect ratio, which allows us to use the set of reduced MHD equations (RMHD). The model is based on a shell technique in the wave vector space orthogonal to the strong magnetic field, while the dependence on the longitudinal coordinate is preserved. Numerical simulations show that injected energy is efficiently stored in the loop where a significant level of magnetic and velocity fluctuations is obtained. Nonlinear interactions give rise to an energy cascade towards smaller scales where energy is dissipated in an intermittent fashion. Due to the strong longitudinal magnetic field, dissipative structures propagate along the loop, with the typical speed of the Alfvén waves. The statistical analysis on the intermittent dissipative events compares well with all observed properties of nanoflare emission statistics. Moreover the recent observations of non thermal velocity measurements during flare occurrence are well described by the numerical results of the simulation model. All these results naturally emerge from the model dynamical evolution without any need of an ad-hoc hypothesis

Numerical modeling of landquakes

International audienceThe Thurwieser landslide that occurred in Italy in 2004 and the seismic waves it generated have been simulated and compared to the seismic signal recorded a few tens of kilometers from the landslide source (i.e., landquake). The main features of the low frequency seismic signal are reproduced by the simulation. Topography effects on the flowing mass have a major impact on the generated seismic signal whereas they weakly affect low-frequency wave propagation. Simulation of the seismic signal makes it possible to discriminate between possible alternative scenarios for flow dynamics and to provide first estimates of the rheological parameters during the flow. As landquakes are continuously recorded by seismic networks, our results provide a new way to collect data on the dynamics and rheology of natural flows

A rockslide-generated tsunami in a Greenland fjord rang the Earth for 9 days

Climate change is increasingly predisposing polar regions to large landslides. Tsunamigenic landslides have occurred recently in Greenland, but none have been reported from the eastern fjords. In September 2023, we detected the start of a 9-day-long global 10.88 mHz (92 s) monochromatic very long-period (VLP) seismic signal, originating from East Greenland. We demonstrate how this event started with a 25×106 m3 glacial thinning-induced rock-ice avalanche plunging into Dickson Fjord, triggering a 200 m high tsunami. Simulations show the tsunami stabilized into a 7 m-high long-duration seiche with a near-identical frequency (11.45 mHz) and slow amplitude decay as the seismic signal. An oscillating, fjord-transverse single-force with a maximum amplitude of 5×1011 N reproduces the seismic amplitudes and their radiation pattern relative to the fjord, demonstrating how a seiche directly caused the 9-day-long seismic signal. Our findings highlight how climate change is causing cascading, hazardous feedbacks between the cryosphere, hydrosphere, and lithosphere

Memory of the Unjamming Transition during Cyclic Tiltings of a Granular Pile

Discrete numerical simulations are performed to study the evolution of the micro-structure and the response of a granular packing during successive loading-unloading cycles, consisting of quasi-static rotations in the gravity field between opposite inclination angles. We show that internal variables, e.g., stress and fabric of the pile, exhibit hysteresis during these cycles due to the exploration of different metastable configurations. Interestingly, the hysteretic behaviour of the pile strongly depends on the maximal inclination of the cycles, giving evidence of the irreversible modifications of the pile state occurring close to the unjamming transition. More specifically, we show that for cycles with maximal inclination larger than the repose angle, the weak contact network carries the memory of the unjamming transition. These results demonstrate the relevance of a two-phases description -strong and weak contact networks- for a granular system, as soon as it has approached the unjamming transition.Comment: 13 pages, 15 figures, soumis \`{a} Phys. Rev.

Relating hygroscopicity and optical properties to chemical composition and structure of secondary organic aerosol particles generated from the ozonolysis of α-pinene

Secondary organic aerosol (SOA) were generated from the ozonolysis of α-pinene in the CESAM (French acronym for Experimental Multiphasic Atmospheric Simulation Chamber) simulation chamber. The SOA formation and aging were studied by following their optical, hygroscopic and chemical properties. The optical properties were investigated by determining the particle complex refractive index (CRI). The hygroscopicity was quantified by measuring the effect of relative humidity (RH) on the particle size (size growth factor, GF) and on the scattering coefficient (scattering growth factor, f(RH)). The oxygen to carbon atomic ratios (O : C) of the particle surface and bulk were used as a sensitive parameter to correlate the changes in hygroscopic and optical properties of the SOA composition during their formation and aging in CESAM. The real CRI at 525 nm wavelength decreased from 1.43–1.60 (±0.02) to 1.32–1.38 (±0.02) during the SOA formation. The decrease in the real CRI correlated to the O : C decrease from 0.68 (±0.20) to 0.55 (±0.16). In contrast, the GF remained roughly constant over the reaction time, with values of 1.02–1.07 (±0.02) at 90% (±4.2%) RH. Simultaneous measurements of O : C of the particle surface revealed that the SOA was not composed of a homogeneous mixture, but contained less oxidised species at the surface which may limit water absorption. In addition, an apparent change in both mobility diameter and scattering coefficient with increasing RH from 0 to 30% was observed for SOA after 14 h of reaction. We postulate that this change could be due to a change in the viscosity of the SOA from a predominantly glassy state to a predominantly liquid state