The Greenland Ice Sheet's surface mass balance in a seasonally sea-ice free Arctic

General circulation models predict a rapid decrease in sea ice extent with concurrent increases in near surface air temperature and precipitation in the Arctic over the 21st century. This has led to suggestions that some Arctic land ice masses may experience an increase in accumulation due to enhanced evaporation from a seasonally sea ice free Arctic Ocean. To investigate the impact of this phenomenon on Greenland ice sheet climate and surface mass balance (SMB) a regional climate model, HadRM3, was used to force an insolation-temperature melt SMB model. A set of experiments designed to investigate the role of sea ice independently from sea surface temperature (SST) forcing are described. In the warmer and wetter SI + SST simulation Greenland experiences a 23% increase in winter SMB but 65% reduced summer SMB, resulting in a net decrease in the annual value. This study shows that sea ice decline contributes to the increased winter balance, causing 25% of the increase in winter accumulation; this is largest in eastern Greenland as the result of increased evaporation in the Greenland Sea. These results indicate that the seasonal cycle of Greenland's SMB will increase dramatically as global temperatures increase, with the largest changes in temperature and precipitation occurring in winter. This demonstrates that the accurate prediction of changes in sea ice cover is important for predicting Greenland SMB and ice sheet evolution.

Diabatic balance model for the equatorial atmosphere

Using an asymptotic expansion, a balance model is derived for the shallow-water equations (SWE) on the equatorial beta-plane that is valid for planetary-scale equatorial dynamics and includes Kelvin waves. In contrast to many theories of tropical dynamics, neither a strict balance between diabatic heating and vertical motion nor a small Froude number is required. Instead, the expansion is based on the smallness of the ratio of meridional to zonal length scales, which can also be interpreted as a separation in time scale. The leading-order model is characterized by a semigeostrophic balance between the zonal wind and meridional pressure gradient, while the meridional wind v vanishes; the model is thus asymptotically nondivergent, and the nonzero correction to v can be found at the next order. Importantly for applications, the diagnostic balance relations are linear for winds when inferring the wind field from mass observations and the winds can be diagnosed without direct observations of diabatic heating. The accuracy of the model is investigated through a set of numerical examples. These examples show that the diagnostic balance relations can remain valid even when the dynamics do not, and the balance dynamics can capture the slow behavior of a rapidly varying solution.

Effect of uncertainty in surface mass balance--elevation feedback on projections of the future sea level contribution of the Greenland ice sheet

We apply a new parameterisation of the Greenland ice sheet (GrIS) feedback between surface mass balance (SMB: the sum of surface accumulation and surface ablation) and surface elevation in the MAR regional climate model (Edwards et al., 2014) to projections of future climate change using five ice sheet models (ISMs). The MAR (Modèle Atmosphérique Régional: Fettweis, 2007) climate projections are for 2000–2199, forced by the ECHAM5 and HadCM3 global climate models (GCMs) under the SRES A1B emissions scenario. The additional sea level contribution due to the SMB– elevation feedback averaged over five ISM projections for ECHAM5 and three for HadCM3 is 4.3% (best estimate; 95% credibility interval 1.8–6.9 %) at 2100, and 9.6% (best estimate; 95% credibility interval 3.6–16.0 %) at 2200. In all results the elevation feedback is significantly positive, amplifying the GrIS sea level contribution relative to the MAR projections in which the ice sheet topography is fixed: the lower bounds of our 95% credibility intervals (CIs) for sea level contributions are larger than the “no feedback” case for all ISMs and GCMs. Our method is novel in sea level projections because we propagate three types of modelling uncertainty – GCM and ISM structural uncertainties, and elevation feedback parameterisation uncertainty – along the causal chain, from SRES scenario to sea level, within a coherent experimental design and statistical framework. The relative contributions to uncertainty depend on the timescale of interest. At 2100, the GCM uncertainty is largest, but by 2200 both the ISM and parameterisation uncertainties are larger. We also perform a perturbed parameter ensemble with one ISM to estimate the shape of the projected sea level probability distribution; our results indicate that the probability density is slightly skewed towards higher sea level contributions.

Mixing-cell boundary conditions and apparent mass balance errors for advective-dispersive solute transport

In many field or laboratory situations, well-mixed reservoirs like, for instance, injection or detection wells and gas distribution or sampling chambers define boundaries of transport domains. Exchange of solutes or gases across such boundaries can occur through advective or diffusive processes. First we analyzed situations, where the inlet region consists of a well-mixed reservoir, in a systematic way by interpreting them in terms of injection type. Second, we discussed the mass balance errors that seem to appear in case of resident injections. Mixing cells (MC) can be coupled mathematically in different ways to a domain where advective-dispersive transport occurs: by assuming a continuous solute flux at the interface (flux injection, MC-FI), or by assuming a continuous resident concentration (resident injection). In the latter case, the flux leaving the mixing cell can be defined in two ways: either as the value when the interface is approached from the mixing-cell side (MC-RT -), or as the value when it is approached from the column side (MC-RT +). Solutions of these injection types with constant or-in one case-distance-dependent transport parameters were compared to each other as well as to a solution of a two-layer system, where the first layer was characterized by a large dispersion coefficient. These solutions differ mainly at small Peclet numbers. For most real situations, the model for resident injection MC-RI + is considered to be relevant. This type of injection was modeled with a constant or with an exponentially varying dispersion coefficient within the porous medium. A constant dispersion coefficient will be appropriate for gases because of the Eulerian nature of the usually dominating gaseous diffusion coefficient, whereas the asymptotically growing dispersion coefficient will be more appropriate for solutes due to the Lagrangian nature of mechanical dispersion, which evolves only with the fluid flow. Assuming a continuous resident concentration at the interface between a mixing cell and a column, as in case of the MC-RI + model, entails a flux discontinuity. This flux discontinuity arises inherently from the definition of a mixing cell: the mixing process is included in the balance equation, but does not appear in the description of the flux through the mixing cell. There, only convection appears because of the homogeneous concentration within the mixing cell. Thus, the solute flux through a mixing cell in close contact with a transport domain is generally underestimated. This leads to (apparent) mass balance errors, which are often reported for similar situations and erroneously used to judge the validity of such models. Finally, the mixing cell model MC-RI + defines a universal basis regarding the type of solute injection at a boundary. Depending on the mixing cell parameters, it represents, in its limits, flux as well as resident injections. (C) 1998 Elsevier Science B.V. All rights reserved.

Table 2: Firn temperature, ram resistance, mass balance and submergence velo city measurements on Colle Gnifetti

Distributional patterns of glaciological parameters at the Colle Gnifetti core drilling site are described and their interrelationships are brietly discussed. Observations within a stake network established in 1980 furnish information about snow accumulation (short term balance), submergence velocity of ice tlow (long term balance), ram hardness (melt layer stratigraphy), and firn temperature. In addition, a numerical model was used to estimate local variations of available radiant energy. Melt layer formation is considerably more intensive on the south facing parts of the firn saddie where incoming radiation is high. These melt layers seem to effectively protect some of the fallen snow from wind erosion. As a result, balance ist up to one order of magnitude larger on south facing slopes. Heat applied to the surface is therefore positively correlated with balance, whereas the relation between solar radiation and firn temperature is less dear. Distributional patterns of submergence velocity confirm that the observed spatial variability of surface balance is representative for longer time periods and greatly intluences the time scale and the stratigraphy of firn and ice cores from Colle Gnifetti.

An enhanced mill extraction model

An enhanced mill extraction model has been developed to calculate mill performance parameters and to predict the extraction performance of a milling unit. The model takes into account the fibre suspended in juice streams and calculates filling ratio, reabsorption factor, imbibition coefficient, and separation efficiency using more complete definitions than those used in previous extraction models. A mass balance model is used to determine the fibre, brix and moisture mass flows between milling units so that a complete milling train, including the return stream from the juice screen, is modelled. Model solutions are presented to determine the effect of different levels of fibre in juice and efficiency of fibre separation in the juice screen on brix extraction. The model provides more accurate results than earlier models leading to better understanding and improvement of the milling process.

A novel population balance model for the dilute acid hydrolysis of hemicellulose

Acid hydrolysis is a popular pretreatment for removing hemicellulose from lignocelluloses in order to produce a digestible substrate for enzymatic saccharification. In this work, a novel model for the dilute acid hydrolysis of hemicellulose within sugarcane bagasse is presented and calibrated against experimental oligomer profiles. The efficacy of mathematical models as hydrolysis yield predictors and as vehicles for investigating the mechanisms of acid hydrolysis is also examined. Experimental xylose, oligomer (degree of polymerisation 2 to 6) and furfural yield profiles were obtained for bagasse under dilute acid hydrolysis conditions at temperatures ranging from 110C to 170C. Population balance kinetics, diffusion and porosity evolution were incorporated into a mathematical model of the acid hydrolysis of sugarcane bagasse. This model was able to produce a good fit to experimental xylose yield data with only three unknown kinetic parameters ka, kb and kd. However, fitting this same model to an expanded data set of oligomeric and furfural yield profiles did not successfully reproduce the experimental results. It was found that a ``hard-to-hydrolyse'' parameter, $\alpha$, was required in the model to ensure reproducibility of the experimental oligomer profiles at 110C, 125C and 140C. The parameters obtained through the fitting exercises at lower temperatures were able to be used to predict the oligomer profiles at 155C and 170C with promising results. The interpretation of kinetic parameters obtained by fitting a model to only a single set of data may be ambiguous. Although these parameters may correctly reproduce the data, they may not be indicative of the actual rate parameters, unless some care has been taken to ensure that the model describes the true mechanisms of acid hydrolysis. It is possible to challenge the robustness of the model by expanding the experimental data set and hence limiting the parameter space for the fitting parameters. The novel combination of ``hard-to-hydrolyse'' and population balance dynamics in the model presented here appears to stand up to such rigorous fitting constraints.

Regolith mass balance inferred from combined mineralogical, geochemical and geophysical studies: Mule Hole gneissic watershed, South India

The aim of this study is to propose a method to assess the long-term chemical weathering mass balance for a regolith developed on a heterogeneous silicate substratum at the small experimental watershed scale by adopting a combined approach of geophysics, geochemistry and mineralogy. We initiated in 2003 a study of the steep climatic gradient and associated geomorphologic features of the edge of the rifted continental passive margin of the Karnataka Plateau, Peninsular India. In the transition sub-humid zone of this climatic gradient we have studied the pristine forested small watershed of Mule Hole (4.3 km(2)) mainly developed on gneissic substratum. Mineralogical, geochemical and geophysical investigations were carried out (i) in characteristic red soil profiles and (ii) in boreholes up to 60 m deep in order to take into account the effect of the weathering mantle roots. In addition, 12 Electrical Resistivity Tomography profiles (ERT), with an investigation depth of 30 m, were generated at the watershed scale to spatially characterize the information gathered in boreholes and soil profiles. The location of the ERT profiles is based on a previous electromagnetic survey, with an investigation depth of about 6 m. The soil cover thickness was inferred from the electromagnetic survey combined with a geological/pedological survey. Taking into account the parent rock heterogeneity, the degree of weathering of each of the regolith samples has been defined using both the mineralogical composition and the geochemical indices (Loss on Ignition, Weathering Index of Parker, Chemical Index of Alteration). Comparing these indices with electrical resistivity logs, it has been found that a value of 400 Ohm m delineates clearly the parent rocks and the weathered materials, Then the 12 inverted ERT profiles were constrained with this value after verifying the uncertainty due to the inversion procedure. Synthetic models based on the field data were used for this purpose. The estimated average regolith thickness at the watershed scale is 17.2 m, including 15.2 m of saprolite and 2 m of soil cover. Finally, using these estimations of the thicknesses, the long-term mass balance is calculated for the average gneiss-derived saprolite and red soil. In the saprolite, the open-system mass-transport function T indicates that all the major elements except Ca are depleted. The chlorite and biotite crystals, the chief sources for Mg (95%), Fe (84%), Mn (86%) and K (57%, biotite only), are the first to undergo weathering and the oligoclase crystals are relatively intact within the saprolite with a loss of only 18%. The Ca accumulation can be attributed to the precipitation of CaCO3 from the percolating solution due to the current and/or the paleoclimatic conditions. Overall, the most important losses occur for Si, Mg and Na with -286 x 10(6) mol/ha (62% of the total mass loss), -67 x 10(6) mol/ha (15% of the total mass loss) and -39 x 10(6) mol/ha (9% of the total mass loss), respectively. Al, Fe and K account for 7%, 4% and 3% of the total mass loss, respectively. In the red soil profiles, the open-system mass-transport functions point out that all major elements except Mn are depleted. Most of the oligoclase crystals have broken down with a loss of 90%. The most important losses occur for Si, Na and Mg with -55 x 10(6) mol/ha (47% of the total mass loss), -22 x 10(6) mol/ha (19% of the total mass loss) and -16 x 10(6) mol/ha (14% of the total mass loss), respectively. Ca, Al, K and Fe account for 8%, 6%, 4% and 2% of the total mass loss, respectively. Overall these findings confirm the immaturity of the saprolite at the watershed scale. The soil profiles are more evolved than saprolite but still contain primary minerals that can further undergo weathering and hence consume atmospheric CO2.

An overview of land application of pig effluent-P using soil P chemistry and mass balance calculations.

Attention is directed at land application of piggery effluent (containing urine, faeces, water, and wasted feed) as a potential source of water resource contamination with phosphorus (P). This paper summarises P-related properties of soil from 0-0.05 m depth at 11 piggery effluent application sites, in order to explore the impact that effluent application has had on the potential for run-off transport of P. The sites investigated were situated on Alfisol, Mollisol, Vertisol, and Spodosol soils in areas that received effluent for 1.5-30 years (estimated effluent-P applications of 100-310000 kg P/ha in total). Total (PT), bicarbonate extractable (PB), and soluble P forms were determined for the soil (0-0.05 m) at paired effluent and no-effluent sites, as well as texture, oxalate-extractable Fe and Al, organic carbon, and pH. All forms of soil P at 0-0.05 m depth increased with effluent application (PB at effluent sites was 1.7-15 times that at no-effluent sites) at 10 of the 11 sites. Increases in PB were strongly related to net P applications (regression analysis of log values for 7 sites with complete data sets: 82.6 % of variance accounted for, p <0.01). Effluent irrigation tended to increase the proportion of soil PT in dilute CaCl2-extractable forms (PTC: effluent average 2.0 %; no-effluent average 0.6%). The proportion of PTC in non-molybdate reactive forms (centrifuged supernatant) decreased (no-effluent average, 46.4 %; effluent average, 13.7 %). Anaerobic lagoon effluent did not reliably acidify soil, since no consistent relationship was observed for pH with effluent application. Soil organic carbon was increased in most of the effluent areas relative to the no-effluent areas. The four effluent areas where organic carbon was reduced had undergone intensive cultivation and cropping. Current effluent management at many of the piggeries failed to maximise the potential for waste P recapture. Ten of the case-study effluent application areas have received effluent-P in excess of crop uptake. While this may not represent a significant risk of leaching where sorption retains P, it has increased the risk of transport of P by run-off. Where such sites are close to surface water, run-off P loads should be managed.

Nutrient mass balance for Lockyer vegetable production

The potential for fertiliser use in Lockyer Valleys intensive vegetable production to impact on the Moreton Bay Waterways (MBW) is not well defined. Notwithstanding nutrient runoff through soil erosion of agricultural lands has been identified as a process that significantly contributes artificial fertiliser to the MBW (SEQ Healthy Waterways Draft Strategy 2006). In order to better understand this issue the present study undertakes a nutrient mass balance to evaluate nitrogen use efficiency in the intensive horticultural industry of the Lockyer Valley.

Hanging drop: An in vitro air toxic exposure model using human lung cells in 2D and 3D structures

Using benzene as a candidate air toxicant and A549 cells as an in vitro cell model, we have developed and validated a hanging drop (HD) air exposure system that mimics an air liquid interface exposure to the lung for periods of 1 h to over 20 days. Dose response curves were highly reproducible for 2D cultures but more variable for 3D cultures. By comparing the HD exposure method with other classically used air exposure systems, we found that the HD exposure method is more sensitive, more reliable and cheaper to run than medium diffusion methods and the CULTEX (R) system. The concentration causing 50% of reduction of cell viability (EC50) for benzene, toluene, p-xylene, m-xylene and o-xylene to A549 cells for 1 h exposure in the HD system were similar to previous in vitro static air exposure. Not only cell viability could be assessed but also sub lethal biological endpoints such as DNA damage and interleukin expressions. An advantage of the HD exposure system is that bioavailability and cell concentrations can be derived from published physicochemical properties using a four compartment mass balance model. The modelled cellular effect concentrations EC50(cell) for 1 h exposure were very similar for benzene, toluene and three xylenes and ranged from 5 to 15 mmol/kg(dry weight) which corresponds to the intracellular concentration of narcotic chemicals in many aquatic species, confirming the high sensitivity of this exposure method. (C) 2013 Elsevier B.V. All rights reserved.

Development of mass balance models for lakes

Yhteenveto: Järvien ainetasemallien kehittäminen.

A spring-mass-lever model, stiffness and inertia maps for single-input, single-output compliant mechanisms

A spring-mass-lever (SML) model is introduced in this paper for a single-input-single-output compliant mechanism to capture its static and dynamic behavior. The SML model is a reduced-order model, and its five parameters provide physical insight and quantify the stiffness and inertia(1) at the input and output ports as well as the transformation of force and displacement between the input and output. The model parameters can be determined with reasonable accuracy without performing dynamic or modal analysis. The paper describes two uses of the SML model: computationally efficient analysis of a system of which the compliant mechanism is a part; and design of compliant mechanisms for the given user-specifications. During design, the SML model enables determining the feasible parameter space of user-specified requirements, assessing the suitability of a compliant mechanism to meet the user-specifications and also selecting and/or re-designing compliant mechanisms from an existing database. Manufacturing constraints, material choice, and other practical considerations are incorporated into this methodology. A micromachined accelerometer and a valve mechanism are used as examples to show the effectiveness of the SML model in analysis and design. (C) 2012 Published by Elsevier Ltd.

Variation of Snowline and Mass Balance of Glaciers of Warwan and Bhut Basins of Western Himalaya Using Remote Sensing Technique

Glaciers are natural reservoirs of fresh water in frozen state and sensitive indicators of climate change. Among all the mountainous glaciated regions, glaciers of Himalayas form one of the largest concentrations of ice outside the Polar Regions. Almost all the major rivers of northern India originate from these glaciers and sustain perennial flow. Therefore, in view of the importance and role of the glaciers in sustaining the life on the Earth, monitoring the health of glaciers is necessary. Glacier's health is monitored in two ways (i) by mapping the change in extent of glaciers (ii) by finding variation in the annual mass balance. This paper has been discussed the later approach for monitoring the health of glaciers of Warwan and Bhut basins. Mass balance of glaciers of these two basins was determined based on the extraction of snow line at the end of ablation season. A series of satellite images of AWiFS sensor were analysed for extraction of snowline on the glaciers for the period of 2005, 2006 and 2007. The snow line at the end of ablation season is used to compute accumulation area ratio (AAR = Accumulation area/Glacier area) for each glacier of basins. An approach based on relationship of AAR to specific mass balance (computed in field) for glaciers of Basapa basin was employed in the present study. Mean of specific mass balance of individual glacier for the year 2005, 2006 and 2007 of Warwan basin was found to be -ve 0.19 m, -ve 0.27 m and -ve 0.2 m respectively. It is 0.05 m, -ve 0.11 m and -ve 0.19 m for Bhut basin. The analysis suggests a loss of 4.3 and 0.83 kmA(3) of glacier in the monitoring period of 3 years for Warwan and Bhut basins respectively. The overall results suggest that the glaciers of Warwan basin and Bhut basins have suffered more loss of ice than gain in the monitoring period of 3 years.