Preliminary spectral characterization of snow in a high altitude tropical glacier and potential effects of impurities in snow on albedo of tropical glaciers

This paper presents some preliminary measurements of snow spectral reflectance on the tropical Bolivian Zongo glacier. Measurements show a correct agreement with theoretical spectral albedo of pure snow in the near infrared region, but lower values in the visible region (by 10–20%) probably due to aerosols contained in snow. Impurity contents ranged from 10 to 100 ppmw in one‐week‐old snow collected from the Zongo glacier, but measurements are scarce. Large amounts of snowfall partly compensate the proximity of dust sources in mid‐latitude glaciers, whereas on outer‐tropical glaciers precipitations are not abundant and are very seasonal, and sources of aerosols are proximate. Copyright © 2007 John Wiley & Sons, Ltd.     

Using a H∞ filter assimilation procedure to estimate root zone soil water content

Root zone soil water content impacts plant water availability, land energy and water balances. Because of unknown hydrological model error, observation errors and the statistical characteristics of the errors, the widely used Kalman filter (KF) and its extensions are challenged to retrieve the root zone soil water content using the surface soil water content. If the soil hydraulic parameters are poorly estimated, the KF and its extensions fail to accurately estimate the root zone soil water. The H‐infinity filter (HF) represents a robust version of the KF. The HF is widely used in data assimilation and is superior to the KF, especially when the performance of the model is not well understood. The objective of this study is to study the impact of uncertain soil hydraulic parameters, initial soil moisture content and observation period on the ability of HF assimilation to predict in situ soil water content. In this article, we study seven cases. The results show that the soil hydraulic parameters hold a critical role in the course of assimilation. When the soil hydraulic parameters are poorly estimated, an accurate estimation of root soil water content cannot be retrieved by the HF assimilation approach. When the estimated soil hydraulic parameters are similar to actual values, the soil water content at various depths can be accurately retrieved by the HF assimilation. The HF assimilation is not very sensitive to the initial soil water content, and the impact of the initial soil water content on the assimilation scheme can be eliminated after about 5–7 days. The observation interval is important for soil water profile distribution retrieval with the HF, and the shorter the observation interval, the shorter the time required to achieve actual soil water content. However, the retrieval results are not very accurate at a depth of 100 cm. Also it is complex to determine the weighting coefficient and the error attenuation parameter in the HF assimilation. In this article, the trial‐and‐error method was used to determine the weighting coefficient and the error attenuation parameter. After the first establishment of limited range of the parameters, ‘the best parameter set’ was selected from the range of values. For the soil conditions investigated, the HF assimilation results are better than the open‐loop results. Copyright © 2010 John Wiley & Sons, Ltd.     

Towards a comprehensive approach to parameter estimation in land surface parameterization schemes

In climate models, the land–atmosphere interactions are described numerically by land surface parameterization (LSP) schemes. The continuing improvement in realism in these schemes comes at the expense of the need to specify a large number of parameters that are either directly measured or estimated. Also, an emerging problem is whether the relationships used in LSPs are universal and globally applicable. One plausible approach to evaluate this is to first minimize uncertainty in model parameters by calibration. In this paper, we conduct a comprehensive analysis of some model diagnostics using a slightly modified version of the Simple Biosphere 3 model for a variety of biomes located mainly in the Amazon. First, the degree of influence of each individual parameter in simulating surface fluxes is identified. Next, we estimate parameters using a multi‐operator genetic algorithm applied in a multi‐objective context and evaluate simulations of energy and carbon fluxes against observations. Compared with the default parameter sets, these parameter estimates improve the partitioning of energy fluxes in forest and cropland sites and provide better simulations of daytime increases in assimilation of net carbon during the dry season at forest sites. Finally, a detailed assessment of the parameter estimation problem was performed by accounting for the decomposition of the mean squared error to the total model uncertainty. Analysis of the total prediction uncertainty reveals that the parameter adjustments significantly improve reproduction of the mean and variability of the flux time series at all sites and generally remove seasonality of the errors but do not improve dynamical properties. Our results demonstrate that error decomposition provides a meaningful and intuitive way to understand differences in model performance. To make further advancements in the knowledge of these models, we encourage the LSP community to adopt similar approaches in the future. Copyright © 2012 John Wiley & Sons, Ltd.     

Impact of grid resolution on the integrated and distributed response of a coupled surface–subsurface hydrological model for the des Anglais catchment,Quebec

Digital elevation models (DEMs) at different resolutions (180, 360, and 720 m) are used to examine the impact of different levels of landscape representation on the hydrological response of a 690‐km² catchment in southern Quebec. Frequency distributions of local slope, plan curvature, and drainage area are calculated for each grid size resolution. This landscape analysis reveals that DEM grid size significantly affects computed topographic attributes, which in turn explains some of the differences in the hydrological simulations. The simulations that are then carried out, using a coupled, process‐based model of surface and subsurface flow, examine the effects of grid size on both the integrated response of the catchment (discharge at the main outlet and at two internal points) and the distributed response (water table depth, surface saturation, and soil water storage). The results indicate that discharge volumes increase as the DEM is coarsened, and that coarser DEMs are also wetter overall in terms of water table depth and soil water storage. The reasons for these trends include an increase in the total drainage area of the catchment for larger DEM cell sizes, due to aggregation effects at the boundary cells of the catchment, and to a decrease in local slope and plan curvature variations, which in turn limits the capacity of the watershed to transmit water downslope and laterally. The results obtained also show that grid resolution effects are less pronounced during dry periods when soil moisture dynamics are mostly controlled by vertical fluxes of evaporation and percolation. Copyright © 2010 John Wiley & Sons, Ltd.     

Simulation of evapotranspiration and carbon dioxide flux in the wheat‐maize rotation croplands of the North China Plain using the Simple Biosphere Model

The North China Plain, which is critical for food production in China, is encountering serious water shortage due to heavy agricultural water requirement. The accurate modelling of carbon dioxide flux and evapotranspiration (ET) in croplands is thus essential for yield prediction and water resources management. The land surface model is powerful in simulating energy and carbon dioxide fluxes between land and atmosphere. Some key processes in the Simple Biosphere Model (Version 2, SiB2) were parameterized based on the observations. The simulated fluxes were tested against the eddy covariance flux measurements over two typical winter wheat/maize double cropping fields. A simple diagnostic parameterisation of soil respiration, not included in SiB2, was added and calibrated using the observations to model the carbon budget. The Ball‐Berry stomatal conductance model was calibrated using observed leaf gas exchange rate, showing that the original SiB2 could significantly underpredict the ET in the wheat field. SiB2 significantly underpredicted soil resistance at the Weishan site, leading to overpredict the ET. Overall, there was a close agreement between the simulated and observed latent heat fluxes and net CO₂ exchange using the re‐parameterized SiB2. These findings are important when the model is used for the regional simulation in the North China Plain. Copyright © 2011 John Wiley & Sons, Ltd.     

Modelling the impacts of subsurface drainage on surface runoff and sediment yield in the Le Sueur Watershed,Minnesota, USA

Abstract

Tile drainage influences infiltration and surface runoff and is thus an important factor in the erosion process. Tile drainage reduces surface runoff, but questions abound on its influence on sediment transport through its dense network and into the stream network. The impact of subsurface tiling on upland erosion rates in the Le Sueur River watershed, USA, was assessed using the Water Erosion Prediction Project (WEPP) model. Six different scenarios of tile drainage with varying drainage coefficient and management type (no till and autumn mulch-till) were evaluated. The mean annual surface runoff depth, soil loss rate and sediment delivery ratio (SDR) for croplands, based on a 30-year simulation for the watershed with untiled autumn mulch-till (Scenario 1), were estimated to be 83.5 mm, 0.27 kg/m² and 86.7%, respectively; on no-till management systems (Scenario 4), the respective results were 72.3 mm, 0.06 kg/m² and 88.2%. Tile drains reduced surface runoff, soil loss and SDR estimates for Scenario 1 by, on average, 14.5, 8.1 and 7.9%, respectively; and for Scenario 4 by an estimated 31.5, 22.1 and 20.2%, respectively. The impact of tile drains on surface runoff, soil loss and SDR was greater under the no-till management system than under the autumn mulch-till management system. Comparison of WEPP outputs with those of the Soil Water Assessment Tool (SWAT) showed differences between the two methods.

Editor Z.W. Kundzewicz

Citation Maalim, F.K. and Melesse. A.M., 2013. Modelling impacts of subsurface drainage on surface runoff and sediment yield in the Le Sueur Watershed in Minnesota, USA. Hydrological Sciences Journal, 58 (3), 570–586.     

Microtopography,water storage and flow patterns in a fine‐textured soil under agricultural use

Agricultural use of soils implies tillage and often compaction and therefore influences processes on soil surface and affects infiltration of water into the subsoil. Although many studies on soil surface processes or flow patterns in soils exist, works relating both are rare in literature. We did two tracer experiments with Brilliant Blue FCF on a tilled and compacted plot and a non‐tilled one to investigate water storage on the soil surface during simulated rainfall and changes of soil microtopography, to analyse the associated flow patterns in the soil and to relate both to tillage and compaction. Our results show that storage was larger on the tilled and compacted plot than on the non‐tilled one. After tillage, transport processes above the plough pan were partly disconnected from those underneath because macropores were disrupted and buried by the tillage operation. However, preferential flow along cracks occurred on both plots and the macropores buried below the tillage pan still functioned as preferential flow paths. Therefore, we conclude that the studied soil is susceptible to deep vertical solute propagation at dry conditions when cracks are open, irrespective of tillage and compaction. Copyright © 2012 John Wiley & Sons, Ltd.     

Age and origin of the lowlands of Mars

One of the many significant findings of the Mars Global Surveyor mission is the presence of hundreds of quasi-circular depressions (QCDs) observed from high-resolution MOLA topography data. Their presence has recently been interpreted to reflect a northern lowlands that archive some of the earliest recorded rocks on Mars, mostly below a veneer of Hesperian and Amazonian materials. Here we analyze these data, coupled with a recent synthesis of geologic, geophysical, geomorphic, topographic, and magnetic information. Such analysis allows us to suggest a potential plate tectonic phase during the recorded Early into Middle Noachian martian history that transitioned into a monoplate world with episodic magmatic-driven activity persisting to present. This working hypothesis is based on: (1) the observation that the basement of the northern plains is younger than the basement of the southern highlands, but older than the material exposures of the cratered highlands, suggesting different formational ages for each one of the three geologic-time units; (2) the observation that parts of the very ancient highland's crust are highly magnetized, thus suggesting that most if not all of the formation of the lowlands basement postdates the shut off of the martian dynamo, some 4 Gyr ago, and so allowing hundreds of millions of years for the shaping of the buried lowlands. Consequently, the role of endogenic processes in the earliest geological evolution of Mars (Early perhaps into Middle Noachian) requires reconsideration, since MOLA topographic and MGS magnetic data afford a temporal window sufficient for very early, non-primordial shaping of the northern lowlands' basement.     

Evaluating sulfur dynamics during storm events for three watersheds in the northeastern USA: a combined hydrological,chemical and isotopic approach

Concerns related to climate change have resulted in an increasing interest in the importance of hydrological events such as droughts in affecting biogeochemical responses of watersheds. The effects of an unusually dry summer in 2002 had a marked impact on the biogeochemistry of three watersheds in the north‐eastern USA. Chemical, isotopic and hydrological responses with particular emphasis on S dynamics were evaluated for Archer Creek (New York), Sleepers River (Vermont) and Cone Pond (New Hampshire) watersheds. From 1 August to 14 September 2002, all three watersheds had very low precipitation (48 to 69 mm) resulting in either very low or no discharge (mean 0·015, 0·15 and 0·000 mm day^?1 for Archer Creek, Sleepers River and Cone Pond, respectively). From 15 September to 31 October 2002, there was a substantial increase in precipitation totals (212, 246 and 198 mm, respectively) with increased discharge. Archer Creek was characterized by a large range of SO₄^2? concentrations (152 to 389 µeq L^?1, mean = 273 µeq L^?1) and also exhibited the greatest range in δ³⁴S values of SO₄^2? (?1·4 to 8·8 ‰ ). Sleepers River's SO₄^2? concentrations ranged from 136 to 243 µeq L^?1 (mean = 167 µeq L^?1) and δ³⁴S values of SO₄^2? ranged from 4·0 to 9·0 ‰ . Cone Pond's SO₄^2? concentrations (126–187 µeq L^?1, mean = 154 µeq L^?1) and δ³⁴S values (2·4 to 4·3 ‰ ) had the smallest ranges of the three watersheds. The range and mean of δ¹⁸O‐SO₄^2? values for Archer Creek and Cone Pond were similar (3·0 to 8·9 ‰ , mean = 4·5 ‰ ; 3·9 to 6·3 ‰ , mean = 4·9 ‰ ; respectively) while δ¹⁸O‐SO₄^2? values for Sleepers River covered a larger range with a lower mean (1·2 to 10·0 ‰ , mean = 2·5). The difference in Sleepers River chemical and isotopic responses was attributed to weathering reactions contributing SO₄^2?. For Archer Creek wetland areas containing previously reduced S compounds that were reoxidized to SO₄^2? probably provided a substantial source of S. Cone Pond had limited internal S sources and less chemical or isotopic response to storms. Differences among the three watersheds in S biogeochemical responses during these storm events were attributed to differences in S mineral weathering contributions, hydrological pathways and landscape features. Further evaluations of differences and similarities in biogeochemical and hydrological responses among watersheds are needed to predict the impacts of climate change. Copyright © 2008 John Wiley & Sons, Ltd.     

Moraine development at a small High‐Arctic valley glacier: Rieperbreen,Svalbard

Ice‐cored lateral and frontal moraine complexes, formed at the margin of the small, land‐based Rieperbreen glacier, central Svalbard, have been investigated through field observations and interpretations of aerial photographs (1936, 1961 and 1990). The main focus has been on the stratigraphical and dynamic development of these moraines as well as the disintegration processes. The glacier has been wasting down since the ‘Little Ice Age’ (LIA) maximum, and between 1936 and 1990 the glacier surface was lowered by 50–60 m and the front retreated by approximately 900 m. As the glacier wasted, three moraine ridges developed at the front, mainly as melting out of sediments from debris‐rich foliation and debris‐bands formed when the glacier was polythermal, probably during the LIA maximum. The disintegration of the moraines is dominated by wastage of buried ice, sediment gravity‐flows, meltwater activity and some frost weathering. A transverse glacier profile with a northward sloping surface has developed owing to the higher insolation along the south‐facing ice margin. This asymmetric geometry also strongly affects the supraglacial drainage pattern. Lateral moraines have formed along both sides of the glacier, although the insolation aspect of the glacier has resulted in the development of a moraine 60 m high along its northern margin. Copyright © 2001 John Wiley & Sons, Ltd.