Retrieval of grassland biophysical parameters using multitemporal optical and radar satellite data

The amount and quality of available biomass is a key factor for the sustainable livestock industry and agricultural management related decision making. Globally 31.5% of land cover is grassland while 80% of Ireland’s agricultural land is grassland. In Ireland, grasslands are intensively managed and provide the cheapest feed source for animals. This dissertation presents a detailed state of the art review of satellite remote sensing of grasslands, and the potential application of optical (Moderate–resolution Imaging Spectroradiometer (MODIS)) and radar (TerraSAR-X) time series imagery to estimate the grassland biomass at two study sites (Moorepark and Grange) in the Republic of Ireland using both statistical and state of the art machine learning algorithms. High quality weather data available from the on-site weather station was also used to calculate the Growing Degree Days (GDD) for Grange to determine the impact of ancillary data on biomass estimation. In situ and satellite data covering 12 years for the Moorepark and 6 years for the Grange study sites were used to predict grassland biomass using multiple linear regression, Neuro Fuzzy Inference Systems (ANFIS) models. The results demonstrate that a dense (8-day composite) MODIS image time series, along with high quality in situ data, can be used to retrieve grassland biomass with high performance (R2 = 0:86; p < 0:05, RMSE = 11.07 for Moorepark). The model for Grange was modified to evaluate the synergistic use of vegetation indices derived from remote sensing time series and accumulated GDD information. As GDD is strongly linked to the plant development, or phonological stage, an improvement in biomass estimation would be expected. It was observed that using the ANFIS model the biomass estimation accuracy increased from R2 = 0:76 (p < 0:05) to R2 = 0:81 (p < 0:05) and the root mean square error was reduced by 2.72%. The work on the application of optical remote sensing was further developed using a TerraSAR-X Staring Spotlight mode time series over the Moorepark study site to explore the extent to which very high resolution Synthetic Aperture Radar (SAR) data of interferometrically coherent paddocks can be exploited to retrieve grassland biophysical parameters. After filtering out the non-coherent plots it is demonstrated that interferometric coherence can be used to retrieve grassland biophysical parameters (i. e., height, biomass), and that it is possible to detect changes due to the grass growth, and grazing and mowing events, when the temporal baseline is short (11 days). However, it not possible to automatically uniquely identify the cause of these changes based only on the SAR backscatter and coherence, due to the ambiguity caused by tall grass laid down due to the wind. Overall, the work presented in this dissertation has demonstrated the potential of dense remote sensing and weather data time series to predict grassland biomass using machine-learning algorithms, where high quality ground data were used for training. At present a major limitation for national scale biomass retrieval is the lack of spatial and temporal ground samples, which can be partially resolved by minor modifications in the existing PastureBaseIreland database by adding the location and extent ofeach grassland paddock in the database. As far as remote sensing data requirements are concerned, MODIS is useful for large scale evaluation but due to its coarse resolution it is not possible to detect the variations within the fields and between the fields at the farm scale. However, this issue will be resolved in terms of spatial resolution by the Sentinel-2 mission, and when both satellites (Sentinel-2A and Sentinel-2B) are operational the revisit time will reduce to 5 days, which together with Landsat-8, should enable sufficient cloud-free data for operational biomass estimation at a national scale. The Synthetic Aperture Radar Interferometry (InSAR) approach is feasible if there are enough coherent interferometric pairs available, however this is difficult to achieve due to the temporal decorrelation of the signal. For repeat-pass InSAR over a vegetated area even an 11 days temporal baseline is too large. In order to achieve better coherence a very high resolution is required at the cost of spatial coverage, which limits its scope for use in an operational context at a national scale. Future InSAR missions with pair acquisition in Tandem mode will minimize the temporal decorrelation over vegetation areas for more focused studies. The proposed approach complements the current paradigm of Big Data in Earth Observation, and illustrates the feasibility of integrating data from multiple sources. In future, this framework can be used to build an operational decision support system for retrieval of grassland biophysical parameters based on data from long term planned optical missions (e. g., Landsat, Sentinel) that will ensure the continuity of data acquisition. Similarly, Spanish X-band PAZ and TerraSAR-X2 missions will ensure the continuity of TerraSAR-X and COSMO-SkyMed.

Spring snowmelt and midwinter thaw and refreeze north of 60°N based on Seawinds QuikScat 2000-2009

The scatterometer SeaWinds on QuikSCAT provided regular measurements at Ku-band from 1999 to 2009. Although it was designed for ocean applications, it has been frequently used for the assessment of seasonal snowmelt patterns aside from other terrestrial applications such as ice cap monitoring, phenology and urban mapping. This paper discusses general data characteristics of SeaWinds and reviews relevant change detection algorithms. Depending on the complexity of the method, parameters such as long-term noise and multiple event analyses were incorporated. Temporal averaging is a commonly accepted preprocessing step with consideration of diurnal, multi-day or seasonal averages.

Wetlands in central Siberia from ASAR WS 2003, link to shapefiles

An area in central Siberia (partial coverage of Turukhansky und Yeniseysky districts) was investigated using satellite data. It covers freshwater ecosystems of non-forested peatlands in boreal forests. The satellite data represent the growing seasons of 2003/2004. Microwave data were acquired by the Advanced Synthetic Aperture Radar (ASAR) instrument onboard ENVISAT. The multi-temporal capabilities and resolution (150mx150m in WS mode) of the ASAR wide swath mode enabled the detection of dynamic features >2ha over this vast area. Scatterometer (QuikScat) data could be employed to distinguish hydro-periods. Wetland types have been identified on the basis of seasonal changes in backscatter. Results for peatlands have been compared with Russian forest inventory data which contain information on wetland distribution.

Methane fluxes and carbonate deposits at eastern Mediterranean cold seeps

High acoustic seafloor-backscatter signals characterize hundreds of patches of methane-derived authigenic carbonates and chemosynthetic communities associated with hydrocarbon seepage on the Nile Deep Sea Fan (NDSF) in the Eastern Mediterranean Sea. During a high-resolution ship-based multibeam survey covering a ~ 225 km**2 large seafloor area in the Central Province of the NDSF we identified 163 high-backscatter patches at water depths between 1500 and 1800 m, and investigated the source, composition, turnover, flux and fate of emitted hydrocarbons. Systematic Parasound single beam echosounder surveys of the water column showed hydroacoustic anomalies (flares), indicative of gas bubble streams, above 8% of the high-backscatter patches. In echosounder records flares disappeared in the water column close to the upper limit of the gas hydrate stability zone located at about 1350 m water depth due to decomposition of gas hydrate skins and subsequent gas dissolution. Visual inspection of three high-backscatter patches demonstrated that sediment cementation has led to the formation of continuous flat pavements of authigenic carbonates typically 100 to 300 m in diameter. Volume estimates, considering results from high-resolution autonomous underwater vehicle (AUV)-based multibeam mapping, were used to calculate the amount of carbonate-bound carbon stored in these slabs. Additionally, the flux of methane bubbles emitted at one high-backscatter patch was estimated (0.23 to 2.3 × 10**6 mol a**-1) by combined AUV flare mapping with visual observations by remotely operated vehicle (ROV). Another high-backscatter patch characterized by single carbonate pieces, which were widely distributed and interspaced with sediments inhabited by thiotrophic, chemosynthetic organisms, was investigated using in situ measurements with a benthic chamber and ex situ sediment core incubation and allowed for estimates of the methane consumption (0.1 to 1 × 10**6 mol a**-1) and dissolved methane flux (2 to 48 × 10**6 mol a**-1). Our comparison of dissolved and gaseous methane fluxes as well as methane-derived carbonate reservoirs demonstrates the need for quantitative assessment of these different methane escape routes and their interaction with the geo-, bio-, and hydrosphere at cold seeps.

Raw data of SCADCP (self-contained Acoustic Doppler Current Profiler) from three moorings in the Eastern Weddell Sea in 2005

Ten-month time series of mean volume backscattering strength (MVBS) and vertical velocity obtained from three moored acoustic Doppler current profilers (ADCPs) deployed from February until December 2005 at 64°S, 66.5°S and 69°S along the Greenwich Meridian were used to analyse the diel vertical zooplankton migration (DVM) and its seasonality and regional variability in the Lazarev Sea. The estimated MVBS exhibited distinct patterns of DVM at all three mooring sites. Between February and October, the timing of the DVM and the residence time of zooplankton at depth were clearly governed by the day-night rhythm. Mean daily cycles of the ADCP-derived vertical velocity were calculated for successive months and showed maximum ascent and descent velocities of 16 and -15 mm/s. However, a change of the MVBS pattern occurred in late spring/early austral summer (October/November), when the zooplankton communities ceased their synchronous vertical migration at all three mooring sites. Elevated MVBS values were then concentrated in the uppermost layers (<50 m) at 66.5°S. This period coincided with the decay of sea ice coverage at 64°S and 66.5°S between early November and mid-December. Elevated chlorophyll concentrations, which were measured at the end of the deployment, extended from 67°S to 65°S and indicated a phytoplankton bloom in the upper 50 m. Thus, we propose that the increased food supply associated with an ice edge bloom caused the zooplankton communities to cease their DVM in favour of feeding.

On the microstructure, growth pattern and original porosity of belemnite rostra: insights from calcitic Jurassic belemnites

Calcitic belemnite rostra are usually employed to perform paleoenvironmental studies based on geochemical data. However, several questions, such as their original porosity and microstructure, remain open, despite they are essential to make accurate interpretations based on geochemical analyses.This paper revisits and enlightens some of these questions. Petrographic data demonstrate that calcite crystals of the rostrum solidum of belemnites grow from spherulites that successively develop along the apical line, resulting in a “regular spherulithic prismatic” microstructure. Radially arranged calcite crystals emerge and diverge from the spherulites: towards the apex, crystals grow until a new spherulite is formed; towards the external walls of the rostrum, the crystals become progressively bigger and prismatic. Adjacent crystals slightly vary in their c-axis orientation, resulting in undulose extinction. Concentric growth layering develops at different scales and is superimposed and traversed by a radial pattern, which results in the micro-fibrous texture that is observed in the calcite crystals in the rostra.Petrographic data demonstrate that single calcite crystals in the rostra have a composite nature, which strongly suggests that the belemnite rostra were originally porous. Single crystals consistently comprise two distinct zones or sectors in optical continuity: 1) the inner zone is fluorescent, has relatively low optical relief under transmitted light (TL) microscopy, a dark-grey color under backscatter electron microscopy (BSEM), a commonly triangular shape, a “patchy” appearance and relatively high Mg and Na contents; 2) the outer sector is non-fluorescent, has relatively high optical relief under TL, a light-grey color under BSEM and low Mg and Na contents. The inner and fluorescent sectors are interpreted to have formed first as a product of biologically controlled mineralization during belemnite skeletal growth and the non-fluorescent outer sectors as overgrowths of the former, filling the intra- and inter-crystalline porosity. This question has important implications for making paleoenvironmental and/or paleoclimatic interpretations based on geochemical analyses of belemnite rostra.Finally, the petrographic features of composite calcite crystals in the rostra also suggest the non-classical crystallization of belemnite rostra, as previously suggested by other authors.

Microstructure characterization of 316L deformed at high strain rates using EBSD

Specimens from split Hopkinson pressure bar experiments, at strain rates between ~ 1000–9000 s− 1 at room temperature and 500 °C, have been studied using electron backscatter diffraction. No significant differences in the microstructures were observed at different strain rates, but were observed for different strains and temperatures. Size distribution for subgrains with boundary misorientations > 2° can be described as a bimodal lognormal area distribution. The distributions were found to change due to deformation. Part of the distribution describing the large subgrains decreased while the distribution for the small subgrains increased. This is in accordance with deformation being heterogeneous and successively spreading into the undeformed part of individual grains. The variation of the average size for the small subgrain distribution varies with strain but not with strain rate in the tested interval. The mean free distance for dislocation slip, interpreted here as the average size of the distribution of small subgrains, displays a variation with plastic strain which is in accordance with the different stages in the stress-strain curves. The rate of deformation hardening in the linear hardening range is accurately calculated using the variation of the small subgrain size with strain.

Validation of AltiKa Matching Pursuit Rain Flag

The major drawback of Ka band, operating frequency of the AltiKa altimeter on board SARAL, is its sensitivity to atmospheric liquid water. Even light rain or heavy clouds can strongly attenuate the signal and distort the signal leading to erroneous geophysical parameters estimates. A good detection of the samples affected by atmospheric liquid water is crucial. As AltiKa operates at a single frequency, a new technique based on the detection by a Matching Pursuit algorithm of short scale variations of the slope of the echo waveform plateau has been developed and implemented prelaunch in the ground segment. As the parameterization of the detection algorithm was defined using Jason-1 data, the parameters were re-estimated during the cal-val phase, during which the algorithm was also updated. The measured sensor signal-to-noise ratio is significantly better than planned, the data loss due to attenuation by rain is significantly smaller than expected (<0.1%). For cycles 2 to 9, the flag detects about 9% of 1Hz data, 5.5% as rainy and 3.5 % as backscatter bloom (or sigma0 bloom). The results of the flagging process are compared to independent rain data from microwave radiometers to evaluate its performances in term of detection and false alarms.

Assessment of SARAL/AltiKa Wave Height Measurements Relative to Buoy, Jason-2, and Cryosat-2 Data

SARAL/AltiKa GDR-T are analyzed to assess the quality of the significant wave height (SWH) measurements. SARAL along-track SWH plots reveal cases of erroneous data, more or less isolated, not detected by the quality flags. The anomalies are often correlated with strong attenuation of the Ka-band backscatter coefficient, sensitive to clouds and rain. A quality test based on the 1Hz standard deviation is proposed to detect such anomalies. From buoy comparison, it is shown that SARAL SWH is more accurate than Jason-2, particularly at low SWH, and globally does not require any correction. Results are better with open ocean than with coastal buoys. The scatter and the number of outliers are much larger for coastal buoys. SARAL is then compared with Jason-2 and Cryosat-2. The altimeter data are extracted from the global altimeter SWH Ifremer data base, including specific corrections to calibrate the various altimeters. The comparison confirms the high quality of SARAL SWH. The 1Hz standard deviation is much less than for Jason-2 and Cryosat-2, particularly at low SWH. Furthermore, results show that the corrections applied to Jason-2 and to Cryosat-2, in the data base, are efficient, improving the global agreement between the three altimeters.

A new technique to automatically quantify microstructures of fine grained carbonate mylonites: two-step etching combined with SEM imaging and image analysis

A two-step etching technique for fine-grained calcite mylonites using 0.37% hydrochloric and 0.1% acetic acid produces a topographic relief which reflects the grain boundary geometry. With this technique, calcite grain boundaries become more intensely dissolved than their grain interiors but second phase minerals like dolomite, quartz, feldspars, apatite, hematite and pyrite are not affected by the acid and therefore form topographic peaks. Based on digital backscatter electron images and element distribution maps acquired on a scanning electron microscope, the geometry of calcite and the second phase minerals can be automatically quantified using image analysis software. For research on fine-grained carbonate rocks (e.g. dolomite calcite mixtures), this low-cost approach is an attractive alternative to the generation of manual grain boundary maps based on photographs from ultra-thin sections or orientation contrast images.

DopSCAT: A mission concept for simultaneous measurements of marine winds and surface currents

A radar scatterometer operates by transmitting a pulse of microwave energy toward the ocean's surface and measuring the normalized (per-unit-surface) radar backscatter coefficient (σ°). The primary application of scatterometry is the measurement of near-surface ocean winds. By combining σ° measurements from different azimuth angles, the 10 m vector wind can be determined through a Geophysical Model Function (GMF), which relates wind and backscatter. This paper proposes a mission concept for the measurement of both oceanic winds and surface currents, which makes full use of earlier C-band radar remote sensing experience. For the determination of ocean currents, in particular, the novel idea of using two chirps of opposite slope is introduced. The fundamental processing steps required to retrieve surface currents are given together with their associated accuracies. A detailed description of the mission proposal and comparisons between real and retrieved surface currents are presented. The proposed ocean Doppler scatterometer can be used to generate global surface ocean current maps with accuracies better than 0.2 m/s at a spatial resolution better than 25 km (i.e., 12.5 km spatial sampling) on a daily basis. These maps will allow gaining some insights on the upper ocean mesoscale dynamics. The work lies at a frontier, given that the present inability to measure ocean currents from space in a consistent and synoptic manner represents one of the greatest weaknesses in ocean remote sensing.

Altimeter estimation of sea surface wind stress for light to moderate winds

Aircraft altimeter and in situ measurements are used to examine relationships between altimeter backscatter and the magnitude of near-surface wind and friction velocities. Comparison of altimeter radar cross section with wind speed is made through the modified Chelton-Wentz algorithm. Improved agreement is found after correcting 10-m winds for both surface current and atmospheric stability. An altimeter friction velocity algorithm is derived based on the wind speed model and an open-ocean drag coefficient. Close agreement between altimeter- and in situ-derived friction velocities is found. For this dataset, quality of the altimeter inversion to surface friction velocity is comparable to that for adjusted winds and clearly better than the inversion to true 10-m wind speed.

Effect of ray and speed perturbations on ionospheric tomography by over-the-horizon radar: A new method

International audience

Effects of microstructure on the variation of the unloading behavior of DP780 steels

The nonlinear unloading behavior of three different commercial dual-phase steels (DP780 grade equivalent) was examined. These steels exhibited small variations in chemical composition (0.07 to 0.10 mass percent carbon) and martensite volume fraction (0.23 to 0.28), and they demonstrated similar hardening behavior. Uniaxial loading-unloading-loading tests were conducted at room temperature and quasi-static strain rates between engineering strains of 0.5 and 8%. Steel microstructures were examined using electron backscatter diffraction and nanoindentation techniques. The microplastic component of the unloading strain exhibited no dependence on the martensite volume fraction or the ferrite grain size within the small range encountered in this investigations. Instead, the magnitude of the microplastic component of the unloading strain increased as the strength ratio between the martensite and ferrite phases increased. Correspondingly, the apparent unloading modulus, or chord modulus, exhibited a greater reduction for equivalent increments of strain hardening as the strength ratio increased. These results suggest that springback can be reduced in structures containing two ductile phases if the strength ratio between the harder and softer phases is reduced.