Tomographic Inversion Research Articles

Acoustic tomographic inversion of 3D temperature fields with mesoscale anomaly in the South China Sea

Acoustic tomographic inversion is based on travel times measured along the transmission paths between all station pairs to reconstruct three-dimensional temperature structures with mesoscale anomalies. In this study, tomographic simulation experiments were designed based on the Hybrid Coordinate Ocean Model (HYCOM) reanalysis data to reconstruct mesoscale phenomena from travel time data obtained from five, seven, and nine stations in the South China Sea over a domain of 100 × 100 km. The travel times for each station pair were calculated in the vertical section using the Bellhop acoustic ray simulation method. Six Empirical orthogonal function (EOF) modes of sound speed along the sound transmission paths in a vertical slice were used to formulate the inversion equations. The horizontal-slice distributions of temperature in the tomography domain were reconstructed using the grid-segmented method for each depth layer. For station-to-station distances greater than 100 km, the performance of inversion was best for the seven-station case rather than for the nine-station case, with the highest horizontal resolution of the three cases. This case study concluded that the seven-station case rather than the nine-station case provided an optimal station number for reconstructing the three-dimensional temperature fields.

On the intricacies of soot volume fraction measurements in counterflow diffusion flames with light extinction: Effects of curtain flow

Counterflow diffusion flame is an attractive platform for fundamental research on kinetics of soot formation. Accurate determination of soot volume fraction in the flame is a prerequisite for in-depth analysis of the sooting characteristics and assessment of predictive soot models. Light extinction has been proven to be an efficient technique for measuring soot volume fraction thanks to its non-intrusiveness and its simple optical setup. Nevertheless, tomographic inversion needs to be performed if spatially resolved soot volume fraction is to be obtained from the measured light extinction data which is essentially a projection along the line-of-sight. In this regard, radial distribution of soot volume fraction would affect the accuracy of the measurement through its influences on the inversion processes. In this work, we show that the curtain flow, which is necessary to avoid the formation of the undesired secondary diffusion flame and to keep the core counterflow from ambient disturbances, has notable effects on spatially resolved soot volume fraction measurements with line-of-sight measurements. In particular, different flow rate settings of the curtain flow can result in different soot distributions at the edges of soot fields: upwards curved, outwards extended, and downwards curved, which may influence the measurement of centerline soot volume fraction distribution. The necessity of tomographic inversion, the minimal region of the projection image necessary for tomographic inversion (when necessary), the quasi-one-dimensional feature of soot distribution, and the sensitivity of measurement to slight flame asymmetry were investigated where possible to determine the most suitable curtain flow configuration for soot volume fraction measurements by light extinction. Recommendations on curtain flow setting are finally made.

Subsurface anatomy of the Irazú–Turrialba volcanic complex, inferred from the integration of local and ambient seismic tomographic methods

SUMMARY Irazú and Turrialba are a twin volcanic complex that marks a distinct stop in volcanism along the Central America volcanic arc. We present a new traveltime velocity model of the crust beneath Irazú and Turrialba volcanoes, Costa Rica, and interpret it considering the results of previous ambient noise tomographic inversions. Data were acquired by a temporary seismic network during a period of low activity of the Irazú–Turrialba volcanic complex in 2018–2019. Beneath the Irazú volcano, we observe low P-wave velocities (VP = 5 km s−1) and low velocity ratios (VP/VS = 1.6). In contrast, below the Turrialba volcano, we observe low S-wave velocities (VS = 3 km s−1) and a high VP/VS (= 1.85) anomaly. We found that locations of low VP and VS anomalies (−15 %) correspond well with shear wave velocity anomalies retrieved from ambient noise tomography. At shallower depths, we observe high VP and VS anomalies (+15 %) located between the summits of the volcanoes. Subvertical velocity anomalies are also observed at greater depths, with high VP and VS anomalies appearing at the lower limits of our models. We propose a complex structure of an intermediate magmatic reservoir, presenting multiphase fluid states of a liquid-to-gas transition beneath Irazú and a juvenile store of magmatic fluid beneath Turrialba, while shallow fluid transport provides evidence of magmatic–hydrothermal interactions.

3D radiated power analysis of JET SPI discharges using the Emis3D forward modeling tool

Precise values for radiated energy in tokamak disruption experiments are needed to validate disruption mitigation techniques for burning plasma tokamaks like ITER and SPARC. Control room analysis of radiated power (P rad) on JET assumes axisymmetry, since fitting 3D radiation structures with limited bolometry coverage is an under-determined problem. In mitigated disruptions, radiation is toroidally asymmetric and 3D, due to fast-growing 3D MHD modes and localized impurity sources. To address this problem, Emis3D adopts a physics motivated forward modeling (‘guess and check’) approach, comparing experimental bolometry data to synthetic data from user-defined radiation structures. Synthetic structures are observed with the Cherab modeling framework and a best fit chosen using a reduced χ 2 statistic. 2D tomographic inversion models are tested, as well as helical flux tubes and 3D MHD simulated structures from JOREK. Two nominally identical pure neon shattered pellet injection (SPI) mitigated discharges in JET are analyzed. 2D tomographic inversions with added toroidal freedom are the best fits in the thermal quench (TQ) and current quench (CQ). In the pre-TQ, 2D reconstructions are statistically the best fits, but are likely over-optimized and do not capture the 3D radiation structure seen in fast camera images. The next-best pre-TQ fits are helical structures that extend towards the high-field side, consistent with an impurity flow under the magnetic nozzle effect also observed in JOREK simulations. Whole-disruption radiated fractions of 0.98+0.03/ −0.29 and 1.01+0.02/−0.17 are found, suggesting that the stored energy may have been fully mitigated by each SPI, although mitigation efficiencies well below ITER and SPARC requirements for high energy pulses are still within the large uncertainties. Emis3D is also used to validate JOREK SPI simulations, and confirms improvements in matching experiment from changes to impurity modeling. Time-dependent toroidal peaking factors are calculated and discussed.

Calorimetry measurement for energy balance and energy distribution in WEST for L-mode plasmas

This paper presents the energy balance of 602 pulses from four different experimental campaigns for the WEST tokamak. Different magnetic configurations have been studied, with lower single null (LSN) and upper single null (USN) configuration with deuterium or helium plasmas. The energy balance is closed with an imbalance of about 5% of the total injected energy for most of the campaigns and for different magnetic configurations. The distribution over the whole machine is shown, with the outer first wall receiving most of the energy due to its large surface area with about 30% of the total heat load, and the divertor with 25% due to the heat loads deposited by the convected power in the scrape-off layer (SOL). Finally, the tomography inversion of the bolometry measurement allows us to disentangle the contribution of the radiated and convected power in the energy absorbed by each type of plasma-facing component. We show that in the USN configuration about 63% of the available energy in the SOL is deposited in the upper divertor (UDIV) through convected heat loads, while in LSN this value is spread over the lower divertor with 45% and the baffle and UDIV with about 10% for both.

Topography effect on seismic waveform tomography: a quantitative study

SUMMARY In seismic tomography practices the Earth's surface is sometimes assumed to be either spherical or flat for convenience in forward modelling calculations. The effect of irregular surface topography on seismic wave propagation is thus ignored, resulting in biases in the phases and amplitudes of synthetic seismograms, which contribute to the residuals that are mapped into velocity structures in tomography inversions. In this study, we conduct a series of inversion experiments based on the adjoint waveform tomography method to quantitatively assess the topography effect on waveform-based inversion results. We first employ models with simplified topography to better highlight and quantify the topography effect. Results show that when topography effect is ignored in the forward modelling, it is mapped into velocity perturbations, leading to spurious velocity anomalies in tomography models. The strength of the spurious velocity anomaly is quasi-linearly related to locally averaged topography gradient. Our inversion experiments demonstrate that in places of strong topography variation, such as the Longmenshan Fault Zone region where the 2008 Mw 7.9 Wenchuan earthquake occurred, topography effect can lead to spurious relative velocity anomalies of up to 10 per cent, which cannot be ignored in waveform-based tomography inversions.

SItomo – A toolbox for splitting intensity tomography and application in the Eastern Alps

The tomographic inversion of shear wave splitting data for upper mantle anisotropy has been a longstanding challenge. This is due to the ray-based approximation of classical approaches and the near-vertical incidence of the core-mantle converted phases such as SKS that are often used. Recent developments include the calculation of finite-frequency sensitivity kernels for SKS splitting intensity observations, which allows us to accurately take into account the sensitivity to anisotropic structure with depth. A requirement of this tomographic technique is a dense station spacing, which results in overlapping sensitivity kernels at depth and allows for the localization of anisotropic structure. This is satisfied by a growing number of temporary seismic deployments, which motivates the desire to image anisotropic complexities with depth. Here, we introduce and make available a toolbox for the MATLAB environment that facilitates the application of finite-frequency splitting intensity tomography to dense seismic arrays. Our implementation includes several key features, including: 1) A forward calculation of splitting intensities and sensitivity kernels for a complex anisotropic model space. 2) Consideration of the dominant period of the wave, allowing for multiple-frequency analysis, as well as the incoming wave’s non-vertical incidence. 3) The inversion can be based on a classical gradient descent, on a form of the conjugate gradient method known as the BFGS algorithm, or on a gradient-informed stochastic reversible jump algorithm, allowing for a data-driven parametrization of the model space. 4) Importing splitting intensity measurements from waveforms processed in SplitRacer allows for fast pre-processing of large data sets due to its fully automatic design. To illustrate our method, we present both synthetic tests and an application to real data. We apply our inversion procedure to data from the Swath-D network, which densely covers the transition of the Central to the Eastern Alps. Previous studies showed evidence for an abrupt lateral change of layered seismic anisotropy that had been attributed to an opening for channeled asthenospheric flow. Using an SKS splitting intensity tomography approach, we can confirm previous inferences while providing additional constraints on the distribution of anisotropy laterally and with depth.

Multi-Sensor Seismic Processing Approach Using Geophones and HWC DAS in the Monitoring of CO2 Storage at the Hellisheiði Geothermal Field in Iceland

Geothermal power production may result in significant CO2 emissions as part of the produced steam. CO2 capture, utilisation, subsurface storage (CCUS) and developments to exploit geothermal resources are focal points for future clean and renewable energy strategies. The Synergetic Utilisation of CO2 Storage Coupled with Geothermal Energy Deployment (SUCCEED) project aims to demonstrate the feasibility of using produced CO2 for re-injection in the geothermal field to improve geothermal performance, while also storing the CO2 as an action for climate change mitigation. Our study has the aim to develop innovative reservoir-monitoring technologies via active-source seismic data acquisition using a novel electric seismic vibrator source and permanently installed helically wound cable (HWC) fibre-optic distributed acoustic sensing (DAS) system. Implemented together with auxiliary multi-component (3C and 2C) geophone receiver arrays, this approach gave us the opportunity to compare and cross-validate the results using wavefields from different acquisition systems. We present the results of the baseline survey of a time-lapse monitoring project at the Hellisheiði geothermal field in Iceland. We perform tomographic inversion and multichannel seismic processing to investigate both the shallower and the deeper basaltic rocks targets. The wavefield analysis is supported by seismic modelling. The HWC DAS and the geophone-stacked sections show good consistency, highlighting the same reflection zones. The comparison of the new DAS technology with the well-known standard geophone acquisition proves the effectiveness and reliability of using broadside sensitivity HWC DAS in surface monitoring applications.

Surface distributed acoustic sensing for high-resolution near-surface characterization: Results from a 3D onshore field experiment

Distributed acoustic sensing technology enables high-density seismic acquisition at a fraction of the cost. When deployed on the surface, surface distributed acoustic sensing (S-DAS) acquisition provides a cost-effective solution for dense high-resolution near-surface characterization through the analysis and inversion of surface waves. This is made possible by the relatively low cost of the fiber and the dense spatial sampling of the realized seismic data. The S-DAS data were collected during the acquisition of a 3D land large-scale field test and processed with a focus on recent advancements in the use of surface wave analysis and inversion. We compare and validate the result from the S-DAS recording with colocated multicomponent (3C) geophones and a conventional high-density surface seismic nodal acquisition. The comparison to 3C geophones demonstrated that for applications such as surface-wave inversion, S-DAS can outperform conventional geophones and shows consistency between electrical resistivity tomography and surface seismic inversion from S-DAS. In addition, continuous passive recording of environmental noise also offers a convenient alternative to active shooting allowing for surface wave inversion from reconstructed virtual shots.

The 2023 Turkey earthquake doublet: Earthquake relocation, seismic tomography, and stress field inversion

The 2023 Turkey earthquake doublet: Earthquake relocation, seismic tomography, and stress field inversion

Bayesian Inversion, Uncertainty Analysis and Interrogation Using Boosting Variational Inference

AbstractGeoscientists use observed data to estimate properties of the Earth's interior. This often requires non‐linear inverse problems to be solved and uncertainties to be estimated. Bayesian inference solves inverse problems under a probabilistic framework, in which uncertainty is represented by a so‐called posterior probability distribution. Recently, variational inference has emerged as an efficient method to estimate Bayesian solutions. By seeking the closest approximation to the posterior distribution within any chosen family of distributions, variational inference yields a fully probabilistic solution. It is important to define expressive variational families so that the posterior distribution can be represented accurately. We introduce boosting variational inference (BVI) as a computationally efficient means to construct a flexible approximating family comprising all possible finite mixtures of simpler component distributions. We use Gaussian mixture components due to their fully parametric nature and the ease with which they can be optimized. We apply BVI to seismic travel time tomography and full waveform inversion, comparing its performance with other methods of solution. The results demonstrate that BVI achieves reasonable efficiency and accuracy while enabling the construction of a fully analytic expression for the posterior distribution. Samples that represent major components of uncertainty in the solution can be obtained analytically from each mixture component. We demonstrate that these samples can be used to solve an interrogation problem: to assess the size of a subsurface target structure. To the best of our knowledge, this is the first method in geophysics that provides both analytic and reasonably accurate probabilistic solutions to fully non‐linear, high‐dimensional Bayesian full waveform inversion problems.

Lithospheric underthrusting beneath the northern Eastern Cordillera plateau, Colombian Andes: Insights from analyzing tomographic inversion using local earthquake travel time and gravity data

The northern Eastern Cordillera in the Colombian Andes is a wide (~200 km) high-elevation (>2.5 km) subduction-related asymmetric plateau. This plateau is currently subjected to shortening and magmatic addition, and yet the nature of the lowermost crust and its transition to the underlying mantle remains poorly constrained. To improve our knowledge about the structure of the deep crust beneath the EC plateau, we conducted a joint inversion of travel times of local earthquakes and gravity data. Gravity contributes to up to 0.5% in dVp and dVs with no affectation of the shape or location of the discussed anomalies. Along the latitudinal profile with the highest resolution (~5.7°N), two anomalies are identified at depths of 40–60 km beneath the plateau. A western low-velocity anomaly is interpreted as crustal material underthrust eastward beneath the northwestern EC. This process is triggering the abrupt change in topography between the adjacent low-elevation basin and the orogenic plateau. Additionally, a high-velocity anomaly beneath the eastern flank of the EC is likely related to mantle metasomatism and westward underthrusting of the foreland lithosphere. This metasomatism is either related to the interaction between mafic magmas and the uppermost mantle or silica enrichment that occurred during a past episode of flat subduction. Evidence for foreland underthrusting includes relocated earthquakes at 33–42 km depth within the thrust system between the EC and the foreland, along with increased shortening in the eastern part of the plateau, an increase in exhumation rates, the eastward migration of deformation towards the foreland, and dominated thick-skinned deformation in the upper crust of the foreland with an eastward vergence of thrusts and related folds. Furthermore, within the central part of the plateau, our results show low velocities between 30 and 40 km depth, consistent with previous constraints, suggesting magmatic underplating beneath the Paipa-Iza volcanic complex.

Hydrologic Impacts of a Strike-Slip Fault Zone: Insights from Joint 3D Body-Wave Tomography of Rock Valley

ABSTRACT The Rock Valley fault zone (RVFZ), an intraplate strike-slip fault zone in the southern Nevada National Security Site (NNSS), hosted a series of very shallow (<3 km) earthquakes in 1993. The RVFZ may also have hydrological significance within the NNSS, potentially playing a role in regional groundwater flow, but there is a lack of local hydrological data. In the Spring of 2021, we collected active-source accelerated weight drop seismic data over part of the RVFZ to better characterize the shallow subsurface. We manually picked ∼17,000 P-wave travel times and over 14,000 S-wave travel times, which were inverted for P-wave velocity (VP), S-wave velocity (VS), and VP/VS ratio in a 3D joint tomographic inversion scheme. Seismic velocities are imaged as deep as ∼700 m in areas and generally align with geologic and structural expectations. VP and VS are relatively reduced near mapped and inferred faults, with the most prominent lower VP and VS zone around the densest collection of faults. We image VP/VS ratios ranging from ∼1.5 to ∼2.4, the extremes of which occur at a depth of ∼100 m and are juxtaposed across a fault. One possible interpretation of the imaged seismic velocities is enhanced fault damage near the densest collection of faults with relatively higher porosity and/or crack density at ∼100 m depth, with patches of semiperched groundwater present in the sedimentary rock in higher VP/VS areas and drier rock in lower VP/VS areas. A relatively higher VP/VS area beneath the densest faults persists at depth, which suggests percolation of groundwater via the fault damage zone to the regionally connected lower carbonate aquifer. Potentially, the presence and movement of groundwater may have played a role in the 1993 earthquake aftershocks.

Anisotropic structure in the mantle wedge beneath southeastern Mexico from shear-wave splitting tomography

Cocos intraslab earthquakes were used to make a 3-D tomographic inversion to define a crystallographic orientation model for the mantle wedge beneath southeastern Mexico. This model provided insights regarding the pattern of the mantle wedge flow and its relationship to the geometry of the subducting slab. The mantle wedge was parametrized as a 3-D block model of crystallographic orientations assuming the elastic constants of olivine and orthopyroxene with orthorhombic symmetry (hexagonal symmetry was also tested). A linearized, damped, and iterative least-squares approach was used to account for the nonlinear behavior of the shear-wave splitting, numerically recalculating partial derivatives after each iteration. The best-fitting model is consistent with two main flow regimes: (1) 2-D corner flow in a mantle wedge core made up of A-type olivine fabric northwest of the Tehuantepec Ridge extension, and (2) 3-D trench-parallel mantle flow in a mantle wedge core made up of A-, C-, or E-type olivine fabric southeast of this geological feature. Around the Tehuantepec Ridge extension, a partially serpentinized mantle wedge tip is inferred since olivine a-axis orientations are trench-parallel regardless of whether a 2-D corner flow or a 3-D trench-parallel flow prevails. Right above the Tehuantepec Ridge extension (beyond the 100 km isodepth contour of the subducting slab), a change of well-resolved olivine a-axis orientations from trench-normal to trench-parallel while going from northwest to southeast is observed. It signals an abrupt change in the mantle flow pattern possibly through a vertical tear in the Cocos slab. 3-D toroidal flow could be driving subslab mantle material into the mantle wedge around the deepest slab segment. Lastly, approximately trench-normal olivine a-axis orientations are observed in the mantle wedge tip near the Mexico and Guatemala border region, and they could be explained by assuming the presence of B-type olivine fabric.

3D shear-wave velocity structure of the Eastern Indian Shield from ambient noise tomography: Ultramafic lower crust underneath the Singhbhum Craton

Singhbhum Craton (SC) in the Eastern Indian Shield (EIS) is known to have preserved Archean crust which might provide insights into the Precambrian tectonics. In this study, ambient noise cross-correlation of about 1 year of continuous data recorded in 15 broadband seismic stations across EIS is performed. Rayleigh wave group velocities at 4–22 s periods obtained by frequency-time analysis of the cross-correlograms were subjected to tomographic inversion to produce group velocity maps. The spatial resolution of these maps was found to be ∼100 km. Inverted 1D shear-wave velocity (VS) models at regular grid points of these tomographic maps were combined to construct the first 3D VSmodel of the SC and adjoining regions within the EIS. The 3D model highlights a high VS(3.6–3.8 km/s) structure in the 4–20 km depth range and an ultra-high VS(∼4.0–4.2 km/s) structure at depths ≥ 20 km within the SC only. The VS of the upper layer and its regionalisation match with the well-studied granite batholith within the SC. The ultra-high velocitylayer in the lower crust, spatially spanning ∼100 km × 150 km, is distinctly imaged for the first time and its VS indicates probable ultramafic composition. The existence of ultramafic composition at ∼800 MPa lithostatic pressure requires ultra-high temperature metamorphism, possible through mantle upwelling/melting. The previous studies of upwarped crust-lithosphere beneath the SC and the lower-crustal ultramafic composition inferred in this study probably confirm the delamination of the lithosphere beneath the SC at some point during the Precambrian eon.

The Role of Subslab Low‐Velocity Anomalies Beneath the Nazca Ridge and Iquique Ridge on the Nazca Plate and Their Possible Contribution to the Subduction Angle

AbstractSubducting the buoyant crustal material of an aseismic oceanic ridge has been regarded as a dominant contributor to flat slab subduction. However, normal‐dip subduction is also observed in some cases where ridges are subducting. In this study, we compare the subduction of two ridges on the Nazca Plate: Nazca Ridge (flat slab) and Iquique Ridge (normal‐dip slab). Anisotropy determined by shear wave splitting observation suggests that the low‐velocity anomalies found beneath the ridges are mapping anisotropic structure into isotropic velocities. After a tomographic inversion incorporating anisotropy models for both ridges, we find that the low‐velocity anomalies found beneath the Nazca Ridge are not anisotropic and therefore likely represent warm mantle, and those beneath the Iquique Ridge are caused by anisotropy. We conclude that subslab mantle buoyancy has a larger impact on the subduction angle than the crustal material of the ridge.

Optimizing Boundary Conditions in GNSS Tomography: A Continuous 7‐Month Case Study in the Amazon

AbstractDifferent estimates of the regional water vapor scale height, taken from ERA5 reanalysis, in situ observations and the direct optimization of retrieved water vapor profiles in GNSS tomography, are found to have major impact in the performance of tomographic inversions, with the better results displaying mean errors comparable to radiosondes. The analysis uses 7 months of GNSS (Global Navigation Satellite Systems) observations in the Amazon Dense GNSS Network near Manaus, Brazil, in 2011–2012, to compute a time series of water vapor profiles, with a tomographic technique capable of producing quasi‐instantaneous inversions with minimal external data or constraints. Results compare very well with 12‐hourly in situ radiosondes, especially in the lower troposphere above 2 km, and its daily‐to‐seasonal variability compares well with WRF (Weather Research and Forecasting) convective‐permitting simulations driven by ERA5 boundary conditions, suggesting that GNSS tomography may be an important source of atmospheric water vapor data for different applications.

Assessment of the ITER divertor bolometer diagnostic performance

The ITER bolometer diagnostic will provide measurements of the total radiation emitted from the plasma, a part of the overall energy balance. It will consist of some 550 lines of sight (LOS) of bolometer detector, bundled in multiple individual cameras, which will be located in the gaps between blanket modules on the vacuum vessel wall, in five divertor cassettes, in two upper port plugs and in one equatorial port plug (Meister et al., 2017). The LOS are optimised as much as possible with design constraints to capture the details of different plasma regions in which the intensity and local distribution of radiation will vary over a large range. This is especially true in the divertor, where up to 70 % of the total thermal exhaust power from the core will have to be radiated during burning plasma operation. This radiation will be tightly concentrated meaning that not only will tomographic reconstruction of the distribution be challenging, but the bolometer cameras themselves need to be properly designed to handle the resulting heat fluxes. This paper first presents an assessment of the photonic heat fluxes falling on bolometers in the divertor region during high performance operation (where the heat fluxes will be highest). These heat fluxes are computed using ray-tracing from radiation distributions obtained with the SOLPS (plasma boundary region) and JINTRAC (plasma core region) codes and will serve as essential input to the thermal design studies of bolometers. The same SOLPS and JINTRAC simulations are used as models for the second focus of the paper which examines how closely the radiation distribution obtained from tomographic inversion of synthetic signals from the bolometer LOS matches the model input. Some effort has also been devoted to a study of how the reconstructions are affected by loss of individual LOS.

CSTNet: A Dual-Branch Convolutional Neural Network for Imaging of Reactive Flows Using Chemical Species Tomography.

Chemical species tomography (CST) has been widely used for in situ imaging of critical parameters, e.g., species concentration and temperature, in reactive flows. However, even with state-of-the-art computational algorithms, the method is limited due to the inherently ill-posed and rank-deficient tomographic data inversion and by high computational cost. These issues hinder its application for real-time flow diagnosis. To address them, we present here a novel convolutional neural network, namely CSTNet, for high-fidelity, rapid, and simultaneous imaging of species concentration and temperature using CST. CSTNet introduces a shared feature extractor that incorporates the CST measurements and sensor layout into the learning network. In addition, a dual-branch decoder with internal crosstalk, which automatically learns the naturally correlated distributions of species concentration and temperature, is proposed for image reconstructions. The proposed CSTNet is validated both with simulated datasets and with measured data from real flames in experiments using an industry-oriented sensor. Superior performance is found relative to previous approaches in terms of reconstruction accuracy and robustness to measurement noise. This is the first time, to the best of our knowledge, that a deep learning-based method for CST has been experimentally validated for simultaneous imaging of multiple critical parameters in reactive flows using a low-complexity optical sensor with a severely limited number of laser beams.

Unifying VLF Transmitter and Sferic Modeling Efforts Via Tomography

AbstractWe demonstrate a methodology for utilizing measurements from very low frequency (VLF, 3−30 kHz) transmitters and lightning emissions to produce 3D lower electron density maps, and apply it to multiple geophysical disturbances. The D‐region lower ionosphere (60−90 km) forms the upper boundary of the Earth‐ionosphere waveguide which allows VLF radio waves to propagate to global distances. Measurements of these signals have, in many prior studies, been used to infer path‐average electron density profiles within the D region. Historically, researchers have focused on either measurements of VLF transmitters or radio atmospherics (sferics) from lightning. In this work, we build on recently published methods for each and present a method to unify the two approaches via tomography. The output of the tomographic inversion produces maps of electron density over a large portion of the United States and Gulf of Mexico. To illustrate the benefits of this unified approach, daytime and nighttime maps are compared between a sferic‐only model and the new approach suggested here. We apply the model to characterize two geophysical disturbances: solar flares and lower ionospheric changes associated with thunderstorms.