Data and code repository for "Subsurface fluvial sediments beneath InSight on Mars from geophysical constraints"

This repository contains the data for subsurface structure inversion, codes for data processing, and inverted structure models. The codes have been tested for Python version 3.8.10 and MATLAB version R2023b with WINDOWS 11 system (22631.3880). The repository is organized in the following: data/data/CVC/ WS1_ver_norm.dat: Vertical ground compliance extracted from convective vortex events, normalized at a wind speed of 4 m/s ([1]). WS2_ver_norm.dat: Similar to 'WS1_ver_norm.dat', normalized at a wind speed of 8 m/s ([1]).data/HVSR/ S1222a_HVSR.dat: Original horizontal-to-vertical spectral ratio (HVSR) under diffused wavefield conditions for event S1222a. S1222a_HVSR_with2.4_cutLSA.dat: Similar to 'S1222a_HVSR.dat', with the Load Shunt Assembly contamination at 9.0 Hz removed. S1222a_HVSR_without2.4_cutLSA.dat: Similar to 'S1222a_HVSR_with2.4_cutLSA.dat', with the HVSR trough caused by the 2.4 Hz signal removed. HVSR_frequency_sample.dat: Frequency sample points of HVSR used for inversion.data/ICHV/ ICHV_S0986c.dat: Horizontal-to-vertical ratio of infrasound compliance, for S0986c impact event.data/S1222a/ Data_format: Data format for sub-directories in 'data/S1222a/'. cal_DFAHV.m: To calculate HVSR in 'data/S1222a/DFAHV/'. cal_PSD.m: To calculate power spectral density (PSD) in 'data/S1222a/PSD/'.data/S1222a/Coher/ Data for cross-component coherence.data/S1222a/DFAHV/ Data for HVSR under diffused wavefield conditions (DFA).data/S1222a/lib/ Related subroutines.data/S1222a/PSD/ Data for power spectral density.data/S1222a/waveform/ Data for waveforms of event S1222a. model/ Data_format: Data format for sub-directories in 'model/'. Model_WITH2.4_HVSR_CVC_ICHV_0.50_0.25_0.25: The adopted representative model for inversion case 'WITH2.4_HVSR_CVC_ICHV_0.50_0.25_0.25'. Model_WITHOUT2.4_HVSR_CVC_ICHV_0.50_0.25_0.25: The adopted representative model for inversion case 'WITHOUT2.4_HVSR_CVC_ICHV_0.50_0.25_0.25'.model/WITH2.4_HVSR_CVC_ICHV_0.50_0.25_0.25/ Models from the inversion case WITH the 2.4 Hz signal in the HVSR curve, with weight coefficients HVSR:CVC:ICHV = 0.50:0.25:0.25.model/WITH2.4_HVSR_CVC_ICHV_0.70_0.15_0.15/ Similar to 'WITH2.4_HVSR_CVC_ICHV_0.50_0.25_0.25/', with weight coefficients HVSR:CVC:ICHV = 0.70:0.15:0.15.model/WITH2.4_HVSR_CVC_ICHV_0.34_0.33_0.33/ Similar to 'WITH2.4_HVSR_CVC_ICHV_0.50_0.25_0.25/', with weight coefficients HVSR:CVC:ICHV = 0.34:0.33:0.33.model/WITH2.4_HVSR_CVC_ICHV_0.10_0.45_0.45/ Similar to 'WITH2.4_HVSR_CVC_ICHV_0.50_0.25_0.25/', with weight coefficients HVSR:CVC:ICHV = 0.10:0.45:0.45. model/WITHOUT2.4_HVSR_CVC_ICHV_0.50_0.25_0.25/ Similar to 'WITH2.4_HVSR_CVC_ICHV_0.50_0.25_0.25/', WITHOUT the 2.4 Hz signal in the HVSR curve. lithology/ fig3.py: To plot velocity fit to inverted middle layers by rock physics models ('Figure 3'). figE1.py: To plot velocity fit to inverted regolith layer by rock physics models ('Extended Figure 1'). figS15.py: To plot velocity fit to inverted middle layers by rock physics models ('Supplementary Figure 15').lithology/lib/ Related subroutines. References[1] Onodera, K. et al. Description of Martian convective vortices observed by InSight and implications for vertical vortex structure and subsurface physical properties. J. Geophys. Res. Planets 128, e2023JE007896 (2023).