Numerical simulations in full space heterogeneous media

This dataset contains the output from numerical simulations of both the wave equation and the radiative transfer equation, performed in two heterogeneous media. The models considered are uniformly heterogeneous, with heterogeneities correlated with a von Kármán autocorrelation function (ACF) and a correlation length of 320 meters. The dataset includes synthetic from two simulation scenarios that differ in the amplitude of velocity fluctuations: one with 10% fluctuations and the other with 25%. SIMULATION OF THE WAVE EQUATION The simulations were performed using the SEM46 spectral-element code (SEM46 documentation), with fourth-order interpolation on cubic elements of 200 m edge length. The source is isotropic, placed at the center of the medium, modeled as a Ricker wavelet centered at 3 Hz. Scenario with 10% fluctuation: files labeled with suffix 10pc_full, scenario with 25% fluctuation: files labeled with suffix 25pc_full Receiver Layout and Sampling 10pc_full files: Contain 38034 seismograms Distributed on concentric spherical shells, centered in the source location, each with 2113 receivers Shells have a minimum radius of 2 km, increasing in 1 km increments 25pc_full files: Contain 29582 seismograms Same spherical distribution as above For both datasets: Displacement: x, y,zStrain components: exx, eyy, ezzRotations: wx, wy, wz Sampling 833Hz Reading the Data The seismograms are stored in binary format and can be read using Python. For example, to load the x-component of displacement from the 10pc_full case: #### import numpy as np data_type = np.dtype('float32').newbyteorder('<') seismo_x = np.fromfile('fsismos_x0000_10pc_full', dtype=data_type) #### Radial Mean Field Additionally, the files: radial_field_10pc_full.npy radial_field_25pc_full.npy Contain the mean radial displacement field recorded for each scenario. These simulations were run with the same media described above, using a total of 90756 receivers distributed across 36 concentric spherical shells (2,521 receivers per shell), starting at a radius of 2 km with 0.5 km increments. Sampling frequency: 500Hz in the 10% case, 833Hz in the 25% case. SIMULATION OF THE RADIATIVE TRANSFER EQUATION In addition to the data obtained from numerical simulations of the wave equation, we include results derived from solving the radiative transfer equation, applied to heterogeneous media with the same properties described previously. The datasets rholSSglob25 and rholSPglob25 represent the S and P-wave energy densities, respectively, for a scenario with 25% fluctuations in material properties. These values are recorded at distances ranging from 0.5 km to 50 km from the seismic source, with increments of 0.5 km. Similarly, rholSSglob10 and rholSPglob10 correspond to the case with 10% fluctuations. The file tglob contains the temporal sampling points at which the energy densities are evaluated.