Constraints implementation for IQML and mode direction-of-arrival estimators

The iterative quadratic maximum likelihood IQML and the method of direction estimation MODE are well known high resolution direction-of-arrival DOA estimation methods. Their solutions lead to an optimization problem with constraints. The usual linear constraint presents a poor performance for certain DOA values. This work proposes a new linear constraint applicable to both DOA methods and compare their performance with two others: unit norm and usual linear constraint. It is shown that the proposed alternative performs better than others constraints. The resulting computational complexity is also investigated.

A search for anisotropy in the arrival directions of ultra high energy cosmic rays recorded at the Pierre Auger Observatory

Observations of cosmic rays arrival directions made with the Pierre Auger Observatory have previously provided evidence of anisotropy at the 99% CL using the correlation of ultra high energy cosmic rays (UHECRs) with objects drawn from the Veron-Cetty Veron catalog. In this paper we report on the use of three catalog independent methods to search for anisotropy. The 2pt-L, 2pt+ and 3pt methods, each giving a different measure of self-clustering in arrival directions, were tested on mock cosmic ray data sets to study the impacts of sample size and magnetic smearing on their results, accounting for both angular and energy resolutions. If the sources of UHECRs follow the same large scale structure as ordinary galaxies in the local Universe and if UHECRs are deflected no more than a few degrees, a study of mock maps suggests that these three method can efficiently respond to the resulting anisotropy with a P-value = 1.0% or smaller with data sets as few as 100 events. using data taken from January 1, 2004 to July 31, 2010 we examined the 20, 30, ... , 110 highest energy events with a corresponding minimum energy threshold of about 49.3 EeV. The minimum P-values found were 13.5% using the 2pt-L method, 1.0% using the 2pt+ method and 1.1% using the 3pt method for the highest 100 energy events. In view of the multiple (correlated) scans performed on the data set, these catalog-independent methods do not yield strong evidence of anisotropy in the highest energy cosmic rays.

LARGE-SCALE DISTRIBUTION OF ARRIVAL DIRECTIONS OF COSMIC RAYS DETECTED ABOVE 10(18) eV AT THE PIERRE AUGER OBSERVATORY

A thorough search for large-scale anisotropies in the distribution of arrival directions of cosmic rays detected above 10(18) eV at the Pierre Auger Observatory is presented. This search is performed as a function of both declination and right ascension in several energy ranges above 10(18) eV, and reported in terms of dipolar and quadrupolar coefficients. Within the systematic uncertainties, no significant deviation from isotropy is revealed. Assuming that any cosmic-ray anisotropy is dominated by dipole and quadrupole moments in this energy range, upper limits on their amplitudes are derived. These upper limits allow us to test the origin of cosmic rays above 10(18) eV from stationary Galactic sources densely distributed in the Galactic disk and predominantly emitting light particles in all directions.

HIGH PERFORMANCE, LOW COST SUBSPACE DECOMPOSITION AND POLYNOMIAL ROOTING FOR REAL TIME DIRECTION OF ARRIVAL ESTIMATION: ANALYSIS AND IMPLEMENTATION

This thesis develops high performance real-time signal processing modules for direction of arrival (DOA) estimation for localization systems. It proposes highly parallel algorithms for performing subspace decomposition and polynomial rooting, which are otherwise traditionally implemented using sequential algorithms. The proposed algorithms address the emerging need for real-time localization for a wide range of applications. As the antenna array size increases, the complexity of signal processing algorithms increases, making it increasingly difficult to satisfy the real-time constraints. This thesis addresses real-time implementation by proposing parallel algorithms, that maintain considerable improvement over traditional algorithms, especially for systems with larger number of antenna array elements. Singular value decomposition (SVD) and polynomial rooting are two computationally complex steps and act as the bottleneck to achieving real-time performance. The proposed algorithms are suitable for implementation on field programmable gated arrays (FPGAs), single instruction multiple data (SIMD) hardware or application specific integrated chips (ASICs), which offer large number of processing elements that can be exploited for parallel processing. The designs proposed in this thesis are modular, easily expandable and easy to implement. Firstly, this thesis proposes a fast converging SVD algorithm. The proposed method reduces the number of iterations it takes to converge to correct singular values, thus achieving closer to real-time performance. A general algorithm and a modular system design are provided making it easy for designers to replicate and extend the design to larger matrix sizes. Moreover, the method is highly parallel, which can be exploited in various hardware platforms mentioned earlier. A fixed point implementation of proposed SVD algorithm is presented. The FPGA design is pipelined to the maximum extent to increase the maximum achievable frequency of operation. The system was developed with the objective of achieving high throughput. Various modern cores available in FPGAs were used to maximize the performance and details of these modules are presented in detail. Finally, a parallel polynomial rooting technique based on Newton’s method applicable exclusively to root-MUSIC polynomials is proposed. Unique characteristics of root-MUSIC polynomial’s complex dynamics were exploited to derive this polynomial rooting method. The technique exhibits parallelism and converges to the desired root within fixed number of iterations, making this suitable for polynomial rooting of large degree polynomials. We believe this is the first time that complex dynamics of root-MUSIC polynomial were analyzed to propose an algorithm. In all, the thesis addresses two major bottlenecks in a direction of arrival estimation system, by providing simple, high throughput, parallel algorithms.

Parametric estimation of pulse arrival time: a robust approach to pulse wave velocity

Pulse wave velocity (PWV) is a surrogate of arterial stiffness and represents a non-invasive marker of cardiovascular risk. The non-invasive measurement of PWV requires tracking the arrival time of pressure pulses recorded in vivo, commonly referred to as pulse arrival time (PAT). In the state of the art, PAT is estimated by identifying a characteristic point of the pressure pulse waveform. This paper demonstrates that for ambulatory scenarios, where signal-to-noise ratios are below 10 dB, the performance in terms of repeatability of PAT measurements through characteristic points identification degrades drastically. Hence, we introduce a novel family of PAT estimators based on the parametric modeling of the anacrotic phase of a pressure pulse. In particular, we propose a parametric PAT estimator (TANH) that depicts high correlation with the Complior(R) characteristic point D1 (CC = 0.99), increases noise robustness and reduces by a five-fold factor the number of heartbeats required to obtain reliable PAT measurements.

Sequence of arrival determines plant-mediated interactions between herbivores

Induced changes in plant quality can mediate indirect interactions between herbivores. Although the sequence of attack by different herbivores has been shown to influence plant responses, little is known about how this affects the herbivores themselves. We therefore investigated how induction by the leaf herbivore Spodoptera frugiperda influences resistance of teosinte (Zea mays mexicana) and cultivated maize (Zea mays mays) against root-feeding larvae of Diabrotica virgifera virgifera. The importance of the sequence of arrival was tested in the field and laboratory. Spodoptera frugiperda infestation had a significant negative effect on colonization by D. virgifera larvae in the field and weight gain in the laboratory, but only when S. frugiperda arrived on the plant before the root herbivore. When S. frugiperda arrived after the root herbivore had established, no negative effects on larval performance were detected. Yet, adult emergence of D. virgifera was reduced even when the root feeder had established first, indicating that the negative effects were not entirely absent in this treatment. The defoliation of the plants was not a decisive factor for the negative effects on root herbivore development, as both minor and major leaf damage resulted in an increase in root resistance and the extent of biomass removal was not correlated with root-herbivore growth. We propose that leaf-herbivore-induced increases in feeding-deterrent and/or toxic secondary metabolites may account for the sequence-specific reduction in root-herbivore performance. Synthesis. Our results demonstrate that the sequence of arrival can be an important determinant of plant-mediated interactions between insect herbivores in both wild and cultivated plants. Arriving early on a plant may be an important strategy of insects to avoid competition with other herbivores. To fully understand plant-mediated interactions between insect herbivores, the sequence of arrival should be taken into account. © 2011 The Authors. Journal of Ecology © 2011 British Ecological Society.

An open letter addressed to Sir Moses Montefiore ... on the day of his arrival in the holy city of Jerusalem

by Meyer Auerbach and Samuel Salant

Path characteristics and colony arrival-departure statistics of king penguins (Aptenodytes patagonicus) on Possession and Kerguelen Island

Orientation based on visual cues can be extremely difficult in crowded bird colonies due to the presence of many individuals. We studied king penguins (Aptenodytes patagonicus) that live in dense colonies and are constantly faced with such problems. Our aims were to describe adult penguin homing paths on land and to test whether visual cues are important for their orientation in the colony. We also tested the hypothesis that older penguins should be better able to cope with limited visual cues due to their greater experience. We collected and examined GPS paths of homing penguins. In addition, we analyzed 8 months of penguin arrivals to and departures from the colony using data from an automatic identification system. We found that birds rearing chicks did not minimize their traveling time on land and did not proceed to their young (located in creches) along straight paths. Moreover, breeding birds' arrivals and departures were affected by the time of day and luminosity levels. Our data suggest that king penguins prefer to move in and out of the colony when visual cues are available. Still, they are capable of navigating even in complete darkness, and this ability seems to develop over the years, with older breeding birds more likely to move through the colony at nighttime luminosity levels. This study is the first step in unveiling the mysteries of king penguin orientation on land.

Fig. S1. Loci of possible event locations based on differential ts - tp arrival times at the single functioning geophone

Geological storage of CO2 that has been captured at large, point source emitters represents a key potential method for reduction of anthropogenic greenhouse gas emissions. However, this technology will only be viable if it can be guaranteed that injected CO2 will remain trapped in the subsurface for thousands of years or more. A signi?cant issue for storage security is the geomechanical response of the reservoir. Concerns have been raised that geomechanical deformation induced by CO2 injection will create or reactivate fracture networks in the sealing caprocks, providing a pathway for CO2 leakage. In this paper, we examine three large-scale sites where CO2 is injected at rates of ab. 1 megatonne/y or more: Sleipner, Weyburn, and In Salah. We compare and contrast the observed geomechanical behavior of each site, with particular focus on the risks to storage security posed by geomechanical deformation. At Sleipner, the large, high-permeability storage aquifer has experienced little pore pressure increase over 15 y of injection, implying little possibility of geomechanical deformation. At Weyburn, 45 y of oil production has depleted pore pressures before increases associated with CO2 injection. The long history of the ?eld has led to complicated, sometimes nonintuitive geomechanical deformation. At In Salah, injection into the water leg of a gas reservoir has increased pore pressures, leading to uplift and substantial microseismic activity. The differences in the geomechanical responses of these sites emphasize the need for systematic geomechanical appraisal before injection in any potential storage site.

Controlled time of arrival feasibility analysis

Previous research studies and operational trials have shown that using the airborne Required Time of Arrival (RTA) function, an aircraft can individually achieve an assigned time to a metering or merge point accurately. This study goes a step further and investigates the application of RTA to a real sequence of arriving aircraft into Melbourne Australia. Assuming that the actual arrival times were Controlled Time of Arrivals (CTAs) assigned to each aircraft, the study examines if the airborne RTA solution would work. Three scenarios were compared: a baseline scenario being the actual flown trajectories in a two hour time-span into Melbourne, a scenario in which the sequential landing slot times of the baseline scenario were assigned as CTAs and a third scenario in which the landing slots could be freely redistributed to the inbound traffic as CTAs. The research found that pressure on the terminal area would sometimes require aircraft to lose more time than possible through the RTA capability. Using linear holding as an additional measure to absorb extensive delays, up to 500NM (5%) of total track reduction and 1300kg (3%) of total fuel consumption could be saved in the scenario with landing slots freely distributed as CTAs, compared to the baseline scenario. Assigning CTAs in an arrival sequence requires the ground system to have an accurate trajectory predictor to propose additional delay measures (path stretching, linear holding) if necessary. Reducing the achievable time window of the aircraft to add control margin to the RTA function, had a negative impact and increased the amount of intervention other than speed control required to solve the sequence. It was concluded that the RTA capability is not a complete solution but merely a tool to assist in managing the increasing complexity of air traffic.

Efficient speed advisories for multi-stage-metering arrival management

This paper presents a methodology and algorithm for Air Traffic Control (ATC) to efficiently achieve schedules arrival times through speed control in the presence of uncertainty. The methodology does not assume the availability of airborne time of arrival control and can therefore be applied to legacy aircraft. The speed advisories are calculated in a manner that allows for sufficient control margin to, if required, adjust the aircraft's trajectory at a later stage to correct for estimated arrival time drift at the lowest impact to efficiency. The methodology is therefore envisioned to prevent major last-minute interventions and instead assists ATC in allowing more continuous descent approaches to be conducted by aircraft leading to more efficient operations.