3D Velocity Research Articles

XookSuut: A BAYESIAN TOOL FOR MODELING CIRCULAR AND NON–CIRCULAR FLOWS ON 2D VELOCITY MAPS

We present XookSuut, a Python implementation of the DiskFit algorithm, optimized to perform robust Bayesian inference on parameters describing models of circular and noncircular rotation in galaxies. XookSuut surges as a Bayesian alternative for kinematic modeling of 2D velocity maps; it implements effcient sampling methods, specifically Markov Chain Monte Carlo (MCMC) and Nested Sampling (NS), to obtain the posteriors and marginalized distributions of kinematic models including circular motions, axisymmetric radial flows, bisymmetric flows, and harmonic decomposition of the LoS velocity. In this way, kinematic models are obtained by pure sampling methods, rather than standard minimization techniques based on the Χ2. All together, XookSuut represents a sophisticated tool for deriving rotational curves and to explore the error distribution and covariance between parameters.

Multichannel Doppler fiber-imager spectrometer for spatiotemporal characterization of high-intensity laser-driven plasmas.

We demonstrate a flexible multichannel fiber-based imaging Doppler spectrometer to characterize plasmas in high intensity (≥1 × 1018 W/cm2) laser-plasma experiments at high repetition rates. This instrument collects data from ×21 different plasma locations combining optical fibers and a single imaging spectrometer. This diagnostic maps the plasma velocity evolution as a function of time with sub-pico-second resolution. Experimental results showing 2D velocity measurements of plasma with 20 μm spatial resolution are presented. Intensities of the order of 1018 W/cm2 were used to generate a plasma, while a much less intense, frequency doubled (400nm), probe beam (1011 W/cm2) was used to measure the Doppler shift from the plasma critical surface. The instrument can be scaled to a larger number of channels (e.g., 100) still using a single spectrometer.

Using time-lapse seismic velocity changes to monitor the Domo de San Pedro Geothermal field, Mexico

The Domo de San Pedro Geothermal field, located in the westernmost part of the trans-Mexican volcanic belt, is part of a silicic volcanic complex. Geothermal energy is produced since 2015 from 9 wells as deep as 3 km. We monitor geothermal operations using a continuous seismic dataset from a temporary network of 20 broadband stations from March 2021 to January 2022. The well logs describing the geothermal operations allow us to investigate the seismic response of the geothermal field to fluid extraction and re-injection across the geothermal wells. We produced a catalog of 217 local earthquakes and derive an optimized 1D velocity model for the geothermal field between 1.2 m above the sea level (masl) and 12 km depth. Apparent velocity changes (dv/v) are measured using seismic ambient noise interferometry from the coda of noise cross-correlations. Next, we identify and locate subsurface velocity anomalies, combining the methods of station collapse sensitivity and ambient noise imaging applied to the continuous seismic noise records. We distinguish three periods of seismic activity where prominent changes in velocity coincide with the re-injection of fluids and significant variation in the production rate of the geothermal wells. We determine the distribution of these velocity transients inside the geothermal field and propose that they are associated with variations in reservoir fluid content and effective stress changes beneath the volcanic system and are responsible for induced microseismicity. Our study shows that implementing dv/v quantification as an additional tool for geothermal seismic monitoring can lead to an improved understanding of the behavior of the geothermal system over time and ultimately de-risk geothermal operations.

Analyzing 50 years of the Lacq induced seismicity (SouthWestern, France) highlights the role of fluid injection

Summary The Lacq area in southwest France has been associated with continuous moderate induced seismic activity since 1969. However, the mechanisms driving this induced seismicity are not fully understood: reservoir depletion has been proposed as the main factor, and more recently wastewater injection has been suggested to play a more important role (Grasso et al., 2021). The interpretation of these mechanisms relies heavily on the quality of earthquake locations, which we prove to be weak due to a lack of local instrumentation for several years. In order to provide the most complete and reliable induced event catalog for the studies of the Lacq induced seismicity mechanisms & seismic hazard, we made an exhaustive compilation, analysis and improvement of all available catalogs. We also provided new earthquake detections & relocations in a 3D velocity model from past and present temporary deployments never used for studying the Lacq area. Important remaining location uncertainties lead us to also carefully sort the events according to their location confidence, defining 3 classes of events (unconstrained location, location constrained within 2-3 km and 1-2 km respectively). This new harmonized catalog and the identification of well-constrained events, covering 50 years of induced seismicity, allow us to propose that wastewater injection is almost certainly the main mechanism driving the seismicity, with (i) most of the constrained events located within the reservoir boundaries and (ii) the released seismic energy variations following variations in injection operations at different scales. In particular, we have also highlighted a change in the injection-seismicity relationship around 2010–2013. From 2013, despite lower injection volumes, seismicity remained persistent and some clusters of earthquakes were detected predominantly in spring, summer, and early autumn, except in winter periods. From 2016, we observed a strong temporal relationship between days with higher rate/volume injections (approximately above 400m3/day) and both clustered events and higher magnitude earthquakes (greater than 2.4).

Mechanism for High-Precision Control of Movement at Maximum Output in the Vertical Jump Task.

Human movements are governed by a tradeoff between speed and accuracy. Previous studies that have investigated the tradeoff relationship in sports movements involving whole-body movements have been limited to examining the relationship from the perspective of competition-specific movements, and the findings on whether the relationship is valid have not been unified. Therefore, this study incorporated a vertical jump task with the introduction of a condition in which landing position control was added to evaluate the essence of a sports movement that requires both speed and accuracy. Accuracy was examined using a method that quantifies the coordinates of the landing and takeoff positions using entropy. The mechanism of that tradeoff was then examined by confirming the phenomenon and analyzing the 3D vector trajectories. An increase in accuracy and a decrease in speed were observed when the landing position was the control target, even in the vertical jumping task normally performed at maximum effort, and the 3D velocity vector was characterized by the following: a reduced scalar and a more vertical direction. While the entropy from the takeoff to the landing position seemed to decrease when the accuracy of the landing position improved, the following noteworthy results were obtained given the characteristics of the vertical jump. Unlike traditional feedback control in the entropy reduction in hand movements, the trajectory is predetermined in a feedforward-like manner by controlling the initial velocity vector at takeoff, which allows the landing point to be adjusted.

Seismic evidence of a plume conduit in a metasomatized lithosphere beneath the Barmer rift in northwestern India

Teleseismic tomography was done to investigate the lithosphere beneath northwestern India which was affected by several magmatic episodes since Neoproterozoic times. Broadband stations were deployed in phases at 29 locations in an area of 300 km × 300 km in northwestern India. About 8000 high-quality P phases corresponding to 934 events were used. Relative residuals estimated through adaptive stacking were inverted using the Fast Marching teleseismic tomography to obtain the 3D velocity perturbations. The tomograms revealed an anomalous ≈ 100 km wide, cylindrical zone of low velocities (−2% to −2.3%) extending from a depth of ≈ 60 km to ≈140 km beneath the Barmer rift. A broader zone of low velocities (−1%) encompasses the upper mantle down to a depth of at least 200 km. The anomaly lies beneath the localized, small-volume, eruptive centers of the earliest phase of Deccan volcanism within and on the flanks of the rift which falls on the projected extension of the Reunion hotspot track along the rift system of northwestern India. Explaining such large anomalies requires a significant increase in temperature, presence of melt and/or compositional changes. Geochemical studies indicate that the polychronous (early Cretaceous to Paleogene) alkaline complexes above the anomalous zone are derived from hydrated mantle peridotite and contain an abundance of phlogopites, suggesting a hydrated and metasomatized lithosphere beneath northwestern India. A thermal source impinging from below upon a metasomatized, weak, and rifted lithosphere may result in the observed signature akin to that of a plume conduit that overprints the older signatures of the Neoproterozoic Malani magmatic event.

The Effects of Layer-Induced Elastic Anisotropy on Microseismic Monitoring Records in Underground Coal Mines

Real-time microseismic observations are extensively used in underground hard rock mines. However, microseismic monitoring in coal mines is not common. In this study, the features of microseismic (MS) event waveforms were discussed under conditions of thin-layered sedimentary media. The possibility of source mechanism evaluation was considered for the MS event recorded in coal mine. Moment tensor (MT) inversion of the waveforms can be challenging due to influence of thin-layered medium structure. Splitting of shear SH- and SV-waves is one of the characteristic features that affect MS records in coal mines. The effect of splitting was demonstrated using real microseismic data collected by an underground seismic array installed in the Polosukhinskaya coal mine, Russia. Velocities of shear SH- and SV-waves were estimated from travel times. The results reveal significant variations in shear wave velocities composed of a radial anisotropy of 75%, thus confirming the importance of a thin-layered model for moment tensor inversion. The thin-layered model was built for geological conditions at the Polosukhinskaya coal mine and synthetic seismograms were generated using QSEIS (Pyrocko OS). A full moment tensor of a single typical event was found employing the layered 1D velocity model.

Understanding earthquake location uncertainties using global sensitivity analysis framework

SUMMARY Earthquake hypocentres are routinely obtained by a common inversion problem of P- and S-phase arrivals observed on a seismological network. Improving our understanding of the uncertainties associated with the hypocentral parameters is crucial for reliable seismological analysis, understanding of tectonic processes and seismic hazard assessment. However, current methods often overlook uncertainties in velocity models and variable trade-offs during inversion. Here, we propose to unravel the effects of the main sources of uncertainty in the location process using techniques derived from the Global Sensitivity Analysis (GSA) framework. These techniques provide a quantification of the effects of selected variables on the variance of the earthquake location using an iterative model that challenges the inversion scheme. Specifically, we consider the main and combined effects of (1) variable network geometry, (2) the presence of errors in the analyst’s observations and (3) errors in velocity parameters from a 1-D velocity model. These multiple sources of uncertainty are described by a dozen of random variables in our model. Using a Monte Carlo sampling approach, we explore the model configurations and analyse the differences between the initial reference location and 100 000 resulting hypocentral locations. The GSA approach using Sobol's variance decomposition allows us to quantify the relative importance of our choice of variables. It highlights the critical importance of the velocity model approximation and provides a new objective and quantitative insight into understanding the sources of uncertainty in the inversion process.

The influence of an anticline structure on ambient noise spectral anomalies at an underground gas storage

SUMMARY The purpose of this study is to investigate the seismic ambient noise spectral anomalies that occur near gas reservoirs. These anomalies involve a significant spectral amplification of the vertical component for frequencies generally between 1.5 and 4 Hz and have been reported at various hydrocarbon sites worldwide. There are differing views on the mechanisms responsible for these anomalies. The guideline for this study is that many hydrocarbon reservoirs share a common geological feature: an anticline structure. It appears to cause site effects that influence the amplitude of the ambient noise wavefield. This research examines a dense real data set of ambient noise recorded at the Chémery underground gas storage site in France. The analysis identifies stable spectral anomalies between 1.2 and 2.4 Hz that are correlated to the position of the anticline structure, which also corresponds to the position of the gas bubble. We use a beamforming technique to study the composition and the origin of the ambient noise, and show that the variations of the spectral anomalies over time are correlated to changes in the source wavefield. Finally, we perform numerical simulations of Rayleigh wave propagation within a realistic 3-D velocity model of the Chémery site, while using source distributions directly extracted from real data analysis. The comparison of the simulated anomalies with real data yields a satisfactory qualitative fit. We conclude that the fundamental-mode Rayleigh wave site effect on the anticline is the main mechanism of the spectral anomaly.

3D Stresses and Velocities Caused by Continental Plateaus: Scaling Analysis and Numerical Calculations With Application to the Tibetan Plateau

AbstractUnderstanding stresses is crucial for geodynamics since they govern rock deformation and metamorphic reactions. However, the magnitudes and distribution of crustal stresses are still uncertain. Here, we use a 3D numerical model in spherical coordinates to investigate stresses and velocities that result from lateral crustal thickness variations around continental plateaus like those observed for the Tibetan plateau. We do not consider any far‐field deformation so that the plateau deforms by horizontal dilatation and vertical thinning. We assume viscous creep, a simplified plateau geometry, and simplified viscosity and density distributions to couple the numerical results with a scaling analysis. Specifically, we study the impact of the viscosity ratio between crust and lithospheric mantle, a rectangular plateau corner, a stress‐dependent power‐law flow law and Earth's double curvature on the crustal stress field and horizontal velocities. Two orders of magnitude variation in crustal and lithospheric mantle viscosities change the maximum crustal differential stress only by a factor of ≈2. We derive simple analytical estimates for the crustal deviatoric stress and horizontal velocity which agree to first order with 3D numerical calculations. We apply these estimates to calculate the average crustal viscosity in the eastern Tibetan plateau as ≈1022 Pa · s. Furthermore, our results show that a corner strongly affects the stress distribution, particularly the shear stresses, while Earth's curvature has a minor impact on the stresses. We further discuss potential implications of our results to strike‐slip faulting and fast exhumation around the Tibetan plateau's syntaxes.

The Mass of the Large Magellanic Cloud from the Three-dimensional Kinematics of Its Globular Clusters

We estimate the mass of the Large Magellanic Cloud (LMC) using the kinematics of 30 LMC globular clusters (GCs). We combine proper motions measured with the Hubble Space Telescope, Gaia, or a combination of the two, from a recent study by Bennet et al. with literature line-of-sight velocities (LOSVs) to give three components of motion. With these, we derive a 3D velocity dispersion anisotropy β=−0.72−1.07+0.62 , consistent with the GCs forming a flattened system with significant azimuthal motion. We then apply a tracer mass estimator and measure an enclosed mass M(<13.2kpc)=2.66−0.36+0.42×1010 M ⊙. This is broadly consistent with results from previous studies of the LOSVs of GCs and other luminous tracers. Assuming a cosmologically constrained Navarro–Frenk–White distribution for the dark matter, this implies a virial mass Mvirial=1.80−0.54+1.05×1011 M ⊙. Despite being an extrapolation by almost an order of magnitude in radius, this result is consistent with published estimates from other methods that are directly sensitive to the LMC’s total mass. Our results support the conclusion that the LMC is approximately 17−6+10 % of the Milky Way’s mass, making it a significant contributor to the Local Group potential.

An investigation of the effect of stratification stability and saltation sand flux on the anisotropy of atmospheric surface layer wall turbulence

In the atmospheric surface layer (ASL), the anisotropic characteristics of turbulence become more complex due to the strong influence of heat flux and particle motion. In this paper, the effects of stratification stability and total saltation sand flux on turbulent anisotropy of ASL wall turbulence are systematically analyzed by using high-frequency 3D velocity, temperature, and saltation sand flux measured in the field. Our results show that the temperature gradient intensifies the spanwise turbulent motion and enhances the vertical motion asymmetry. Although the buoyancy-driven turbulence enhances the spanwise variation, the momentum transfer mode between the streamwise–vertical and spanwise–vertical does not change. With the increase in instability, the large-scale motions (LSM) of the spanwise velocity may be destroyed and more spanwise small-scale motions (SSM) formed. In addition, saltation increases the sharpness, randomness, and extreme frequency of spanwise velocity distribution, but does not change the momentum transfer mode. Saltation may enlarge the fluctuation range and scale of the spanwise velocity, but may not be the main factor affecting the anisotropy of SSM. These results can deepen people's understanding of wind-blown sand movement and turbulence properties in ASL, and provides empirical data and insights that can significantly contribute to the development and refinement of turbulence models, particularly those that need to account for the complex interplay between stratification stability, saltation sand flux, and turbulence anisotropy.

On the Effect of 3D Wave Propagation on 2D Regional‐Scale Velocity Model Building

AbstractActive seismic surveys are routinely employed by academia to study geological structure of the crust and upper mantle. Wavefields generated during these surveys are sampled at the receiver locations, but the wave‐paths traveled from a source to a sensor remains unknown. Although seismic acquisition layouts designed to investigate complex crustal‐scale environments are often two‐dimensional, the seismogram recorded at the receiver location represents information gathered along the three‐dimensional wavepaths that might offset from the 2D source/receiver profile along its transverse direction. This so‐called 3D‐effect distorts the results of 2D seismic imaging, which is unable to handle the out‐of‐plane propagation. Despite the numerous 2D seismic imaging case studies, the assessment of this issue is often overlooked. However, the problem exists ‐ especially for crustal‐scale profiles, where seismic energy propagates over distances of hundreds of kilometers and probes different crustal units. In this work we investigate the impact of 3D‐effect on the results of 2D velocity model building from the academic ocean‐bottom seismometer data. We show with polarization analysis how the 3D‐effect can manifest itself in the data domain. Using various scenarios of acquisition we evaluate the imprint of the out‐of‐plane propagation on the data and the results of full‐waveform inversion. We show that 2D velocity model building from the seismic profiles acquired in the complex geological setting can lead to wrong solution. Looking for the remedy to this issue we couple different configurations of acquisition geometries with 3D full‐waveform inversion that allow to handle the 3D effect and provide correct model reconstruction.

Angle-dependent 4D seismic time-strain inversion for estimating subsurface thickness and velocity changes

In time-lapse (also called 4D) seismic analysis, the time shift of a certain seismic reflection event is caused by changes in the seismic velocity and the depth of the event. An interpretation of 4D time shifts is normally simplified by neglecting displacement changes (strains) or assuming linear relations between thickness and velocity strains. Here, we go beyond these assumptions and propose a least-squares optimization method to simultaneously estimate the thickness and velocity strains in vertical transverse isotropic media from angle-dependent 4D seismic time strains. Through examples from synthetic and field data, we show that 4D thickness strains, velocity strains, and anisotropic parameter changes can be estimated simultaneously without prior knowledge about the geomechanics of the survey area. Our time-strain inversion method can be applied to any other 4D seismic data set in which angle-stack images are available. We see that our method has high potential in many other applications, because thickness and velocity strains are the fundamental components of most physical properties used in 4D seismic and geomechanics applications.

A convolutional neural network for creating near-surface 2D velocity images from GPR antenna measurements

Abstract In this research, our focus lies in exploring the effectiveness of a frequency-velocity convolutional neural network (CNN) in the efficient and non-intrusive acquisition of 2D wave velocity visuals of near-surface geological substances, accomplished through the analysis of data from ground-penetrating radar (GPR). To learn complex correlations between antenna readings and subsurface velocities, the proposed CNN model makes use of the spatial features present in the GPR data. By employing a network architecture capable of accurately detecting both local and global patterns within the data, it becomes feasible to efficiently extract valuable velocity information from GPR readings. The CNN model is trained and validated using a substantial dataset consisting of GPR readings along with corresponding ground truth velocity images. Diverse subsurface settings, encompassing different soil types and geological characteristics, are employed to gather the GPR measurements. In the supervised learning approach employed to train the CNN model, the GPR measurements serve as input, while the associated ground truth velocity images are utilized as target outputs. The model is trained using backpropagation and optimized using a suitable loss function to reduce the difference between the predicted velocity images and the actual images. The experimental results demonstrate the effectiveness of the proposed CNN method in accurately deriving 2D velocity images of near-surface materials from GPR antenna observations. Compared to traditional techniques, the CNN model exhibits superior velocity calculation precision and achieves high levels of accuracy. Moreover, when applied to unseen GPR data, the trained model exhibits promising generalization abilities, highlighting its potential for practical subsurface imaging applications.

Nonergodic Ground-Motion Models for Subduction Zone and Crustal Earthquakes in Japan

ABSTRACT We investigate the nonergodic behavior of ground motions from subduction zone earthquakes and crustal earthquakes in Japan using the National Research Institute for Earth Science and Disaster Resilience strong-motion flatfile together with various reference ergodic ground-motion models (GMMs). For the nonergodic path effects, our nonergodic GMM has two path effects terms: a cell-specific linear-distance scaling, denoted as δP2PQ, that mimics the effects of a 3D Q structure, and a source- and site-specific term related to the effects of the 3D velocity structure, denoted δP2PV. The resulting model provides spatially varying nonergodic source, site, and path terms. The GMM smoothly interpolates and extrapolates the nonergodic terms in space so that the GMM can be applied to any combination of source and site locations in the region. In regions where data are sparse, the median nonergodic terms will approach zero but with large epistemic uncertainty. Over the period range of 0–10 s, the ranges of the standard deviations of the nonergodic source, path, and site terms are 0.2–0.65, 0.25–0.40, and 0.15–0.40 natural log units, respectively. The fully nonergodic models lead to a 40%–55% decrease in the aleatory standard deviation compared to the reference ergodic GMMs. This large reduction in the aleatory standard deviation combined with the change in the median given by the nonergodic terms can exert a significant impact on the computed seismic hazard for the Japan region.

Quality Control of DAS VSP Data in Desert Environment Using Simulations and Matching Filters.

The unconsolidated near surface and large, daily temperature variations in the desert environment degrade the vertical seismic profiling (VSP) data, posing the need for rigorous quality control. Distributed acoustic sensing (DAS) VSP data are often benchmarked using geophone surveys as a gold standard. This study showcases a new simulation-based way to assess the quality of DAS VSP acquired in the desert without geophone data. The depth uncertainty of the DAS channels in the wellbore is assessed by calibrating against formation depth based on the concept of conservation of the energy flux. Using the 1D velocity model derived from checkshot data, we simulate both DAS and geophone VSP data via an elastic pseudo-spectral finite difference method, and estimate the source and receiver signatures using matching filters. These field geophone data show high amplitude variations between channels that cannot be replicated in the simulation. In contrast, the DAS simulation shows a high visual similarity with the field DAS first arrival waveforms. The simulated source and receiver signatures are visually indistinguishable from the field DAS data in this study. Since under perfect conditions, the receiver signatures should be invariant with depth, we propose a new DAS data quality control metric based on local variations of the receiver signatures which does not require geophone measurements.

Spatio-spectral 4D coherent ranging using a flutter-wavelength-swept laser

Coherent light detection and ranging (LiDAR), particularly the frequency-modulated continuous-wave LiDAR, is a robust optical imaging technology for measuring long-range distance and velocity in three dimensions (3D). We propose a spatio-spectral coherent LiDAR based on a unique wavelength-swept laser to enable both axial coherent ranging and lateral spatio-spectral beam scanning simultaneously. Instead of the conventional unidirectional wavelength-swept laser, a flutter-wavelength-swept laser (FWSL) successfully decoupled bidirectional wavelength modulation and continuous wavelength sweep, which overcame the measurable distance limited by the sampling process. The decoupled operation in FWSL enabled sequential sampling of flutter-wavelength modulation across its wide spectral bandwidth of 160 nm and, thus, allowed simultaneous distance and velocity measurement over an extended measurable distance. Herein, complete four-dimensional (4D) imaging, combining real-time 3D distance and velocity measurements, was implemented by solid-state beam scanning. An acousto-optic scanner was synchronized to facilitate the other lateral beam scanning, resulting in an optimized solid-state coherent LiDAR system. The proposed spatio-spectral coherent LiDAR system achieved high-resolution coherent ranging over long distances and real-time 4D imaging with a frame rate of 10 Hz, even in challenging environments.

4D time lapse tomography for monitoring cave propagation and stress distribution in Deep Mill Level Zone (DMLZ) PT Freeport Indonesia

Caving activity results in an increased induced seismicity which should be monitored to avoid massive and uncontrolled rock damage. This research was conducted at the Deep Mill Level Zone (DMLZ) underground mine, the deepest underground mine in Indonesia operated by PT Freeport Indonesia. This research aims to monitor cave propagation by using 4D tomography with a catalogue of microseismic for 57 days with a total of 14,821 events recorded by 84 stations consisting of 176,265 P phases and 133,472 S phases. The data is divided into four subsets to see the velocity evolution related to cave progress. Checkerboard Resolution Test (CRT) and Derivative Weight Sum (DWS) are used to assess the resolution of the inversion. 3D initial velocity model is constructed based on geological information and coring data. We have succeeded in identifying the interpreted cave propagation of a 60 m extension to the NW at around 100 m above undercut level based on 4D changes in velocity tomogram validate by Time Domain Reflectometry data. The decrease of Vp and Vs in subset 3 is interpreted due to the fracturing processes as the cave progresses. Furthermore, we observe a stress redistribution along with the progress of the cave, which is characterized by high velocities (Vp and Vs) due to compensation for changes in low velocity values ​​in the area in front of the cave, which is starting to collapse. We suggest that a considerable change in the velocity tomogram as an indicator of impending caving.

Evidence of strong wave turbulence and of Bolgiano temperature spectra in katabatic winds on steep slopes

The katabatic winds on steep slopes investigated in the present study reveal a novel spectral behavior, observed in the outer part of the jet. At low wavenumbers, the one-dimensional (1D) velocity spectra show evidence of a kx−1 range for the three components of the velocity vector: Eu(kx),Ev(kx),Ew(kx)∝kx−1 [as well as for the 1D temperature spectrum Eθ(kx)∝kx−1]. This suggests the existence of strong wave turbulence. A necessary condition for strong wave turbulence to be manifest is that the flow direction wavenumber, kx, extends to much lower values than the slope normal one, kz. This is satisfied in the present field experiment where wave energy is injected at wavenumber kx=kN=(Na sin α)/uj¯, while kz∼1/Δz, with Na the ambient stratification, α the slope angle, uj¯ the maximum wind velocity, and Δz the shear layer thickness of the jet. In the inertial range, the velocity spectra exhibit a power law kx−5/3 over two decades, whereas the temperature-buoyancy spectra show evidence of a −7/5 slope in the buoyancy sub-range, followed by a −5/3 slope. The change in spectral slopes occurs at the Bolgiano scale LB that is close to the Dougherty–Ozmidov scale LOZ. The high Reynolds number based on the Taylor micro-scale, Reλ∼103, allows clear identification of the spectral laws.