scholarly journals Evaluation of different methods to model near-surface turbulent fluxes for an alpine glacier in the Cariboo Mountains, BC, Canada

Author(s):  
Valentina Radić ◽  
Brian Menounos ◽  
Joseph Shea ◽  
Noel Fitzpatrick ◽  
Mekdes A. Tessema ◽  
...  

Abstract. As part of surface energy balance models used to simulate glacier melting, choosing parameterizations to adequately estimate turbulent heat fluxes is extremely challenging. This study aims to evaluate a set of bulk methods commonly used to estimate turbulent heat fluxes for a sloped glacier surface. The methods differ in their parameterizations of the bulk exchange coefficient that relates the fluxes to the mean meteorological variables measured 2 m above a glacier surface. The performance of 23 bulk approaches in simulating 30-min sensible and latent heat fluxes is evaluated against the measured fluxes from an open path eddy-covariance (OPEC) method. The evaluation is performed at a point scale of an alpine glacier, using one-level meteorological and OPEC observations from a multi-day periods in the 2010 and 2012 summer season. The analysis of the two independent seasons yielded similar findings, listed as following. The bulk method, with or without the commonly used Monin–Obukhov (M–O) stability functions, overestimates the turbulent heat fluxes over the observational period, mainly due to an overestimation of the momentum flux. In the absence of OPEC-derived M–O stability parameter, no method can successfully predict this parameter, which results in poor performances of the M–O stability corrections and consequently the bulk method. The OPEC-derived 30-min momentum flux is linearly related to the measured wind speed, contrary to the proposed quadratic relation by the commonly used bulk methods. An approach based on a katabatic flow model, which assumes a linear relation between the shear stress and the wind speed, outperforms any other bulk approach that we tested in simulating the momentum flux. In agreement with the katabatic flow model, we show that in a more stable atmosphere the bulk exchange coefficient for momentum is smaller. The sensible heat flux can be more successfully modeled if the bulk exchange coefficients for momentum and heat are allowed to follow different parametrization schemes, rather than assuming equal schemes as is the case in the common bulk methods. Further data from different glaciers are needed to investigate any universality of these findings.

2017 ◽  
Vol 11 (6) ◽  
pp. 2897-2918 ◽  
Author(s):  
Valentina Radić ◽  
Brian Menounos ◽  
Joseph Shea ◽  
Noel Fitzpatrick ◽  
Mekdes A. Tessema ◽  
...  

Abstract. As part of surface energy balance models used to simulate glacier melting, choosing parameterizations to adequately estimate turbulent heat fluxes is extremely challenging. This study aims to evaluate a set of four aerodynamic bulk methods (labeled as C methods), commonly used to estimate turbulent heat fluxes for a sloped glacier surface, and two less commonly used bulk methods developed from katabatic flow models. The C methods differ in their parameterizations of the bulk exchange coefficient that relates the fluxes to the near-surface measurements of mean wind speed, air temperature, and humidity. The methods' performance in simulating 30 min sensible- and latent-heat fluxes is evaluated against the measured fluxes from an open-path eddy-covariance (OPEC) method. The evaluation is performed at a point scale of a mountain glacier, using one-level meteorological and OPEC observations from multi-day periods in the 2010 and 2012 summer seasons. The analysis of the two independent seasons yielded the same key findings, which include the following: first, the bulk method, with or without the commonly used Monin–Obukhov (M–O) stability functions, overestimates the turbulent heat fluxes over the observational period, mainly due to a substantial overestimation of the friction velocity. This overestimation is most pronounced during the katabatic flow conditions, corroborating the previous findings that the M–O theory works poorly in the presence of a low wind speed maximum. Second, the method based on a katabatic flow model (labeled as the KInt method) outperforms any C method in simulating the friction velocity; however, the C methods outperform the KInt method in simulating the sensible-heat fluxes. Third, the best overall performance is given by a hybrid method, which combines the KInt approach with the C method; i.e., it parameterizes eddy viscosity differently than eddy diffusivity. An error analysis reveals that the uncertainties in the measured meteorological variables and the roughness lengths produce errors in the modeled fluxes that are smaller than the differences between the modeled and observed fluxes. This implies that further advances will require improvement to model theory rather than better measurements of input variables. Further data from different glaciers are needed to investigate any universality of these findings.


2020 ◽  
Vol 37 (4) ◽  
pp. 589-603 ◽  
Author(s):  
Xiangzhou Song

AbstractSea surface currents are commonly neglected when estimating the air–sea turbulent heat fluxes in bulk formulas. Using buoy observations in the Bohai Sea, this paper investigated the effects of near-coast multiscale currents on the quantification of turbulent heat fluxes, namely, latent heat flux (LH) and sensible heat flux (SH). The maximum current reached 1 m s−1 in magnitude, and a steady northeastward current of 0.16 m s−1 appeared in the southern Bohai Strait. The predominant tidal signal was the semidiurnal current, followed by diurnal components. The mean absolute surface wind was from the northeast with a speed of approximately 3 m s−1. The surface winds at a height of 11 m were dominated by the East Asian monsoon. As a result of upwind flow, the monthly mean differences in LH and SH between the estimates with and without surface currents ranged from 1 to 2 W m−2 in July (stable boundary layer) and November (unstable boundary layer). The hourly differences were on average 10 W m−2 and ranged from 0 to 24 W m−2 due to changes in the relative wind speed by high-frequency rotating surface tidal currents. The diurnal variability in LH/SH was demonstrated under stable and unstable boundary conditions. Observations provided an accurate benchmark for flux comparisons. The newly updated atmospheric reanalysis products MERRA-2 and ERA5 were superior to the 1° OAFlux data at this buoy location. However, future efforts in heat flux computation are still needed to, for example, consider surface currents and resolve diurnal variations.


2000 ◽  
Vol 46 (154) ◽  
pp. 445-452 ◽  
Author(s):  
Bruce Denby ◽  
W. Greuell

AbstractA one-dimensional second-order closure model and in situ observations on a melting glacier surface are used to investigate the suitability of bulk and profile methods for determining turbulent fluxes in the presence of the katabatic wind-speed maximum associated with glacier winds. The results show that profile methods severely underestimate turbulent fluxes when a wind-speed maximum is present. The bulk method, on the other hand, only slightly overestimates the turbulent heat flux in the entire region below the wind-speed maximum and is thus much more appropriate for use on sloping glacier surfaces where katabatic winds dominate and wind-speed maxima are just a few meters above the surface.


2013 ◽  
Vol 54 (63) ◽  
pp. 41-51 ◽  
Author(s):  
J.P. Conway ◽  
N.J. Cullen

AbstractThe turbulent sensible and latent heat fluxes are important components of the surface energy balance over glaciers in the Southern Alps of New Zealand, contributing over half the energy available for ablation during large melt events. To calculate these terms confidently in glacier mass-balance models it is essential to use appropriate parameterizations for surface roughness and atmospheric stability. Eddy covariance measurements at Brewster Glacier were obtained over an ice surface to help facilitate an assessment of the calculation of the turbulent heat fluxes. The roughness length for momentum was found to be 3.6 x 10−3m, while the roughness lengths for temperature and humidity were two orders of magnitude smaller, in agreement with surface renewal theory. A Monte Carlo approach was used to assess the uncertainty in turbulent heat fluxes calculated using the bulk aerodynamic method. It was found that input-data and roughness-length uncertainty could not explain underestimates of observed sensible heat fluxes during periods with low wind speed and large temperature gradients. During these periods a katabatic wind speed maximum alters the formulation of the turbulent exchange coefficient to that typically observed in a neutral atmosphere and this has implications for glacier mass-balance sensitivity.


2005 ◽  
Vol 42 ◽  
pp. 367-372 ◽  
Author(s):  
C. Rolstad ◽  
J. Oerlemans

AbstractThe surface energy balance of glaciers is studied to determine their sensitivity to climate variations. It is known that the turbulent heat fluxes are sensitive to increases in temperature. Automatic weather station data from ablation regions are used to measure melt rates, radiative fluxes and the meteorological data required to determine turbulent heat fluxes using bulk formulas. The turbulent exchange coefficient must be determined for closure of the energy budget. The available methods are the eddy correlation method, the profile method and the residual method, which is applied and tested here. In the residual method the coefficient is determined by fitting a calculated melt curve to an observed melt curve. The coefficients are estimated for three sites: for Vatnajökull, Iceland, Ch = (1.3 ± 0.55) × 10–3 (1998) and Ch = (2.5 ±1.1) × 10–3 (1999); for Morteratschgletscher, Switzerland, Ch = (2.1 ±0.55) × 10–3 (1998); and for West Greenland, Ch = (2.0 ±0.52) × 10–3 (1998-2000). It is found that the coefficient can be determined to within 26% uncertainty under the following conditions: all terms in the energy balance are measured, there is no differential melt on the glacier surface, the melt curves are fitted when the entire snow layer has melted, and the measurement period is several weeks.


1999 ◽  
Vol 11 (1) ◽  
pp. 93-99 ◽  
Author(s):  
S. Argentini ◽  
G. Mastrantonio ◽  
A. Viola

Simultaneous acoustic Doppler sodar and tethersonde measurements were used to study some of the characteristics of the unstable boundary layer at Dumont d'Urville, Adélie Land, East Antarctica during the summer 1993–94. A description of the convective boundary layer and its behaviour in connection with the wind regime is given along with the frequency distribution of free convection episodes. The surface heat flux has been evaluated using the vertical velocity variance derived from sodar measurements. The turbulent exchange coefficients, estimated by coupling sodar and tethered balloon measurements, are in strong agreement with those present in literature for the Antarctic regions.


2021 ◽  
Author(s):  
Andrew Bennett ◽  
Bart Nijssen

<p>Machine learning (ML), and particularly deep learning (DL), for geophysical research has shown dramatic successes in recent years. However, these models are primarily geared towards better predictive capabilities, and are generally treated as black box models, limiting researchers’ ability to interpret and understand how these predictions are made. As these models are incorporated into larger models and pushed to be used in more areas it will be important to build methods that allow us to reason about how these models operate. This will have implications for scientific discovery that will ensure that these models are robust and reliable for their respective applications. Recent work in explainable artificial intelligence (XAI) has been used to interpret and explain the behavior of machine learned models.</p><p>Here, we apply new tools from the field of XAI to provide physical interpretations of a system that couples a deep-learning based parameterization for turbulent heat fluxes to a process based hydrologic model. To develop this coupling we have trained a neural network to predict turbulent heat fluxes using FluxNet data from a large number of hydroclimatically diverse sites. This neural network is coupled to the SUMMA hydrologic model, taking imodel derived states as additional inputs to improve predictions. We have shown that this coupled system provides highly accurate simulations of turbulent heat fluxes at 30 minute timesteps, accurately predicts the long-term observed water balance, and reproduces other signatures such as the phase lag with shortwave radiation. Because of these features, it seems this coupled system is learning physically accurate relationships between inputs and outputs. </p><p>We probe the relative importance of which input features are used to make predictions during wet and dry conditions to better understand what the neural network has learned. Further, we conduct controlled experiments to understand how the neural networks are able to learn to regionalize between different hydroclimates. By understanding how these neural networks make their predictions as well as how they learn to make predictions we can gain scientific insights and use them to further improve our models of the Earth system.</p>


2021 ◽  
Vol 22 (10) ◽  
pp. 2547-2564
Author(s):  
Georg Lackner ◽  
Daniel F. Nadeau ◽  
Florent Domine ◽  
Annie-Claude Parent ◽  
Gonzalo Leonardini ◽  
...  

AbstractRising temperatures in the southern Arctic region are leading to shrub expansion and permafrost degradation. The objective of this study is to analyze the surface energy budget (SEB) of a subarctic shrub tundra site that is subject to these changes, on the east coast of Hudson Bay in eastern Canada. We focus on the turbulent heat fluxes, as they have been poorly quantified in this region. This study is based on data collected by a flux tower using the eddy covariance approach and focused on snow-free periods. Furthermore, we compare our results with those from six Fluxnet sites in the Arctic region and analyze the performance of two land surface models, SVS and ISBA, in simulating soil moisture and turbulent heat fluxes. We found that 23% of the net radiation was converted into latent heat flux at our site, 35% was used for sensible heat flux, and about 15% for ground heat flux. These results were surprising considering our site was by far the wettest site among those studied, and most of the net radiation at the other Arctic sites was consumed by the latent heat flux. We attribute this behavior to the high hydraulic conductivity of the soil (littoral and intertidal sediments), typical of what is found in the coastal regions of the eastern Canadian Arctic. Land surface models overestimated the surface water content of those soils but were able to accurately simulate the turbulent heat flux, particularly the sensible heat flux and, to a lesser extent, the latent heat flux.


2015 ◽  
Vol 42 (6) ◽  
pp. 1856-1862 ◽  
Author(s):  
A. B. Villas Bôas ◽  
O. T. Sato ◽  
A. Chaigneau ◽  
G. P. Castelão

Sign in / Sign up

Export Citation Format

Share Document