Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “seismic anisotropy”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Seismic monitoring of CO 2 geosequestration using multi-well 4D DAS VSP: Stage 3 of the CO2CRC Otway project

An important part of any CO 2 geosequestration project is to ensure CO 2 containment and conformance in the subsurface. This is generally done by implementing a comprehensive, risk-based Measurement, Monitoring and Verification plan, a key element of which is active time-lapse seismic monitoring. However, high cost and environmental impact of the standard surface seismic monitoring dictate the need for a cost-effective and environmentally friendly alternative. An opportunity to develop such method emerges with advances in distributed acoustic sensing (DAS) technology, which turns an optical fibre into a seismic sensor with dense spatial sampling. DAS can be permanently deployed in multiple wells across the geosequestration site providing a robust and non-intrusive network of seismic receivers. This approach was developed and tested in the CO2CRC Otway project, where injection of 15 kt of CO 2 at 1.5 km depth was monitored with a 4D vertical seismic profiling (VSP) using five borehole DAS arrays and mobile vibroseis sources. The 4D DAS VSP in each of the five wells provides broadly consistent images of the CO 2 plume with some differences due to different illumination of the target horizon, lateral variation of velocities, and seismic anisotropy. When the newly injected CO 2 reaches a CO 2 plume created as a result of an earlier injection into the same formation ~600 m updip, 4D DAS VSP shows a change in reflectivity in that area and beyond. Furthermore, this shows a potential of 4D DAS VSP for monitoring gas injection into gas-saturated reservoirs.

58 GEOSCIENCES↗
DOE OSTI2021

Seismic monitoring of a small CO 2 injection using a multi-well DAS array: Operations and initial results of Stage 3 of the CO 2 CRC Otway project

Active time-lapse seismic is widely employed for monitoring CO 2 geosequestration due to its ability to track the distribution of fluids in space and time. However, standard 4D seismic monitoring suffers from several challenges, including high cost, disruption to other land uses, and, consequently, relatively large intervals between monitor surveys. Some of these challenges can be mitigated using permanently installed sources and receivers. Such an approach was tested at the CO 2 CRC Otway site by continuous offset VSP monitoring of 15,000 t of supercritical CO 2 injected into an aquifer 1,500 m deep with nine permanent seismic sources (surface orbital vibrators or SOVs) and five downhole fibre-optic receivers. This continuous monitoring is complemented by multi-well 4D VSP using a mobile vibroseis source and the same DAS receivers, which included one baseline and two monitor surveys after injection of 4,000 and 12,000 t of CO 2 . The continuous DAS-SOV monitoring detected an abrupt increase of travel times below the injection interval on the second day of injection (after injection of 300 t of CO 2 ) and tracked the growth of the areal CO 2 plume by mapping changes of reflection amplitudes. The plume is also detected by time-lapse changes of reflection amplitudes in multi-well 4D VSPs. The plume images obtained from continuous offset VSP and 4D VSP are broadly consistent with each other but with some differences due to differences in illumination, lateral variations of velocities and seismic anisotropy. Furthermore, these differences also serve as a measure of uncertainty of 4D VSP images.

4D VSP↗
DOE OSTI2021

Seismic Visibility of Eclogite in the Earth’s Upper Mantle—Implications From High Pressure‐Temperature Single‐Crystal Elastic Properties of Omphacite

Abstract Identifying and locating the geochemical and geophysical heterogeneities in the Earth’s interior is one of the most important and challenging tasks for the deep Earth scientists. Subducted oceanic crust metamorphizes into the dense eclogite in the upper mantle and is considered as a major cause of geochemical and geophysical heterogeneities in the deep Earth. In order to detect eclogitic materials inside the Earth, precise measurements of the high pressure‐temperature single‐crystal elasticity of major minerals in eclogite are thus exceedingly important. Omphacite, a Na,Al‐bearing clinopyroxene, constitutes up to 75 vol% of eclogite. In the present study, we performed the first high pressure‐temperature single‐crystal elasticity measurements of omphacite using Brillouin spectroscopy. Utilizing the finite‐strain approach, we obtained the following thermoelastic parameters for omphacite: K S0 ’ = 4.5(1), G 0 ’ = 1.53(5), ∂K S0 /∂ T = −0.029(5) GPa/K, ∂G 0 /∂ T = −0.013(5) GPa/K, with K S0 = 123(3) GPa, G 0 = 74(2) GPa, and ρ 0 = 3.34(1) g/cm 3 . We found that the seismic velocities of undeformed eclogite are similar to pyrolite at the depths of 200–300 and 410–500 km, thus eclogite is seismically invisible at these depths. Combined with the lattice‐preferred orientations of the omphacite in naturally deformed eclogites, we also modeled seismic anisotropy of eclogite at various pressure‐temperature conditions. A 10 km thick subducted eclogitic crust can result in ∼0.2 s shear wave splitting in the Earth’s upper mantle.

Hao, Ming↗
DOE OSTI2020

Global adjoint tomography—model GLAD-M25

SUMMARY Building on global adjoint tomography model GLAD-M15, we present transversely isotropic global model GLAD-M25, which is the result of 10 quasi-Newton tomographic iterations with an earthquake database consisting of 1480 events in the magnitude range 5.5 ≤ Mw ≤ 7.2, an almost sixfold increase over the first-generation model. We calculated fully 3-D synthetic seismograms with a shortest period of 17 s based on a GPU-accelerated spectral-element wave propagation solver which accommodates effects due to 3-D anelastic crust and mantle structure, topography and bathymetry, the ocean load, ellipticity, rotation and self-gravitation. We used an adjoint-state method to calculate Fréchet derivatives in 3-D anelastic Earth models facilitated by a parsimonious storage algorithm. The simulations were performed on the Cray XK7 ‘Titan’ and the IBM Power 9 ‘Summit’ at the Oak Ridge Leadership Computing Facility. We quantitatively evaluated GLAD-M25 by assessing misfit reductions and traveltime anomaly histograms in 12 measurement categories. We performed similar assessments for a held-out data set consisting of 360 earthquakes, with results comparable to the actual inversion. We highlight the new model for a variety of plumes and subduction zones.

58 GEOSCIENCES↗
DOE OSTI2021

Characterization of damage processes in Montney siltstone under triaxial compression using acoustic emission and diagnostic imaging

SUMMARY Crack nucleation and rock failure processes in a fine-grained siltstone (Montney Formation) under triaxial compression are investigated using combined diagnostic techniques, including ultrasonic-wave measurement, acoustic-emission (AE) monitoring, computed tomography (CT) scanning, and thin-section imaging. The sample displays a weak-to-moderate inherent seismic anisotropy and noticeable stress-induced anisotropy prior to failure. No AE event was detected until the applied axial stress reached 95 per cent of the peak value. The signal-to-noise ratio is relatively low, however, and detectable AE events are more diffuse than those observed in highly brittle rocks. The AE locations correlate with a shear fracture zone imaged by CT scanning. AE moment-tensor analysis reveals that events with larger relative magnitudes are characterized by high volumetric (tensile or compressive) components, and the initiation of the failure zone is dominated by combined shear-tensile (ST) failure. Stress inversion of the AE events with high tensile components is in good agreement with the known applied stress. Microscopic imaging of thin sections from the failed sample shows that the failure zone is an en echelon structure consisting of a major fracture with branching micro and minor cracks. This failure mechanism is consistent with a ST source mechanism and is interpreted to be associated with the fine granular structure and mineral composition of Montney siltstone.

Jia, Suzie Q. (ORCID:0000000326496190)↗
DOE OSTI2023

Crystallographic texture formation in Fe-9wt%Si alloy during deformation and phase transition at high pressure

SUMMARY The seismic anisotropy of the Earth's solid inner core has been the topic of much research. It could be explained by the crystallographic preferred orientation (CPO) developing during convection. The likely phase is hexagonal close-packed iron (hcp), alloyed with nickel and some lighter elements. Here we use high energy synchrotron X-rays to study CPO in Fe-9wt%Si, uniaxially compressed in a diamond anvil cell in radial geometry. The experiments reveal that strong preferred orientation forms in the low-pressure body-centred cubic (bcc) phase that appears to be softer than pure iron. CPO is attributed to dominant {110}<111> slip. The onset of the bcc→hcp transition occurs at a pressure of ≈15 GPa, and the alloy remains in a two phase bcc + hcp state up to 40 GPa. The hcp phase forms first with a distinct {11$\bar{2}$0} maximum perpendicular to compression. Modelling shows that this is a transformation texture, which can be described by Burgers orientation relationship with variant selection. Experimental results suggest that bcc grains oriented with <100> parallel to compression transform into hcp first. The CPO of the hcp changes only slowly during further pressure and deviatoric stress increase at ambient temperature. After heating to 1600 K, a change in the hcp CPO is observed with alignment of (0001) planes perpendicular to compression that can be interpreted as dominant (0001)<11$\bar{2}$0> slip, combined with {10$\bar{1}$2}<$\bar{1}$011> mechanical twinning, which is similar to the deformation modes suggested previously for pure hcp iron at inner core conditions.

Vasin, Roman N.↗
NASA NTRS1990

The nature of the earth's core

The properties of the earth's core are overviewed with emphasis on seismologically determined regions and pressures and seismologically measured density, elastic wave velocities, and gravitational acceleration. Attention is given to solid-state convection of the inner core, and it is noted that though seismological results do not conclusively prove that the inner core is convective, the occurrence and magnitude of seismic anisotropy are explained by the effects of solid-state convection. Igneous petrology and geochemistry of the inner core, a layer at the base of the mantle and contact metasomatism at the core-mantle boundary, and evolution of the core-mantle system are discussed. It is pointed out that high-pressure melting experiments indicate that the temperature of the core is ranging from 4500 to 6500 K, and a major implication of such high temperature is that the tectonics and convection of the mantle, as well as the resulting geological processes observed at the surface, are powered by heat from the core. As a result of the high temperatures, along with the compositional contrast between silicates and iron alloy, the core-mantle boundary is considered to be most chemically active region of the earth.

Jeanloz, Raymond↗
NASA NTRS2019

Water Storage and Transport Processes During Impact Experiments on Nominally Anhydrous Minerals

Water, in the form of structurally bound hydrogen in the crystal lattice of nominally anhydrous minerals (NAMs), strongly influences many important physical processes on terrestrial planets and planetary objects. Water enhances the rates of plastic deformation and controls the degree of partial melting in silicate rocks, which influences the generation of melt and therefore the nature of planetary volcanism. Water has also been experimentally demonstrated to influence the nature of lattice preferred orientation in deformed aggregates, and thus may be important in the interpretation of seismic anisotropy data collected from planetary bodies, such as from the current InSight mission on Mars. Therefore, much attention has been focused on characterizing the distribution and concentration of water in the planets and rocky bodies of our solar system.

Tielke, J.↗
DOE OSTI2021

SPiRaL: a multiresolution global tomography model of seismic wave speeds and radial anisotropy variations in the crust and mantle

SUMMARY SPiRaL is a joint global-scale model of wave speeds (P and S) and anisotropy (vertical transverse isotropy, VTI) variations in the crust and mantle. The model is comprised of >2.1 million nodes with five parameters at each node that capture velocity variations for P- and S-waves travelling at arbitrary directions in transversely isotropic media with a vertical symmetry axis (VTI). The crust (including ice, water, sediments and crystalline layers) is directly incorporated into the model. The default node spacing is approximately 2° in the lower mantle and 1° in the crust and upper mantle. The grid is refined with ∼0.25° minimum node spacing in highly sampled regions of the crust and upper mantle throughout North America and Eurasia. The data considered in the construction of SPiRaL includes millions of body wave traveltimes (crustal, regional and teleseismic phases with multiples) and surface wave (Rayleigh and Love) dispersion. A multiresolution inversion approach is employed to capture long-wavelength heterogeneities commonly depicted in global-scale tomography images as well as more localized details that are typically resolved in more focused regional-scale studies. Our previous work has demonstrated that such global-scale models with regional-scale detail can accurately predict both teleseismic and regional body wave traveltimes, which is necessary for more accurate location of small seismic events that may have limited signal at teleseismic distances. SPiRaL was constructed to predict traveltimes for event location and long-period waveform dispersion for seismic source inversion applications in regions without sufficiently tuned models. SPiRaL may also serve as a starting model for full-waveform inversion (FWI) with the goal of fitting waves with periods 10–50 s over multiple broad regions (thousands of kilometres) and potentially the globe. To gain insight to this possibility, we simulated waveforms for a small set of events using SPiRaL and independent waveform-based models for comparison. For the events tested, the performance of the traveltime-based SPiRaL model is shown to be generally on par with regional 3-D waveform-based models in three regions (western United States, Middle East, Korean Peninsula) suggesting SPiRaL may serve as a starting model for FWI over broad regions.

58 GEOSCIENCES↗
DOE OSTI2023

Seismic Characterization of the Blue Mountain Geothermal Field

Subsurface characterization is crucial for geothermal energy exploration and production. Yet hydrothermal reservoirs usually reside in highly fractured and faulted zones where accurate characterization is very challenging because of low signal-to-noise ratios of land seismic data and lack of coherent reflection signals. We perform an active-source seismic characterization for the Blue Mountain geothermal field in Nevada using active seismic data to reveal the elastic medium property complexity and fault distribution at this field. We first employ an unsupervised machine learning method to attenuate groundroll and near-surface guided-wave noise and enhance coherent reflection and scattering signals from noisy seismic data. We then build a smooth initial P-wave velocity model based on an existing magnetotellurics survey result, and use 3D first-arrival traveltime tomography to refine the initial velocity model. We then derive a set of elastic wave velocities and anisotropic parameters using elastic full-waveform inversion, and obtain PP and PS images using elastic reverse-time migration. We identify major faults by analyzing the variations of seismic velocities and anisotropy parameters, and reveal mid- to small-scale faults by applying a supervised machine learning method to the seismic migration images. Our characterization reveals complex velocity heterogeneities and anisotropies, as well as faults, with a high spatial resolution. These results can provide valuable information for optimal placement of future injection and production wells to increase geothermal energy production at the Blue Mountain geothermal power plant.

58 GEOSCIENCES↗
DOE OSTI2024

Heterogeneous mantle effects on the behavior of SmKS waves and outermost core imaging

Seismic traveltime anomalies of waves that traverse the uppermost 100–200 km of the outer core have been interpreted as evidence of reduced seismic velocities (relative to radial reference models) just below the core–mantle boundary (CMB). These studies typically investigate differential traveltimes of SmKS waves, which propagate as P waves through the shallowest outer core and reflect from the underside of the CMB m times. The use of SmKS and S(m-1)KS differential traveltimes for core imaging are often assumed to suppress contributions from earthquake location errors and unknown and unmodelled seismic velocity heterogeneity in the mantle. The goal of this study is to understand the extent to which differential SmKS traveltimes are, in fact, affected by anomalous mantle structure, potentially including both velocity heterogeneity and anisotropy. Velocity variations affect not only a wave's traveltime, but also the path of a wave, which can be observed in deviations of the wave's incoming direction. Since radial velocity variations in the outer core will only minimally affect the wave path, in contrast to other potential effects, measuring the incoming direction of SmKS waves provides an additional diagnostic as to the origin of traveltime anomalies. Here we use arrays of seismometers to measure traveltime and direction anomalies of SmKS waves that sample the uppermost outer core. We form subarrays of EarthScope's regional Transportable Array stations, thus measuring local variations in traveltime and direction. We observe systematic lateral variations in both traveltime and incoming wave direction, which cannot be explained by changes to the radial seismic velocity profile of the outer core. Moreover, we find a correlation between incoming wave direction and traveltime anomaly, suggesting that observed traveltime anomalies may be caused, at least in part, by changes to the wave path and not solely by perturbations in outer core velocity. Modelling of 1-D ray and 3-D wave propagation in global 3-D tomographic models of mantle velocity anomalies match the trend of the observed traveltime anomalies. Overall, we demonstrate that observed SmKS traveltime anomalies may have a significant contribution from 3-D mantle structure, and not solely from outer core structure.

58 GEOSCIENCES↗
DOE OSTI2024

Heterogeneous mantle effects on the behavior of SmKS waves and outermost core imaging

Seismic traveltime anomalies of waves that traverse the uppermost 100–200 km of the outer core have been interpreted as evidence of reduced seismic velocities (relative to radial reference models) just below the core–mantle boundary (CMB). These studies typically investigate differential traveltimes of SmKS waves, which propagate as P waves through the shallowest outer core and reflect from the underside of the CMB m times. The use of SmKS and S(m- 1)KS differential traveltimes for core imaging are often assumed to suppress contributions from earthquake location errors and unknown and unmodelled seismic velocity heterogeneity in the mantle. The goal of this study is to understand the extent to which differential SmKS traveltimes are, in fact, affected by anomalous mantle structure, potentially including both velocity heterogeneity and anisotropy. Velocity variations affect not only a wave’s traveltime, but also the path of a wave, which can be observed in deviations of the wave’s incoming direction. Since radial velocity variations in the outer core will only minimally affect the wave path, in contrast to other potential effects, measuring the incoming direction of SmKS waves provides an additional diagnostic as to the origin of traveltime anomalies. Here we use arrays of seismometers to measure traveltime and direction anomalies of SmKS waves that sample the uppermost outer core. We form subarrays of EarthScope’s regional Transportable Array stations, thus measuring local variations in traveltime and direction. We observe systematic lateral variations in both traveltime and incoming wave direction, which cannot be explained by changes to the radial seismic velocity profile of the outer core. Moreover, we find a correlation between incoming wave direction and traveltime anomaly, suggesting that observed traveltime anomalies may be caused, at least in part, by changes to the wave path and not solely by perturbations in outer core velocity. Modelling of 1-D ray and 3-D wave propagation in global 3-D tomographic models of mantle velocity anomalies match the trend of the observed traveltime anomalies. Overall, we demonstrate that observed SmKS traveltime anomalies may have a significant contribution from 3-D mantle structure, and not solely from outer core structure.

58 GEOSCIENCES↗
DOE OSTI2021

High-Quality Revision of the Israeli Seismic Bulletin

Seismic bulletins, with trustworthy phase picks, origin times, and source locations are key for regional seismic studies, such as travel-time (TT) tomography, attenuation tomography, and anisotropy studies. To lay the groundwork for such studies in Israel, we revised the seismic bulletin of Israel and the surrounding area and obtained a trustworthy TT data set. From the earthquake and explosion bulletins of the Geophysical Institute of Israel, we compiled a starting data set of about 123,000 earthquakes and explosions that occurred during the past 40 yr. After screening out the poorly recorded events, we were left with a data set of ~38,000 well-recorded events. We then revised the remaining data set in two consecutive steps. In the first, we reviewed and updated station metadata, including changes in station metadata parameters over time. In the second step, we jointly relocated a list of selected seismic events, using the Bayesian hierarchical location software package (BayesLoc) of Myers et al. (2007) that performs joint relocation of multiple events. We observed striking dissimilarities between the spatial distributions of the newly relocated catalog and the initial locations. Although the depth distribution of the starting catalog is trimodal with peaks at 0, 5, and 10 km, the distribution in this study is unimodal, with a broad peak between 7.5 and 12.5 km. By differencing the observed arrival times and the origin times obtained through relocation with BayesLoc, we obtained a revised TT database that consists of 261,336 Pg, 132,876 Pn, 114,816 Sg, and 60,394 Sn arrivals, from a set of 30,458 jointly relocated seismic sources. In this work, we compared prerevision and postrevision TTs as a function of epicentral distance and concluded that the revised data set contains far fewer outliers and inconsistencies than the original data set. The revised TT data set may be used for seismic studies, such as TT tomography, attenuation tomography, and anisotropy studies.

58 GEOSCIENCES↗
DOE OSTI2024

Determining Stress Orientation in Rock Valley, Nevada, Using Ambient Seismic Noise

The stress field and the mechanical properties of rocks are important to consider for nuclear explosion monitoring due to their effect on seismic wave radiation from earthquakes and explosions. At the Rock Valley Direct Comparison site, the regional orientation of the maximum horizontal compressive stress (SH max ) is well constrained, but it is unknown whether there are local heterogeneities. Here, I show that stress-induced anisotropy in nonlinear elasticity can be used to estimate the orientation of SH max . Rocks have compliant internal contacts, such as fractures and mineral grain boundaries, that respond to applied strains more strongly than individual mineral crystals. This strain response is asymmetric between compressive and dilatational strains and is affected by anisotropy in the ambient stress field. Traditional seismic velocity measurements are of strain-averaged velocity, which is less sensitive than nonlinear elasticity to the behavior of compliant internal contacts and therefore the stress field and fracture behavior. My results show that the orientation of SH max measured using ambient seismic noise aligns with regional estimates, with some heterogeneity that can also be explained by the limitations of the seismic array. I demonstrate the potential of a passive technique for monitoring the stress field in places that currently lack local measurements.

58 GEOSCIENCES↗
NASA NTRS1982

Upper mantle anisotropy - Evidence from free oscillations

There is a growing body of evidence that much of the upper mantle may be anisotropic to the propagation of seismic waves. The present investigation shows that a small amount of anisotropy completely changes the nature of the surface wave and normal mode problem. In particular, the apparent lack of sensitivity of many of the spheroidal modes to the compressional velocity structure is due to the degeneracy in the isotropic case. The normal mode data set appears to be adequate to resolve the five elastic constants of a transversely isotropic upper mantle. Dziewonski and Anderson (1981) have shown that these data can be fitted with anisotropy restricted to the upper 200 km of the mantle.

Anderson, D. L.↗
DOE OSTI2020

Numerical Modeling of Anisotropic Elastic-Wave Sensitivity Propagation for Optimal Design of Time-Lapse Seismic Surveys

Reliable subsurface time-lapse seismic monitoring is crucial for many geophysical applications, such as enhanced geothermal system characterization, geologic carbon utilization and storage, and conventional and unconventional oil/gas reservoir characterization, etc. In this work, we develop an elastic-wave sensitivity propagation method for optimal design of cost-effective time-lapse seismic surveys considering the fact that most of subsurface geologic layers and fractured reservoirs are anisotropic instead of isotropic. For anisotropic media, we define monitoring criteria using qP- and qS-wave sensitivity energies after decomposing qP- and qS-wave components from the total elastic-wave sensitivity wavefield using a hybrid time- and frequency-domain approach. Geophones should therefore be placed at locations with significant qP- and qS-wave sensitivity energies for cost-effective time-lapse seismic monitoring in an anisotropic geology setting. Our numerical modeling results for a modified anisotropic Hess model demonstrate that, compared with the isotropic case, subsurface anisotropy changes the spatial distributions of elastic-wave sensitivity energies. Consequently, it is necessary to consider subsurface anisotropies when designing the spatial distribution of geophones for cost-effective time-lapse seismic monitoring. This finding suggests that it is essential to use our new anisotropic elastic-wave sensitivity modeling method for optimal design of time-lapse seismic surveys to reliably monitor the changes in subsurface reservoirs, fracture zones or target monitoring regions.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

$S$ Hmax orientation in the Alpine region from observations of stress-induced anisotropy of nonlinear elasticity

The orientation of $S$ Hmax is commonly estimated from in situ borehole breakouts and earthquake focal mechanisms. Borehole measurements are expensive, and therefore sparse, and earthquake measurements can only be made in regions with many well-characterized earthquakes. Here, we derive the stress-field orientation using stress-induced anisotropy in nonlinear elasticity. In this method, we measure the strain derivative of velocity as a function of azimuth. We use a natural pump-probe (NPP) approach which consists of measuring elastic wave speed using empirical Green’s functions (probe) at different points of the earth tidal strain cycle (pump). The approach is validated using a larger data set in the Northern Alpine Foreland region where the orientation of maximum horizontal compressive stress is known from borehole breakouts and drilling-induced fractures. The technique resolves NNW-SSW to N-S directed $S$ Hmax which is in good agreement with conventional methods and the recent crustal stress model. We confirm that the NPP method can be applied to dense large-scale seismic arrays. The technique is then applied to the Southern Alps to understand the contemporary stress pattern associated with the ongoing deformation due to counterclockwise rotation of the Adriatic plate with respect to the European plate. Our results explain why the two major faults in Northeastern Italy, the Giudicarie Fault and the Periadriatic Line (Pustertal–Gailtal Fault) are currently inactive, while the currently acting stress field allows faults in Slovenia to deform actively. We have demonstrated that the pump-probe method has the potential to fill in the measurement gap left by conventional approaches, both in terms of regional coverage and in depth.

58 GEOSCIENCES↗