A Spectrum-Independent Procedure for Correcting Eddy Fluxes Measured with Separated Sensors

We investigate flux underestimates in eddy correlation measurements that are caused by horizontal separation of the sensors. A common eddy correlation setup consists of a sonic anemometer and a humidity sensor which, because of its bulk, must be placed some distance away from the sonic path, leading to a flux loss (of latent heat). Utilizing an additional fast temperature sensor placed near the humidity sensor, we develop a procedure for correcting for this loss. The procedure simultaneously corrects the sensible heat flux for the difference between true temperature and sonic temperature. Our correction procedure, which does not depend on the shape of the cospectrum, is then compared to the widely-used procedure following Moore (1986), which assumes a cospectral model ('Kansas Model). Both correction methods are applied to data collected within the internal boundary layer over a rice paddy, downwind of arid land. Under conditions of good fetch, they were found to agree well. Under poor fetch conditions, the model-based correction tended to be too small, while the spectrum-independent combined correction was robust. The latter is thus recommended for situations where the cospectral shape can be expected to deviate from the 'Kansas shape.     

Rare earth element geochemistry of potassic lavas from the Birunga and Toro-Ankole regions of Uganda,Africa

Neutron activation determination of La, Ce, Sm, Eu, Tb, Yb, Lu, Ta, Hf, Sc, Co and Th in potassic lavas from the Birunga and Toro-Ankole regions show that the rocks are characterized by high rare earth element (REE) contents (161–754 ppm) and form two groups based upon differing La/Yb ratios. One group is made up of katungite, ugandite and mafurite with La/Yb =146–312, and the other of rocks of the leucitite and phonolitic tephrite series, La/Yb =30–56. The trace element content of the ugandite group is similar to that of kimberlites. The data do not indicate any trends of differentiation or simple relationships between the two groups of rocks, although katungite is unlikely to be parental to rocks of lower La/Yb ratios. It is unlikely that in terms of La/Yb ratios that partial melting of mica-garnet-lherzolite mantle can form katungite because of the very small amounts of partial melting required (0.2%), although the La/Yb ratios of 150–200 (ugandites, mafurites) and 30–60 (leucitites, phonolitic tephrites) can be accounted for by 0.3–1.5% and 1–9% melting respectively, if the REE are then concentrated without further La and Yb fractionation. Partial melting of mantle which has been metasomatized by alkaline earths and REE bearing fluids or mixing of carbonatite and nephelenite are also compatable with the observed geochemistry of the lavas. It is considered that gas transfer processes which selectively enrich the light REE may have obscured REE evidence pertaining to early partial melting and/or differentiation processes and therefore that REE geochemistry is of little use in determining the petrogenetic processes involved in the formation of potassic lavas.     

Carbon and sulfur in the Swaziland sequence

Sedimentary and igneous rocks of the Swaziland Sequence of South Africa (>3,200 my. old) contain total carbon (0.09–8.28%), organic carbon (0.04–1.95%) and sulfur (9–53,000 ppm). Amounts of organic carbon are generally high in the lower part of the Sequence (based on two samples from the Theespruit Formation). Sulfur content varies widely with no obvious correlations or trends. These rocks may contain clues to the early carbon and sulfur cycles on earth.     

A discrete particle representation of hillslope hydrology: hypothesis testing in reproducing a tracer experiment at Gårdsjön,Sweden

Despite the long history of the continuum equation approach in hydrology, it is not a necessary approach to the formulation of a physically based representation of hillslope hydrology. The Multiple Interacting Pathways (MIPs) model is a discrete realization that allows hillslope response and transport to be simultaneously explored in a way that reflects the potential occurrence of preferential flows and lengths of pathways. The MIPs model uses random particle tracking methods to represent the flow of water within the subsurface alongside velocity distributions that acknowledge preferential flows and transition probability matrices, which control flow pathways. An initial realization of this model is presented here in application to a tracer experiment carried out in Gårdsjön, Sweden. The model is used as an exploratory tool, testing several hypotheses in relation to this experiment. Copyright © 2011 John Wiley & Sons, Ltd.     

Environmental controls on phytoplankton community composition in the Thames plume, U.K.

The aim of this study was to investigate controls on the phytoplankton community composition and biogeochemistry of the estuarine plume zone of the River Thames, U.K. using an instrumented moored buoy for in situ measurements and preserved sample collection, and laboratory-based measurements from samples collected at the same site. Instrumentation on the moored buoy enabled high frequency measurements of a suite of environmental variables including in situ chlorophyll, water-column integrated irradiance, macronutrients throughout an annual cycle for 2001 e.g. nitrate and silicate, and phytoplankton biomass and species composition. The Thames plume region acts as a conduit for fluvial nutrients into the wider southern North Sea with typical winter concentrations of 45 μM nitrate, 17 μM silicate and 2 μM phosphate measured. The spring bloom resulted from water-column integrated irradiance increasing above 60 W h m^− 2 d^− 1 and was initially dominated by a diatom bloom mainly composed of Nitzschia sp. and Odontella sinesis. The spring bloom then switched after  30 days to become dominated by the flagellate Phaeocystis reaching a maximum chlorophyll concentration of 37.8 μg L^− 1. During the spring bloom there were high numbers of the heterotrophic dinoflagellates Gyrodinium spirale and Katodinium glaucum that potentially grazed the phytoplankton bloom. This diatom–flagellate switch was predicted to be due to a combination of further increasing water-column integrated irradiance > 100 W h m^− 2 d^− 1 and/or silicate reaching potentially limiting concentrations (< 1 μM). Post spring bloom, diatom dominance of the lower continuous summer phytoplankton biomass occurred despite the low silicate concentrations (Av. 0.7 μM from June–August). Summer diatom dominance, generally due to Guinardia delicatula, was expected to be as a result of microzooplankton grazing, dominated by the heterotrophic dinoflagellate Noctiluca scintillans, controlling 0.7–5.0 μm ‘flagellate’ fraction of the phytoplankton community with grazing rates up to 178% of ‘flagellate’ growth rate. The Thames plume region was therefore shown to be an active region of nutrient and phytoplankton processing and transport to the southern North Sea. The use of a combination of moorings and ship-based sampling was essential in understanding the factors influencing nutrient transport, phytoplankton biomass and species composition in this shelf sea plume region.     

Examining the sensitivity of spatial scale in cellular automata Markov chain simulation of land use change

Understanding the spatial scale sensitivity of cellular automata is crucial for improving the accuracy of land use change simulation. We propose a framework based on a response surface method to comprehensively explore spatial scale sensitivity of the cellular automata Markov chain (CA-Markov) model, and present a hybrid evaluation model for expressing simulation accuracy that merges the strengths of the Kappa coefficient and of Contagion index. Three Landsat-Thematic Mapper remote sensing images of Wuhan in 1987, 1996, and 2005 were used to extract land use information. The results demonstrate that the spatial scale sensitivity of the CA-Markov model resulting from individual components and their combinations are both worthy of attention. The utility of our proposed hybrid evaluation model and response surface method to investigate the sensitivity has proven to be more accurate than the single Kappa coefficient method and more efficient than traditional methods. The findings also show that the CA-Markov model is more sensitive to neighborhood size than to cell size or neighborhood type considering individual component effects. Particularly, the bilateral and trilateral interactions between neighborhood and cell size result in a more remarkable scale effect than that of a single cell size.     

Tomography: A window on the role of sulfur in the structure of micrometeorites

Abstract– To determine the role played by sulfides in the formation of vesicles and FeNi metal beads, we mapped the locations and tabulated the numbers of sulfides, metal beads, and vesicles in 1583 sectioned micrometeorites (MMs) using conventional microscopy and in 190 whole MMs using synchrotron computed microtomography (SCMT). Both the section and the SCMT images show that sulfides melt, coalesce, and migrate to the MMs surface. The decomposition of sulfides may occur during all these stages. Given the sulfide morphologies and compositions that we see in section, we think the breakdown of Ni sulfides produces the FeNi beads. The SCMT images show that metal beads are common in melted MMs, >50% have them. Vesicles in porphyritic and scoriaceous MMs are also probably formed as sulfides decompose. Not only do sulfides abut the vesicles but also the temperatures at which sulfides decompose overlap those at which MM surfaces first melt and temporarily seal, suggesting that S gases could produce most of these vesicles. As the vesicle shapes and patterns of distribution differ among MM classes, tomography can be used to nondestructively screen for specific types of MMs. Tomography is a powerful tool for visualizing the three‐dimensional distribution of metal beads, sulfides, mean densities, and vesicles in MMs.     

Toxicity in benthic freshwater cyanobacteria (blue‐green algae): First observations in New Zealand

Toxic cyanobacteria have been linked to dog deaths in the Waikanae River (1998) and Mataura River (1999 and 2000), New Zealand. These were Oscillatoria‐like species which formed benthic mats. This is the first time in New Zealand that animal deaths have been linked to toxins from benthic cyanobacteria.     

Mapping block-and-ash flow hazards based on Titan 2D simulations: a case study from Mt. Taranaki,NZ

Numerical models for simulation of mass flows are typically focussed upon accurately predicting the paths, travel times and inundation from a single flow or collapse event. When considering catchment-based hazards from a volcano, this is complicated by often being faced with several possible scenarios. Over the last 800 years at Mt. Taranaki/Egmont, a number of dome growth and collapse events have resulted in the genesis and emplacement of block-and-ash flows (BAFs). Each BAF was directed northwestward by a breach in the crater rim. The latest dome collapse events in the AD 1880s and AD 1755 inundated the northwestern flank and had run-out lengths 10 km from source. Future activity of this type could have a devastating effect on the Taranaki region’s communities, infrastructure and economy. Hazard planning has involved constructing volcanic hazard maps based upon the areas inundated by past volcanic flows, with little consideration of present-day topography. Here, a numerical geophysical mass flow modelling approach is used to forecast the hazards of future comparable BAF events on NW Mt. Taranaki. The Titan2D programme encompasses a “shallow water”, continuum solution-based, granular flow model. Flow mechanical properties needed for this approach include estimates of internal and basal friction as well as the physical dimensions of the initial collapse. Before this model can be applied to Taranaki BAFs, the input parameters must be calibrated by simulating a range of past collapse events. By using AD 1860 and AD 1755 scenarios, initial collapse volumes can be well constrained and internal and basal friction angles can be evaluated through an iterative approach from previous run-out lengths. A range of possible input parameters was, therefore, determined to produce a suite of potentially inundated areas under present-day terrain. A suite of 10 forecasts from a uniformly distributed range were combined to create a map of relative probabilities of inundation by future BAF events. These results were combined in a GIS package to produce hazard zones related to user-specified hazard thresholds. Using these input parameter constraints, future hazard forecasts for this scale and type of event can also take into account changing summit and topographic configurations following future eruptive or collapse events.