Second-order Structure Function Research Articles

Numerical simulations of a stochastic dynamics leading to cascades and loss of regularity: Applications to fluid turbulence and generation of fractional Gaussian fields

Motivated by the modeling of the spatial structure of the velocity field of three-dimensional turbulent flows, and the phenomenology of cascade phenomena, a linear dynamics was recently proposed that can generate high velocity gradients from a smooth-in-space forcing term. It is based on a linear partial differential equation stirred by an additive random forcing term that is δ-correlated in time. The underlying proposed deterministic mechanism corresponds to a transport in Fourier space that aims to transfer energy injected at large scales towards small scales. The key role of the random forcing is to realize these transfers in a statistically homogeneous way. Whereas at finite times and positive viscosity the solutions are smooth, a loss of regularity is observed for the statistically stationary state in the inviscid limit. We present here simulations, based on finite volume methods in the Fourier domain and a splitting method in time, which are more accurate than the pseudospectral simulations. We show that our algorithm is able to reproduce accurately the expected local and statistical structure of the predicted solutions. We conduct numerical simulations in one, two, and three spatial dimensions, and we display the solutions both in physical and Fourier space. We additionally display key statistical quantities such as second-order structure functions and power spectral densities at various viscosities. Published by the American Physical Society 2024

Toward mapping turbulence in the intra-cluster medium

Context. Future X-ray observatories with high spectral resolution and imaging capabilities will enable measurements and mappings of emission line shifts in the intracluster medium (ICM). Such direct measurements can serve as unique probes of turbulent motions in the ICM. Determining the level and scales of turbulence will improve our understanding of the galaxy cluster dynamical evolution and assembly, together with a more precise evaluation of the non thermal support pressure budget. This will allow for more accurate constraints to be placed on the masses of galaxy clusters, among other potential benefits. Aims. In this view, we implemented the methods presented in the previous instalments of our work to characterising the turbulence in the intra-cluster medium in a feasibility study with the X-ray Integral Field Unit (X-IFU) on board the future European X-ray observatory, Athena. Methods. From idealized mock observations of a toy model cluster, we reconstructed the second-order structure function built with the observed velocity field to constrain the turbulence. We carefully accounted for the various sources of errors to derive the most realistic and comprehensive error budget within the limits of our approach. With prior assumptions on the dissipation scale and power spectrum slope, we constrained the parameters of the turbulent power spectrum model through the use of Markov chain Monte Carlo (MCMC) sampling. Results. With a very long exposure time, a favourable configuration, and a prior assumption of the dissipation scale, we were able to retrieve the injection scale, velocity dispersion, and power spectrum slope, with 1σ uncertainties for better than ∼15% of the input values. We demonstrated the efficiency of our carefully set framework to constrain the turbulence in the ICM from high-resolution X-ray spectroscopic observations, paving the way for more in-depth investigation of the optimal required observing strategy within a more restrictive observational setup with the future Athena/X-IFU instrument.

Temperature statistics of settling particles in homogeneous isotropic turbulence

Heat transfer of settling particles in homogeneous isotropic turbulence is investigated in one-way coupled Eulerian–Lagrangian direct numerical simulations. In the absence of gravity, our observations highlight that there was a correlation between particles concentration and temperature fronts, which are characterized by high-temperature gradients. The particle clustering is mainly influenced by the preferential concentration and the non-local effect, the latter being the memory of their interaction with the flow along their trajectories. Nevertheless, gravity significantly affects the clustering of particles, which further influence the heat transfer of particles with the fluid. We find that the heat transfer of particles is intricately related to particle dynamic inertia and thermal inertia, which can be quantified by the particle Stokes number St and particle thermal Stokes number Stθ, respectively. In this study, we observe that the variance of particle temperature rate of change is independent of Stθ at small thermal inertia but inversely proportional to Stθ2 for large thermal inertia. In the presence of gravity, the variance converges to Stθ−1 for particles with negligible thermal inertia. Our findings reveal that the second-order structure function of the particle temperature, indicating the enhancement of Lagrangian scalar intermittency, adheres to a power-law behavior with limited thermal inertial particles in the dissipation region, denoted as r2. However, as Stθ increases, the power law of the structure function of the particle temperature is disrupted leading to ‘thermal caustics’. The present findings of the heat transfer behavior of settling particles advance the understanding of gravity effect, and are valuable for the modeling of heat transfer in particle-laden turbulent flows.

Finite Reynolds Number Effect on Small-Scale Statistics in Decaying Grid Turbulence

Since about 1997, the realisation that the finite Reynolds number (FRN) effect needs to be carefully taken into account when assessing the behaviour of small-scale statistics came to the fore. The FRN effect can be analysed either in the real domain or in the spectral domain via the scale-by-scale energy budget equation or the transport equation for the energy spectrum. This analysis indicates that the inertial range (IR) is established only when the Taylor microscale Reynolds number Reλ is infinitely large, thus raising doubts about published power-law exponents at finite values of Reλ, for either the second-order velocity structure function (δu)2¯ or the energy spectrum. Here, we focus on the transport equation of (δu)2¯ in decaying grid turbulence, which represents a close approximation to homogeneous isotropic turbulence. The effect on the small-scales of the large-scale forcing term associated with the streamwise advection decreases as Reλ increases and finally disappears when Reλ is sufficiently large. An approach based on the dual scaling of (δu)2¯, i.e., a scaling based on the Kolmogorov scales (when the separation r is small) and another based on the integral scales (when r is large), yields (δu)2¯∼r2/3 when Reλ is infinitely large. This approach also yields (δu)n¯∼rn/3 when Reλ is infinitely large. These results seem to be supported by the trend, as Reλ increases, of available experimental data. Overall, the results for decaying grid turbulence strongly suggest that a tendency towards the predictions of K41 cannot be dismissed at least at Reynolds numbers which are currently beyond the reach of experiments and direct numerical simulations.

Adopting molecular tagging velocimetry for high-resolution measurements of oscillating grid turbulence

AbstractSingle-component molecular tagging velocimetry (1c-MTV) experiments are conducted in an oscillating grid turbulence facility to systematically clarify the trade-offs between maximizing measurement dynamic range and minimizing measurement uncertainty. The primary aim is to obtain reliable turbulence data with the maximum possible vector resolution ($$\varDelta {x}_{1}$$ Δ x 1 ). Four optical magnifications ($$M{_{0}}$$ M 0 =0.11, 0.22, 0.45, 1.11) with four interframe time delay ($$\varDelta {t}$$ Δ t ) values, for each optical magnification case, ranging from 3 to 10 ms, are investigated. The grid oscillates with a frequency (f) of 3.2 and 4.1 Hz, and a single fixed stroke length (S) of 55 mm. The measurement quality is quantified using three important turbulence descriptors — the second-order transverse velocity correlation function $$(\langle {g(r)}\rangle )$$ ( ⟨ g ( r ) ⟩ ) , the second-order transverse velocity structure function $$([\Delta {_{g}u}]^2)$$ ( [ Δ g u ] 2 ) , and the turbulence energy dissipation rate ($$\langle \epsilon \rangle$$ ⟨ ϵ ⟩ ). Other descriptors such as the longitudinal integral length scale (l) and the longitudinal Taylor microscale ($${\lambda }$$ λ ) are also used to compare with reference values reported in literature. The degree of agreement with the reference values from literature, and insensitivity of the descriptor estimates to different MTV design parameters are used to determine the most robust means of obtaining high-resolution turbulence data from 1c-MTV. The experimental design parameters are contextualized using the Kolmogorov length ($${\eta }$$ η ), time ($${\tau _{\eta }}$$ τ η ), and velocity ($${u_{\eta }}$$ u η ) scales that are based on the most reliable $$\langle {\epsilon }\rangle$$ ⟨ ϵ ⟩ estimates from the current experiments. The pixel-displacement corresponding to $${u_{\eta }}$$ u η $$({\varDelta {x}{_{2}}{_{\eta }}})$$ ( Δ x 2 η ) describes the interdependencies between $${M_{0}}$$ M 0 and $$\varDelta {t}$$ Δ t . The data set with $${M_{0}}=0.45$$ M 0 = 0.45 ($${\eta =}18{\varDelta {x}_{1}}$$ η = 18 Δ x 1 for f=3.2 Hz, and $${\eta =}13{\varDelta {x}_{1}}$$ η = 13 Δ x 1 for f=4.1 Hz) represents the most reliable and the highest resolution case. This conclusion was deduced by taking the performance of all the turbulence descriptors into account. As far as the interframe time delay is concerned, $$\varDelta {t} \ge$$ Δ t ≥ 4 ms ($${\varDelta {x}{_{2}}{_{\eta }}}\ge 0.5 \text {pixels}$$ Δ x 2 η ≥ 0.5 pixels ) works well for all $${{M}_{0}}$$ M 0 . To the authors’ knowledge, the present study documents the first instance of a series of 1c-MTV experiments conducted for a systematic clarification of the nature of balancing the available dynamic range and the measurement uncertainty. Additionally, since the study involves turbulent flow with negligible mean-flow, it serves as an adequate representation of 1c-MTV performance in a three-dimensional flow.

Robust relation of streamwise velocity autocorrelation in atmospheric surface layers based on an autoregressive moving average model

We construct an autoregressive moving average (ARMA) model consisting of the history and random effects for the streamwise velocity fluctuation in boundary-layer turbulence. The distance to the wall and the boundary-layer thickness determine the time step and the order of the ARMA model, respectively. Based on the autocorrelation's analytical expression of the ARMA model, we obtain a global analytical expression for the second-order structure function, which asymptotically captures the inertial, dynamic and large-scale ranges. Specifically, the exponential autocorrelation of the ARMA model arises from the autoregressive coefficients and is modified to logarithmic behaviour by the moving-average coefficients. The asymptotic expressions enable us to determine model coefficients by existing parameters, such as the Kolmogorov and the Townsend–Perry constants. A consequent double-log expression for the characteristic length scale is derived and is justified by direct numerical simulation data with $Re_\tau \approx 5200$ and field-measured neutral atmospheric surface layer data with $Re_\tau \sim O(10^6)$ from the Qingtu Lake Observation Array site. This relation is robust because it applies to $Re_\tau$ from $O(10^4)$ to $O(10^6)$ , and even when the statistics of natural ASL deviate from those of canonical boundary-layer turbulence, e.g. in the case of imbalance in energy production and dissipation, and when the Townsend–Perry constant deviates from traditional values.

Effects of streamwise rotation on helicity and vortex in channel turbulence

Helicity plays a key role in the evolution of vortex structures and turbulent dynamics. The helicity dynamics and vortex structures in streamwise-rotating channel turbulence are discussed in this paper using the helicity budget equation and the differentiated second-order structure function equation of helicity. Generally, rotation and Reynolds numbers exhibit opposing effects on the interscale helicity dynamics and the vortices. Under the buffer layer, the positions of the helicity peaks are proportional to the ratio between the Reynolds and rotation numbers. The mechanism is related to the opposing effects of convection and rotation. Rotation directly affects the helicity balance through the Coriolis term and corresponding pressure term. In the buffer layer, the scale helicity is negative at small scales but positive at large scales, which is mainly induced by the spatial effects (the production and the spatial turbulent convection) but reduced by interscale cascades. Examination of structures reveals the close association between scale helicity and streaks, with streak lift angles exhibiting an increase with rotation and a decrease with Reynolds numbers. In the log-law layer, the Coriolis terms and corresponding pressure terms are proportional to the rotation numbers but remain independent of the Reynolds numbers. The negative scale helicity is forward cascaded towards small scales. Generally, spanwise vortices in the log-law layer are related to sweep events and forward cascades. Our findings indicate that these spanwise vortices are suppressed by rotation but recover with increasing Reynolds numbers, aligning with the effects observed in the scale helicity balance.

Characterization of physical properties of a coastal upwelling filament with evidence of enhanced submesoscale activity and transition from balanced to unbalanced motions in the Benguela upwelling region

Abstract. We combine high-resolution in situ data (acoustic Doppler current profiler (ADCP), Scanfish, and surface drifters) and remote sensing to investigate the physical characteristics of a major filament observed in the Benguela upwelling region. The 30–50 km wide and about 400 km long filament persisted for at least 40 d. Mixed-layer depths were less than 40 m in the filament and over 60 m outside of it. Observations of the Rossby number Ro from the various platforms provide the spatial distribution of Ro for different resolutions. Remote sensing focuses on geostrophic motions of the region related to the mesoscale eddies that drive the filament formation and thereby reveals |Ro|<0.1. Ship-based measurements in the surface mixed layer reveal 0.5<|Ro|<1, indicating the presence of unbalanced, ageostrophic motions. Time series of Ro from triplets of surface drifters trapped within the filament confirm these relatively large Ro values and show a high variability along the filament. A scale-dependent analysis of Ro, which relies on the second-order velocity structure function, was applied to the latter drifter group and to another drifter group released in the upwelling zone. The two releases explored the area nearly distinctly and simultaneously and reveal that at small scales (<15 km) Ro values are twice as large in the filament in comparison to its environment with Ro>1 for scales smaller than ∼500 m. This suggests that filaments are hotspots of ageostrophic dynamics, pointing to the presence of a forward energy cascade. The different dynamics indicated by our Ro analysis are confirmed by horizontal kinetic energy wavenumber spectra, which exhibit a power law k−α with α∼5/3 for wavelengths 2π/k smaller than a transition scale of 15 km, supporting significant submesoscale energy at scales smaller than the first baroclinic Rossby radius (Ro1∼30 km). The detected transition scale is smaller than those found in regions with less mesoscale eddy energy, consistent with previous studies. We found evidence for the processes which drive the energy transfer to turbulent scales. Positive Rossby numbers 𝒪(1) associated with cyclonic motion inhibit the occurrence of positive Ertel potential vorticity (EPV) and stabilize the water column. However, where the baroclinic component of EPV dominates, submesoscale instability analysis suggests that mostly gravitational instabilities occur and that symmetric instabilities may be important at the filament edges.

Scaling law of the second-order structure function over the entire sandstorm process

This study reveals the competitive evolutionary process of the main driving factors in the early, middle and late stages of sandstorms, as shear turbulence becomes dominant and is then suppressed by enhanced thermal stability, based on quadrant analysis of the sand-laden turbulent wind field acquired from field observations over the entire sandstorm process. Moreover, the self-organized state of multiscale structures in the energy-containing region of the sand-laden turbulence is found to change significantly as the sandstorm develops. The logarithmic scaling law that governs the cumulative turbulent kinetic energy for the non-stationary flow in the early and late stages of the sandstorm is different from the existing theoretical formula. The corresponding rate of increase in the cumulative kinetic energy with increasing scale is much higher in these stages than in the middle stage of the sandstorm with steady flow. The change in self-organized state of turbulence is responsible for the flow acceleration and the thermal superimposed effect, rather than the addition of sand particles.

Two-dimensional inverse energy cascade in a laboratory surf zone for varying wave directional spread

Surfzone eddies enhance the dispersion and transport of contaminants, bacteria, and larvae across the nearshore, altering coastal water quality and ecosystem health. During directionally spread wave conditions, vertical vortices (horizontal eddies) are injected near the ends of breaking crests. Energy associated with these eddies may be transferred to larger-scale, low-frequency rotational motions through an inverse energy cascade, consistent with two-dimensional turbulence. However, our understanding of the relationships between the wave conditions and the dynamics and energetics of low-frequency surfzone eddies are largely based on numerical modeling. Here, we test these relationships with remotely sensed and in situ observations from large-scale directional wave basin experiments with varying wave conditions over alongshore-uniform barred bathymetry. Surface velocities derived with particle image velocimetry were employed to assess the spatial scales of low-frequency surfzone eddies and compute structure functions with alongshore velocities. Second-order structure functions for directionally spread waves (σθ≥10°) are consistent with energy flux to larger or smaller length scales, while normally incident, unidirectional waves do not display this behavior. Third-order structure functions suggest that the surfzone flows exhibit a bidirectional energy cascade—a direct cascade to smaller and inverse cascade to larger length scales—during large directional spreads waves (σθ≥18°). However, there is not decisive evidence of an inverse energy cascade for moderate directional spreads (σθ=10°). Energy flux varies by cross-shore location and increases with increasing directional spread and wave height. Eddy decorrelation length scales weakly depend on wave directional spread. These findings advance our understanding of the dynamics linking wave breaking to large-scale rotational motions that enhance mixing and lead to rip currents, important conduits for cross-shore material exchange.

Curvature effects on the structure of near-wall turbulence

The interaction between near-wall turbulence and wall curvature is described for the incompressible flow in a plane channel with a small concave–convex–concave bump on the bottom wall, with height comparable to the wall-normal location of the main turbulent structures. The analysis starts from a database generated by a direct numerical simulation and hinges upon the anisotropic generalised Kolmogorov equations, i.e. the exact budget equations for the second-order structure function tensor. The influence of the bump on the wall cycle and on the energy production, redistribution and transfers is described in the physical and scale spaces. Over the upstream side of the bump, the energy drained from the mean flow to sustain the streamwise fluctuations decreases, and the streaks of high and low streamwise velocity weaken and are stretched spanwise. After the bump tip, instead, the production of streamwise fluctuations grows and the streaks intensify, progressively recovering their characteristic spanwise scale. The wall-normal fluctuations, and thus the quasi-streamwise vortices, are sustained by the mean flow over the upstream side of the bump, while energy flows from the vertical fluctuations to the mean field over the downstream side. On the concave portion of the upstream side, the near-wall fluctuations form structures of spanwise velocity which are consistent with Taylor–Görtler vortices at an early stage of development. Their evolution is described by analysing the scale-space pressure–strain term. A schematic description of the bump flow is presented, in which various regions are identified according to the signs of curvature and streamwise pressure gradient.

Characteristics of clustered particle relative velocity in homogeneous and isotropic turbulence

Particle collisions are mainly governed by the preferential concentration of inertia particles and the formation of fold caustics. By fold caustics, we mean that relative velocity does not go smoothly to zero when the particle separations decrease due to inertia. Despite the importance of the second-order relative velocity structure function, there has been relatively little experimental research on the formation of caustics due to the high accuracy requirements for the particle relative velocity measurements. In the dissipation range, an obvious departure between the second-order structure function of particles normalized by the square of the Kolmogorov velocity and the Kolmogorov turbulent scaling of r2 was observed for all four experimental conditions. In the inertial range, the second-order structure function normalized by the square of the Kolmogorov velocity was consistent with the scale form of r2/3 in all cases. The conditional second-order relative velocity structure function of particles in clusters differs from that of the arithmetically averaged particle statistics in both the dissipation range and the inertial subrange. In the dissipation range, caustics are present for all categories of particles. In the inertial range, the scaling of the second-order velocity structure function is almost identical for all categorized particle types.

Turbulent cascade in fully developed turbulent channel flow

We show that Kolmogorov scale-by-scale equilibrium in the intermediate layer of a fully developed turbulent channel flow is only achieved asymptotically around the Taylor length and, therefore, not in an inertial range. Furthermore, we analyse scale-by-scale turbulence production and interscale turbulence energy transfer in terms of alignments/anti-alignments of fluctuating velocities, straining/compressive relative motions, forward/inverse interscale transfer/cascade and homogeneous/non-homogeneous interscale transfer rate contributions. We also propose leading order scalings for second- and third-order two-point statistics, including the extremum interscale turbulence energy transfer rate and a second-order anisotropic structure function, which acts as a scale-by-scale Reynolds shear stress and determines the scale-by-scale (two-point) turbulence production rate.

Turbulence in compact to giant H ii regions

ABSTRACT Radial velocity fluctuations on the plane of the sky are a powerful tool for studying the turbulent dynamics of emission line regions. We conduct a systematic statistical analysis of the H α velocity field for a diverse sample of nine ${\rm H\, \small{II}}$ regions, spanning two orders of magnitude in size and luminosity, located in the Milky Way and other Local Group galaxies. By fitting a simple model to the second-order spatial structure function of velocity fluctuations, we extract three fundamental parameters: the velocity dispersion, the correlation length, and the power-law slope. We determine credibility limits for these parameters in each region, accounting for observational limitations of noise, atmospheric seeing, and the finite map size. The plane-of-sky velocity dispersion is found to be a better diagnostic of turbulent motions than the line width, especially for lower luminosity regions where the turbulence is subsonic. The correlation length of velocity fluctuations is found to be always roughly 2 per cent of the ${\rm H\, \small{II}}$ region diameter, implying that turbulence is driven on relatively small scales. No evidence is found for any steepening of the structure function in the transition from subsonic to supersonic turbulence, possibly due to the countervailing effect of projection smoothing. Ionized density fluctuations are too large to be explained by the action of the turbulence in any but the highest luminosity sources. A variety of behaviours are seen on scales larger than the correlation length, with only a minority of sources showing evidence for homogeneity on the largest scales.

The influence of coarse particle abundance and spatial distribution on sediment transport and cluster evolution in steep channels under sediment-starved conditions

Large particles strongly influence flow resistance, energy dissipation, bed stability, and channel morphology in mountain streams. We conducted flume experiments to evaluate effects of the density and spatial distribution of coarse particles (i.e., keystones) on bed stability, surface texture and sediment transport in steep channels. Keystones were placed on the channel bed surface at the beginning of each experiment, both in clustered and in uniformly spaced (i.e., anti-clustered) configurations and in different numbers (i.e., for different values of keystone density). The flow rate was increased by 20 % every hour and each experiment continued until the bed was completely scoured. Sediment transport was measured with a sediment trap located at the flume outlet. Topographic data from high-resolution surface imagery and second-order structure functions (SSF) of bed elevations, were used to characterize bed surface structuring while the Ripley’s K function was used to assess the degree of clustering of keystones relative to random spatial distributions. Our results indicate that (1) the density and spatial distribution of keystones exert a minor control on the bed surface grain-size distribution (GSD), (2) for increasing flow rates, the spatial distribution of keystones naturally evolves towards a random distribution, regardless of the initial spatial arrangment, and (3) sediment transport has a higher correlation with the proportion of dislodged keystones than with flow discharge.

An analysis of the turbulence in the central region of the Orion Nebula (M42) – II. Homogeneity and power-spectrum analyses

ABSTRACT In this second paper, we continue our analysis of the turbulence in the Huygens region of the Orion Nebula (M42). We calculate the associated transverse structure functions up to eighth order and find that the higher-order transverse structure functions are almost proportional to the second-order transverse structure function. We find that, after proper normalization, the higher-order transverse structure functions only differ by very small deviations from the second-order transverse structure function in a subinterval of the inertial range. We demonstrate that this implies that the turbulence in the Huygens region is quasi-log–homogeneous or, to a better degree of approximation, binomially weighted log–homogeneous in the statistical sense. This implies that there is some type of invariant statistical structure in the velocity field of the Huygens region. We also obtain and analyse the power spectrum of the turbulent field. We find that it displays a large tail that follows, very approximately, two power laws: one of the form E(k) ∝ k−2.7 for the initial side of the tail, and one of the form E(k) ∝ k−1 for the end of the tail. We find that the power law with exponent β ∼ −2.7 corresponds to spatial scales of 0.0301–0.6450 pc. We find that the exponent of the first power law β ∼ −2.7 is related to the exponent α2 of the second-order structure function in the inertial range. We interpret the second power law with exponent β ∼ −1 as an indicator of viscous-dissipative processes occurring at scales of δr = 1–5 pixels, which correspond to spatial scales of 0.00043–0.00215 pc.

Structure function tensor equations with triple decomposition

Exact budget equations are derived for the coherent and stochastic contributions to the second-order structure function tensor. They extend the anisotropic generalised Kolmogorov equations (AGKE) by considering the coherent and stochastic parts of the Reynolds stress tensor, and are useful for the statistical description of turbulent flows with periodic or quasi-periodic features, like, for example, the alternate shedding after a bluff body. While the original AGKE describe production, transport, inter-component redistribution and dissipation of the Reynolds stresses in the combined space of scales and positions, the new equations, called $\varphi$ AGKE, contain the phase $\varphi$ as an additional independent variable, and describe the interplay among the mean, coherent and stochastic fields at the various phases. The newly derived $\varphi$ AGKE are then applied to a case where an exactly periodic external forcing drives the flow: a turbulent plane channel flow modified by harmonic spanwise oscillations of the wall to reduce drag. The phase-by-phase action of the oscillating transversal Stokes layer generated by the forcing on the near-wall turbulent structures is observed, and a detailed description of the scale-space interaction among mean, coherent and stochastic fields is provided thanks to the $\varphi$ AGKE.

Scaling law of the second-order structure function over the entire sandstorm process

Scaling law of the second-order structure function over the entire sandstorm process

Statistical Characteristics of the Multiscale SST Fractal Structure over the Kuroshio Extension Region Using VIIRS Data

The ocean behaves as a typical multiscale fractal structure, whose dynamic and thermal variabilities extend over a wide range of spatial scales, r, spanning from 10−3 to 107 m. Studying the statistical characteristics of multiscale fractal structures is crucial to understanding the interactions and energy cascade processes between different spatial scales. Remote sensing data are one of the best choices for revealing these statistical characteristics. This work analyzes the multiscale (1–1000 km) fractal structures of sea surface temperature (SST) from the Level-2+ (L2P) satellite orbit Visible Infrared Imaging Radiometer Suite (VIIRS) products over the Kuroshio Extension (KE) region (145°E–160°W, 20°N–50°N), using a conventional method (second-order structure function, D(r)) and a newly developed statistical method (spatial variance, V(r)). The results show that both the power-law distribution slopes of D(r) and V(r) are close to 2/3, which is equivalent to the −5/3 wavenumber spectrum. V(r) is found to be more robust when depicting the fractal structure and variance density, V’(r), compared with D(r). V’(r) is slightly larger at the mesoscale (50–150 km) than at the large scale (higher than 150 km) and is much smaller than that at the submesoscale (smaller than 50 km). Additionally, V’(r) has an indiscernible diurnal variation but remarkable seasonal and latitudinal variations. For the seasonal variability, the maximum V’(r) appears in winter at the large scale and mesoscale, and gradually shifts towards spring at the submesoscale, which implies that a forward energy cascade process may occur during this period. The maximum of the latitude-dependent V’(r) appears around 40°N for all the scales. It is consistent with the latitude of the strongest background SST gradient, indicating that the background SST front is the main source of the strong SST multiscale spatial variabilities over the KE region. This work benefits the application of other high-resolution remote sensing data in this research field, including the forthcoming Surface Water Ocean Topography (SWOT) satellite product.

Evaluation and analysis of Voyager 1 48-s resolution magnetic field data

The aim of this paper is to provide information on the radial dependence of turbulence quantities as input for ab initio numerical cosmic-ray modulation models. We present an evaluation and analysis of 48-s resolution Voyager 1 data for the normal component of the magnetic field. These data were preprocessed to fill in data gaps and filter out spurious peaks. The effectiveness of this process is illustrated for a sample data set, for which we show that the underlying spectrum is indeed Kolmogorov. We use a second-order structure function to find correlation times, and also the spectral index and spectral level in the inertial range from expressions derived by Huang et al. (2010). The usefulness of the second-order structure function for data with significant gaps is illustrated by using synthetic data. Gaps up to 68% cause changes of less than 5% in the quantities that we calculate. We show that beyond about 20 au there are significant changes in the behaviour of the magnetic field magnitude and the inertial range spectral index, neither of which has a reasonable explanation. We find that the ratio δBN/BN remains almost constant out to 20 au, that the correlation times is consistent with results from other studies, and that the spectral index in the inertial range is consistent with the Kolmogorov value. A new result is that we could infer that the break between the energy- and the inertial range has a radial dependence that is consistent with that of the correlation scale.