-
Notifications
You must be signed in to change notification settings - Fork 3
/
MimiBRICK.jl
168 lines (131 loc) · 8.24 KB
/
MimiBRICK.jl
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
module MimiBRICK
# Load required packages
using CSVFiles, DataFrames, Distributions, Mimi
# Load BRICK Mimi components.
include(joinpath("components", "antarctic_icesheet_component.jl"))
include(joinpath("components", "antarctic_ocean_component.jl"))
include(joinpath("components", "glaciers_small_icecaps_component.jl"))
include(joinpath("components", "global_sea_level_component.jl"))
include(joinpath("components", "greenland_icesheet_component.jl"))
include(joinpath("components", "landwater_storage_component.jl"))
include(joinpath("components", "thermal_expansion_component.jl"))
# include creation functions for two other model variants
include(joinpath("create_models", "SNEASY_BRICK.jl"))
include(joinpath("create_models", "BRICK_DOECLIM.jl"))
# include other helper functions
include(joinpath("downscale.jl"))
# include calibration functions
include(joinpath("calibration/main_calibration.jl"))
# function to create 'Building blocks for Relevant Ice and Climate Knowledge' (BRICK) model.
"""
get_model(;rcp_scenario::String="RCP85", start_year::Int=1850, end_year::Int=2020)
Return a MimiBRICK model instance that can be modified and run.
Function Arguments:
rcp_scenario = RCP scenario for exogenous forcing
start_year = initial year of the simulation period
end_year = ending year of the simulation period
"""
function get_model(;rcp_scenario::String="RCP85", start_year::Int=1850, end_year::Int=2020)
#-----------------------#
# ----- Load Data ----- #
#-----------------------#
# Load exogenous time-series for global surface temperature and ocean heat content (output from SNEASY under RCP8.5).
# NOTE: for now, only `rcp_scenario = "RCP85"` is supported
temperature_scenario = DataFrame(load(joinpath(@__DIR__, "..", "data", "model_data", "sneasy_temperature_"*rcp_scenario*"_1850_2300.csv")))
oceanheat_scenario = DataFrame(load(joinpath(@__DIR__, "..", "data", "model_data", "sneasy_oceanheat_"*rcp_scenario*"_1850_2300.csv")))
#-------------------------#
# ----- Build BRICK ----- #
#-------------------------#
# Create a Mimi-model instance.
brick = Model()
# Set model time horizon, defaults to 1850-2300.
set_dimension!(brick, :time, start_year:end_year)
# Add in BRICK components.
add_comp!(brick, antarctic_ocean)
add_comp!(brick, antarctic_icesheet)
add_comp!(brick, glaciers_small_icecaps)
add_comp!(brick, greenland_icesheet)
add_comp!(brick, thermal_expansion)
add_comp!(brick, landwater_storage)
add_comp!(brick, global_sea_level)
#-------------------------------------#
# ----- Assign Model Parameters ----- #
#-------------------------------------#
# ----- Antarctic Ocean ----- #
update_param!(brick, :antarctic_ocean, :anto_α, 0.28)
update_param!(brick, :antarctic_ocean, :anto_β, 0.95)
# ----- Antarctic Ice Sheet ----- #
update_param!(brick, :antarctic_icesheet, :ais_ρ_ice, 917.0)
update_param!(brick, :antarctic_icesheet, :ais_ρ_seawater, 1030.0)
update_param!(brick, :antarctic_icesheet, :ais_ρ_rock, 4000.0)
update_param!(brick, :antarctic_icesheet, :ais_sea_level₀, 0.0)
update_param!(brick, :antarctic_icesheet, :ais_ocean_temperature₀, 0.72)
update_param!(brick, :antarctic_icesheet, :ais_radius₀, 1.864e6)
update_param!(brick, :antarctic_icesheet, :ais_bedheight₀, 781.0)
update_param!(brick, :antarctic_icesheet, :ais_slope, 0.0006)
update_param!(brick, :antarctic_icesheet, :ais_μ, 11.0)
update_param!(brick, :antarctic_icesheet, :ais_runoffline_snowheight₀, 1400.0)
update_param!(brick, :antarctic_icesheet, :ais_c, 100.0)
update_param!(brick, :antarctic_icesheet, :ais_precipitation₀, 0.37)
update_param!(brick, :antarctic_icesheet, :ais_κ, 0.062)
update_param!(brick, :antarctic_icesheet, :ais_ν, 0.0086)
update_param!(brick, :antarctic_icesheet, :ais_iceflow₀, 1.2)
update_param!(brick, :antarctic_icesheet, :ais_γ, 2.9)
update_param!(brick, :antarctic_icesheet, :ais_α, 0.23)
update_param!(brick, :antarctic_icesheet, :ais_temperature_coefficient, 0.8365)
update_param!(brick, :antarctic_icesheet, :ais_temperature_intercept, 15.42)
update_param!(brick, :antarctic_icesheet, :ais_local_fingerprint, -1.18)
update_param!(brick, :antarctic_icesheet, :ocean_surface_area, 3.619e14)
update_param!(brick, :antarctic_icesheet, :temperature_threshold, -15.0)
update_param!(brick, :antarctic_icesheet, :λ, 0.0093)
update_param!(brick, :antarctic_icesheet, :include_ais_DSL, true)
# ----- Glaciers & Small Ice Caps ----- #
update_param!(brick, :glaciers_small_icecaps, :gsic_β₀, 0.0013)
update_param!(brick, :glaciers_small_icecaps, :gsic_v₀, 0.376)
update_param!(brick, :glaciers_small_icecaps, :gsic_s₀, -0.0138)
update_param!(brick, :glaciers_small_icecaps, :gsic_n, 0.847)
update_param!(brick, :glaciers_small_icecaps, :gsic_teq, -0.15)
# ----- Greenland Ice Sheet ----- #
update_param!(brick, :greenland_icesheet, :greenland_a, -1.37)
update_param!(brick, :greenland_icesheet, :greenland_b, 8.06)
update_param!(brick, :greenland_icesheet, :greenland_α, 0.0008)
update_param!(brick, :greenland_icesheet, :greenland_β, 0.00009)
update_param!(brick, :greenland_icesheet, :greenland_v₀, 7.52)
# ----- Thermal Expansion ----- #
update_param!(brick, :thermal_expansion, :te_A, 3.619e14)
update_param!(brick, :thermal_expansion, :te_C, 3991.86795711963)
update_param!(brick, :thermal_expansion, :te_ρ, 1027.0)
update_param!(brick, :thermal_expansion, :te_α, 0.16)
update_param!(brick, :thermal_expansion, :te_s₀, 0.0)
update_param!(brick, :thermal_expansion, :ocean_heat_mixed, zeros(length(start_year:end_year)))
oceanheat_idx = findall((in)(start_year:end_year), oceanheat_scenario[!,:Year])
update_param!(brick, :thermal_expansion, :ocean_heat_interior, oceanheat_scenario[oceanheat_idx, :"MAP Ocean Heat"])
# ----- Landwater Storage ----- #
update_param!(brick, :landwater_storage, :lws₀, 0.0)
update_param!(brick, :landwater_storage, :first_projection_year, 2018)
update_param!(brick, :landwater_storage, :lws_random_sample, rand(Normal(0.0003, 0.00018), length(start_year:end_year)))
# ----- Set Parameters With Common Values Across Components ----- #
add_shared_param!(brick, :model_seawater_freeze, -1.8)
connect_param!(brick, :antarctic_icesheet, :seawater_freeze, :model_seawater_freeze)
connect_param!(brick, :antarctic_ocean, :seawater_freeze, :model_seawater_freeze)
temperature_idx = findall((in)(start_year:end_year), temperature_scenario[!,:Year])
add_shared_param!(brick, :model_global_surface_temperature, temperature_scenario[temperature_idx,:"MAP Temperature"], dims = [:time])
connect_param!(brick, :antarctic_icesheet, :global_surface_temperature, :model_global_surface_temperature)
connect_param!(brick, :antarctic_ocean, :global_surface_temperature, :model_global_surface_temperature)
connect_param!(brick, :glaciers_small_icecaps, :global_surface_temperature, :model_global_surface_temperature)
connect_param!(brick, :greenland_icesheet, :global_surface_temperature, :model_global_surface_temperature)
#-----------------------------------------#
#----- Create Component Connections ----- #
#-----------------------------------------#
# Set parameter connections (:component => :parameter).
connect_param!(brick, :global_sea_level => :slr_glaciers_small_ice_caps, :glaciers_small_icecaps => :gsic_sea_level)
connect_param!(brick, :global_sea_level => :slr_greeland_icesheet, :greenland_icesheet => :greenland_sea_level)
connect_param!(brick, :global_sea_level => :slr_antartic_icesheet, :antarctic_icesheet => :ais_sea_level)
connect_param!(brick, :global_sea_level => :slr_thermal_expansion, :thermal_expansion => :te_sea_level)
connect_param!(brick, :global_sea_level => :slr_landwater_storage, :landwater_storage => :lws_sea_level)
connect_param!(brick, :antarctic_icesheet => :antarctic_ocean_temperature, :antarctic_ocean => :anto_temperature)
connect_param!(brick, :antarctic_icesheet => :global_sea_level, :global_sea_level => :sea_level_rise)
# Return BRICK model.
return brick
end
end # Module