Archive of the iSALE simulation datasets for Jackson et al. 2022, "Impact generation of holes in the early lunar crust I: scaling relations"

Archive of the iSALE simulation datasets used in Jackson et al. 2022, "Impact generation of holes in the early lunar crust I: scaling relations". This archive contains all necessary information to re-run any of the simulations in the dataset along with plots that can be used to check the outcomes of each simulation. In addition, it also contains a full compilation of all of the derived data used in Jackson et al. 2022 as well as the analysis scripts used to extract that data from the simulations and the scripts that use the derived data to generate the figures in the manuscript. The archive is provided in a standard compressed (.zip) format. <strong>Detailed description of the contents</strong> Note that this description can also be found in the ReadMe file provided with the archive <strong>1) Simulation folders</strong> The folder '<em>simsmain</em>' contains files for the iSALE2D simulations that comprise the primary simulation database. For each simulation we provide: asteroid.inp - the initialisation file for the simulation that would allow it to be re-run [].dat - a text file containing data derived from the final simulation output. These values from all simulations are compiled in alldata.csv. A description of the columns is provided below for alldata.csv []-Crthickprof.pdf - a figure showing the radial crust thickness profile at the end of the simulation. Light pink shows the cell-by-cell profile, solid red a moving average. []-Damprof.pdf - a figure showing the radial profile of the mean damage factor (averaged over the depth of the crust at each radial cell) at the end of the simulation. Light grey shows the cell-by-cell profile, solid black a moving average. []-TmpDam.dpdf - a figure showing the temperature and damage structure of the entire high-resolution zone at the end of the simulation at 1:1 aspect ratio (i.e. square cells are square in the image). The boundary of the crust is outlined in black. []-TmpDam-zoom.pdf - as above, but zoomed in on the central region of interest in the simulation. This image is stretched to highlight the region of interest (i.e. square cells are not square in the image). []-temporal.dat - a text file containing data derived from the full simulation time series. These values from all simulations are compiled in alldata_temporal.csv. A description of the columns is provided below for alldata_temporal.csv []-crdepthtime.png - a figure showing the change in the maximum depression of the upper surface of the simulation over time. This records the excavation and collapse of the initial transient cavity. []-holeradtime.png - a figure recording the radius of the hole (as measured by the locations at which the crust thickness passes 90 and 95 percent of the initial value) over time. A final plateau is indicative of the simulation having settled. []-holevoltime.png - a figure recording the volume of the hole within the 95% crust thickness radius over time. []-trcavtime.png - a figure recording the volume of the surface depression over time. This can be quite sensitive to small fluctuations at the edges of the surface depression and does not always plateau. Within <em>simsmain</em> each sub-folder is named in the format dW_vX_cY_gZ, where d is the diameter of the impactor (in km), v is the impact velocity (in m/s), c is the crust thickness (in km) and g is the grid resolution of the high-resolution zone (in km). W, X, Y and Z are all in exponential notation with two significant figures (i.e. N.NE+NN). The folder '<em>simsres</em>' contains files for the additional iSALE2D simulations run for resolution testing purposes. These files have the same formats and naming conventions as <em>simsmain</em>. The folder '<em>material</em>' contains the file material.inp, the initialisation file for the material properties which is common to all of the simulations. <strong>2) Analysis scripts</strong> The folder '<em>analysisscripts</em>' contains the Python scripts that were used to generate the files in <em>simsmain</em> and <em>simsres</em>. There are two scripts, iSALEanalysis_batch.py which generates the end-of-simulation plots and data, and iSALEanalysis_batch_temporal.py which generates the time series plots and data. The temporal script typically takes around 400 times longer to run than the end-of-simulation script since the former must examine all (usually 800) timesteps, whereas the end-of-simulation script examines only the first and last timestep. <strong>3) Figures</strong> The folder 'figures' contains CSV files summarising all of the data derived from the iSALE2D simulations along with plotting scripts that use this derived data to generate the plots in Jackson et al. 2022. -----------<br> alldata.csv is a compilation of the data in the [].dat files for all 252 simulations in the primary simulation database. There are 14 columns which in order are: dimp - impactor diameter in km vimp - impact velocity in km/s crd - initial crust thickness in km tdat - simulation time (in seconds) at which the data was taken gres - grid resolution in the high-resolution zone in km rhires - radius of the high-resolution zone in km r_thincr0_95- radial distance (in km) at which the crust thickness exceeds 95% of the initial value r_thincr0_9 - radial distance (in km) at which the crust thickness exceeds 90% of the initial value r_magcr - radius (in km) within which there is direct exposure of magma to the surface (i.e. crust thickness is zero) r_frac09 - radial distance (in km) within which the vertically averaged damage factor in the crust exceeds 0.9 r_frac08 - radial distance (in km) within which the vertically averaged damage factor in the crust exceeds 0.8 r_frac05 - radial distance (in km) within which the vertically averaged damage factor in the crust exceeds 0.5 holevolume - volume (in km^3) of crust excavated from within r_thincr0_95 comment2 - impact regime classification, 'partial'= partial penetration, 'complete' = complete penetration, 'cfd' = cratering with full-depth fracturing, 'classical' = classical crater. -----------<br> alldata_temporal.csv is a compilation of the data in the []-temporal.dat files for all 252 simulations in the primary simulation database. There are 20 columns which in order are: dimp - impactor diameter in km vimp - impact velocity in m/s (Note the difference from alldata.csv) crd - initial crust thickness in km tdat - simulation time (in seconds) at which the data was taken gres - grid resolution in the high-resolution zone in km rhires - radius of the high-resolution zone in km r_thincr0_95- radial distance (in km) at which the crust thickness exceeds 95% of the initial value r_thincr0_9 - radial distance (in km) at which the crust thickness exceeds 90% of the initial value r_magcr - radius (in km) within which there is direct exposure of magma to the surface (i.e. crust thickness is zero) r_frac09 - radial distance (in km) within which the vertically averaged damage factor in the crust exceeds 0.9 r_frac08 - radial distance (in km) within which the vertically averaged damage factor in the crust exceeds 0.8 r_frac05 - radial distance (in km) within which the vertically averaged damage factor in the crust exceeds 0.5 holevolume - volume (in km^3) of crust excavated from within r_thincr0_95 mcrd - maximum depth of the transient cavity in km t_mcrd - simulation time (in seconds) at which transient cavity reaches maximum depth trcav_mcrd - volume (in km^3) of the transient cavity at the time of maximum depth mtrcav - maximum volume of the transient cavity t_mtrcav - simulation time (in seconds) at which transient cavity reaches maximum volume crd_mtrcav - depth of the transient cavity (in km) at the time of maximum volume comment2 - impact regime classification, 'partial'= partial penetration, 'complete' = complete penetration, 'cfd' = cratering with full-depth fracturing, 'classical' = classical crater. -----------<br> restest-rthin.csv and restest-frac.csv are in the same format as alldata.csv and contain the data for the sets of simulations used to test the resolution dependence of r_thin and r_frac respectively. ----------- There are 7 Python plotting scripts: plotter_main.py - the main plotting script, produces plots that don't fall under another heading plotter_KE.py - produces plots as a function of impact kinetic energy plotter_temporal.py - produces plots using time series data from alldata_temporal plotter_cox.py - produces plots comparing our results to those of Bray et al. 2014, Cox et al. 2008, Cox &amp; Bauer 2015 and Miljkovic et al. 2015. Note that we do not include the CSV files containing the data from those papers since it is not ours to distribute, however they can be easily extracted from the relevant publications. plotter_res.py - produces resolution comparison plots plotter_pi.py - produces Pi-scaling plots kplot-sims.py - produces a plot of the 'k' factor (hole size/impactor mass) for comparison with Perera et al. 2018 All of the plots and fitting summaries produced by the plotting scripts are also included本存档为Jackson等人2022年发表于论文《早期月球地壳中孔洞的撞击成因I：尺度关系》的iSALE模拟数据集存档。该存档包含复现数据集中所有模拟所需的全部信息，以及可用于校验各模拟结果的绘图文件。此外，存档还收录了Jackson等人2022年研究中使用的全部衍生数据汇编，以及从模拟结果中提取该衍生数据的分析脚本，和使用衍生数据生成论文配图的脚本。本存档采用标准压缩（.zip）格式。**内容详细说明** 本说明也可在存档附带的ReadMe文件中查阅。 **1) 模拟文件夹** 文件夹`simsmain`包含构成主模拟数据库的iSALE2D模拟相关文件。每个模拟均提供以下文件： - `asteroid.inp`：模拟初始化文件，可用于复现该模拟。 - `[].dat`：包含从模拟最终输出结果提取的文本数据。所有模拟的该类数据均汇编至`alldata.csv`中，下文将对`alldata.csv`的列项进行说明。 - `[]-Crthickprof.pdf`：展示模拟结束时径向地壳厚度剖面的图像。浅粉色表示逐单元格剖面，实线红色为移动平均曲线。 - `[]-Damprof.pdf`：展示模拟结束时平均损伤因子的径向剖面（按每个径向单元格的地壳深度平均）。浅灰色表示逐单元格剖面，实线黑色为移动平均曲线。 - `[]-TmpDam.dpdf`：展示模拟结束时整个高分辨率区域的温度与损伤结构的图像，采用1:1纵横比（即图像中单元格为正方形）。地壳边界以黑色线条勾勒。 - `[]-TmpDam-zoom.pdf`：与上述图像内容一致，但对模拟中核心感兴趣区域进行了放大。该图像经过拉伸以突出目标区域（即图像中单元格非正方形）。 - `[]-temporal.dat`：包含从完整模拟时间序列提取的文本数据。所有模拟的该类数据均汇编至`alldata_temporal.csv`中，下文将对`alldata_temporal.csv`的列项进行说明。 - `[]-crdepthtime.png`：展示模拟过程中地壳上表面最大凹陷随时间变化的图像，记录了初始瞬态空腔的挖取与坍塌过程。 - `[]-holeradtime.png`：记录孔洞半径随时间变化的图像（孔洞半径以地壳厚度恢复至初始值90%和95%的位置进行测量）。最终趋于平稳的曲线表明模拟已进入稳定阶段。 - `[]-holevoltime.png`：记录以95%地壳厚度半径为边界的孔洞体积随时间变化的图像。 - `[]-trcavtime.png`：记录地表凹陷体积随时间变化的图像。该指标对地表凹陷边缘的微小波动较为敏感，不一定会趋于平稳。 `simsmain`内的每个子文件夹均采用`dW_vX_cY_gZ`的命名格式，其中：d为撞击体直径（单位：km），v为撞击速度（单位：m/s），c为地壳厚度（单位：km），g为高分辨率区域的网格分辨率（单位：km）。W、X、Y、Z均采用两位有效数字的指数计数法（即`N.NE+NN`格式）。 文件夹`simsres`包含用于分辨率测试的额外iSALE2D模拟文件，其文件格式与命名约定与`simsmain`一致。 文件夹`material`包含`material.inp`文件，即所有模拟通用的材料属性初始化文件。 **2) 分析脚本** 文件夹`analysisscripts`包含用于生成`simsmain`和`simsres`中文件的Python脚本。包含两个脚本： - `iSALEanalysis_batch.py`：用于生成模拟结束后的绘图与数据文件。 - `iSALEanalysis_batch_temporal.py`：用于生成时间序列绘图与数据文件。时间序列脚本的运行时长通常约为结束后脚本的400倍，因为前者需要处理所有（通常为800个）时间步，而后者仅需处理首个和最后一个时间步。 **3) 配图文件** 文件夹`figures`包含汇总所有iSALE2D模拟衍生数据的CSV文件，以及使用该衍生数据生成Jackson等人2022年论文配图的绘图脚本。 --- `alldata.csv`为主模拟数据库中全部252个模拟的`[].dat`文件数据汇编，共包含14个列项，按顺序分别为： 1. `dimp`：撞击体直径（单位：km） 2. `vimp`：撞击速度（单位：km/s） 3. `crd`：初始地壳厚度（单位：km） 4. `tdat`：数据采集时刻的模拟时间（单位：s） 5. `gres`：高分辨率区域的网格分辨率（单位：km） 6. `rhires`：高分辨率区域的半径（单位：km） 7. `r_thincr0_95`：地壳厚度恢复至初始值95%以上的径向距离（单位：km） 8. `r_thincr0_9`：地壳厚度恢复至初始值90%以上的径向距离（单位：km） 9. `r_magcr`：岩浆直接暴露于地表的径向范围半径（单位：km，即地壳厚度为0的区域） 10. `r_frac09`：地壳垂直平均损伤因子超过0.9的径向范围半径（单位：km） 11. `r_frac08`：地壳垂直平均损伤因子超过0.8的径向范围半径（单位：km） 12. `r_frac05`：地壳垂直平均损伤因子超过0.5的径向范围半径（单位：km） 13. `holevolume`：以`r_thincr0_95`为边界的挖取地壳体积（单位：km³） 14. `comment2`：撞击体制分类别，"partial"为部分穿透，"complete"为完全穿透，"cfd"为全深度断裂成坑，"classical"为经典撞击坑。 --- `alldata_temporal.csv`为主模拟数据库中全部252个模拟的`[]-temporal.dat`文件数据汇编，共包含20个列项，按顺序分别为： 1. `dimp`：撞击体直径（单位：km） 2. `vimp`：撞击速度（单位：m/s，注意与`alldata.csv`的单位差异） 3. `crd`：初始地壳厚度（单位：km） 4. `tdat`：数据采集时刻的模拟时间（单位：s） 5. `gres`：高分辨率区域的网格分辨率（单位：km） 6. `rhires`：高分辨率区域的半径（单位：km） 7. `r_thincr0_95`：地壳厚度恢复至初始值95%以上的径向距离（单位：km） 8. `r_thincr0_9`：地壳厚度恢复至初始值90%以上的径向距离（单位：km） 9. `r_magcr`：岩浆直接暴露于地表的径向范围半径（单位：km，即地壳厚度为0的区域） 10. `r_frac09`：地壳垂直平均损伤因子超过0.9的径向范围半径（单位：km） 11. `r_frac08`：地壳垂直平均损伤因子超过0.8的径向范围半径（单位：km） 12. `r_frac05`：地壳垂直平均损伤因子超过0.5的径向范围半径（单位：km） 13. `holevolume`：以`r_thincr0_95`为边界的挖取地壳体积（单位：km³） 14. `mcrd`：瞬态空腔的最大深度（单位：km） 15. `t_mcrd`：瞬态空腔达到最大深度的模拟时刻（单位：s） 16. `trcav_mcrd`：瞬态空腔达到最大深度时的体积（单位：km³） 17. `mtrcav`：瞬态空腔的最大体积（单位：km³） 18. `t_mtrcav`：瞬态空腔达到最大体积的模拟时刻（单位：s） 19. `crd_mtrcav`：瞬态空腔达到最大体积时的深度（单位：km） 20. `comment2`：撞击体制分类别，"partial"为部分穿透，"complete"为完全穿透，"cfd"为全深度断裂成坑，"classical"为经典撞击坑。 --- `restest-rthin.csv`与`restest-frac.csv`格式与`alldata.csv`一致，分别包含用于测试`r_thin`与`r_frac`分辨率依赖性的模拟数据集。 --- 本存档共包含7个Python绘图脚本： 1. `plotter_main.py`：主绘图脚本，生成未归入其他类别的绘图结果 2. `plotter_KE.py`：生成以撞击动能为自变量的绘图结果 3. `plotter_temporal.py`：使用`alldata_temporal`中的时间序列数据生成绘图结果 4. `plotter_cox.py`：生成本研究结果与Bray等人2014年、Cox等人2008年、Cox & Bauer 2015年以及Miljkovic等人2015年研究结果的对比绘图。需注意，本存档未包含上述研究的CSV数据文件，因该数据不属于本研究的分发范围，但可从相关公开出版物中直接获取。 5. `plotter_res.py`：生成分辨率对比绘图 6. `plotter_pi.py`：生成Pi尺度分析绘图 7. `kplot-sims.py`：生成"k"因子（孔洞尺寸/撞击体质量）的绘图，用于与Perera等人2018年的研究结果进行对比。 所有绘图脚本生成的配图与拟合汇总结果均已包含在存档中。