Quantifying uncertainty in remotely sensed land cover maps

Remotely sensed land cover maps are increasingly used as inputs into environmental simulation models whose outputs inform decisions and policy-making. Risks associated with these decisions are dependent on model output uncertainty, which is in turn affected by the uncertainty of land cover inputs. This article presents a method of quantifying the uncertainty that results from potential mis-classification in remotely sensed land cover maps. In addition to quantifying uncertainty in the classification of individual pixels in the map, we also address the important case where land cover maps have been upscaled to a coarser grid to suit the users’ needs and are reported as proportions of land cover type. The approach is Bayesian and incorporates several layers of modelling but is straightforward to implement. First, we incorporate data in the confusion matrix derived from an independent field survey, and discuss the appropriate way to model such data. Second, we account for spatial correlation in the true land cover map, using the remotely sensed map as a prior. Third, spatial correlation in the mis-classification characteristics is induced by modelling their variance. The result is that we are able to simulate posterior means and variances for individual sites and the entire map using a simple Monte Carlo algorithm. The method is applied to the Land Cover Map 2000 for the region of England and Wales, a map used as an input into a current dynamic carbon flux model.

Continuous sky view factor maps from high resolution urban digital elevation models

This paper presents a new method to calculate sky view factors (SVFs) from high resolution urban digital elevation models using a shadow casting algorithm. By utilizing weighted annuli to derive SVF from hemispherical images, the distance light source positions can be predefined and uniformly spread over the whole hemisphere, whereas another method applies a random set of light source positions with a cosine-weighted distribution of sun altitude angles. The 2 methods have similar results based on a large number of SVF images. However, when comparing variations at pixel level between an image generated using the new method presented in this paper with the image from the random method, anisotropic patterns occur. The absolute mean difference between the 2 methods is 0.002 ranging up to 0.040. The maximum difference can be as much as 0.122. Since SVF is a geometrically derived parameter, the anisotropic errors created by the random method must be considered as significant.

Impact of US Brown Swiss genetics on milk quality from low-input herds in Switzerland: interactions with grazing intake and pasture type

This study investigated the effect of, and interactions between, contrasting crossbreed genetics (US Brown Swiss [BS] × Improved Braunvieh [BV] × Original Braunvieh [OB]) and feeding regimes (especially grazing intake and pasture type) on milk fatty acid (FA) profiles. Concentrations of total polyunsaturated FAs, total omega-3 FAs and trans palmitoleic, vaccenic, α-linolenic, eicosapentaenoic and docosapentaenoic acids were higher in cows with a low proportion of BS genetics. Highest concentrations of the nutritionally desirable FAs, trans palmitoleic, vaccenic and eicosapentaenoic acids were found for cows with a low proportion of BS genetics (0-24% and/or 25-49%) on high grazing intake (75-100% of dry matter intake) diets. Multivariate analysis indicated that the proportion of OB genetics is a positive driver for nutritionally desirable monounsaturated and polyunsaturated FAs while BS genetics proportion was positive driver for total and undesirable individual saturated FAs. Significant genetics × feeding regime interactions were also detected for a range of FAs.

On the structure of infinity-harmonic maps

Let H ∈ C 2(ℝ N×n ), H ≥ 0. The PDE system arises as the Euler-Lagrange PDE of vectorial variational problems for the functional E ∞(u, Ω) = ‖H(Du)‖ L ∞(Ω) defined on maps u: Ω ⊆ ℝ n → ℝ N . (1) first appeared in the author's recent work. The scalar case though has a long history initiated by Aronsson. Herein we study the solutions of (1) with emphasis on the case of n = 2 ≤ N with H the Euclidean norm on ℝ N×n , which we call the “∞-Laplacian”. By establishing a rigidity theorem for rank-one maps of independent interest, we analyse a phenomenon of separation of the solutions to phases with qualitatively different behaviour. As a corollary, we extend to N ≥ 2 the Aronsson-Evans-Yu theorem regarding non existence of zeros of |Du| and prove a maximum principle. We further characterise all H for which (1) is elliptic and also study the initial value problem for the ODE system arising for n = 1 but with H(·, u, u′) depending on all the arguments.

'The debatable territory where geology and archaeology meet': reassessing the early archaeobotanical work of Clement Reid and Arthur Lyell at Roman Silchester

The first large-scale archaeobotanical study in Britain, conducted from 1899 to 1909 by Clement Reid and Arthur Lyell at Silchester, provided the first evidence for the introduction of Roman plant foods to Britain, yet the findings have thus far remained unverified. This paper presents a reassessment of these archaeobotanical remains, now stored as part of the Silchester Collection in Reading Museum. The documentary evidence for the Silchester study is summarised, before the results are presented for over a 1000 plant remains including an assessment of preservation, identification and modern contamination. The dataset includes both evidence for the presence of nationally rare plant foods, such as medlar, and several archaeophytes. The methodologies and original interpretations of Reid and Lyell’s study are reassessed in light of current archaeobotanical knowledge. Spatial and contextual patterns in the distribution of plant foods and ornamental taxa are also explored. Finally, the legacy of the study for the development of archaeobotany in the 20th century is evaluated.

Daily thin-ice thickness maps from modis thermal infrared imagery

Accurate knowledge of ice-production rates within the marginal ice zones of the Arctic Ocean requires monitoring of the thin-ice distribution within polynyas. The thickness of the ice layer controls the heat loss and hence the new-ice formation. An established thinice algorithm using high-resolution MODIS data allows deriving the ice-thickness distribution within polynyas. The average uncertainty is ±4.7 cm for ice thicknesses below 0.2 m. In this study, the ice-thickness distributions within the Laptev Sea polynya for the two winter seasons 2007/08 and 2008/09 are calculated. Then, a new method is applied to determine a daily MODIS thin-ice product.

Triple isotope (δD, δ17O, δ18O) study on precipitation, drip water and speleothem fluid inclusions for a central European cave (NW Switzerland)

Deuterium (dD) and oxygen (d18O) isotopes are powerful tracers of the hydrological cycle and have been extensively used for paleoclimate reconstructions as they can provide information on past precipitation, temperature and atmospheric circulation. More recently, the use of 17Oexcess derived from precise measurement of d17O and d18O gives new and additional insights in tracing the hydrological cycle whereas uncertainties surround this proxy. However, 17Oexcess could provide additional information on the atmospheric conditions at the moisture source as well as about fractionations associated with transport and site processes. In this paper we trace water stable isotopes (dD, d17O and d18O) along their path from precipitation to cave drip water and finally to speleothem fluid inclusions for Milandre cave in northwestern Switzerland. A two year-long daily resolved precipitation isotope record close to the cave site is compared to collected cave drip water (3 months average resolution) and fluid inclusions of modern and Holocene stalagmites. Amount weighted mean dD, d18O and d17O are �71.0‰, �9.9‰, �5.2‰ for precipitation, �60.3‰, �8.7‰, �4.6‰ for cave drip water and �61.3‰, �8.3‰, �4.7‰ for recent fluid inclusions respectively. Second order parameters have also been derived in precipitation and drip water and present similar values with 18 per meg for 17Oexcess whereas d-excess is 1.5‰ more negative in drip water. Furthermore, the atmospheric signal is shifted towards enriched values in the drip water and fluid inclusions (D of ~ þ 10‰ for dD). The isotopic composition of cave drip water exhibits a weak seasonal signal which is shifted by around 8e10 months (groundwater residence time) when compared to the precipitation. Moreover, we carried out the first d17O measurement in speleothem fluid inclusions, as well as the first comparison of the d17O behaviour from the meteoric water to the fluid inclusions entrapment in speleothems. This study on precipitation, drip water and fluid inclusions will be used as a speleothem proxy calibration for Milandre cave in order to reconstruct paleotemperatures and moisture source variations for Western Central Europe.

Timing and nature of the penultimate deglaciation in a high alpine stalagmite from Switzerland

The timing and nature of the penultimate deglaciation, also known as Termination II (T-II), is subject of controversial discussions due to the scarcity of precisely-dated palaeoclimate records. Here we present a new precisely-dated and highly-resolved multi-proxy stalagmite record covering T-II from the high alpine Schafsloch Cave in Switzerland, an area where climate is governed by the North Atlantic. The inception of stalagmite growth at 137.4 ± 1.4 kyr before present (BP) indicates the presence of drip water and cave air temperatures of above 0°C, and is related to a climate-induced change in the thermal state (from cold- to warm-based) of the glacier above the cave. The cessation of stalagmite growth between 133.1 ± 0.7 and 131.9 ± 0.6 kyr BP is most likely related to distinct drop in temperature associated with Heinrich stadial 11. The resumption of stalagmite growth at 131.9 ± 0.6 kyr BP is accompanied by an abrupt increase in temperature and precipitation as indicated by distinct shifts in the oxygen and carbon isotopic composition as well as in trace element concentrations. The mid-point of T-II is around 131.8 ± 0.6 kyr BP in the Schafsloch Cave record is significantly earlier compared to the age of 129.1 ± 0.1 kyr BP in the Sanbao Cave record from China. The different ages between both records can be best explained by the competing effects of insolation and glacial boundary forcing on seasonality and snow cover extent in Eurasia.