Digital beamforming and moving target indication for future spaceborne synthetic aperture radar

During the last decade the interest on space-borne Synthetic Aperture Radars (SAR) for remote sensing applications has grown as testified by the number of recent and forthcoming missions as TerraSAR-X, RADARSAT-2, COSMO-kyMed, TanDEM-X and the Spanish SEOSAR/PAZ. In this sense, this thesis proposes to study and analyze the performance of the state-of-the-Art space-borne SAR systems, with modes able to provide Moving Target Indication capabilities (MTI), i.e. moving object detection and estimation. The research will focus on the MTI processing techniques as well as the architecture and/ or configuration of the SAR instrument, setting the limitations of the current systems with MTI capabilities, and proposing efficient solutions for the future missions. Two European projects, to which the Universitat Politècnica de Catalunya provides support, are an excellent framework for the research activities suggested in this thesis. NEWA project proposes a potential European space-borne radar system with MTI capabilities in order to fulfill the upcoming European security policies. This thesis will critically review the state-of-the-Art MTI processing techniques as well as the readiness and maturity level of the developed capabilities. For each one of the techniques a performance analysis will be carried out based on the available technologies, deriving a roadmap and identifying the different technological gaps. In line with this study a simulator tool will be developed in order to validate and evaluate different MTI techniques in the basis of a flexible space-borne radar configuration. The calibration of a SAR system is mandatory for the accurate formation of the SAR images and turns to be critical in the advanced operation modes as MTI. In this sense, the SEOSAR/PAZ project proposes the study and estimation of the radiometric budget. This thesis will also focus on an exhaustive analysis of the radiometric budget considering the current calibration concepts and their possible limitations. In the framework of this project a key point will be the study of the Dual Receive Antenna (DRA) mode, which provides MTI capabilities to the mission. An additional aspect under study is the applicability of the Digital Beamforming on multichannel and/or multistatic radar platforms, which conform potential solutions for the NEWA project with the aim to fully exploit its capability jointly with MTI techniques.

Estimativa do teor de água no solo utilizando dados de radar de abertura sintética

Em virtude da crescente demanda mundial por alimentos, um monitoramento eficaz e em larga escala da umidade do solo constitui fator de grande importância para a previsão de safras. Este trabalho teve por objetivo apresentar uma técnica para o cálculo do teor de água no solo, utilizando modelos preditivos de umidade do solo, baseados em dados de radar de abertura sintética (SAR). Foram utilizados dados do SAR a bordo do JERS-1 ("Japanese Earth Resources Satellite") e dois modelos empíricos. O primeiro relaciona o coeficiente de retroespalhamento com a permissividade complexa (modelo de Dubois), e o segundo relaciona a permissividade complexa com o teor de água do solo e algumas de suas características físico-hídricas, tais como percentagem de areia e argila (modelo de Hallikainen). Inicialmente, os dados do SAR/JERS-1 foram calibrados e, por meio do modelo de Dubois, foram calculados os valores de permissividade complexa. Para tanto, foi necessário inserir níveis estimados de rugosidade do solo. A partir destes resultados, utilizou-se o modelo de Hallikainen para calcular a umidade volumétrica. A análise geral dos resultados indica que a técnica de estimação de umidade do solo a partir de imagens de radar de abertura sintética, utilizada neste estudo, mostrou-se física e matematicamente exeqüível. No entanto, apresentou uma precisão moderada, não sendo ainda recomendada para o uso operacional no mapeamento de umidade do solo. A análise dos resultados revelou também que a precisão dos dados é bastante influenciada pela precisão dos valores de rugosidade introduzidos.

Estimation of surface roughness in a semiarid region from C-band ERS-1 synthetic aperture radar data

In this study, we investigated the feasibility of using the C-band European Remote Sensing Satellite (ERS-1) synthetic aperture radar (SAR) data to estimate surface soil roughness in a semiarid rangeland. Radar backscattering coefficients were extracted from a dry and a wet season SAR image and were compared with 47 in situ soil roughness measurements obtained in the rocky soils of the Walnut Gulch Experimental Watershed, southeastern Arizona, USA. Both the dry and the wet season SAR data showed exponential relationships with root mean square (RMS) height measurements. The dry C-band ERS-1 SAR data were strongly correlated (R² = 0.80), while the wet season SAR data have somewhat higher secondary variation (R² = 0.59). This lower correlation was probably provoked by the stronger influence of soil moisture, which may not be negligible in the wet season SAR data. We concluded that the single configuration C-band SAR data is useful to estimate surface roughness of rocky soils in a semiarid rangeland.

Uso do radar penetrante no solo (GPR) na investigação dos solos dos tabuleiros costeiros no litoral norte do estado da Bahia

Foram estudados nove perfis ao longo de uma toposseqüência sobre os sedimentos do Grupo Barreiras, na Fazenda Rio Negro, município de Entre Rios (BA), utilizando a prospecção eletromagnética por meio do Radar Penetrante no Solo - "Ground-penetrating radar - GPR", objetivando analisar a utilização dessa ferramenta na aquisição de informações sobre as feições que ocorrem no solo, mediante a comparação entre os radargramas obtidos e a descrição pedológica. O equipamento utilizado foi um Geophysical Survey System modelo GPR SR system-2, com antena de 80 MHz. A análise radargramétrica confirmou o aparecimento dos fragipãs e duripãs em profundidade, que ocorrem sempre acompanhados de um processo de transformação dos solos do tipo Latossolo Amarelo e Argissolo Amarelo em Espodossolo. Os padrões de reflexão mostram claramente os domínios dos solos argilosos e dos solos arenosos, com e sem a presença dos horizontes endurecidos.

Radar analysis of the life cycle of mesoscale convective systems during the 10 June 2000 event

The 10 June 2000 event was the largest flash flood event that occurred in the Northeast of Spain in the late 20th century, both as regards its meteorological features and its considerable social impact. This paper focuses on analysis of the structures that produced the heavy rainfalls, especially from the point of view of meteorological radar. Due to the fact that this case is a good example of a Mediterranean flash flood event, a final objective of this paper is to undertake a description of the evolution of the rainfall structure that would be sufficiently clear to be understood at an interdisciplinary forum. Then, it could be useful not only to improve conceptual meteorological models, but also for application in downscaling models. The main precipitation structure was a Mesoscale Convective System (MCS) that crossed the region and that developed as a consequence of the merging of two previous squall lines. The paper analyses the main meteorological features that led to the development and triggering of the heavy rainfalls, with special emphasis on the features of this MCS, its life cycle and its dynamic features. To this end, 2-D and 3-D algorithms were applied to the imagery recorded over the complete life cycle of the structures, which lasted approximately 18 h. Mesoscale and synoptic information were also considered. Results show that it was an NS-MCS, quasi-stationary during its stage of maturity as a consequence of the formation of a convective train, the different displacement directions of the 2-D structures and the 3-D structures, including the propagation of new cells, and the slow movement of the convergence line associated with the Mediterranean mesoscale low.

Modelling of weather radar echoes from anomalous propagation using a hybrid parabolic equation method and NWP model data

Contamination of weather radar echoes by anomalous propagation (anaprop) mechanisms remains a serious issue in quality control of radar precipitation estimates. Although significant progress has been made identifying clutter due to anaprop there is no unique method that solves the question of data reliability without removing genuine data. The work described here relates to the development of a software application that uses a numerical weather prediction (NWP) model to obtain the temperature, humidity and pressure fields to calculate the three dimensional structure of the atmospheric refractive index structure, from which a physically based prediction of the incidence of clutter can be made. This technique can be used in conjunction with existing methods for clutter removal by modifying parameters of detectors or filters according to the physical evidence for anomalous propagation conditions. The parabolic equation method (PEM) is a well established technique for solving the equations for beam propagation in a non-uniformly stratified atmosphere, but although intrinsically very efficient, is not sufficiently fast to be practicable for near real-time modelling of clutter over the entire area observed by a typical weather radar. We demonstrate a fast hybrid PEM technique that is capable of providing acceptable results in conjunction with a high-resolution terrain elevation model, using a standard desktop personal computer. We discuss the performance of the method and approaches for the improvement of the model profiles in the lowest levels of the troposphere.

Mass conservative three-dimensional water tracer distribution from Markov chain Monte Carlo inversion of time-lapse ground-penetrating radar data

Time-lapse geophysical measurements are widely used to monitor the movement of water and solutes through the subsurface. Yet commonly used deterministic least squares inversions typically suffer from relatively poor mass recovery, spread overestimation, and limited ability to appropriately estimate nonlinear model uncertainty. We describe herein a novel inversion methodology designed to reconstruct the three-dimensional distribution of a tracer anomaly from geophysical data and provide consistent uncertainty estimates using Markov chain Monte Carlo simulation. Posterior sampling is made tractable by using a lower-dimensional model space related both to the Legendre moments of the plume and to predefined morphological constraints. Benchmark results using cross-hole ground-penetrating radar travel times measurements during two synthetic water tracer application experiments involving increasingly complex plume geometries show that the proposed method not only conserves mass but also provides better estimates of plume morphology and posterior model uncertainty than deterministic inversion results.

Integration of ground-penetrating radar and microgravimetric methods to map shallow caves

Ground-penetrating radar (GPR) and microgravimetric surveys have been conducted in the southern Jura mountains of western Switzerland in order to map subsurface karstic features. The study site, La Grande Rolaz cave, is an extensive system in which many portions have been mapped. By using small station spacing and careful processing for the geophysical data, and by modeling these data with topographic information from within the cave, accurate interpretations have been achieved. The constraints on the interpreted geologic models are better when combining the geophysical methods than when using only one of the methods, despite the general limitations of two-dimensional (2D) profiling. For example, microgravimetry can complement GPR methods for accurately delineating a shallow cave section approximately 10 X 10 mt in size. Conversely, GPR methods can be complementary in determining cavity depths and in verifying the presence of off-line features and numerous areas of small cavities and fractures, which may be difficult to resolve in microgravimetric data.

Uso de sistema GPR (Ground Penetring Radar) na avaliação de atributos de um solo sob plantio de cana-de-açúcar

O GPR utiliza técnica indireta para investigação das estruturas de subsuperfície, sendo caracterizado como método não-invasivo, o qual permite extrair informações ao longo do perfil de solo sem perfurar, sondar ou escavar. Este trabalho teve o objetivo de avaliar a utilização do GPR no estabelecimento de atributos físicos do solo em área de cultivo com cana-de-açúcar, além de estabelecer sua correlação com o método do cilindro volumétrico. A área de estudo está localizada no município de Bom Sucesso - PR, sendo utilizada com cultura da cana há mais de 20 anos. Foram obtidas amostras de solos indeformadas com cilindro volumétrico nas profundidades de 0-0,1, 0,1-0,2, 0,2-0,3, 0,3-0,4 e 0,4-0,5 m, ao longo de três transeptos paralelos, sendo o central referente à linha de plantio (transepto 2) e os outros dois correspondentes às entrelinhas laterais (transeptos 1 e 3). Os valores de densidade e porosidade total, obtidos por meio das amostras físicas, foram inicialmente interpretados para verificar o estado do solo nas linhas e nas entrelinhas de cana, nas cinco profundidades. Posteriormente, foram correlacionados com os valores das constantes dielétricas (K) adquiridos pelo georradar, porém apenas para a camada de 0,4-0,5 m, uma vez que as camadas superficiais até a profundidade de 0,4 m não puderam ser avaliadas devido às limitações técnicas do equipamento, aliadas às condições físicas e químicas do solo, que acabaram por prejudicar as reflexões. A densidade mostrou-se menor nas linhas, quando comparada às entrelinhas, para as camadas de 0-0,1 e 0,2-0,3 m. A porosidade total foi menor nas entrelinhas do que nas linhas na camada de 0-0,1 m. Os valores de K foram maiores nas linhas, onde a umidade e a porosidade total foram maiores, e menores nas entrelinhas, onde as densidades encontravam se próximas ao limite crítico. A correlação de K com a densidade e porosidade total não foi significativa. Pequenas modificações na porosidade do solo em função do aumento da densidade foram suficientes para diminuir o valor de K nas entrelinhas, devido à redução na quantidade de água nos poros. Esse fato pode ser confirmado pela correlação de K com a umidade. Assim, pode-se dizer que o georradar foi extremamente sensível à variação de umidade no solo, apresentando melhores resultados para esta variável do que para a densidade e porosidade total determinada por meio do cilindro volumétrico.

Seismic reflectivity of the sediment-covered seafloor: Effects of velocity gradients and fine-scale layering

Knowledge of the reflectivity of the sediment-covered seabed is of significant importance to marine seismic data acquisition and interpretation as it governs the generation of reverberations in the water layer. In this context pertinent, but largely unresolved, questions concern the importance of the typically very prominent vertical seismic velocity gradients as well as the potential presence and magnitude of anisotropy in soft surficial seabed sediments. To address these issues, we explore the seismic properties of granulometric end-member-type clastic sedimentary seabed models consisting of sand, silt, and clay as well as scale-invariant stochastic layer sequences of these components characterized by realistic vertical gradients of the P- and S-wave velocities. Using effective media theory, we then assess the nature and magnitude of seismic anisotropy associated with these models. Our results indicate that anisotropy is rather benign for P-waves, and that the S-wave velocities in the axial directions differ only slightly. Because of the very high P- to S-wave velocity ratios in the vicinity of the seabed our models nevertheless suggest that S-wave triplications may occur at very small incidence angles. To numerically evaluate the P-wave reflection coefficient of our seabed models, we apply a frequency-slowness technique to the corresponding synthetic seismic wavefields. Comparison with analytical plane-wave reflection coefficients calculated for corresponding isotropic elastic half-space models shows that the differences tend to be most pronounced in the vicinity of the elastic equivalent of the critical angle as well as in the post-critical range. We also find that the presence of intrinsic anisotropy in the clay component of our layered models tends to dramatically reduce the overall magnitude of the P-wave reflection coefficient as well as its variation with incidence angle.

Evaluation of Bond Retainage in Portland Cement Concrete Overlays Using Infrared Thermography and Ground Penetrating Radar, HR-537, Final Report, 1989

When concrete deterioration begins to occur in highway pavement, repairs become necessary to assure the rider safety, extend its useful life and restore its riding qualities. One rehabilitation technique used to restore the pavement to acceptable highway standards is to apply a thin portland cement concrete (PCC) overlay to the existing pavement. First, any necessary repairs are made to the existing pavement, the surface is then prepared, and the PCC overlay is applied. Brice Petrides-Donohue, Inc. (Donohue) was retained by the Iowa Department of Transportation (IDOT) to evaluate the present condition with respect to debonding of the PCC overlay at fifteen sites on Interstate 80 and State Highway 141 throughout the State of Iowa. This was accomplished by conducting an infrared thermographic and ground penetrating radar survey of these sites which were selected by the Iowa Department of Transportation. The fifteen selected sites were all two lanes wide and one-tenth of a mile long, for a total of three lane miles or 190,080 square feet. The selected sites are as follows: On Interstate 80 Eastbound, from milepost 35.25 to 35.35, milepost 36.00 to 36.10, milepost 37.00 to 37.10, milepost 38.00 to 38.10 and milepost 39.00 to 39.10, on State Highway 141 from milepost 134.00 to 134.10, milepost 134.90 to milepost 135.00, milepost 135.90 to 136.00, milepost 137.00 to 137.10 and milepost 138.00 to 138.10, and on Interstate 80 Westbound from milepost 184.00 to 184.10, milepost 185.00 to 185.10, milepost 186.00 to 186.10, milepost 187.00 to 187.10, and from milepost 188.00 to 188.10.

Analysis of an iterative deconvolution approach for estimating the source wavelet during waveform inversion of crosshole ground-penetrating radar data

A major issue in the application of waveform inversion methods to crosshole georadar data is the accurate estimation of the source wavelet. Here, we explore the viability and robustness of incorporating this step into a time-domain waveform inversion procedure through an iterative deconvolution approach. Our results indicate that, at least in non-dispersive electrical environments, such an approach provides remarkably accurate and robust estimates of the source wavelet even in the presence of strong heterogeneity in both the dielectric permittivity and electrical conductivity. Our results also indicate that the proposed source wavelet estimation approach is relatively insensitive to ambient noise and to the phase characteristics of the starting wavelet. Finally, there appears to be little-to-no trade-off between the wavelet estimation and the tomographic imaging procedures.