Estimation of regional evapotranspiration over the North China Plain using geostationary satellite data

Data from the first operational Chinese geostationary satellite Fengyun-2C (FY-2C) satellite are applied in combination with Moderate Resolution Imaging Spectroradiometer (MODIS) satellite products for the assessment of regional evapotranspiration over the North China Plain. The approach is based on the improved triangle method, where the temperature–vegetation index space includes thermal inertia. Two thermal infrared channels from FY-2C are used to estimate surface temperature (Ts) based on a split window algorithm originally proposed for the MSG-SEVIRI sensor. Subsequently the high temporal resolution of FY-2C data is exploited to give the morning rise in Ts. Combined with the 16 days composite MODIS vegetation indices product (MOD13) at a spatial resolution of 5 km, evaporative fraction (EF) is estimated by interpolation in the ΔTs–NDVI triangular-shaped scatter space. Finally, regional actual evapotranspiration (ET) is derived from the evaporative fraction and available energy estimated from MODIS surface albedo products MCD43. Spatial variations of estimated surface variables (Ts, EF and ET) corresponded well to land cover patterns and farmland management practices. Estimated ET and EF also compared well to lysimeter data collected for the period June 2005–September 2007. The improved triangle method was also applied to MODIS products for comparison. Estimates based on FY-2C products proved to provide slightly better results than those based on MODIS products. The consistency of the estimated spatial variation with other spatial data supports the use of FY-2C data for ET estimation using the improved triangle method. Of particular value is the high temporal frequency of image acquisitions from FY-2C which improves the likelihood of obtaining cloud free image acquisitions as compared to polar orbiting sensors like MODIS.     

Interpretations and observations of ocean wave spectra

The paper starts with a discussion of the linear stochastic theory of ocean waves and its various nonlinear extensions. The directional spectrum, with its unique dispersion relation connecting frequency (ω) and wavenumber (k), is no longer valid for nonlinear waves, and examples of $\left( \mathbf{k},\omega\right) The paper starts with a discussion of the linear stochastic theory of ocean waves and its various nonlinear extensions. The directional spectrum, with its unique dispersion relation connecting frequency (ω) and wavenumber (k), is no longer valid for nonlinear waves, and examples of ( k,w)\left( \mathbf{k},\omega\right) -spectra based on analytical expressions and computer simulations of nonlinear waves are presented. Simulations of the dynamic nonlinear evolution of unidirectional free waves using the nonlinear Schr?dinger equation and its generalizations show that components above the spectral peak have larger phase and group velocities than anticipated by linear theory. Moreover, the spectrum does not maintain a thin well-defined dispersion surface, but rather develops into a continuous distribution in ( k,w)\left( \mathbf{k,}\omega\right) -space. The majority of existing measurement systems rely on linear theory for the interpretation of their data, and no measurement systems are currently able to measure the full spectrum in the open ocean with high accuracy. Nevertheless, there exist a few low-resolution systems where data may be interpreted within a minimal assumption of a non-restricted ( k,w)\left( \mathbf{k,}\omega\right) -spectrum. The theory is reviewed, and analyses based on conventional spectral analysis as well as a directional wavelet analysis are carried out on data from a compact laser array at the Ekofisk field in the North Sea. The investigation confirms the strong impact of the second order spectrum below the spectral peak, but is non-conclusive about the off-set in the support of the first order spectrum seen in the dynamical simulations.     

Post-deformational annealing of calcite rocks

The evolution of microstructure and crystallographic preferred orientation (CPO) during post-deformational annealing was studied on three calcite rock types differing in purity and grain size: Carrara marble (98% calcite, mean grain size of 115 μm), Solnhofen limestone (96%, 5 μm) and synthetic calcite aggregates (99%, 7 μm). Samples were first deformed in torsion at 727 °C at a shear strain rate of 3 × 10^{− 4} s^{− 1} to a shear strain of 5 and subsequently heat-treated at 727 °C for various durations between 0 and 24 h. Microstructures and CPOs were analysed by optical microscopy, image analysis and electron backscatter diffraction (EBSD).All rock types deformed in the dislocation creep field at the same applied conditions, but their microstructures and CPOs after deformation and after annealing differed depending on starting grain size and material composition. In Carrara marble and in the synthetic calcite aggregate, a strong CPO developed during deformation accompanied by dynamic recrystallisation with significant changes in grain size. During annealing, widespread grain growth and subtle changes of CPO occurred, and equilibrated foam microstructures were approached after long annealing times. The CPO is the only feature in annealed samples indicating an earlier deformation phase, although it is not always identical to the CPO formed during deformation. In the more impure Solnhofen limestone, secondary phases on grain boundaries suppressed grain boundary mobility and prevented both the formation of a recrystallisation CPO during deformation and grain size modification during deformation and annealing.     

Multi-Constrained Catchment Scale Optimization of Groundwater Abstraction Using Linear Programming

Due to increasing water demands globally, freshwater ecosystems are under constant pressure. Groundwater resources, as the main source of accessible freshwater, are crucially important for irrigation worldwide. Over-abstraction of groundwater leads to declines in groundwater levels; consequently, the groundwater inflow to streams decreases. The reduction in baseflow and alteration of the streamflow regime can potentially have an adverse effect on groundwater-dependent ecosystems. A spatially distributed, coupled groundwater–surface water model can simulate the impacts of groundwater abstraction on aquatic ecosystems. A constrained optimization algorithm and a simulation model in combination can provide an objective tool for the water practitioner to evaluate the interplay between economic benefits of groundwater abstractions and requirements to environmental flow. In this study, a holistic catchment-scale groundwater abstraction optimization framework has been developed that allows for a spatially explicit optimization of groundwater abstraction, while fulfilling a predefined maximum allowed reduction of streamflow (baseflow [Q95] or median flow [Q50]) as constraint criteria for 1484 stream locations across the catchment. A balanced K-Means clustering method was implemented to reduce the computational burden of the optimization. The model parameters and observation uncertainties calculated based on Bayesian linear theory allow for a risk assessment on the optimized groundwater abstraction values. The results from different optimization scenarios indicated that using the linear programming optimization algorithm in conjunction with integrated models provides valuable information for guiding the water practitioners in designing an effective groundwater abstraction plan with the consideration of environmental flow criteria important for the ecological status of the entire system.     

Corrections for near surface effects: Decomposition of the magnetotelluric impedance tensor and scaling corrections for regional resistivities: A tutorial

This paper primarily examines the effects of small-scale or near-surface conductivity inhomogeneities on the magnetotelluric (MT) impedance tensor. These effects cause three different types of distortion results. (1) The well-knownstatic shifts of sounding curves. (2) When the underlying regional setting is two-dimensional then the two regional impedances are mixed in an arbitrary coordinate system. Thus the level and shape of each sounding curve is distorted as are the phases. (3) At sufficiently high frequencies these effects generate anomalous magnetic fields that in turn alter the background phases.This tutorial first explores the usefulness of various MT tensor analysis techniques to overcome the problem of phase mixing and to recover regional information in the presence of local geological noise. Synthetic and experimental data are considered. A sequence ofa priori models of increasing complexity are described. The use of appropriate decompositions of the MT tensor each with an increasing number of parameters is emphasised. In a second part, phase mixing and static shifts are examined from a synoptic view. Some static shift removal techniques that can be used in conjunction with the decomposition are discussed. This paper is not a review but rather an investigation of a few methods that the authors have found useful with field data.     

Invasion genetics of Pacific oyster Crassostrea gigas shaped by aquaculture stocking practices

As a result of aquaculture activities Pacific oysters Crassostrea gigas (Thunberg, 1793) have invaded the European Wadden Sea. Using a variable noncoding mitochondrial marker, we show that the invaded range is the result of two independent invasions. Haplotype frequencies point towards two separate groups, one in the southern and the other in the northern Wadden Sea. We found virtually no genetic differentiation throughout the southern range and the putative source from British Columbia, Canada, suggesting that the Southern region can be considered as a closed population. In the North, mismatch distributions, haplotype ordination and isolation-by-distance analysis suggest a stronger, persistent impact of aquaculture on invasive populations. Due to the ongoing supply of new genetic material from hatchery production the northern invasive populations can therefore be considered as an open population highlighting the importance of aquaculture practice on the genetics of this keystone invader in the Wadden Sea.     

Petrogenesis of ultramafic and mafic xenoliths from Mesozoic basanites in southern Sweden: constraints from mineral chemistry

Jurassic basanite necks occurring at the junction of two major fault zones in Scania contain ultramafic (peridotites, pyroxenites) and mafic xenoliths, which together indicate a diversity of upper mantle and lower crustal assemblages beneath this region. The peridotites can be subdivided into lherzolites, dunites and harzburgites. Most lherzolites are porphyroclastic, containing orthopyroxene and olivine porphyroclasts. They consist of Mg-rich silicates (Mg# = Mg/(Mg + Fe_tot) × 100; 88–94) and vermicular spinel. Calculated equilibration temperatures are lower in porphyroclastic lherzolites (975–1,007°C) than in equigranular lherzolite (1,079°C), indicating an origin from different parts of the upper mantle. According to the spinel composition the lherzolites represent residues of 8–13% fractional melting. They are similar in texture, mineralogy and major element composition to mantle xenoliths from Cenozoic Central European volcanic fields. Dunitic and harzburgitic peridotites are equigranular and only slightly deformed. Silicate minerals have lower to similar Mg# (83–92) as lherzolites and lack primary spinel. Resorbed patches in dunite and harzburgite xenoliths might be the remnants of metasomatic processes that changed the upper mantle composition. Pyroxenites are coarse, undeformed and have silicate minerals with partly lower Mg# than peridotites (70–91). Pyroxenitic oxides are pleonaste spinels. According to two-pyroxene thermometry pyroxenites show a large range of equilibration temperatures (919–1,280°C). In contrast, mafic xenoliths, which are mostly layered gabbronorites with pyroxene- and plagioclase-rich layers, have a narrow range of equilibration temperatures (828–890°C). These temperature ranges, together with geochemical evidence, indicate that pyroxenites and gabbroic xenoliths represent mafic intrusions within the Fennoscandian crust.     

Documentation of Istanbul Historic Peninsula by Kinematic Terrestrial Laser Scanning

Although the Historic Peninsula of old Istanbul was added to the UNESCO World Heritage List in 1985, complete documentation of this important area has not yet been carried out. In 2006 the Istanbul Metropolitan Municipality's Historical Environment Protection Directorate initiated the "Historic Peninsula Project", which covers an area of 1500 ha and includes approximately 48 000 buildings in crowded and frequently narrow streets. BİMTAŞ, a company owned by the Municipality, immediately started the documentation of all buildings in the project area using terrestrial laser scanning. This created the challenge of building up an efficient production environment with new high-end technology to fulfil the requirements of this project in a very short timeframe of 2 years. This paper describes the entire production environment for documentation of all buildings, detailing the frequent adaptations of the process resulting from learning on the job. Although the data acquisition and mapping environment was established in the course of production, it was always essential to optimise the technical solutions in order to meet the requirements for data quality and delivery deadlines. Only 80 ha of the required 1500 ha was completed using static scanning during the first 6 months, thus requiring a change from static to mobile terrestrial laser scanning in order to accelerate the work and to conclude the scanning phase for the remaining major area within 3 months.