A compact Earth system model.
Download a release. Read The Fine Manual.
OSCAR v3 is developed in Python 3.9 and run preferentially through IPython. It makes heavy use of the xarray
package (v0.20.1), and netCDF4
for saving data (v1.5.7). It also relies on other common scientific packages: numpy
(v1.23.3) and scipy
(v1.9.3). Although the latest versions of the packages used to run OSCAR v3 are given, older versions are likely to work.
The source code is provided firstly for transparency, and only secondly for dissemination. This means that it is provided as is, and no support of any kind is guaranteed (but feel free to ask). Feedbacks and contributions are welcome.
get_final=True
to request a second output that is the final state of the called simulation.mod_process
that could lead to NaN
values under some uncommon combination of regional aggregations and land carbon cycle parameterizations.scipy.odr
), to be used for future module calibrations.pCO2
functional form.D_CO2
. The latter can be turned off using adapt_nt=False
when calling the model.x_OH_NOX
, x_OH_CO
and x_OH_VOC
when mod_Foh_fct == log
. The sensitivity of the OH sink to ozone precursors would be of the opposite sign in half the parameterizations.p_shift
in the bookkeeping module. It would previously overwrite some minor transitions with NaN
, leading to a slight underestimation of land-use change emissions.fct_process_alt
in the next version.Houghton_2017
regions that had Colombia and Guinea swapped. This had minimal impact on the land C cycle.This version is exactly the one used in the IPCC AR6. A number of functionalities were added, mostly for convenience.
fct_process
was renamed mod_process
for consistency, which may affect existing run scripts.fct_process_alt
to alter the structure of the bookkeeping module. The default structure is now defined over only one bio_land
dimension instead of two (bio_from
and bio_to
) to reduce the model's memory demand. Mathematically, this is equivalent to the full structure, although extremely small differences can be seen because of the numerical solving. This can be tested using the provided check_LUC_structure
script.adjust_config
function in fct_genMC
.e_ohu
parameter representing the ocean heat uptake efficacy in the climate module (see e.g. here). For consistency, this parameter is set to one by default for all configurations.D_Flasc
, d_CO2
, CFF
, d_OHC
) in mod_process
to help diagnose the model ex-post.get_IPCC_AR6_parameters
script, to generate a set of prior parameters identical to those used in the IPCC AR6.D_ewet
(wetlands areal emissions), to account for the new flexible structure of the land carbon cycle. This reduces the speed at which CH4 emissions respond to a change in net primary productivity.The physical equations and parameter values of this version are exactly the same as in v2.4. A few notable changes are worth mentioning here.
This is a complete revamp of OSCAR, for which the code was rewritten from scratch and moved to Python 3. OSCAR v3 relies heavily on the xarray
package (and therefore netCDF saving format), to better structure the input and output data, to more easily deal with the many dimensions of internal data, and to parallelize simulations so that the run time of a typical Monte Carlo ensemble has been significantly reduced. It also uses its own object classes (Model
and Process
) that make it easy to change, extend or tune the model's structure and/or equations. Because of this complete overhaul, manipulation of the model has to be learned from scratch, however.
This version is the last update of OSCAR v2. It is meant to be as close to OSCAR v3.0 as possible.
mod_biomeV3
option to force biome aggregation to that of v3.0.mod_EHWPspeed
option that introduces further variations in the decay time of harvested wood products. In addition to the normal
configuration, the fast
one scales the value of tau_HWP
so that 20% of the pool remains after 80% of the initial decay time (a rescale by ~0.5), and the slow
one scales the value of tau_HWP
so that 30% of the pool remains after 150% of the initial decay time (a rescale by ~1.25).mod_EHWPfct
option, so that only exponential decay of harvested wood products is now possible. This slightly increases CO2 emissions from land-use change in a Monte Carlo run.p_wet
on the mod_LSNKcover
option. It is now the average of all possible configurations, which has very little effect on the simulated wetlands CH4 emissions.p_HWP1_BB
quantifying how much of the harvested wood products in pool 1 are actually burnt in the open, and thus accounted for in non-CO2 anthropogenic biomass burning emissions. It is set to 0.5
to roughly match present-day estimates.GISS-E2-R-TOMAS
to the mod_O3Tradeff
option.Laube-HL
configuration of the mod_O3Sfracrel
option, since it makes little physical sense to have parameters specific to high latitudes in a global model like OSCAR.Daniel2010-lin
configuration of the mod_O3Snitrous
option, since a non-saturating effect of N2O onto stratospheric O3 lacked physical ground.r_HWP
(the response function for harvested wood products). The previous discretization caused delayed and therefore too low emissions.k_BC_adjust
for the option mod_BCadjust == CSIRO
. This has very little impact on a Monte Carlo run.radeff_O3t
for the option mod_O3Tradeff == mean_ACCMIP
. This has no impact on a Monte Carlo run.beta_npp0
and CO2_comp
parameters when isolating the urban biome, to prevent NaN
from appearing during the simulation.NaN
values no longer appear in the alpha_BB
parameters when isolating the urban biome.AeroChem_ACCMIP
input data for the CSIRO-Mk360
configuration. This was causing biased atmospheric lifetimes for POA and BC aerosols under this configuration (and slightly biased ones under the average mean-ACCMIP
configuration).f_pCO2
functions, causing a too efficient ocean carbon sink under high warming and high atmospheric CO2.f_pH
function, causing unrealistic surface ocean pH changes.Initial release on GitHub. Exact model used by Gasser et al. (2017).
v3.1 (partial) | Gasser, T., L. Crepin, Y. Quilcaille, R. A. Houghton, P. Ciais & M. Obersteiner. "Historical CO2 emissions from land-use and land-cover change and their uncertainty." Biogeosciences 17: 4075–4101 (2020). doi:10.5194/bg-17-4075-2020
v2.3 (partial) | Gasser, T., M. Kechiar, P. Ciais, E. J. Burke, T. Kleinen, D. Zhu, Y. Huang, A. Ekici & M. Obersteiner. "Path-dependent reductions in CO2 emission budgets caused by permafrost carbon release." Nature Geoscience 11: 830-835 (2018). doi:10.1038/s41561-018-0227-0
v2.2 (full) | Gasser, T., P. Ciais, O. Boucher, Y. Quilcaille, M. Tortora, L. Bopp & D. Hauglustaine. "The compact Earth system model OSCAR v2.2: description and first results." Geoscientific Model Development 10: 271-319 (2017). doi:10.5194/gmd-10-271-2017
RetroSearch is an open source project built by @garambo | Open a GitHub Issue
Search and Browse the WWW like it's 1997 | Search results from DuckDuckGo
HTML:
3.2
| Encoding:
UTF-8
| Version:
0.7.3