skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Unraveling the Connection between Subsurface Stress and Geomorphic Features: Dataset
This repository stores data using for the manuscript: Unraveling the Connection between Subsurface Stress and Geomorphic Features The data file used in this study is 'Input_stress_fault_river_BK_091525.csv'. The code used to reproduce all figures in the manuscript is 'Kuhasubpasin_et_al_2025.ipynb' The file contain these following data: Column unit range description lat degree (-90, 90) Latitude lon degree (-180, 180) Longitude azi_R degree (0, 180)* Interpolated azimuth of river network (interpolate without considering river order) azi_r1 degree (0, 180)* Interpolated azimuth of 1'-order river azi_r2 degree (0, 180)* Interpolated azimuth of 2'-order river azi_r3 degree (0, 180)* Interpolated azimuth of 3'-order river azi_r4 degree (0, 180)* Interpolated azimuth of 4'-order river azi_r5 degree (0, 180)* Interpolated azimuth of 5'-order river Drainage_area cell - Drainage area river_order order (1, 7) Majority of the order river in grid cell elev km (0, 5.1375) Elevation TcstDens g/cm^3 (2.7439,2.962) Average crustal density from CRUST 1.0 TcstThk km (5.0731 73.517) Total crustal thickness from CRUST 1.0 crust_type     Crustal type from ECM1 Te km (1,200) Effective elastic thickness MI - (-1,1) Mantle influence index azi_Z degree (0, 180)* Topographic aspect azi_F degree (0, 180)* Interpolated azimuth of faults reg_F - (0, 1) Regime of F azi_SO degree (0, 180)* Interpolated azimuth of feature πœŽπ‘‚ from WSM reg_SO - (0, 1) Regime of πœŽπ‘‚ azi_SO_010 degree (0, 180)* Interpolated azimuth of πœŽπ‘‚ measured between 0-10 km azi_SO_1020 degree (0, 180)* Interpolated azimuth of πœŽπ‘‚ measured between 10-20 km azi_SO_2030 degree (0, 180)* Interpolated azimuth of πœŽπ‘‚ measured between 20-30 km azi_SO_3040 degree (0, 180)* Interpolated azimuth of πœŽπ‘‚ measured between 30-40 km azi_SO_nofm degree (0, 180)* Interpolated azimuth of πœŽπ‘‚ measured from focal mechanism azi_SO_fm degree (0, 180)* Interpolated azimuth of πœŽπ‘‚ measured from other techniques azi_SL degree (0, 180)* Interpolated azimuth of 𝜎𝐿 reg_SL - (0, 1) Regime of 𝜎𝐿 sp1_SL Pa - Magnitude of principal stress 1 for 𝜎𝐿 sp2_SL Pa - Magnitude of principal stress 2 for 𝜎𝐿 azi_SM degree (0, 180)* Interpolated azimuth of feature πœŽπ‘€ reg_SM - (0, 1) Regime of πœŽπ‘€ sp1_SM Pa - Magnitude of principal stress 1 for πœŽπ‘€ sp2_SM Pa - Magnitude of principal stress 2 for πœŽπ‘€ azi_ST degree (0, 180)* Interpolated azimuth of feature πœŽπ‘‡ reg_ST - (0, 1) Regime of πœŽπ‘‡ sp1_ST Pa - Magnitude of principal stress 1 for πœŽπ‘‡ sp2_ST Pa - Magnitude of principal stress 2 for πœŽπ‘‡ azi_SB degree (0, 180)* Interpolated azimuth of feature 𝜎𝐡 delta_SO_F degree (0, 90) Ξ”πœŽπ‘‚βˆ’πΉ delta_SL_F degree (0, 90) Ξ”πœŽπΏβˆ’πΉ delta_SM_F degree (0, 90) Ξ”πœŽπ‘€βˆ’πΉ delta_ST_F degree (0, 90) Ξ”πœŽπ‘‡βˆ’πΉ delta_SB_F degree (0, 90) Ξ”πœŽπ΅βˆ’πΉ delta_SO_R1 degree (0, 90) Ξ”πœŽπ‘‚βˆ’π‘…1 :1' order river delta_SL_R1 degree (0, 90) Ξ”πœŽπΏβˆ’π‘…1 delta_SM_R1 degree (0, 90) Ξ”πœŽπ‘€βˆ’π‘…1 delta_ST_R1 degree (0, 90) Ξ”πœŽπ‘‡βˆ’π‘…1 delta_SB_R1 degree (0, 90) Ξ”πœŽπ΅βˆ’π‘…1 delta_F_R1 degree (0, 90) Ξ”πΉβˆ’π‘…1 delta_SO_R2 degree (0, 90) Ξ”πœŽπ‘‚βˆ’π‘…2 :2' order river delta_SL_R2 degree (0, 90) Ξ”πœŽπΏβˆ’π‘…2 delta_SM_R2 degree (0, 90) Ξ”πœŽπ‘€βˆ’π‘…2 delta_ST_R2 degree (0, 90) Ξ”πœŽπ‘‡βˆ’π‘…2 delta_SB_R2 degree (0, 90) Ξ”πœŽπ΅βˆ’π‘…2 delta_F_R2 degree (0, 90) Ξ”πΉβˆ’π‘…2 delta_SO_R3 degree (0, 90) Ξ”πœŽπ‘‚βˆ’π‘…3 :3' order river delta_SL_R3 degree (0, 90) Ξ”πœŽπΏβˆ’π‘…3 delta_SM_R3 degree (0, 90) Ξ”πœŽπ‘€βˆ’π‘…3 delta_ST_R3 degree (0, 90) Ξ”πœŽπ‘‡βˆ’π‘…3 delta_SB_R3 degree (0, 90) Ξ”πœŽπ΅βˆ’π‘…3 delta_F_R3 degree (0, 90) Ξ”πΉβˆ’π‘…3 delta_SO_R4 degree (0, 90) Ξ”πœŽπ‘‚βˆ’π‘…4 :4' order river delta_SL_R4 degree (0, 90) Ξ”πœŽπΏβˆ’π‘…4 delta_SM_R4 degree (0, 90) Ξ”πœŽπ‘€βˆ’π‘…4 delta_ST_R4 degree (0, 90) Ξ”πœŽπ‘‡βˆ’π‘…4 delta_SB_R4 degree (0, 90) Ξ”πœŽπ΅βˆ’π‘…4 delta_F_R4 degree (0, 90) Ξ”πΉβˆ’π‘…4 delta_SO_R5 degree (0, 90) Ξ”πœŽπ‘‚βˆ’π‘…5 :5' order river delta_SL_R5 degree (0, 90) Ξ”πœŽπΏβˆ’π‘…5 delta_SM_R5 degree (0, 90) Ξ”πœŽπ‘€βˆ’π‘…5 delta_ST_R5 degree (0, 90) Ξ”πœŽπ‘‡βˆ’π‘…5 delta_SB_R5 degree (0, 90) Ξ”πœŽπ΅βˆ’π‘…5 delta_F_R5 degree (0, 90) Ξ”πΉβˆ’π‘…5 delta_SO_R>1 degree (0, 90) Ξ”πœŽπ‘‚βˆ’π‘…>1 :>1' order river delta_SL_R>1 degree (0, 90) Ξ”πœŽπΏβˆ’π‘…>1 delta_SM_R>1 degree (0, 90) Ξ”πœŽπ‘€βˆ’π‘…>1 delta_ST_R>1 degree (0, 90) Ξ”πœŽπ‘‡βˆ’π‘…>1 delta_SB_R>1 degree (0, 90) Ξ”πœŽπ΅βˆ’π‘…>1 delta_F_R>1 degree (0, 90) Ξ”πΉβˆ’π‘…>1 delta_SO_Z degree (0, 90) Ξ”πœŽπ‘‚βˆ’π‘ delta_SL_Z degree (0, 90) Ξ”πœŽπΏβˆ’π‘ delta_SM_Z degree (0, 90) Ξ”πœŽπ‘€βˆ’π‘ delta_ST_Z degree (0, 90) Ξ”πœŽπ‘‡βˆ’π‘ delta_SB_Z degree (0, 90) Ξ”πœŽπ΅βˆ’π‘ delta_F_Z degree (0, 90) Ξ”πΉβˆ’π‘ delta_Z_R1 degree (0, 90) Ξ”π‘βˆ’π‘…1 :1' order river delta_Z_R2 degree (0, 90) Ξ”π‘βˆ’π‘…2 :2' order river delta_Z_R3 degree (0, 90) Ξ”π‘βˆ’π‘…3 :3' order river delta_Z_R4 degree (0, 90) Ξ”π‘βˆ’π‘…4 :4' order river delta_Z_R5 degree (0, 90) Ξ”π‘βˆ’π‘…5 :5' order river delta_Z_R>1 degree (0, 90) Ξ”π‘βˆ’π‘…>1 :>1' order river *The range is not (0,360) because we only consider azimuth not direction  more » « less
Award ID(s):
2012073 1945431
PAR ID:
10641523
Author(s) / Creator(s):
; ;
Publisher / Repository:
Zenodo
Date Published:
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract Marine multichannel and wide‐angle seismic data constrain the distribution of seamounts, sediment cover sequence and crustal structure along a 460 km margin‐parallel transect of the Hikurangi Plateau. Seismic reflection data reveals five seamount up‐to 4.5 km high and 35–75 km wide, with heterogeneous internal velocity structure. Sediment cover decreases south‐to‐north from ∼4.5 km to ∼1–2 km. The Hikurangi Plateau crust (VP5.5–7.5 km/s) is 11 Β± 1 km thick in the south, but thins by 3–4 km further north (∼7–8 km). Gravity models constructed along two seismic lines show the reduction in crustal thickness persists further east, coinciding with a bathymetric scarp. Gravity data suggest the transition in crustal thickness may reflect spatial variability in deformation and lithospheric extension associated with plateau breakup. Variability in the thickness of subducting crust may contribute to differences in megathrust geometry, upper‐plate stress state and high‐rates of contraction and uplift along the southern Hikurangi margin. 
    more » « less
  2. Abstract The extremely oblique Indo‐Burma subduction zone exhibits dextral strike‐slip faulting along the Sagaing, Kabaw, and Churachandpur‐Mao Faults as well as east‐west shortening between the Sagaing Fault and Bengal Basin. Through regional stress analysis, considering areas from central Tibet, around the eastern Himalaya Syntaxis, to Burma, it has been determined that the principal compressive stress directions align with the principal strain rates. The northeast‐southwest oriented compressive stress direction from the western Shan Plateau continues into Burma. Notably, P axes align with the topographic gradients, and T axes are sub‐parallel to the topographic contours in the Shan Plateau region south of 27Β°N. These stress patterns are consistent with a gravitational potential energy induced crustal and mantle flow. The alignment of the fast shear wave with the maximum strain rate and the colinear NW‐SE to E‐W fast direction of the SKS wave and T axis determined from focal mechanisms in the Shan Plateau suggest that the mantle lithosphere deforms in concert with the crust. We suggest crust and mantle flow south of the Red River Fault has resulted in widening of the lithosphere in the Shan Plateau in an east‐west direction. Therefore, the Sagaing Fault has bowed approximately 50–100 km westward if we assume that the Sagaing Fault was originally straight. Our results of regional stress inversion are consistent with late Miocene to present E‐W shortening in the Indo‐Burma subduction zone resulting from the release of gravitational potential energy from the central Tibetan Plateau. 
    more » « less
  3. Abstract Processes controlling the formation of continental whole-lithosphere shear zones are debated, but their existence requires that the lithosphere is mechanically coupled from base to top. We document the formation of a dextral, whole-lithosphere shear zone in the Death Valley region (DVR), southwest United States. Dextral deflections of depth gradients in the lithosphere-asthenosphere boundary and Moho are stacked vertically, defining a 20–50-km-wide, lower lithospheric shear zone with ~60 km of shear. These deflections underlie an upper-crustal fault zone that accrued ~60 km of dextral slip since ca. 8–7 Ma, when we infer that whole-lithosphere shear began. This dextral offset is less than net dextral offset on the upper-crustal fault zone (~90 km, ca. 13–0 Ma) and total upper-crustal extension (~250 km, ca. 16–0 Ma). We show that, before ca. 8–7 Ma, weak middle crust decoupled upper-crustal deformation from deformation in the lower crust and mantle lithosphere. Between 16 and 7 Ma, detachment slip thinned, uplifted, cooled, and thus strengthened the middle crust, which is exposed in metamorphic core complexes collocated with the whole-lithosphere shear zone. Midcrustal strengthening coupled the layered lithosphere vertically and therefore enabled whole-lithosphere dextral shear. Where thick crust exists (as in pre–16 Ma DVR), midcrustal strengthening is probably a necessary condition for whole-lithosphere shear. 
    more » « less
  4. Archived data for the manuscript β€œDifferentiating Between Simultaneous Loss Drivers in Earth’s Outer Radiation Belt: Multi-Dimensional Phase Space Density Analysis”  Staples et al., submitted to Geophysical Research Letters 2022.</p> These files contain multi-mission phase space density measurements including, Van Allen Probes, GOES 13, 14, 15, GPS, MMS, and THEMIS, computed in adiabatic coordinates. All data is from September 2017. For detailed description of the method used in the computation of this data, see section 2.1 of the submitted manuscript. The THEMIS, Van Allen Probe, MMS, and GOES data used in computations is publicly available via http://cdaweb.gsfc.nasa.gov  The GPS data is available via https://www.ngdc.noaa.gov/stp/space-weather/satellite-data/satellite-systems/gps/</p> FILES:</p> 'psd_intp_T89_20170901-20170930_allsc.cdf'</p> 'psd_intp_T89_20170901-20170930_rbsp-b.cdf'</p> DATA OWNER: Adam Kellerman DATA PREPERATION: Frances Staples</p> CONTACT: Adam Kellerman: akellerman@epss.ucla.edu Frances Staples: frances.staples@ucl.ac.uk</p> FS was supported by NASA grants 80NSSC20K1402 and NSF grant 2149782. ACK acknowledges support from NASA grants 80NSSC20K1402 and 80NSSC20K1281, and NSF grant 2149782. 
    more » « less
  5. {"Abstract":["This is an extracted data product for radar bed reflectivity from Whillans Ice Plain, West Antarctica. The original data are hosted by the Center for Remote Sensing of Ice Sheets (CReSIS; see associated citation below). The files here can be recalculate and are meant to be used within a set of computational notebooks here:https://doi.org/10.5281/zenodo.10859135\n\nThere are two csv files included here, each structured as a Pandas dataframe. You can load them in Python like:df = pd.read_csv('./Picked_Bed_Power.csv')\n\nThe first file, 'Picked_Bed_Power.csv' is the raw, uncorrected power from the radar image at the bed pick provided by CReSIS. There are also other useful variables for georeferencing, flight attributes, etc.\n\nThe second file, 'Processed_Reflectivity.csv' is processed from the first file. Processing includes: 1) a spreading correction; 2) an attenuation correction; and, 3) a power adjustment flight days based on compared power at crossover points. This file also has identifiers for regions including "grounded ice", "ungrounded ice", and "subglacial lakes"."]} 
    more » « less