Simulated Mass Balance Research Articles

The Effect of Physically Based Ice Radiative Processes on Greenland Ice Sheet Albedo and Surface Mass Balance in E3SM

AbstractA significant portion of surface melt on the Greenland Ice Sheet (GrIS) is due to dark ice regions in the ablation zone, where solar absorption is influenced by the physical properties of the ice, light absorbing constituents (LACs), and the overlying crustal surface or melt ponds. Earth system models (ESMs) typically prescribe the albedo of ice surfaces as a constant value in the visible and near‐infrared spectral regions. This work advances ESM ice radiative transfer modeling by (a) incorporating a physically based radiative transfer model (SNow, ICe and Aerosol Radiation model Adding‐Doubling Version 4; SNICAR‐ADv4) into the Energy Exascale Earth System Model (E3SM), (b) determining spatially and temporally varying bare ice physical properties over the GrIS ablation zone from satellite observations to inform SNICAR‐ADv4, and (c) assessing the impacts on simulated GrIS albedo and surface mass balance associated with modeling of more realistic bare ice albedo. GrIS‐wide bare ice albedo in E3SMv2 is overestimated by ∼4% in the visible and ∼7% in the near‐infrared wavelengths compared to the Moderate Resolution Imaging Spectroradiometer. Our bare ice physical property retrieval method found that LACs, ice crustal surfaces, and melt ponds reduce visible albedo by 30% in the bare ice region of the GrIS ablation zone. The realistic bare ice albedo reduces surface mass balance by ∼145 Gt, or 0.4 mm of sea‐level equivalent between 2000 and 2021 compared to the default E3SM. This work highlights the importance of simulating bare ice albedo accurately and realistically to improve our ability to quantify changes in the GrIS surface mass and radiative energy budgets.

Groundwater flow modeling and recharge estimation of heterogeneous aquifer: Applied to Matmata aquifer, southeastern, Tunisia

Groundwater is an important source of water supply for several socioeconomic sectors. In particular, heterogeneous carbonate aquifers are an important potential reservoir of groundwater resources. The Matmata carbonate aquifer, located in southern Tunisia, suffers from over abstraction and severe anthropogenic pressure. This study aims to estimate diffusive recharge rates using climate and soil datasets, and proposes a combination of CMB techniques and numerical simulation methods. Groundwater chloride mass balance and numerical simulations were best practices for achieving these goals. The Hydrus-1D model was applied, calibrated and validated with soil water content measured at different depths and using long-term, diurnal-scale field and climate data. The results highlight that Richard's retention parameters give a close correlation with the observed saturated moisture content with RMSE of 0.07. The calculated evaporation is around 22.5 mm for the entire study period, and the contribution of rainfall to infiltration, which was estimated at 9.8 mm as 8.4%. The CMB model shows a recent recharge rate 4.83 mm. yr−1, corresponding to a groundwater recharge rate of 5.84 % of precipitation. Groundwater recharge estimation is one of the most important concepts of water balance, and its assessment is crucial for groundwater management and protection.

Simulation of Furfural and Levulinic Acid Production from Lignocellulosic Biomass

Furfural and levulinic acid are chemicals that can be produced from lignocellulosic biomass. In this research, mass balance simulation of furfural and levulinic acid production at scale of 100 tons of biomass/year was examined. This research using various models of kinetic reaction and biomass to solvent ratio from lignocellulosic biomass raw materials consisting of cellulose, hemicellulose, lignin, ash, and water composition of 35%, 30%, 25%, 3%, and 7%. The kinetic models used in this research were separate kinetic models and simultaneous model (model-3). The separate kinetic models were divided into two models which separate kinetic model 1 (model-1) and separate kinetic model 2 (model-2). SuperPro Designer 9.0 software was used to calculate mass balance simulation. From the research, it was found that variations in kinetic reaction model affected furfural and levulinic acid production. Higher biomass to solvent ratio produced higher furfural and levulinic acid production. The highest furfural produced from simulation process was kinetic reaction model-2 with a biomass to solvent ratio of 1:30 and 0.67 liter furfural/hour. While the highest levulinic acid produced from the simulation process was kinetic reaction model-1 with a biomass to solvent ratio of 1:30 and 2.37 liter levulinic acid/hour.

Meltwater runoff and glacier mass balance in the high Arctic: 1991–2022 simulations for Svalbard

Abstract. The Arctic is undergoing increased warming compared to the global mean, which has major implications for freshwater runoff into the oceans from seasonal snow and glaciers. Here, we present high-resolution (2.5 km) simulations of glacier mass balance, runoff, and snow conditions on Svalbard from 1991–2022, one of the fastest warming regions in the world. The simulations are created using the CryoGrid community model forced by Copernicus Arctic Regional ReAnalysis (CARRA) (1991–2021) and AROME-ARCTIC forecasts (2016–2022). Updates to the water percolation and runoff schemes are implemented in the CryoGrid model for the simulations. In situ observations available for Svalbard, including automatic weather station data, stake measurements, and discharge observations, are used to carefully evaluate the quality of the simulations and model forcing. We find a slightly negative climatic mass balance (CMB) over the simulation period of −0.08 mw.e.yr-1 but with no statistically significant trend. The most negative annual CMB is found for Nordenskiöldland (−0.73 mw.e.yr-1), with a significant negative trend of −0.27 mw.e. per decade for the region. Although there is no trend in the annual CMB, we do find a significant increasing trend in the runoff from glaciers of 0.14 mw.e. per decade. The average runoff was found to be 0.8 mw.e.yr-1. We also find a significant negative trend in the refreezing of −0.13 mw.e. per decade. Using AROME-ARCTIC forcing, we find that 2021/22 has the most negative CMB and highest runoff over the 1991–2022 simulation period investigated in this study. We find the simulated climatic mass balance and runoff using CARRA and AROME-ARCTIC forcing are similar and differ by only 0.1 mw.e.yr-1 in climatic mass balance and by 0.2 mw.e.yr-1 in glacier runoff when averaged over all of Svalbard. There is, however, a clear difference over Nordenskiöldland, where AROME-ARCTIC simulates significantly higher mass balance and significantly lower runoff. This indicates that AROME-ARCTIC may provide similar high-quality predictions of the total mass balance of Svalbard as CARRA, but regional uncertainties should be taken into consideration. The simulations produced for this study are made publicly available at a daily and monthly resolution, and these high-resolution simulations may be re-used in a wide range of applications including studies on glacial runoff, ocean currents, and ecosystems.

A wind-driven snow redistribution module for Alpine3D v3.3.0: adaptations designed for downscaling ice sheet surface mass balance

Abstract. Ice sheet surface mass balance describes the net snow accumulation at the ice sheet surface. On the Antarctic ice sheet, winds redistribute snow, resulting in a surface mass balance that is variable in both space and time. Representing wind-driven snow redistribution processes in models is critical for local assessments of surface mass balance, repeat altimetry studies, and interpretation of ice core accumulation records. To this end, we have adapted Alpine3D, an existing distributed snow modeling framework, to downscale Antarctic surface mass balance to horizontal resolutions up to 1 km. In particular, we have introduced a new two-dimensional advection-based wind-driven snow redistribution module that is driven by an offline coupling between WindNinja, a wind downscaling model, and Alpine3D. We then show that large accumulation variability can be at least partially explained by terrain-induced wind speed variations which subsequently redistribute snow around rolling topography. By comparing Alpine3D to airborne-derived snow accumulation measurements within a testing domain over Pine Island Glacier in West Antarctica, we demonstrate that our Alpine3D downscaling approach improves surface mass balance estimates when compared to the Modern-Era Retrospective analysis for Research and Applications, Version 2 (MERRA-2), a global atmospheric reanalysis which we use as atmospheric forcing. In particular, when compared to MERRA-2, Alpine3D reduces simulated surface mass balance root mean squared error by 23.4 mmw.e.yr-1 (13 %) and increases variance explained by 24 %. Despite these improvements, our results demonstrate that considerable uncertainty stems from the employed saltation model, confounding simulations of surface mass balance variability.

Brief communication: Non-linear sensitivity of glacier mass balance to climate attested by temperature-index models

Abstract. Temperature-index models have been widely used for glacier-mass projections spanning the 21st century. The ability of temperature-index models to capture non-linear responses of glacier surface mass balance (SMB) to high deviations in air temperature and solid precipitation was recently discussed in the context of mass-balance simulations employing advanced machine-learning techniques. Here, we performed numerical experiments with a classic temperature-index model and confirmed that such models are capable of detecting non-linear responses of glacier SMB to temperature and precipitation changes. Non-linearities derive from the change in the degree-day factor over the ablation season and from the lengthening of the ablation season.

Tourmaline Breccias from the Río Blanco-Los Bronces Porphyry Copper District, Chile: Constraints on the Fluid Source and the Utility of Tourmaline Composition for Exploration

Abstract Tourmaline-cemented breccia bodies host much of the ore in the Río Blanco-Los Bronces porphyry Cu-Mo deposits. We determined the chemical and B isotope composition of tourmaline as well as S isotope ratios of anhydrite and sulfide minerals to shed light on the composition and origin of mineralizing fluids. Also, the utility of tourmaline as an indicator mineral was tested by comparing mineralized and barren breccias. Tourmaline in mineralized samples has a narrow Mg range (1.5–2 apfu) and variable, generally low Al contents (4–6.5 apfu). A strong negative correlation of Al with Fe indicates monovalent substitution of Al and Fe3+, implying relatively oxidizing fluids. In contrast, tourmaline from barren breccias has a narrower Al range (6–7 apfu), lower and more variable Mg (0.2–2.5 apfu), and a strong negative Mg-Fe correlation, suggesting more reduced fluids with a dominance of Fe2+. These features and the implications of redox contrast may have exploration significance. Tourmaline from all breccia samples yielded δ11B values from 1.8 to 7.9‰. A magmatic source of boron is concluded from the identical B isotope values of granite-hosted tourmaline in the district (1.2–7.7‰) and from the similar range of regional volcanic and porphyry rocks in the Central Andes. The δ34S values of coexisting anhydrite (11.6–14.5‰) and chalcopyrite (–1.5 to –0.2‰) in mineralized breccia give S isotope exchange temperatures of 377° to 437°C, consistent with fluid inclusion temperatures. Total sulfur δ34Sfluid estimates between 1.4 ± 3.9 and 8.8 ± 1.3‰ are broadly consistent with a magmatic source but not well constrained. However, published O and H isotope ratios of quartz and tourmaline from the Río Blanco-Los Bronces breccias have a clear magmatic signature, so this is the preferred scenario. Mass balance simulations of the boron budget show that typical magma flux rates, water contents, and boron concentration for the Central Andes could produce the estimated 107 tons of boron in the Río Blanco-Los Bronces breccias within the 4-m.y. duration of porphyry intrusions if (1) magma accumulated and evolved at midcrustal levels before emplacement and (2) boron partitioned strongly to the fluid phase (DBfluid/melt > 3).

Temporospatial Variation, Chemical Composition, and Source Resolution of PM2.5 in the Southeastern Taiwan Island

This study investigated the temporospatial variation, chemical composition, and source resolution of fine particles (PM2.5) in the southeastern seas of the Taiwan Island. 24-hr PM2.5 was sampled simultaneously at two remote sites, the Green Island (West Pacific Ocean; WPO) and the Kenting Peninsula (northern Bashi Channel; BC), in four seasons. After sampling, the chemical fingerprints of PM2.5 were characterized and further applied to resolve the potential sources of PM2.5 and their contribution by using a receptor model on the basis of chemical mass balance (CMB), enrichment factor (EF), and backward trajectory simulation. It showed that PM2.5 concentrations in winter (10.8 µg m–3) and spring (12.0 µg m–3) (i.e., during the period of Asian Northeastern Monsoons; ANMs) were higher than those in summer (4.0 µg m–3) and fall (6.6 µg m–3). In terms of chemical composition of PM2.5, secondary inorganic aerosols (SIAs = NO3–, SO42–, and NH4+) (56.7–67.2%) were the dominant component of water-soluble ions (WSIs) in PM2.5, while crustal elements (Mg, Al, Ca, Fe, and K) (44.0–61.2%) dominated the metallic contents in PM2.5. High EF values (> 10) showed that V, Mn, Ni, Cu, and Zn were potentially contributed from anthropogenic sources. Moreover, organic carbon (OC) (0.6 µg m–3) was superior to elemental carbon (EC) (0.3 µg m–3) in PM2.5. The OC/EC ratios higher than 2.0 showed the potential chemical formation of secondary organic aerosols (SOAs) in the atmosphere in winter and spring. Trajectory simulation indicated that high PM2.5 concentrations were mostly originated from North and Central China, Japan islands, and Korea Peninsula. Major sources of PM2.5 resolved by CMB receptor modeling were ordered as: sea salts (19.9%) > fugitive dust (19.8%) > industrial boilers (oil-fired) (10.8%) > secondary sulfate (9.8%) > mobile sources (8.0%).

Direct, geodetic and simulated mass balance studies of the Kolahoi Glacier in the Kashmir Himalaya, India

Knowledge of glacier mass balance (MB) is critical for understanding glacier-climate interactions and projecting future water availability in the data-scarce Himalayan region. In this study, we present the glaciological (2014–2019), geodetic (2000–2014) and reconstructed (1980–2019) MB of the Kolahoi Glacier (KG) in the Kashmir Himalaya, India. During the five-year observation period, the average in situ MB of the Glacier was −0.83 ± 0.34 m w.e. a−1, with significant inter-annual variation. The mean annual geodetic MB of the KG between 2000 and 2014 was −0.90 ± 0.09 m w.e. The reconstructed MB revealed four distinct glacier recession episodes since 1980. The Period I (1980–1990) showed a reduced mass loss of −0.27 ± 0.42 m w.e. a−1, with four years of positive mass balance. The Period II (2000–2010), on the other hand, experienced the highest mass loss of −1.18 m w.e. a−1 throughout the entire reconstructed MB series. All three MB measurement- in situ, geodetic and simulated agreed well. The MB showed an excellent correspondence with the observed temperature and precipitation. Analysis of the data revealed that the MB is sensitive to temperature at the rate of −0.65 m w.e. a−1 °C−1, but the sensitivity is only 0.13 m w.e. a−1 for a 10% change in precipitation. Furthermore, there is a good correlation between the simulated glacier mass loss and depleting glacier-melt in the autumn streamflow. The mass loss of the KG is expected to exacerbate in future as a result of the projected climate change, and thus further diminishing the streamflow of the transboundary rivers emanating from the region. It is envisaged that the findings will contribute to a better understanding of glacier-climate interactions and their implications for future water availability in the Kashmir Himalaya.

Energy and glacier mass balance of Fürkeleferner, Italy: past, present, and future

The energy and mass balance of mountain glaciers translate into volume changes that play out as area changes over time. From this, together with former moraines during maximum advances, information on past climate conditions and the climatic drivers behind during glacier advances can be obtained. Here, we use the distributed COupled Snowpack and Ice surface energy and mass balance model in PYthon (COSIPY) to simulate the present state of an Italian glacier, named Fürkeleferner, for the mass balance years 2013–2017. Next, we investigate the local climate during the time of the last “Little Ice Age” (LIA) maximum glacier advance using COSIPY together with the LIA glacier outline retrieved from moraine mapping and a digital elevation model (DEM) adapted for the glacier’s geometry at the time of the LIA as a benchmark. Furthermore, the glacier’s sensitivity to future air temperature increase of +1 K and +2 K is investigated using the same model. For all simulations, meteorological data of closely located climate stations are used to force the model. We show the individual monthly contribution of individual energy and mass balance components. Refreezing during the summer months is an important component of the energy and mass balance, on average about 9 % relative to total annual ablation. The results from simulating past climate show a 2.8 times larger glacier area for Fürkeleferner during the LIA than today. This further implies a 2.5 K colder climate, assuming that the amount of precipitation was 10 %–20 % in excess of today’s value. Concerning further temperature increase of 2 K, the glacier would only consist of the ablation area implying sustained mass loss and eventual total mass loss. Even under current climatic conditions, the glacier area would have to decrease to 17 % of its current area to be in a steady state. We discuss the reliability of the results by comparing simulated present mass balance to measured mass balances of neighboring glaciers in the European Alps and with short-term measurements on Fürkeleferner itself. In conclusion, we are able to show how the glacier responds to past and future climate change and determine the climatic drivers behind.

New modelling approach to optimize rainwater harvesting system for non-potable uses and groundwater recharge: A case study from Israel

• A rainwater harvesting system from a building in Israel was used as a case study. • Precipitation, runoff coefficient and non-potable consumption were measured. • Measured data used to model and assess the efficiency of the system. • Results are stochastic, allowing to assess variability of the system performance. Countries with dry climates must optimize the use of natural water sources. Rainwater Harvesting (RWH) is a viable water source for non-potable uses. Besides, rainwater can also be used for managed aquifer recharge (MAR). This work discusses the combined use of RWH for toilet flushing and MAR. A mass-balance simulation was used to estimate efficiency of the combined system in a building in Israel equipped with a RWH system used for MAR only. This enabled the use of actual water consumption patterns and the development of a runoff function using rainfall and runoff measurements. By the novel use and the integration of measured dynamic data into the model, the results are stochastic and the variability and uncertainty of the performance of the RWH system could be assessed. Non-potable water savings were ∼20% for the expected cost-effective tank size, and in this case ∼40% of the rainwater was used for MAR. The stochastic approach of the combined system proved the importance of the non-potable consumption value as well as precipitation series. Finally, it is shown that the potable water savings from using rainwater for toilet flushing would represent an annual CO 2 emissions reduction of 267 kg.

Modelling glacier mass balance and climate sensitivity in the context of sparse observations: application to Saskatchewan Glacier, western Canada

Abstract. Glacier mass balance models are needed at sites with scarce long-term observations to reconstruct past glacier mass balance and assess its sensitivity to future climate change. In this study, North American Regional Reanalysis (NARR) data were used to force a physically based, distributed glacier mass balance model of Saskatchewan Glacier for the historical period 1979–2016 and assess its sensitivity to climate change. A 2-year record (2014–2016) from an on-glacier automatic weather station (AWS) and historical precipitation records from nearby permanent weather stations were used to downscale air temperature, relative humidity, wind speed, incoming solar radiation and precipitation from the NARR to the station sites. The model was run with fixed (1979, 2010) and time-varying (dynamic) geometry using a multitemporal digital elevation model dataset. The model showed a good performance against recent (2012–2016) direct glaciological mass balance observations as well as with cumulative geodetic mass balance estimates. The simulated mass balance was not very sensitive to the NARR spatial interpolation method, as long as station data were used for bias correction. The simulated mass balance was however sensitive to the biases in NARR precipitation and air temperature, as well as to the prescribed precipitation lapse rate and ice aerodynamic roughness lengths, showing the importance of constraining these two parameters with ancillary data. The glacier-wide simulated energy balance regime showed a large contribution (57 %) of turbulent (sensible and latent) heat fluxes to melting in summer, higher than typical mid-latitude glaciers in continental climates, which reflects the local humid “icefield weather” of the Columbia Icefield. The static mass balance sensitivity to climate was assessed for prescribed changes in regional mean air temperature between 0 and 7 ∘C and precipitation between −20 % and +20 %, which comprise the spread of ensemble Representative Concentration Pathway (RCP) climate scenarios for the mid (2041–2070) and late (2071–2100) 21st century. The climate sensitivity experiments showed that future changes in precipitation would have a small impact on glacier mass balance, while the temperature sensitivity increases with warming, from −0.65 to −0.93 m w.e. a−1 ∘C−1. The mass balance response to warming was driven by a positive albedo feedback (44 %), followed by direct atmospheric warming impacts (24 %), a positive air humidity feedback (22 %) and a positive precipitation phase feedback (10 %). Our study underlines the key role of albedo and air humidity in modulating the response of winter-accumulation type mountain glaciers and upland icefield-outlet glacier settings to climate.

High- and low-Mg adakitic rocks in southern Tibet: Implication for the crustal thickening and geodynamic process in the late Cretaceous

The history of crustal thickening is critical for comprehending the evolution process of the Tibetan plateau. The genetic relationship and spatiotemporal distribution of adakitic rocks with high Sr/Y, (La/Yb) N ratios may impose fundamental constraints on specific geodynamic and crustal thickening processes. The Songka intrusive suite consists of adakitic rocks with high-Mg adakite (monzonites), low-Mg adakitic rocks (quartz monzonites), and gabbroic-diorites. Their geochemical characteristics, petrogenesis, and geodynamic implication is discussed in this paper. Songka high-Mg adakites formed when subducting Neo-Tethyan oceanic slab melts were assimilated by mantle peridotite at elevated temperatures (average Ti-in-zircon temperatures of sample SK1319 are 830 °C). Low-Mg adakitic rocks exhibit markedly different geochemical characteristics than those formed in the juvenile thickened lower crust. The least-squares mass balance and trace element simulations indicate that Songka low-Mg adakitic samples were formed at a degree of 60–65% fractional crystallization of melts resembling Songka high-Mg adakites. Songka gabbroic diorites are formed in the metasomatic mantle wedge. The genetic relationship between these rocks indicates that the crust of the southern Tibet did not appear to thicken obviously during the early stage of late Cretaceous (100–90 Ma). The intensive emplacement of late Cretaceous high-Mg adakites (Mg # > 50) and mafic rocks (SiO 2 < 52%) in the southern Lhasa terrane around 95–90 Ma, together with the slow convergence velocity of the Indian-Asian continent, can be explained by the response of rapid crustal thickening and Neo-Tethyan slab rollback at 95–90 Ma. • High-Mg adakitic rocks were formed from partial melting of Neo-Tethyan slab. • Coeval low-Mg adakitic rocks were produced by fractional crystallization. • Rapid crustal growth and Neo-Tethyan slab rollback occurred at 95–90 Ma.

Circular economy and life cycle assessment of alumina production: Simulation-based comparison of Pedersen and Bayer processes

The excessive production of bauxite residue (red mud) in the Bayer process is one of the major challenges amongst alumina producers. The Pedersen process is known as a combination of smelting reduction of bauxite and leaching treatment of the produced slag for alumina production, and the process also produces an inert bauxite residue (grey mud), which is more suitable for further processing compared to the red mud. This paper describes the life cycle assessment for two alumina production processes, the current commercialized Bayer process and the alternative Pedersen process, based on the complete mass and energy balance simulation. Process simulation and environmental software are applied to quantify the metallurgical processes and their related environmental effects. By assuming the same treatments for the bauxite residues, it is found that, with the production of pig iron, the Pedersen process shows the best performance in mineral resource scarcity but requires a higher energy demand and generates twice the amount of bauxite residue compared to the Bayer process. The final, relative performance in climate change impact depends on the carbon intensity of energy and the bauxite ore composition.

The internal melting of landfast sea ice in Prydz Bay, East Antarctica

Summertime internal melting of Antarctic sea ice is common due to the penetration of solar radiation below the snow and ice surface. We focus on the role of internal melting and heat conduction in generating gap layers within the ice. These often occur approximately 0.1 m below the ice surface. In a small-scale survey over land-fast sea ice in Prydz Bay, East Antarctica, we observed, for the first time, gap layers 0.6–1.0 m below the surface for both first-year ice and multi-year ice. A 1D snow/ice thermodynamic model successfully simulated snow and ice mass balance and the evolution of the gap layers. Their spatial distribution was largely controlled by snow thickness and ice thickness. A C-shaped ice temperature profile with the lowest values in the middle of the ice layer resulted in heat flux convergence causing downward progression of the internal melt layer. Multidecadal (1979–2019) seasonal simulations showed decreasing air temperature favored a postposed internal melting onset, reduced total internal melt, and delayed potential ice breakup, which indicated a higher chance for local coastal ice to be shifted from first-year ice to multi-year ice.

Albedo Parametrizations for the Laohugou Glacier No.12 in the Qilian Mountains—Previous Models and an Alternative Approach

Accurate estimates of albedo can be crucial for energy balance models of glaciers. A number of algorithms exist which are often site dependent and rely on accurate measurements or estimates of snow depth. Using the well-established COSIMA model we simulate the energy and mass balance of the Laohugou Glacier No.12 in the Qilian Mountains, on the northern fringe of the Qinghai-Tibetan Plateau, a glacier that has been well studied in the past. Using energy flux and mass balance measurements between 2010 and 2015 we were able to validate the model over multiple seasons. Using the original albedo parametrization, the model fails to reproduce the observed mass balance. We show that this is due to the failure to estimate snow depth accurately. We therefore applied two alternative albedo algorithms, one well established example and one new parametrization only dependent on temperature and time since last snow fall. As a result, mass balance simulations improve considerably from a RMSE of 0.53 m w.e. for the original parametrization to 0.39 and 0.19 m w.e. for the uncalibrated established and the new calibrated model respectively. Modelled albedo during the ablation period (NSE = 0.05, R2 = 0.33) is more accurate than during the accumulation period (NSE = −0.37, R2 = 0.04). Testing the new model at another glacier on the Tibetan Plateau shows that a local recalibration of parameters remains necessary to achieve satisfying results. Investigations into the effect of impurities in snow, regional moisture sources and changing surface characteristics with rising temperatures will be crucial for accurate projections into the future.

Changes in the Surface Elevation of the Laohugou Glacier No. 12 in Western Qilian Mountains

As the largest valley glacier in the Qilian Mountains, the Laohugou glacier No. 12 (LHG12) has shrunk significantly since 1957. In this study, two topographic maps and a WorldView-2 satellite stereopair image data were used to assess the volume and cumulative mass balance of LHG12 located at the western Qilian Mountains during 1957–2015. During the study period, the LHG12 exhibited changes in two processes: slightly ablation and stability in a brief period during 1957–1989 and strong melting and accelerated ablation during 1989–2015. During 1957–2015, the volume of LHG12 decreased by 0.38 km3, the average thickness decreased by 17.23 m, the cumulative mass balance (MB) was −14.69 ± 3.00 m w. e., and ablation was found glacier-wide. By comparing the previous MB simulation and digital elevation model (DEM) differencing results, it was found that the MB simulation results underestimated the strong melting trend of LHG12 since the 1990s. Temperature rose, especially in autumn and winter, and could cause the ice temperature of LHG12 to increase, and LHG12 may become more sensitive to climate change.

Dynamic energy and mass balance model for an industrial alkaline water electrolyzer plant process

This paper proposes a parameter adjustable dynamic mass and energy balance simulation model for an industrial alkaline water electrolyzer plant that enables cost and energy efficiency optimization by means of system dimensioning and control. Thus, the simulation model is based on mathematical models and white box coding, and it uses a practicable number of fixed parameters. Zero-dimensional energy and mass balances of each unit operation of a 3 MW, and 16 bar plant process were solved in MATLAB functions connected via a Simulink environment. Verification of the model was accomplished using an analogous industrial plant of the same power and pressure range having the same operational systems design. The electrochemical, mass flow and thermal behavior of the simulation and the industrial plant were compared to ascertain the accuracy of the model and to enable modification and detailed representation of real case scenarios so that the model is suitable for use in future plant optimization studies. The thermal model dynamically predicted the real case with 98.7 % accuracy. Shunt currents were the main contributor to relative low Faraday efficiency of 86 % at nominal load and steady-state operation and heat loss to ambient from stack was only 2.6 % of the total power loss.

The surface energy balance of Austre Lovénbreen, Svalbard, during the ablation period in 2014

The ability to simulate the surface energy balance is key to studying land–atmosphere interactions; however, it remains a weakness in Arctic polar sciences. Based on the analysis of meteorological data from 1 June to 30 September 2014 from an automatic weather station on the glacier Austre Lovénbreen, near Ny–Ålesund, Svalbard, we established a surface energy balance model to simulate surface melt. The results reveal that the net shortwave radiation accounts for 87% (39 W m–2) of the energy sources, and is controlled by cloud cover and surface albedo. The sensible heat equals 6 W m–2 and is a continuous energy source at the glacier surface. Net longwave radiation and latent heat account for 31% and 5% of heat sinks, respectively. The simulated summer mass balance equals –793 mm w.e., agreeing well with the observation by an ultrasonic ranger.

On the importance of subsurface heat flux for estimating the mass balance of alpine glaciers

Glaciers, as massive freshwater reservoirs, support the planet's living systems and have an impact on our daily lives, even for communities living far away. Ongoing and future climate change is predicted to have strong impacts on the mass balance of alpine glacier around the world. To understand the relationship between climate and glacier dynamics, a range of mass balance models are currently used. Most of these models however, ignore subsurface heat fluxes as a component of glacier mass balance. Here, we set out to investigate the importance of subsurface heat flux for the mass balance of an alpine glacier using a surface energy mass balance model (SEM) coupled with a multilayer subsurface heat conduction model (MSHCM) that resolves the subsurface glacier temperature. As a case study, we investigate the Urumqi Glacier No.1 in the Tianshan Mountains (NW China), which has a long and continuous time series of surface and subsurface glacier temperature measurements. We evaluate the results of both glacier temperature models (SEM and MSHCM) using these in situ observations and investigate the sensitivity of mass balance to five meteorological factors: air temperature, precipitation, incoming shortwave radiation, relative humidity, and wind speed. The mass balance of the glacier was simulated first by including the influence of subsurface heat flux, and second, the subsurface heat flux was neglected. Observed and simulated mass balance and the englacial temperature were found to be reasonably close in both cases. Furthermore, the mass balance was simulated with a zero surface temperature assumption, which resulted in a 6% overestimation of the summer ablation. We concluded that the mass balance of Urumqi Glacier No.1 was most sensitive to variations in temperature, followed by precipitation. Furthermore, our results show that subsurface heat flux in the ablation area can generally be neglected in estimating the mass balance of alpine glaciers during ablation season.