If you are interested in plotting the results of the inversions or would like to use the models for your own research, feel free to download the jupyter notebook or python script here to get started. The inversion models are saved as HDF5 files and can be read into python using the h5py package. Details about the inversion parameters can be found in the config files (YAML files). Additional information, like the station locations and cross-correlations that were used for the inversion, are in CSV files. The following code shows how to read in and access the inversion models. These can be obtained by downloading the files for an inversion here.
import h5py file_path = "path to a HDF5 model file here" file_h5 = h5py.File(file_path,'r') model = file_h5['model'] coordinates = file_h5['coordinates'] # check other keys print(file_h5.keys())
Within the downloaded inversion folder you will also find various different plots in addition to the noise source maps. All of these can be reproduced and adapted to your needs with the following juyter notebook. The following python packages are needed:
The same code is also available as a python script. Note that only the initial model will be plotted and saved in both, the jupyter notebook and the script. The iteration number can be changed within the code. Download Python Script
The full inversion code that is used to perform the inversions can be found on github: https://github.com/jigel/noisi_inv. It includes a Jupyter notebook (and python script) to setup the inversion configuration file and get you started. In theory, this allows you to reproduce the daily inversions as the corresponding configuration file is included in the downloadable output folder. However, this requires an AxiSEM wavefield with a dominant period of 5 seconds and HPC facilities to run it on multiple cores. Pre-computed wavefields can be downloaded from http://ds.iris.edu/ds/products/syngine/ or you can contact us and we can provide you with the wavefield we use. Note that before using a configuration file from the daily inversions the "download_data_date" needs to be changed from "yesterday" to the date you wish to run the inversion for and multiple paths need to be changed; most importantly the "wavefield_type" to "instaseis" and the "wavefield_path" to your local AxiSEM wavefield.
We run two daily inversions: one for a global station distribution and one for stations surrounding the North Atlantic. The stations in those lists are checked for available data every day. Usually up to 70% of the stations have data available. You can download the CSV files with the network, station, latitude, longitude here:
