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.     

Effects of soil heterogeneity on martian ground-ice stability and orbital estimates of ice table depth

Data from the Mars Odyssey Gamma-Ray Spectrometer (GRS) instrument suite and results from numerical simulations of subsurface ground-ice stability have been used to estimate the depth of martian ground-ice. Geographic correlation between these estimates is remarkable; the relative ice table depth distributions also agree well. However, GRS-based estimates of ice table depth are generally deeper than predictions based on ground-ice stability simulations. This discrepancy may be related to heterogeneities in the martian surface such as rocks, dust, and albedo variations. We develop a multi-dimensional numerical model of ground-ice stability in a heterogeneous martian subsurface and use it to place the first quantitative constraints on the response of the ice table to meter-scale heterogeneities. We find that heterogeneities produce significant undulations/topography in the ice table at horizontal length scales of a few meters. Decimeter scale rocks create localized areas of deep ice, producing a vertical depression of 10-60 cm in the ice table over a horizontal range of 1-2 rock radii. Decimeter scale dust lenses produce locally shallow ice; however the magnitude of the vertical deflection of the ice table is small (1-4 cm). The effects of decimeter scale albedo variations on the ice table are nearly negligible, although albedo very weakly enhances the effects of dark rocks and bright dust on the ice table. Additionally, we investigate the role played by rocks in estimates of ice table depth based on orbital data. Surface rocks can account for more than half of the discrepancy between ice table depths inferred from GRS data and those predicted by theoretical ice-stability simulations that utilize thermophysical observations. Our results have considerable relevance to the up-coming Mars Scout Mission, Phoenix, because they indicate that the uncertainty in the ice table depth of a given region is greater than differences between current depth estimates. Likewise, small-scale depth variability due to heterogeneities at the eventual landing site is potentially greater than differences between current depth estimates.     

Palaeosurfaces in central West Greenland as reference for identification of tectonic movements and estimation of erosion

Landform analysis of basement rocks has been undertaken with the aid of digital elevation data, aerial photographs and field observations in central West Greenland (69°15′N–66°00′N). Palaeosurfaces have been identified, dated relatively to each other, used to quantify uplift and fault movements and also used to estimate differential erosion. Two types of palaeosurfaces were mapped across the Precambrian basement: a surface at low elevation with distinct hills (hilly relief), and two planation surfaces formed across different types of basement rocks. The hilly relief surface emerges as an inclined surface from Cretaceous cover rocks in Disko Bugt and is interpreted as a stripped late Mesozoic etch surface. This surface is cut off towards the south by a less inclined planation surface, which is younger and thus of Cenozoic age. It is similar to the post-Eocene (Miocene?) planation surfaces identified on Disko and Nuussuaq in other studies. The planation surface splits in two southwards towards high areas around Nordre Isortoq and Sukkertoppen Ice Cap. The upper planation surface forms near-summit areas of tectonic blocks dipping in different directions and with different tilts. The uplift centres define the crests of two mega blocks, separated by the ‘Sisimiut Line’ which coincides with the Precambrian Ikertôq thrust zone. A partially developed lower planation surface indicates a first uplift of maximum 500 m followed by a second uplift of maximum 1000 m. We infer that these uplift events occurred during the late Neogene based on correlation with similar surfaces on Nuussuaq and the timing of exhumational events estimated from apatite fission track analyses of samples from a deep borehole on Nuussuaq (reported elsewhere). The difference between a reconstruction of the upper planation surface across the entire area and the present topography was used as an estimate of erosion of basement rock since the formation of the upper planation surface. The erosion is unevenly distributed and varies from almost none on the well-preserved planation surfaces to 800–1300 m along valleys, and even more in the fjords. Erosion is less within areas of gneiss in granulite facies, than in areas of gneiss in amphibolite facies.     

A modified land surface temperature split window retrieval algorithm and its applications over China

Due to the difficulties in correcting the influences of the atmosphere absorbability and the Earth surface emissivity diversification, the retrieval of LST (land surface temperature) from satellite data is a challenging task. In this paper, a modified Becker's split window LST inversion algorithm is developed for retrieving LST from the NOAA-16/17 AVHRR data. A new set of parameters for Becker's LST algorithm is proposed. The algorithm is developed from a surface brightness temperature dataset generated from the MODTRAN program, which uses a range of surface parameters and atmospheric quantities as inputs. The 10-day composites of the channels 4 and 5 brightness temperature data of NOAA-17 AVHRR (1-km resolution) are used to generate the clear-sky LST. As a validation of the algorithm, the retrieved LST is compared with MODIS LST of same period and area. The two LST products are found to be consistent, with the absolute difference being about 2.5 K for most areas. The NOAA retrieved LST is also compared with in-situ ground surface 0-cm measurements taken from 257 meteorological stations, which cover overall China area for the three periods of satellite observations. The comparison shows that the correlation between the retrieved LST and in-situ measurements is over 0.90 and the RMSE (root mean square error) is about 3.4 K.     

Layering stratigraphy of eastern Coprates and northern Capri Chasmata, Mars

Distinct competent layers are observed in the slopes of eastern Coprates Chasma, part of the Valles Marineris system on Mars. Our observations indicate that the stratigraphy of Coprates Chasma consists of alternating thin strong layers and thicker sequences of relatively weak layers. The strong, competent layers maintain steeper slopes and play a major role in controlling the overall shape and geomorphology of the chasmata slopes. The topmost competent layer in this area is well preserved and easy to identify in outcrops on the northern rim of Coprates Chasma less than 100 m below the southern Ophir Planum surface. The volume of the topmost emplaced layer is at least 70 km³ and may be greater than 2100 km³ if the unit underlies most of Ophir Planum. The broad extent of this layer allows us to measure elevation offsets within the north rim of the chasma and in a freestanding massif within Coprates Chasma where the layer is also observed. Rim outcrop morphology and elevation differences between Ophir and Aurorae Plana may be indicative of the easternmost extent of the topmost competent layer. These observations allow an insight into the depositional processes that formed the stratigraphic stack into which this portion of the Valles Marineris is carved, and they present a picture of some of the last volcanic activity in this area. Furthermore, the elevation offsets within the layer are evidence of significant subsidence of the massif and surrounding material.