Orbital dynamics control monsoon-driven peat carbon sink during the late Paleozoic Ice Age

This repository provides the data products and fully reproducible analysis scripts for the precession-phase orbital experiments used in this study. The dataset is organized by orbital scenario and precession phase, with case directories such as Emax_Omax_90, Emax_Omean_180, and Emean_Omin_270. Each case contains atmospheric output (folder atm/) and, where applicable, land output (folder lnd/) in netCDF format. The atmospheric files include near-surface temperature (TS), precipitation (PRECC+PRECL), and surface water flux/evaporation (QFLX), while the land files provide grid-cell area (area) used for land-only and land–monsoon area-weighted diagnostics. Seasonal-split files (e.g., *.cam.h0.05-09.nc and *.cam.h0.11-03.nc) are included to support monsoon-region definitions based on seasonal precipitation contrasts. In addition, the repository includes MATLAB and Python scripts that reproduce all key diagnostics and figures, including (i) regional annual-mean TS/P/E calculations for the study area, (ii) land–monsoon masking and area-weighted climatologies, (iii) a multi-layer peat-column burial model driven by orbital (eccentricity–obliquity) scenarios with interpolation in (E, O) space to estimate cumulative burial over 1.2 Myr (68 precession cycles), and (iv) Sobol-based sensitivity analyses. The file sobol_param_50_for_monsoon.npy provides a representative subset of Sobol parameter vectors used for robustness tests and cross-scenario burial comparisons.