R/DTRFunctions.R
diurnal_temp_variation_sineexp.RdThe function estimates temperature across hours using a diurnal temperature variation function incorporating sine and exponential components (Parton and Logan 1981) .
diurnal_temp_variation_sineexp(
T_max,
T_min,
t,
t_r,
t_s,
alpha = 2.59,
beta = 1.55,
gamma = 2.2
)numeric maximum and minimum daily temperatures (C).
numeric time for temperature estimate (hour).
numeric times of sunrise and sunset (hour).
numeric time difference between t_x (time of maximum temperature) and noon (hour).
numeric time difference between t_x and sunrise (hour).
numeric decay parameter for rate of t change from sunset to t_n (time of minimum temp).
numeric temperature (C) at a specified hour.
Default alpha, beta, and gamma values are the average of 5 North Carolina sites (Wann et al. 1985)
.
Other alpha, beta, and gamma parameterizations include values for Denver, Colorado from Parton and Logan (1981)
:
150 cm air temperature: alpha = 1.86, beta = 2.20, gamma = -0.17
10 cm air temperature: alpha = 1.52, beta = 2.00, gamma = -0.18
soil surface temperature: alpha = 0.50, beta = 1.81, gamma = 0.49
10cm soil temperature: alpha = 0.45, beta = 2.28, gamma = 1.83
Parton WJ, Logan JA (1981).
“A model for diurnal variation in soil and air temperature.”
Agricultural Meteorology, 23, 205-216.
doi:10.1016/0002-1571(81)90105-9
, https://www.sciencedirect.com/science/article/abs/pii/0002157181901059.
Wann M, Yen D, Gold HJ (1985).
“Evaluation and calibration of three models for daily cycle of air temperature.”
Agricultural and Forest Meteorology, 34, 121-128.
https://doi.org/10.1016/0168-1923(85)90013-9.
Other microclimate functions:
air_temp_profile_neutral(),
air_temp_profile_segment(),
air_temp_profile(),
degree_days(),
direct_solar_radiation(),
diurnal_radiation_variation(),
diurnal_temp_variation_sinesqrt(),
diurnal_temp_variation_sine(),
monthly_solar_radiation(),
partition_solar_radiation(),
proportion_diffuse_solar_radiation(),
solar_radiation(),
surface_roughness(),
wind_speed_profile_neutral(),
wind_speed_profile_segment()
diurnal_temp_variation_sineexp(T_max = 30,
T_min = 10,
t = 11,
t_r = 6,
t_s = 18,
alpha = 2.59,
beta = 1.55,
gamma = 2.2)
#> [1] 23.91954