You signed in with another tab or window. Reload to refresh your session.You signed out in another tab or window. Reload to refresh your session.You switched accounts on another tab or window. Reload to refresh your session.Dismiss alert
Copy file name to clipboardExpand all lines: paper/paper.md
+20-30Lines changed: 20 additions & 30 deletions
Display the source diff
Display the rich diff
Original file line number
Diff line number
Diff line change
@@ -47,36 +47,26 @@ Lastly, production rates produced by monteProd() can be projected into deep time
47
47
Overall, LithoGas provides a straightforward workflow for modelling multiple samples from a single structured input dataframe through to publication-ready summaries and plots that inform exploration. Functions include: 1) performing the Monte Carlo modelling via monteProd(), ingesting data as a structured dataframe (Table 1); 2) summarising Monte Carlo results via monteSum() and plotting novel source-volume-scaling plots via monteH2Plot() and monteHePlot() (Figure 1). Example datasets include dataframes structured for Monte Carlo model input with and without known rock properties (structuredDF) and the summarised lithology distributions from The Canadian Rock Physical Property Database (CRPPData).
48
48
49
49
50
-
#### Table 1: Example layout of structured dataframe required for input into Radiolysis monteProd() function. The function data from the input dataframe by column name (e.g., $uMin, $rockDenMax), as such column names must match exactly.
|uMin |Numeric|Minimum of modelled uranium concentration distribution (ppm)|
56
-
|uMax |Numeric|Maximum of modelled uranium concentration distribution (ppm)|
57
-
|uMean|Numeric|Mean of modelled uranium concentration distribution (ppm)|
58
-
|uSD|Numeric|Standard deviation of modelled uranium concentration distribution (ppm)|
59
-
|thMin |Numeric|Minimum of modelled thorium concentration distribution (ppm)|
60
-
|thMax |Numeric|Maximum of modelled thorium concentration distribution (ppm)|
61
-
|thMean|Numeric|Mean of modelled thorium concentration distribution (ppm)|
62
-
|thSD|Numeric|Standard deviation of modelled thorium concentration distribution (ppm)|
63
-
|kMin |Numeric|Minimum of modelled potassium concentration distribution (wt%)|
64
-
|kMax |Numeric|Maximum of modelled potassium concentration distribution (wt%)|
65
-
|kMean|Numeric|Mean of modelled potassium concentration distribution (wt%)|
66
-
|kSD|Numeric|Standard deviation of modelled potassium concentration distribution (wt5)|
67
-
|fluDenMin |Numeric|Minimum of modelled fluid density distribution (g/cm3)|
68
-
|fluDenMax |Numeric|Maximum of modelled fluid density distribution (g/cm3)|
69
-
|fluDenMean|Numeric|Mean of modelled fluid density distribution (g/cm3)|
70
-
|fluDenSD|Numeric|Standard deviation of modelled fluid density distribution (g/cm3)|
71
-
lithLith|String|Lithologiy/category of desired rock properties from Canadian Rock Physical Property Database|
72
-
|rockDenMin |Numeric|Minimum of modelled rock density distribution (g/cm3)|
73
-
|rockDenMax |Numeric|Maximum of modelled rock density distribution (g/cm3)|
74
-
|rockDenMean|Numeric|Mean of modelled rock density distribution (g/cm3)|
75
-
|rockDenSD|Numeric|Standard deviation of modelled rock density distribution (g/cm3)|
76
-
|porMin |Numeric|Minimum of modelled grain porosity distribution (decimal fraction)|
77
-
|porMax |Numeric|Maximum of modelled grain porosity distribution (decimal fraction)|
78
-
|porMean|Numeric|Mean of modelled grain porosity distribution (decimal fraction)|
79
-
|porSD|Numeric|Standard deviation of grain porosity distribution (decimal fraction)|
50
+
#### Table 1: Columns required for different models within monteProd() function. monteProd() gathers input data by column name (e.g., $uMin, $rockDenMax), as such column names must match exactly. Parameters listed with "..." end of column the name (e.g., "u...") require a min, max, mean, and standard deviation columns (e.g., "uMin", "uMax", "uMean", "uSD").
51
+
52
+
| Parameter (...Min/...Max/...Mean/...SD) | Type | Description | Model |
| Sample | Text | Sample name – not necessarily named “Sample” but some sort of unique sample ID | All models |
55
+
| Age_Ma | Numeric | The age of the rock in millions of year (Ma) | All models |
56
+
| AgeUnc2S_Ma | Numeric | Uncertainty as of age (2S /2σ) in absolute terms (Ma) | All models |
57
+
| u... | Numeric | Minimum, maximum, mean, standard deviation of modelled distribution of uranium concentration - in ppm. Each uMin, uMean, uMax, uSD, are separate columns. | Radiolysis – monteRad() |
58
+
| th... | Numeric | Minimum, maximum, mean, standard deviation of modelled distribution of thorium concentration - in ppm. Each thMin, thMean, thMax, thSD, are separate columns. | Radiolysis – monteRad() |
59
+
| k... | Numeric | Minimum, maximum, mean, standard deviation of modelled distribution of potassium concentration - in weight percent potassium (wt% K). Each kMin, kMean, kMax, kSD, are separate columns. | Radiolysis – monteRad() |
60
+
| fluDen... | Numeric | Minimum, maximum, mean, standard deviation of modelled distribution of fluid density in g/cm3. Each fluDenMin, fluDenMax, fluDenMean, fluDenSD, are separate columns | Radiolysis – monteRad() |
61
+
| Fe3FeTRatInit... | Numeric | Minimum, maximum, mean, standard deviation of modelled distribution of initial Fe3+/FeT ratio at time of rock formation/crystallisation. This ratio is unitless. | Serpentinization, both monteSerpFeSpecies() and monteSerpFeTotal() |
62
+
| Fe2O3... | Numeric | Minimum, maximum, mean, standard deviation of modelled distribution of measured concentration of Fe2O3, as per iron speciation measurements (iron titration, Mossbauer). In units of weight percent Fe2O3 (wt% Fe2O3). | Serpentinization - monteSerpFeSpecies() |
63
+
| FeO... | Numeric | Minimum, maximum, mean, standard deviation of modelled distribution of measured concentration of FeO, as per iron speciation measurements (iron titration, Mossbauer). In units of weight percent FeO (wt% FeO). | Serpentinization - monteSerpFeSpecies() |
64
+
| Fe2O3T... | Numeric | Minimum, maximum, mean, standard deviation of modelled distribution of measured total iron represented as Fe2O3. In units of weight percent Fe2O3 (wt% Fe2O3T). | Serpentinization - monteSerpFeTotal() |
65
+
| Fe3FeTRatCur... | Numeric | Minimum, maximum, mean, standard deviation of modelled distribution of initial Fe3+/FeT ratio at time of measurement. This ratio is unitless. | Serpentinization - monteSerpFeTotal() |
66
+
| litLith | Categorical selection | Lithologic groups from Enkin et al. (2018), inclusion of a categorial rock type will draw rock density and porosity from the known distribution and ignore porosity and density distributions input. If the column is left blank, or the column is not included, the model will pull from the rock density and porosity columns | All models, if not assigned by rockDen and por Min/Max/Mean/SD |
67
+
| rockDen... | Numeric | Minimum, maximum, mean, standard deviation of modelled distribution, in units of g/cm3 | All models, if not assigned by litLith field and therefore joinLitProp function |
68
+
| por... | Numeric | Minimum, maximum, mean, standard deviation of modelled distribution of rock porosity (grain porosity) in fractional porosity (0-1) | Radiolysis, if not assigned by litLith field and therefore joinLitProp function |
69
+
80
70
81
71
# Figures
82
72

0 commit comments