On independent component analysis applied to unmixing hyperspectral data

One of the most challenging task underlying many hyperspectral imagery applications is the spectral unmixing, which decomposes a mixed pixel into a collection of reectance spectra, called endmember signatures, and their corresponding fractional abundances. Independent Component Analysis (ICA) have recently been proposed as a tool to unmix hyperspectral data. The basic goal of ICA is to nd a linear transformation to recover independent sources (abundance fractions) given only sensor observations that are unknown linear mixtures of the unobserved independent sources. In hyperspectral imagery the sum of abundance fractions associated to each pixel is constant due to physical constraints in the data acquisition process. Thus, sources cannot be independent. This paper address hyperspectral data source dependence and its impact on ICA performance. The study consider simulated and real data. In simulated scenarios hyperspectral observations are described by a generative model that takes into account the degradation mechanisms normally found in hyperspectral applications. We conclude that ICA does not unmix correctly all sources. This conclusion is based on the a study of the mutual information. Nevertheless, some sources might be well separated mainly if the number of sources is large and the signal-to-noise ratio (SNR) is high.

Circumsolar energetic particle distribution on 2011 November 3

Late on 2011 November 3, STEREO-A, STEREO-B, MESSENGER, and near-Earth spacecraft observed an energetic particle flux enhancement. Based on the analysis of in situ plasma and particle observations, their correlation with remote sensing observations, and an interplanetary transport model, we conclude that the particle increases observed at multiple locations had a common single source active region and the energetic particles filled a very broad region around the Sun. The active region was located at the solar backside (as seen from Earth) and was the source of a large flare, a fast and wide coronal mass ejection, and an EIT wave, accompanied by type II and type III radio-emission. In contrast to previous solar energetic particle events showing broad longitudinal spread, this event showed clear particle anisotropies at three widely separated observation points at 1AU, suggesting direct particle injection close to the magnetic footpoint of each spacecraft, lasting for several hours.We discuss these observations and the possible scenarios explaining the extremely broad particle spread for this event.

Multisensor Modeling Underwater with Uncertain Information

This thesis develops an approach to the construction of multidimensional stochastic models for intelligent systems exploring an underwater environment. It describes methods for building models by a three- dimensional spatial decomposition of stochastic, multisensor feature vectors. New sensor information is incrementally incorporated into the model by stochastic backprojection. Error and ambiguity are explicitly accounted for by blurring a spatial projection of remote sensor data before incorporation. The stochastic models can be used to derive surface maps or other representations of the environment. The methods are demonstrated on data sets from multibeam bathymetric surveying, towed sidescan bathymetry, towed sidescan acoustic imagery, and high-resolution scanning sonar aboard a remotely operated vehicle.

gvSIG Sensores

El protocolo SOS (Sensor Observation Service) es una especificación OGC dentro de la iniciativa Sensor Web Enablement (SWE), que permite acceder a las observaciones y datos de sensores heterogéneos de una manera estándar. En el proyecto gvSIG se ha abierto una línea de investigación entorno a la SWE, existiendo en la actualidad dos prototipos de clientes SOS para gvSIG y gvSIG Mobile. La especificación utilizada para describir las medidas proporcionadas por sensores es Observation & Measurement (O&M) y la descripción de los metadatos de los sensores (localización. ID, fenómenos medidos, procesamiento de los datos, etc) se obtiene a partir del esquema Sensor ML. Se ha implementado el siguiente conjunto de operaciones: GetCapabilities para la descripción del servicio; DescribeSensor para acceder a los metadatos del sensor y el GetObservation para recibir las observaciones. En el caso del prototipo para gvSIG escritorio se puede acceder a los datos procedentes de los distintos grupos de sensores “offerings” añadiéndolos en el mapa como nuevas capas. Los procedimientos o sensores que están incluidos en un “offering” son presentados como elementos de la capa que se pueden cartografiar en el mapa. Se puede acceder a las observaciones (GetObservation) de estos sensores filtrando los datos por intervalo de tiempo y propiedad del fenómeno observado. La información puede ser representada sobre el mapa mediante gráficas para una mejor comprensión con la posibilidad de comparar datos de distintos sensores. En el caso del prototipo para el cliente móvil gvSIG Mobile, se ha utilizado la misma filosofía que para el cliente de escritorio, siendo cada “offering” una nueva capa. Las observaciones de los sensores pueden ser visualizadas en la pantalla del dispositivo móvil y se pueden obtener mapas temáticos,con el objetivo de facilitar la interpretación de los datos

Can desert dust explain the outgoing longwave radiation anomaly over the Sahara during July 2003?

Measurements of the top‐of‐the‐atmosphere outgoing longwave radiation (OLR) for July 2003 from Meteosat‐7 are used to assess the performance of the numerical weather prediction version of the Met Office Unified Model. A significant difference is found over desert regions of northern Africa where the model emits too much OLR by up to 35 Wm−2 in the monthly mean. By cloud‐screening the data we find an error of up to 50 Wm−2 associated with cloud‐free areas, which suggests an error in the model surface temperature, surface emissivity, or atmospheric transmission. By building up a physical model of the radiative properties of mineral dust based on in situ, and surface‐based and satellite remote sensing observations we show that the most plausible explanation for the discrepancy in OLR is due to the neglect of mineral dust in the model. The calculations suggest that mineral dust can exert a longwave radiative forcing by as much as 50 Wm−2 in the monthly mean for 1200 UTC in cloud‐free regions, which accounts for the discrepancy between the model and the Meteosat‐7 observations. This suggests that inclusion of the radiative effects of mineral dust will lead to a significant improvement in the radiation balance of numerical weather prediction models with subsequent improvements in performance.

Characteristic magnetic field and speed properties of interplanetary coronal mass ejections and their sheath regions

Prediction of the solar wind conditions in near-Earth space, arising from both quasi-steady and transient structures, is essential for space weather forecasting. To achieve forecast lead times of a day or more, such predictions must be made on the basis of remote solar observations. A number of empirical prediction schemes have been proposed to forecast the transit time and speed of coronal mass ejections (CMEs) at 1 AU. However, the current lack of magnetic field measurements in the corona severely limits our ability to forecast the 1 AU magnetic field strengths resulting from interplanetary CMEs (ICMEs). In this study we investigate the relation between the characteristic magnetic field strengths and speeds of both magnetic cloud and noncloud ICMEs at 1 AU. Correlation between field and speed is found to be significant only in the sheath region ahead of magnetic clouds, not within the clouds themselves. The lack of such a relation in the sheaths ahead of noncloud ICMEs is consistent with such ICMEs being skimming encounters of magnetic clouds, though other explanations are also put forward. Linear fits to the radial speed profiles of ejecta reveal that faster-traveling ICMEs are also expanding more at 1 AU. We combine these empirical relations to form a prediction scheme for the magnetic field strength in the sheaths ahead of magnetic clouds and also suggest a method for predicting the radial speed profile through an ICME on the basis of upstream measurements.

From the sun to the Earth: The 13 May 2005 coronal mass ejection

We report the results of a multi-instrument, multi-technique, coordinated study of the solar eruptive event of 13 May 2005. We discuss the resultant Earth-directed (halo) coronal mass ejection (CME), and the effects on the terrestrial space environment and upper Earth atmosphere. The interplanetary CME (ICME) impacted the Earth’s magnetosphere and caused the most-intense geomagnetic storm of 2005 with a Disturbed Storm Time (Dst) index reaching −263 nT at its peak. The terrestrial environment responded to the storm on a global scale. We have combined observations and measurements from coronal and interplanetary remote-sensing instruments, interplanetary and near-Earth in-situ measurements, remote-sensing observations and in-situ measurements of the terrestrial magnetosphere and ionosphere, along with coronal and heliospheric modelling. These analyses are used to trace the origin, development, propagation, terrestrial impact, and subsequent consequences of this event to obtain the most comprehensive view of a geo-effective solar eruption to date. This particular event is also part of a NASA-sponsored Living With a Star (LWS) study and an on-going US NSF-sponsored Solar, Heliospheric, and INterplanetary Environment (SHINE) community investigation.

Predictability of biomass burning in response to climate changes

Climate is an important control on biomass burning, but the sensitivity of fire to changes in temperature and moisture balance has not been quantified. We analyze sedimentary charcoal records to show that the changes in fire regime over the past 21,000 yrs are predictable from changes in regional climates. Analyses of paleo- fire data show that fire increases monotonically with changes in temperature and peaks at intermediate moisture levels, and that temperature is quantitatively the most important driver of changes in biomass burning over the past 21,000 yrs. Given that a similar relationship between climate drivers and fire emerges from analyses of the interannual variability in biomass burning shown by remote-sensing observations of month-by-month burnt area between 1996 and 2008, our results signal a serious cause for concern in the face of continuing global warming.

Information-based data selection for ensemble data assimilation

Ensemble-based data assimilation is rapidly proving itself as a computationally-efficient and skilful assimilation method for numerical weather prediction, which can provide a viable alternative to more established variational assimilation techniques. However, a fundamental shortcoming of ensemble techniques is that the resulting analysis increments can only span a limited subspace of the state space, whose dimension is less than the ensemble size. This limits the amount of observational information that can effectively constrain the analysis. In this paper, a data selection strategy that aims to assimilate only the observational components that matter most and that can be used with both stochastic and deterministic ensemble filters is presented. This avoids unnecessary computations, reduces round-off errors and minimizes the risk of importing observation bias in the analysis. When an ensemble-based assimilation technique is used to assimilate high-density observations, the data-selection procedure allows the use of larger localization domains that may lead to a more balanced analysis. Results from the use of this data selection technique with a two-dimensional linear and a nonlinear advection model using both in situ and remote sounding observations are discussed.

Data assimilation with correlated observation errors: experiments with a 1-D shallow water model

Remote sensing observations often have correlated errors, but the correlations are typically ignored in data assimilation for numerical weather prediction. The assumption of zero correlations is often used with data thinning methods, resulting in a loss of information. As operational centres move towards higher-resolution forecasting, there is a requirement to retain data providing detail on appropriate scales. Thus an alternative approach to dealing with observation error correlations is needed. In this article, we consider several approaches to approximating observation error correlation matrices: diagonal approximations, eigendecomposition approximations and Markov matrices. These approximations are applied in incremental variational assimilation experiments with a 1-D shallow water model using synthetic observations. Our experiments quantify analysis accuracy in comparison with a reference or ‘truth’ trajectory, as well as with analyses using the ‘true’ observation error covariance matrix. We show that it is often better to include an approximate correlation structure in the observation error covariance matrix than to incorrectly assume error independence. Furthermore, by choosing a suitable matrix approximation, it is feasible and computationally cheap to include error correlation structure in a variational data assimilation algorithm.

Opportunities for magnetospheric research using EISCAT/ESR and CLUSTER

The four Cluster spacecraft offer a unique opportunity to study structure and dynamics in the magnetosphere and we discuss four general ways in which ground-based remote-sensing observations of the ionosphere can be used to support the in-situ measurements. The ionosphere over the Svalbard islands will be studied in particular detail, not only by the ESR and EISCAT incoherent scatter radars, but also by optical instruments, magnetometers, imaging riometers and the CUTLASS bistatic HF radar. We present an on-line procedure to plan coordinated measurements by the Cluster spacecraft with these combined ground-based systems. We illustrate the philosophy of the method, using two important examples of the many possible configurations between the Cluster satellites and the ground-based instruments.

Non-Maxwellian ion velocity distributions observed using EISCAT

Recent observations from the EISCAT incoherent scatter radar have revealed bursts of poleward ion flow in the dayside auroral ionosphere which are consistent with the ionospheric signature of flux transfer events at the magnetopause. These bursts frequently contain ion drifts which exceed the neutral thermal speed and, because the neutral thermospheric wind is incapable of responding sufficiently rapidly, toroidal, non-Maxwellian ion velocity distributions are expected. The EISCAT observations are made with high time resolution (15 seconds) and at a large angle to the geomagnetic field (73.5°), allowing the non-Maxwellian nature of the distribution to be observed remotely for the first time. The observed features are also strongly suggestive of a toroidal distribution: characteristic spectral shape, increased scattered power (both consistent with reduced Landau damping and enhanced electric field fluctuations) and excessively high line-of-sight ion temperatures deduced if a Maxwellian distribution is assumed. These remote sensing observations allow the evolution of the distributions to be observed. They are found to be non-Maxwellian whenever the ion drift exceeds the neutral thermal speed, indicating that such distributions can exist over the time scale of the flow burst events (several minutes).

ERA-Interim/Land: a global land surface reanalysis data set

ERA-Interim/Land is a global land surface reanalysis data set covering the period 1979–2010. It describes the evolution of soil moisture, soil temperature and snowpack. ERA-Interim/Land is the result of a single 32-year simulation with the latest ECMWF (European Centre for Medium-Range Weather Forecasts) land surface model driven by meteorological forcing from the ERA-Interim atmospheric reanalysis and precipitation adjustments based on monthly GPCP v2.1 (Global Precipitation Climatology Project). The horizontal resolution is about 80 km and the time frequency is 3-hourly. ERA-Interim/Land includes a number of parameterization improvements in the land surface scheme with respect to the original ERA-Interim data set, which makes it more suitable for climate studies involving land water resources. The quality of ERA-Interim/Land is assessed by comparing with ground-based and remote sensing observations. In particular, estimates of soil moisture, snow depth, surface albedo, turbulent latent and sensible fluxes, and river discharges are verified against a large number of site measurements. ERA-Interim/Land provides a global integrated and coherent estimate of soil moisture and snow water equivalent, which can also be used for the initialization of numerical weather prediction and climate models.

Evolução tectono-estrutural do Campo de Xaréu (Sub-bacia de Mundaú, Bacia do Ceará - NE do Brasil: abordagem multiescala e pluriferramental

The Xaréu Oil Field, located in the center-southern portion of the Mundaú Sub-Basin (eastern portion of the Ceará Basin), is characterized by a main Iramework of NW-trending and NE-dipping faults. The faults in the Xaréu Oil Field, among which the Xaréu Fautt stands out, are arranged according to an extensional-listriclan, rooted on a detachment surface corresponding to the Mundaú Fault, the border fautt of Mundaú Sub-Basin. During the tectonic-structural evolution of the Xaréu Oil Field and the Mundaú Sub-Basin, the Mundaú Fault played a crucial role on the control of the geometry of both compartments. The main carbonatic unit in the Xaréu Oil Field, named the Trairí Member(Paracuru Formation of Late Aptian to Early Albian age), contains the largest oil volume in the field, concentrated in structurally-controlled accumulations. The Trairí Member is composed by a variety of carbonatic rocks (massive, bedded or laminated calcilutites, ostracodites, calcarenites and carbonatic rudites, all of them presenting variable degrees of dolomitization). The carbonatic rocks are interbedded into thick packages of black shales and marls, besides local beds of siliciclastic conglomerates, sandstones, siltnes and argillites. From the spatial association and the genetic relationships between the carbonatic and siliciclastic units, it is possible to group them in three lithofacies associations (Marginal Plain, Ramp and Lacustrine Interior) that, together, were developed in a lacustrine system associated to a marginal sabkha. Structural studies based on drill coresthat sample the Trairí Member in the Xaréu Oil Field allowed to characterize two generations of meso- to microscale structures: the D1 group presents a typical hydroplastic character, being characterized by intra/interstratal to oblique-bedding shear zones. The hydroplastic character related to these structures allowed to infer their development at an early-lithilication stage of the Trairí Member, leading to infer an Early Cretaceous age to them. The second group of structures identified in the drill cores, nominated D2 and ascribed to a Neogene age, presents a strictly brttle character, being typilied by normal faults and slickenfibers of re-crystallized clayminerals, ali olthem displaying variable orientations. Although the present faults in the Xaréu Oil Field (and, consequently, in the Mundaú Sub-Basin) were classically relerred as struetures of essentially normal displacement, the kinematics analysis of the meso-to microscaie D1 struetures in the drill cores led to deline oblique displacements (normal with a clockwise strike-slip component) to these faults, indicating a main tectonic transport to ENE. These oblique movements would be responsible for the installation of a transtensive context in the Mundaú Sub-Basin, as part of the transcurrent to translormant opening of the Atlantic Equatorial Margin. The balancing of four struetural cross-sections ofthe Xaréu Oil Field indicates that the Mundaú Fault was responsible for more than 50% of the total stretching (ß factor) registered during the Early Aptian. At the initial stages of the "rifting", during Early Aptianuntil the Holocene, the Mundaú Sub-Basin (and consequently the Xaréu Oil Fleld) accumulated a total stretching between 1.21 and 1.23; in other words, the crust in this segment of the Atlantic Equatorial Margin was subjeeted to an elongation of about 20%. From estimates of oblique displacements related to the faults, it ws possible to construct diagrams that allow the determination of stretching factors related to these displacements. Using these diagrams and assuming the sense 01 dominant teetonictransport towards ENE, it was possible to calculate the real stretching lactors related to the oblique movement 0 of the faults in the Mundaú Sub-Basin. which reached actual values between 1.28 and 1.42. ln addnion to the tectonic-structural studies in the Xaréu Oil Field, the interpretation of remote sensing products, coupled wnh characterization of terrain analogues in seleeted areas along the northern Ceará State (continental margins of the Ceará and Potiguar basins), provided addnional data and constraints about the teetonic-structural evolution of the oil lield. The work at the analogue sites was particularly effective in the recognition and mapping, in semidetail scale, several generations of struetures originated under a brittle regime. Ali the obtained information (from the Xaréu Oil Field, the remote sensor data and the terrain analogues) were jointly interpreted, culminating with the proposnion of an evolutionary model lor this segment of the Atlantic Equatorial Margin; this model that can be applied to the whole Margin, as well. This segmentof the Atlantic Equatorial Margin was delormedin an early E-W (when considered lhe present-day position of the South American Plate) transcurrent to transform regime with dextral kinematics, started Irom, at least, the Early Aptian, which left its record in several outcrops along the continental margin of the Ceará State and specilically in the Xaréu off Field. The continuous operation of the regime, through the Albian and later periods, led to the definitive separation between the South American and African plates, with the formation of oceanic lithosphere between the two continental blocks, due to the emplacement off spreading centers. This process involved the subsequent transition of the transcurrent to a translorm dextral regime, creating lhe Equatorial Atlantic Oceano With the separation between the South American and African plates already completed and the increasing separation between lhe continental masses, other tecton ic mechanisms began to act during the Cenozoic (even though the Cretaceous tectonic regime lasted until the Neogene), like an E-W compressive stress líeld (related to the spreading olthe oceanic floor along lhe M id-Atlantic Ridge and to the compression of the Andean Chain) effective Irom the Late Cretaceous, and a state of general extension olthe horizontal surface (due to the thermal uplift ofthe central portion of Borborema Province), effective during the Neogene. The overlap of these mechanisms during the Cenozoic led to the imprint of a complex tectonic framework, which apparently influenced the migration and entrapment 01 hydrocarbon in the Ceará Basin

Satellite and in situ data integrated analysis to study the upper ocean an coastal environment of the Italian seas

The thesis objectives are to develop new methodologies for study of the space and time variability of Italian upper ocean ecosystem through the combined use of multi-sensors satellite data and in situ observations and to identify the capability and limits of remote sensing observations to monitor the marine state at short and long time scales. Three oceanographic basins have been selected and subjected to different types of analyses. The first region is the Tyrrhenian Sea where a comparative analysis of altimetry and lagrangian measurements was carried out to study the surface circulation. The results allowed to deepen the knowledge of the Tyrrhenian Sea surface dynamics and its variability and to defined the limitations of satellite altimetry measurements to detect small scale marine circulation features. Channel of Sicily study aimed to identify the spatial-temporal variability of phytoplankton biomass and to understand the impact of the upper ocean circulation on the marine ecosystem. An combined analysis of the satellite of long term time series of chlorophyll, Sea Surface Temperature and Sea Level field data was applied. The results allowed to identify the key role of the Atlantic water inflow in modulating the seasonal variability of the phytoplankton biomass in the region. Finally, Italian coastal marine system was studied with the objective to explore the potential capability of Ocean Color data in detecting chlorophyll trend in coastal areas. The most appropriated methodology to detect long term environmental changes was defined through intercomparison of chlorophyll trends detected by in situ and satellite. Then, Italian coastal areas subject to eutrophication problems were identified. This work has demonstrated that satellites data constitute an unique opportunity to define the features and forcing influencing the upper ocean ecosystems dynamics and can be used also to monitor environmental variables capable of influencing phytoplankton productivity.