Description Usage Arguments Details Value Note Author(s) References See Also Examples
A simple model to calculate snow melt based on empirical melting factors.
1 2 3 4 5 6 7 | snowMelt(airT, netRad, glacierMask, snowMask, tUnit = "K",
snowTMF = 45*10^-4, disSnowTMF = stack(),
snowRMF = 0.53*10^-4, disSnowRMF = stack(),
tuningFacAirT = 1, disTuningFacAirT = stack(),
decimalPlaces = 4, outType = "mean", writeOutput = FALSE,
outputName = "snowMelt", tmpCreate = FALSE,
tmpDir = "", outDir = "", ... )
|
airT |
An object of class |
netRad |
An object of class |
glacierMask |
An object of class |
snowMask |
An object of class |
tUnit |
An object of class |
snowTMF |
An object of class |
disSnowTMF |
An object of class |
snowRMF |
An object of class |
disSnowRMF |
An object of class |
tuningFacAirT |
An object of class |
disTuningFacAirT |
An object of class |
decimalPlaces |
An object of class |
outType |
An object of class |
writeOutput |
An object of class |
outputName |
An object of class |
tmpCreate |
An object of class |
tmpDir |
An object of class |
outDir |
An object of class |
... |
Further arguments. |
An enhanced degree-day model (e.g. Hock, 2003, 2005; Pellicciotti & alii, 2005) is applied to quantify glacier mass loss ascribed to melted snow and/or firn using empirical temperature (TMF) and radiative melting factors (RMF). Fore more information please refer to the examples below or the original publication (Groos et al., submitted, Equation 10)
An object of class 'RasterLayer'
returning the calculated spatial distribution of snow melt (e.g. in m d-1, depending on 'tmpRes'
).
The following input variables are the requested minimum to run the model:
'airT'
(for every time step)
'netRad'
(for every time step)
'glacierMask'
(stationary or for every time step)
'snowMask'
(stationary or for every time step)
A default value (constant in space and time) is given for each additional argument like 'snowTMF'
or 'snowRMF'
. If desired, the default parameters can be modified. Furthermore, there is the option to pass distributed values (stationary or for every time step) instead of general values using the related 'dis*'-arguments like 'disSnowTMF'
or 'disSnowRMF'
. In this case, the general parameter is ignored.
File format of written ouput: GeoTIFF.
Alexander R. Groos (alexander.groos@giub.unibe.ch)
Groos, A.R., Mayer, C., Smiraglia, C., Diolaiuti, G., and Lambrecht A. (submitted). A first attempt to model region-wide glacier surface mass balances in the Karakoram: findings and future challenges. Geografia Fisica e Dinamica Quaternaria.
Hock, R. (2003). Temperature index melt modelling in mountain areas. Journal of Hydrology 282, 104-115.
Hock, R. (2005). Glacier melt: a review of processes and their modelling. Progress in Physical Geography 29, 362-391.
Pellicciotti F., Brock B., Strasser U., Burlando P., Funk M. and Corripio J. (2005). An enhanced temperature-index glacier melt model including the shortwave radiation balance: development and testing for Haut Glacier d'Arolla, Switzerland. Journal of Glaciology, 51, 573-587.
glacialMelt
, iceMelt
, debrisCoveredIceMelt
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 | # Load the provided RasterLayer objects as exemplary
# input for the function
data(airTemperature_30m_daily, netRad_30m_hourly,
glacierMask_30m, firnMask_30m, package = "glacierSMBM")
# Individual RasterLayer objects should be loaded or
# created using the function raster()
# Include RasterLayer in RasterStack
AirTemperature_30m_daily <- stack(airTemperature_30m_daily * 1.01)
NetRad_30m_hourly <- stack(netRad_30m_hourly)
GlacierMask_30m <- stack(glacierMask_30m)
FirnMask_30m <- stack(firnMask_30m)
# Calculate snow melt using standard settings
output <- snowMelt(airT = AirTemperature_30m_daily,
netRad = NetRad_30m_hourly, glacierMask = GlacierMask_30m,
snowMask = FirnMask_30m)
# Plot output
plot(output, main = "snow melt",
legend.args=list(text='Snow melt (m d-1)', side=3, line=1.5))
# Calculate snow melt using modified setting (e.g. air
# temperature in degree Celsius instead of Kelvin; changes
# melting factors)
# Therefore exemplarily convert temperature from kelvin to
# degree Celsius
airTcelsius <- subset(AirTemperature_30m_daily, 1) - 273.15
# Include RasterLayer in RasterStack
airTcelsius <- stack(airTcelsius)
output <- snowMelt(airT = airTcelsius,
netRad = NetRad_30m_hourly, glacierMask = GlacierMask_30m,
tUnit = "C", snowMask = FirnMask_30m, snowTMF = 75*10^-4,
snowRMF = 1.2*10^-4)
# Plot output
plot(output, main = "snow melt",
legend.args=list(text='Snow melt (m d-1)', side=3, line=1.5))
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.