Layer-Valley Hall Effect under Inversion and Time-Reversal Symmetries

Hall effects have been the central paradigms in modern physics,materials science and practical applications,and have led to many exciting breakthroughs,including the discovery of topological Chern invariants and the revolution of metrological resistance standard.To date,the Hall effects have mainly focused on a single degree of freedom(Do F),and most of them require the breaking of spatial-inversion and/or time-reversal symmetries.Here we demonstrate a new type of Hall effect,i.e.,layer-valley Hall effect,based on a combined layer-valley Do F characterized by the product of layer and valley indices.The layer-valley Hall effect has a quantum origin arising from the layer-valley contrasting Berry curvature,and can occur in nonmagnetic centrosymmetric crystals with both spatial-inversion and time-reversal symmetries,transcending the symmetry constraints of single Do F Hall effect based on the constituent layer or valley index.Moreover,the layer-valley Hall effect is highly tunable and shows a W-shaped pattern in response to the out-of-plane electric fields.Additionally,we discuss the potential detection approaches and material-specific design principles of layer-valley Hall effect.Our results demonstrate novel Hall physics and open up exotic paradigms for new research direction of layer-valleytronics that exploits the quantum nature of the coupled layer-valley DoF.

Linearized waveform inversion for vertical transversely isotropic elastic media:Methodology and multi-parameter crosstalk analysis

Seismic migration and inversion are closely related techniques to portray subsurface images and identify hydrocarbon reservoirs.Seismic migration aims at obtaining structural images of subsurface geologic discontinuities.More specifically,seismic migration estimates the reflectivity function(stacked average reflectivity or pre-stack angle-dependent reflectivity)from seismic reflection data.On the other hand,seismic inversion quantitatively estimates the intrinsic rock properties of subsurface formulations.Such seismic inversion methods are applicable to detect hydrocarbon reservoirs that may exhibit lateral variations in the inverted parameters.Although there exist many differences,pre-stack seismic migration is similar with the first iteration of the general linearized seismic inversion.Usually,seismic migration and inversion techniques assume an acoustic or isotropic elastic medium.Unconventional reservoirs such as shale and tight sand formation have notable anisotropic property.We present a linearized waveform inversion(LWI)scheme for weakly anisotropic elastic media with vertical transversely isotropic(VTI)symmetry.It is based on two-way anisotropic elastic wave equation and simultaneously inverts for the localized perturbations(ΔVp_(0)/Vp_(0)/Vs_(0)/Vs_(0)/,Δ∈,Δδ)from the long-wavelength reference model.Our proposed VTI-elastic LWI is an iterative method that requires a forward and an adjoint operator acting on vectors in each iteration.We derive the forward Born approximation operator by perturbation theory and adjoint operator via adjoint-state method.The inversion has improved the quality of the images and reduces the multi-parameter crosstalk comparing with the adjoint-based images.We have observed that the multi-parameter crosstalk problem is more prominent in the inversion images for Thomsen anisotropy parameters.Especially,the Thomsen parameter is the most difficult to resolve.We also analyze the multi-parameter crosstalk using scattering radiation patterns.The linearized waveform inversion for VTI-elastic media presented in this article provides quantitative information of the rock properties that has the potential to help identify hydrocarbon reservoirs.

True-temperature inversion algorithm for a multi-wavelength pyrometer based on fractional-order particle-swarm optimization

Herein,a method of true-temperature inversion for a multi-wavelength pyrometer based on fractional-order particle-swarm optimization is proposed for difficult inversion problems with unknown emissivity.Fractional-order calculus has the inherent advantage of easily jumping out of local extreme values;here,it is introduced into the particle-swarm algorithm to invert the true temperature.An improved adaptive-adjustment mechanism is applied to automatically adjust the current velocity order of the particles and update their velocity and position values,increasing the accuracy of the true temperature values.The results of simulations using the proposed algorithm were compared with three algorithms using typical emissivity models:the internal penalty function algorithm,the optimization function(fmincon)algorithm,and the conventional particle-swarm optimization algorithm.The results show that the proposed algorithm has good accuracy for true-temperature inversion.Actual experimental results from a rocket-motor plume were used to demonstrate that the true-temperature inversion results of this algorithm are in good agreement with the theoretical true-temperature values.

Surface wave inversion with unknown number of soil layers based on a hybrid learning procedure of deep learning and genetic algorithm

Surface wave inversion is a key step in the application of surface waves to soil velocity profiling.Currently,a common practice for the process of inversion is that the number of soil layers is assumed to be known before using heuristic search algorithms to compute the shear wave velocity profile or the number of soil layers is considered as an optimization variable.However,an improper selection of the number of layers may lead to an incorrect shear wave velocity profile.In this study,a deep learning and genetic algorithm hybrid learning procedure is proposed to perform the surface wave inversion without the need to assume the number of soil layers.First,a deep neural network is adapted to learn from a large number of synthetic dispersion curves for inferring the layer number.Then,the shear-wave velocity profile is determined by a genetic algorithm with the known layer number.By applying this procedure to both simulated and real-world cases,the results indicate that the proposed method is reliable and efficient for surface wave inversion.

A novel efficient energy absorber with free inversion of a metal foam-filled circular tube

In this paper, a novel efficient energy absorber with free inversion of a metal foam-filled circular tube(MFFCT) is designed, and the axial compressive behavior of the MFFCT under free inversion is studied analytically and numerically. The theoretical analysis reveals that the energy is mainly dissipated through the radial bending of the metal circular tube, the circumferential expansion of the metal circular tube, and the metal filled-foam compression. The principle of energy conservation is used to derive the theoretical formula for the minimum compressive force of the MFFCT over free inversion under axial loading. Furthermore, the free inversion deformation characteristics of the MFFCT are analyzed numerically. The theoretical steady values are found to be in good agreement with the results of the finite element(FE) analysis. The effects of the average diameter of the metal tube, the wall thickness of the metal tube, and the filled-foam strength on the free inversion deformation of the MFFCT are considered. It is observed that in the steady deformation stage, the load-carrying and energy-absorbing capacities of the MFFCT increase with the increase in the average diameter of the metal tube, the wall thickness of the metal tube, or the filled-foam strength. The specific energy absorption(SEA) of free inversion of the MFFCT is significantly higher than that of the metal tube alone.

Multi-scale data joint inversion of minerals and porosity in altered igneous reservoirs—A case study in the South China Sea

There are abundant igneous gas reservoirs in the South China Sea with significant value of research,and lithology classification,mineral analysis and porosity inversion are important links in reservoir evaluation.However,affected by the diverse lithology,complicated mineral and widespread alteration,conventional logging lithology classification and mineral inversion become considerably difficult.At the same time,owing to the limitation of the wireline log response equation,the quantity and accuracy of minerals can hardly meet the exploration requirements of igneous formations.To overcome those issues,this study takes the South China Sea as an example,and combines multi-scale data such as micro rock slices,petrophysical experiments,wireline log and element cutting log to establish a set of joint inversion methods for minerals and porosity of altered igneous rocks.Specifically,we define the lithology and mineral characteristics through core slices and mineral data,and establish an igneous multi-mineral volumetric model.Then we determine element cutting log correction method based on core element data,and combine wireline log and corrected element cutting log to perform the lithology classification and joint inversion of minerals and porosity.However,it is always difficult to determine the elemental eigenvalues of different minerals in inversion.This paper uses multiple linear regression methods to solve this problem.Finally,an integrated inversion technique for altered igneous formations was developed.The results show that the corrected element cutting log are in good agreement with the core element data,and the mineral and porosity results obtained from the joint inversion based on the wireline log and corrected element cutting log are also in good agreement with the core data from X-ray diffraction.The results demonstrate that the inversion technique is applicable and this study provides a new direction for the mineral inversion research of altered igneous formations.

Three-dimensional constrained gravity inversion of Moho depth and crustal structural characteristics at Mozambique continental margin

Mozambique's continental margin in East Africa was formed during the break-off stage of the east and west Gondwana lands. Studying the geological structure and division of continent-ocean boundary(COB) in Mozambique's continental margin is considered of great significance to rebuild Gondwana land and understand its movement mode. Along these lines, in this work, the initial Moho was fit using the known Moho depth from reflection seismic profiles, and a 3D multi-point constrained gravity inversion was carried out. Thus, highaccuracy Moho depth and crustal thickness in the study area were acquired. According to the crustal structure distribution based on the inversion results, the continental crust at the narrowest position of the Mozambique Channel was detected. According to the analysis of the crustal thickness, the Mozambique ridge is generally oceanic crust and the COB of the whole Mozambique continental margin is divided.

Crustal and uppermost mantle structure of the northeastern Qinghai-Xizang Plateau from joint inversion of surface wave dispersions and receiver functions with P velocity constraints

Lithospheric structure beneath the northeastern Qinghai-Xizang Plateau is of vital significance for studying the geodynamic processes of crustal thickening and expansion of the Qinghai-Xizang Plateau. We conducted a joint inversion of receiver functions and surface wave dispersions with P-wave velocity constraints using data from the Chin Array Ⅱ temporary stations deployed across the Qinghai-Xizang Plateau. Prior to joint inversion, we applied the H-κ-c method(Li JT et al., 2019) to the receiver function data in order to correct for the back-azimuthal variations in the arrival times of Ps phases and crustal multiples caused by crustal anisotropy and dipping interfaces. High-resolution images of vS, crustal thickness, and vP/vSstructures in the Qinghai-Xizang Plateau were simultaneously derived from the joint inversion. The seismic images reveal that crustal thickness decreases outward from the Qinghai-Xizang Plateau. The stable interiors of the Ordos and Alxa blocks exhibited higher velocities and lower crustal vP/vSratios. While, lower velocities and higher vP/vSratios were observed beneath the Qilian Orogen and Songpan-Ganzi terrane(SPGZ), which are geologically active and mechanically weak, especially in the mid-lower crust.Delamination or thermal erosion of the lithosphere triggered by hot asthenospheric flow contributes to the observed uppermost mantle low-velocity zones(LVZs) in the SPGZ. The crustal thickness, vS, and vP/vSratios suggest that whole lithospheric shortening is a plausible mechanism for crustal thickening in the Qinghai-Xizang Plateau, supporting the idea of coupled lithospheric-scale deformation in this region.

Statistical Inversion Based on Nonlinear Weighted Anisotropic Total Variational Model and Its Application in Electrical Impedance Tomography

Electrical impedance tomography (EIT) aims to reconstruct the conductivity distribution using the boundary measured voltage potential. Traditional regularization based method would suffer from error propagation due to the iteration process. The statistical inverse problem method uses statistical inference to estimate unknown parameters. In this article, we develop a nonlinear weighted anisotropic total variation (NWATV) prior density function based on the recently proposed NWATV regularization method. We calculate the corresponding posterior density function, i.e., the solution of the EIT inverse problem in the statistical sense, via a modified Markov chain Monte Carlo (MCMC) sampling. We do numerical experiment to validate the proposed approach.

Rapid source inversions of the 2023 SE Türkiye earthquakes with teleseismic and strong-motion data

We conducted rapid inversions of rupture process for the 2023 earthquake doublet occurred in SE Türkiye,the first with a magnitude of M_(W)7.8 and the second with a magnitude of M_(W)7.6,using teleseismic and strong-motion data.The teleseismic rupture models of the both events were obtained approximately 88 and 55 minutes after their occurrences,respectively.The rupture models indicated that the first event was an asymmetric bilateral event with ruptures mainly propagating to the northeast,while the second one was a unilateral event with ruptures propagating to the west.This information could be useful in locating the meizoseismal areas.Compared with teleseismic models,the strong-motion models showed relatively higher resolution.A noticeable difference was found for the M_(W)7.6 earthquake,for which the strong-motion models shows a bilateral event,rather than a unilateral event,but the dominant rupture direction is still westward.Nevertheless,all strong-motion models are consistent with the teleseismic models in terms of magnitudes,durations,and dominant rupture directions.This suggests that both teleseismic and strong-motion data can be used for fast determination of major source characteristics.In contrast,the strong-motion data would be preferable in future emergency responses since they are recorded earlier and have a better resolution ability on the source ruptures.

GBO algorithm for seismic source parameters inversion

The use of geodetic observation data for seismic fault parameters inversion is the research hotspot of geodetic inversion, and it is also the focus of studying the mechanism of earthquake occurrence. Seismic fault parameters inversion has nonlinear characteristics, and the gradient-based optimizer(GBO) has the characteristics of fast convergence speed and falling into local optimum hardly. This paper applies GBO algorithm to simulated earthquakes and real LuShan earthquakes in the nonlinear inversion of the Okada model to obtain the source parameters. The simulated earthquake experiment results show that the algorithm is stable, and the seismic source parameters obtained by GBO are slightly closer to the true value than the multi peak particle swarm optimization(MPSO). In the 2013 LuShan earthquake experiment, the root mean square error between the deformation after forwarding of fault parameters obtained by the introduced GBO algorithm and the surface observation deformation was 3.703 mm, slightly better than 3.708 mm calculated by the MPSO. Moreover, the inversion result of GBO algorithm is better than MPSO algorithm in stability. The above results show that the introduced GBO algorithm has a certain practical application value in seismic fault source parameters inversion.

Moment tensor and stress inversion solutions of acoustic emissions during compression and tensile fracturing in crystalline rocks

We investigate the accuracy and robustness of moment tensor(MT)and stress inversion solutions derived from acoustic emissions(AEs)during the laboratory fracturing of prismatic Barre granite specimens.Pre-cut flaws in the specimens introduce a complex stress field,resulting in a spatial and temporal variation of focal mechanisms.Specifically,we consider two experimental setups:(1)where the rock is loaded in compression to generate primarily shear-type fractures and(2)where the material is loaded in indirect tension to generate predominantly tensile-type fractures.In each test,we first decompose AE moment tensors into double-couple(DC)and non-DC terms and then derive unambiguous normal and slip vectors using k-means clustering and an unstructured damped stress inversion algorithm.We explore temporal and spatial distributions of DC and non-DC events at different loading levels.The majority of the DC and the tensile non-DC events cluster around the pre-cut flaws,where macro-cracks later develop.Results of stress inversion are verified against the stress field from finite element(FE)modeling.A good agreement is found between the experimentally derived and numerically simulated stress orientations.To the best of the authors’knowledge,this work presents the first case where stress inversion methodologies are validated by numerical simulations at laboratory scale and under highly heterogeneous stress distributions.

Fast pre-stack multi-channel inversion constrained by seismic reflection features

Classical multi-channel technology can significantly reduce the pre-stack seismic inversion uncertainty, especially for complex geology such as high dipping structures. However, due to the consideration of complex structure or reflection features, the existing multi-channel inversion methods have to adopt the highly time-consuming strategy of arranging seismic data trace-by-trace, limiting its wide application in pre-stack inversion. A fast pre-stack multi-channel inversion constrained by seismic reflection features has been proposed to address this issue. The key to our method is to re-characterize the reflection features to directly constrain the pre-stack inversion through a Hadamard product operator without rearranging the seismic data. The seismic reflection features can reflect the distribution of the stratum reflection interface, and we obtained them from the post-stack profile by searching the shortest local Euclidean distance between adjacent seismic traces. Instead of directly constructing a large-size reflection features constraint operator advocated by the conventional methods, through decomposing the reflection features along the vertical and horizontal direction at a particular sampling point, we have constructed a computationally well-behaved constraint operator represented by the vertical and horizontal partial derivatives. Based on the Alternating Direction Method of Multipliers (ADMM) optimization, we have derived a fast algorithm for solving the objective function, including Hadamard product operators. Compared with the conventional reflection features constrained inversion, the proposed method is more efficient and accurate, proved on the Overthrust model and a field data set.

Rock physics and seismic reflectivity parameterization and amplitude variation with offsets inversion in terms of total organic carbon indicator

Total organic carbon (TOC) prediction with elastic parameter inversions has been widely used in the identification and evaluation of source rocks. However, the elastic parameters used to predict TOC are not only determined by TOC but also depend on the other physical properties of source rocks. Besides, the TOC prediction with the elastic parameters inversion is an indirect method based on the statistical relationship obtained from well logs and experiment data. Therefore, we propose a rock physics model and define a TOC indicator mainly affected by TOC to predict TOC directly. The proposed rock physics model makes the equivalent elastic moduli of source rocks parameterized by the TOC indicator. Combining the equivalent elastic moduli of source rocks and Gray’s approximation leads to a novel linearized approximation of the P-wave reflection coefficient incorporating the TOC indicator. Model examples illustrate that the novel reflectivity approximation well agrees with the exact Zoeppritz equation until incident angles reach 40°. Convoluting the novel P-wave reflection approximation with seismic wavelets as the forward solver, an AVO inversion method based on the Bayesian theory is proposed to invert the TOC indicator with seismic data. The synthetic examples and field tests validate the feasibility and stability of the proposed AVO inversion approach. Using the inversion results of the TOC indicator, TOC is directly and accurately estimated in the target area.

Data inversion of multi-dimensional magnetic resonance in porous media

Since its inception in the 1970s,multi-dimensional magnetic resonance(MR)has emerged as a powerful tool for non-invasive investigations of structures and molecular interactions.MR spectroscopy beyond one dimension allows the study of the correlation,exchange processes,and separation of overlapping spectral information.The multi-dimensional concept has been re-implemented over the last two decades to explore molecular motion and spin dynamics in porous media.Apart from Fourier transform,methods have been developed for processing the multi-dimensional time-domain data,identifying the fluid components,and estimating pore surface permeability via joint relaxation and diffusion spectra.Through the resolution of spectroscopic signals with spatial encoding gradients,multi-dimensional MR imaging has been widely used to investigate the microscopic environment of living tissues and distinguish diseases.Signals in each voxel are usually expressed as multi-exponential decay,representing microstructures or environments along multiple pore scales.The separation of contributions from different environments is a common ill-posed problem,which can be resolved numerically.Moreover,the inversion methods and experimental parameters determine the resolution of multi-dimensional spectra.This paper reviews the algorithms that have been proposed to process multidimensional MR datasets in different scenarios.Detailed information at the microscopic level,such as tissue components,fluid types and food structures in multi-disciplinary sciences,could be revealed through multi-dimensional MR.

Research on inversion method for complex source-term distributions based on deep neural networks

This study proposes a source distribution inversion convolutional neural network (SDICNN), which is deep neural network model for the inversion of complex source distributions, to solve inversion problems involving fixed-source distributions. A function is developed to obtain the distribution information of complex source terms from radiation parameters at individual sampling points in space. The SDICNN comprises two components:a fully connected network and a convolutional neural network. The fully connected network mainly extracts the parameter measurement information from the sampling points,whereas the convolutional neural network mainly completes the fine inversion of the source-term distribution. Finally, the SDICNN obtains a high-resolution source-term distribution image. In this study, the proposed source-term inversion method is evaluated based on typical geometric scenarios. The results show that, unlike the conventional fully connected neural network, the SDICNN model can extract the two-dimensional distribution features of the source terms, and its inversion results are better. In addition, the effects of the shielding mechanism and number of sampling points on the inversion process are examined. In summary, the result of this study can facilitate the accurate assessment of dose distributions in nuclear facilities.

The applicability and underlying factors of frequency-dependent amplitude-versus-offset(AVO)inversion

Recently,the great potential of seismic dispersion attributes in oil and gas exploration has attracted extensive attention.The frequency-dependent amplitude versus offset(FAVO)technology,with dispersion gradient as a hydrocarbon indicator,has developed rapidly.Based on the classical AVO theory,the technology works on the assumption that elastic parameters are frequency-dependent,and implements FAVO inversion using spectral decomposition methods,so that it can take dispersive effects into account and effectively overcome the limitations of the classical AVO.However,the factors that affect FAVO are complicated.To this end,we construct a unified equation for FAVO inversion by combining several Zoeppritz approximations.We study and compare two strategies respectively with(strategy 1)and without(strategy 2)velocity as inversion input data.Using theoretical models,we investigate the influence of various factors,such as the Zoeppritz approximation used,P-and S-wave velocity dispersion,inversion input data,the strong reflection caused by non-reservoir interfaces,and the noise level of the seismic data.Our results show that FAVO inversion based on different Zoeppritz approximations gives similar results.In addition,the inversion results of strategy 2 are generally equivalent to that of strategy 1,which means that strategy 2 can be used to obtain dispersion attributes even if the velocity is not available.We also found that the existence of non-reservoir strong reflection interface may cause significant false dispersion.Therefore,logging,geological,and other relevant data should be fully used to prevent this undesirable consequence.Both the P-and S-wave related dispersion obtained from FAVO can be used as good indicators of a hydrocarbon reservoir,but the P-wave dispersion is more reliable.In fact,due to the mutual coupling of P-and S-wave dispersion terms,the P-wave dispersion gradient inverted from PP reflection seismic data has a stronger hydrocarbon detection ability than the S-wave dispersion gradient.Moreover,there is little difference in using post-stack data or pre-stack angle gathers as inversion input when only the P-wave dispersion is desired.The real application examples further demonstrate that dispersion attributes can not only indicate the location of a hydrocarbon reservoir,but also,to a certain extent,reveal the physical properties of reservoirs.

Petrophysical parameters inversion for heavy oil reservoir based on a laboratory-calibrated frequency-variant rock-physics model

Heavy oil has high density and viscosity, and exhibits viscoelasticity. Gassmann's theory is not suitable for materials saturated with viscoelastic fluids. Directly applying such model leads to unreliable results for seismic inversion of heavy oil reservoir. To describe the viscoelastic behavior of heavy oil, we modeled the elastic properties of heavy oil with varying viscosity and frequency using the Cole-Cole-Maxwell (CCM) model. Then, we used a CCoherent Potential Approximation (CPA) instead of the Gassmann equations to account for the fluid effect, by extending the single-phase fluid condition to two-phase fluid (heavy oil and water) condition, so that partial saturation of heavy oil can be considered. This rock physics model establishes the relationship between the elastic modulus of reservoir rock and viscosity, frequency and saturation. The viscosity of the heavy oil and the elastic moduli and porosity of typical reservoir rock samples were measured in laboratory, which were used for calibration of the rock physics model. The well-calibrated frequency-variant CPA model was applied to the prediction of the P- and S-wave velocities in the seismic frequency range (1–100 Hz) and the inversion of petrophysical parameters for a heavy oil reservoir. The pre-stack inversion results of elastic parameters are improved compared with those results using the CPA model in the sonic logging frequency (∼10 kHz), or conventional rock physics model such as the Xu-Payne model. In addition, the inversion of the porosity of the reservoir was conducted with the simulated annealing method, and the result fits reasonably well with the logging curve and depicts the location of the heavy oil reservoir on the time slice. The application of the laboratory-calibrated CPA model provides better results with the velocity dispersion correction, suggesting the important role of accurate frequency dependent rock physics models in the seismic prediction of heavy oil reservoirs.

Fast 3D joint inversion of gravity and magnetic data based on cross gradient constraint

The gravity and magnetic data can be adopted to interpret the internal structure of the Earth.To improve the calculation efficiency during the inversion process and the accuracy and reliability of the reconstructed physical property models,the triple strategy is adopted in this paper to develop a fast cross-gradient joint inversion for gravity and magnetic data.The cross-gradient constraint contains solving the gradients of the physical property models and performing the cross-product calculation of their gradients.The sparse matrices are first obtained by calculating the gradients of the physical property models derived from the first-order finite difference.Then,the triple method is applied to optimize the storages and the calculations related to the gradients of the physical property models.Therefore,the storage compression amount of the calculations related to the gradients of the physical property models and the cross-gradient constraint are reduced to one-fold of the number of grid cells at least,and the compression ratio increases with the increase of the number of grid cells.The test results from the synthetic data and field data prove that the structural coupling is achieved by using the fast cross-gradient joint inversion method to effectively reduce the multiplicity of solutions and improve the computing efficiency.

A deep-learning-based approach for seismic surface-wave dispersion inversion(SfNet)with application to the Chinese mainland

Surface-wave tomography is an important and widely used method for imaging the crust and upper mantle velocity structure of the Earth.In this study,we proposed a deep learning(DL)method based on convolutional neural network(CNN),named SfNet,to derive the vS model from the Rayleigh wave phase and group velocity dispersion curves.Training a network model usually requires large amount of training datasets,which is labor-intensive and expensive to acquire.Here we relied on synthetics generated automatically from various spline-based vS models instead of directly using the existing vS models of an area to build the training dataset,which enhances the generalization of the DL method.In addition,we used a random sampling strategy of the dispersion periods in the training dataset,which alleviates the problem that the real data used must be sampled strictly according to the periods of training dataset.Tests using synthetic data demonstrate that the proposed method is much faster,and the results for the vS model are more accurate and robust than those of conventional methods.We applied our method to a dataset for the Chinese mainland and obtained a new reference velocity model of the Chinese continent(ChinaVs-DL1.0),which has smaller dispersion misfits than those from the traditional method.The high accuracy and efficiency of our DL approach makes it an important method for vS model inversions from large amounts of surface-wave dispersion data.