3D continental thermal model for the western United States and East Asia and the associated thermal model Python software package

This repository contains the netCDF4 datasets of the 3D continental thermal model for the western United States and East Asia, and a Python-based thermal modeling function for "Lee et al., Building a thermal model of the continental lithosphere: Implications for solid Earth processes and beyond" that is in review. ThermalModel_WUS.nc includes the 3D crustal thermal model for the western United States and all the temperature proxies used to constrain the thermal model, including surface heat flow (Williams et al., 2007; Mordensky and DeAngelo, 2023), seismogenic thickness derived from relocated earthquake catalog (EQ_catalog_WUS.csv), seismogenic thickness uncertainty, crustal thickness (Buehler and Shearer, 2017), and Moho temperature datasets (Schutt et al., 2018). Furthermore, it also includes the modeling results of all the thermal parameters (e.g., crustal, mantle, and surface heat flow, thermal conductivity), temperature at 20 km depth, crustal-scale geothermal gradient, and upper crustal geothermal gradient. ThermalModel_EastAsia.nc includes the lithospheric thermal model for East Asia and all the temperature proxies used to constrain the thermal structure, including surface heat flow (Wu et al., 2013; Sun et al., 2022), Curie depth (Li et al., 2017), and lithospheric thickness (An and Shi, 2006). It contains the same type of results as the ThermalModel_WUS.nc. EQ_catalog_WUS.csv contains the merged earthquake catalog in California and Nevada, USA. This catalog is merged from datasets from northern (Waldhauser and Schaff, 2008; Waldhauser, 2009) and southern California (Hauksson et al., 2012), and Nevada (Trugman, 2023, 2024). Earthquake data range from 1980 to 2024 and include events that has magnitude larger than one. Hypocentral depths that are deeper than the crustal thickness estimates from Buehler and Shearer (2017) were removed from the catalog. The hypocenter datum is referenced to surface elevation. Please refer to the supplemental material of Lee et al. et for the details of merging earthquake catalog. ThermalModelFunction.zip contains the Python function 'ContinentalThermalModel.py', configuration text file 'ContinentalThermalModel_config.txt', and model documentation and instruction Jupyter source file 'ThermalModel_Documentation.ipynb'. Dependent Python functions are also included in this folder, such as 'D95.py', 'Get_time.py', 'make_output_config.py', 'Model_extent.py', 'Parse_config.py', 'Read_config.py', 'resample_tif.py', and 'Run_ThermalModel.py'.Input temperature proxies for the western United States and East Asia preferred model are also located in the folder 'Temperature_proxies_WUS' and 'Temperature_proxies_East_Asia', including earthquake catalog 'EQ_catalog_WUS.csv', surface heat flow (Williams et al., 2007; Mordensky and DeAngelo, 2023) 'heat_flow.tif', crustal thickness (Buehler and Shearer, 2017) 'Moho_depth.tif', and Moho temperature 'Moho_T.tif' and its uncertainty 'Moho_T_Uncer.tif' datasets (Schutt et al., 2018), surface heat flow (Asia_heatflow.tif; Wu et al., 2013; Sun et al., 2022), Curie depth (Curie_Li_et_al_2017.tif; Li et al., 2017), and lithospheric thickness (LAB_An_Shi_2006.tif; An and Shi 2006).