Multi-directional wave spectra from marine X-band radar

The signal measured by heave–pitch–roll directional wave buoys yields the first four coefficients of a Fourier series. Data adaptive methods must be employed to estimate a directional wave spectrum. Marine X-band radars (MRs) have the advantage over buoys that they can measure “model-free” two-dimensional (2D) wave spectra. This study presents the first comprehensive validation of MR-derived multi-directional wave characteristics. It is based on wave data from the 2010 Impact of Typhoons on the Ocean in the Pacific (ITOP) experiment in the Philippine Sea, namely MR measurements from R/V Roger Revelle, Extreme Air–Sea Interaction (EASI) buoy measurements, as well as WAVEWATCH-III (WW3) modeling results. Buoy measurements of mean direction and spreading as function of frequency, which do not require data adaptive methods, are used to validate the WW3 wave spectra. An advanced MR wave retrieval technique is introduced that addresses various shortcomings of existing methods. Spectral partitioning techniques, applied to MR and WW3 results, reveal that multimodal seas are frequently present. Both data sets are in excellent agreement, tracking the evolution of up to 4 simultaneous wave systems over extended time periods. This study demonstrates MR’s and WW3’s strength at measuring and predicting 2D wave spectra in swell-dominated seas.     

Multi-hydro hydrological modelling of a complex peri-urban catchment with storage basins comparing C-band and X-band radar rainfall data

ABSTRACT

The spread of impervious surfaces in urban areas combined with the rise in the intensity of rainfall events as a result of climate change has led to dangerous increases in storm water flows. This paper discusses a new implementation of the fully distributed hydrological model Multi-Hydro (developed at École des Ponts ParisTech), when operating storage basins, and its ability to deal with high-resolution radar rainfall data. The peri-urban area of Massy (south of Paris, France) was selected as a case study for having six of these drainage facilities, used extensively in flood control. Two radar rainfall datasets with different spatiotemporal resolutions were used: Météo-France’s PANTHER rainfall product (C-band) and ENPC’s X-band DPSRI. The rainfall spatiotemporal variability was analysed statistically using Universal Multifractals (UM). Finally, to validate the application, the water level simulations were compared with local measurements in the Cora storage basin located next to the catchment’s single outlet.     

Stochastic estimation of facies using ground penetrating radar data

Explicitly defining large-scale heterogeneity is a necessary step of groundwater model calibration if accurate estimates of flow and transport are to be made. In this work, neural networks are used to estimate radar facies probabilities from ground penetrating radar (GPR) images, yielding stochastic facies-based models that honour the large-scale architecture of the subsurface. For synthetic GPR images, a neural network was able to correctly identify radar facies with an accuracy of approximately 90%. Manual interpretation of a set of 450 MHz GPR field data from the Borden aquifer resulted in the identification of four radar facies. Of these, a neural network was able to identify two facies with an accuracy of near 80% and one with an accuracy of 44%. The neural network was not able to identify the fourth facies, likely due to the choice of defining facies characteristics. Sequential indicator simulation was used to generate facies realizations conditioned to the radar facies probabilities. Numerical simulations indicate that significant improvements in the prediction of solute transport are possible when GPR is used to constrain the facies model compared to using well data alone, especially when data are sparse.This work was supported by funding to R. Knight under Grant No. DE-FG07–00ER15118-A000, Environmental Management Science Program, Office of Science and Technology, Office of Environment Management, United States Department of Energy (DOE). However, any opinions, findings, conclusions, or recommendations expressed herein are those of the authors and do not necessarily reflect the views of DOE. Further support was provided by a Stanford Graduate Fellowship to S. Moysey. The authors would also like to thank James Irving for his assistance with processing of the radar data.     

The generation ofK indices from digitally recorded magnetic data

TheK index was developed by Bartels and is a measure of geomagnetic activity. The important feature of this index is that it is a measure of solar activity superimposed on the regular daily variation.It is the most used estimate of activity being of particular interest to researchers in ionospheric physics and telecommunications.In response to increased availability of digital data from magnetic observatories and their demand by the scientific community it becomes desireable to check whetherK indices produced from magnetograms plotted from digital data are consistent with those handscaled from La Cour charts.This paper presents results of these comparisons and also reports on several methods tested to produce theK index directly from the digital data.Handscaled indices from digital plots are suitable for producing theK index. Computations using digital data do not reproduceK as defined by Bartels precisely but do provide an equivalent index which is satisfactory for most user and research purposes.     

Isostatic response of the Laccadive Ridge from admittance analysis of gravity and bathymetry data

The Laccadive Ridge (L-R), trending roughly parallel to the west coast of India, is an intriguing segment of the northernmost Chagos-Laccadive Ridge (C-L-R) system. Although crustal nature and isostatic response of the southern C-L-R is well known, there are no similar studies on the L-R. In the present study, the isostatic response of the lithosphere beneath the L-R is estimated so as to characterize its crustal nature, total crustal as well as effective elastic plate thickness and mode of compensation. Twelve gravity and bathymetry profiles across the ridge were analyzed using linear transfer function and forward model techniques. The observed admittance function within the diagnostic waveband of 250 < λ > 80 km (0.025 < k > 0.080 km⁻¹) fits well with (i) the Airy model whose average crustal thickness (T_c) and density are 17 ± 2 km and 2.7 × 10³ kg m⁻³, respectively, and (ii) the thin plate flexure model of isostasy with an effective elastic plate thickness (T_e) of 2–3 km. The estimated average crustal thickness and density are in good agreement with published seismic refraction results over the ridge. The results of the present study support an Airy model of isostasy for the L-R. The low T_e value, in view of other published results in the study area, suggests stretched and loaded continental lithosphere of the L-R during the evolution of the western continental margin of India.     

Mapping of a major paleodrainage system in eastern Libya using orbital imaging radar: The Kufrah River

Over the last few decades, remote sensing has revealed buried river channels in a number of regions worldwide, in many cases providing evidence of dramatic paleoenvironmental changes over Cenozoic time scales. Using orbital radar satellite imagery, we mapped a major paleodrainage system in eastern Libya, that could have linked the Kufrah Basin to the Mediterranean coast through the Sirt Basin, possibly as far back as the middle Miocene. Synthetic Aperture Radar images from the PALSAR sensor clearly reveal a 900 km-long river system, which starts with three main tributaries (north-eastern Tibesti, northern Uweinat and western Gilf Kebir/Abu Ras) that connect in the Kufrah oasis region. The river system then flows north through the Jebel Dalmah, and forms a large alluvial fan in the Sarir Dalmah. The sand dunes of the Calanscio Sand Sea prevent deep orbital radar penetration and preclude detailed reconstruction of any possible connection to the Mediterranean Sea, but a 300 km-long link to the Gulf of Sirt through the Wadi Sahabi paleochannel is likely. If this connection is confirmed, and its Miocene antiquity is established, then the Kufrah River, comparable in length to the Egyptian Nile, will have important implications for the understanding of the past environments and climates of northern Africa from the middle Miocene to the Holocene.     

Determination of noise spectra from strong motion data recorded in Greece

A large number of strong ground motions inGreece have been recorded by analogaccelerographs. The processing of analogstrong motion recordings and conversion toa digital form introduces noise in thesignal, due to the digitization andprocessing whichsignificantly affects the record. In thepresent paper a unified processing anddetermination of digitization andprocessing noise for the Greek strongmotion records is presented. Moreover,appropriate relations are proposed for thelower cut-off frequency with respect tothe epicentral distance and earthquakemagnitude for record filtering.     

Extracting rotational components of earthquake ground motion using data recorded at multiple stations

Rotational components of earthquake ground motion have not been considered for seismic analysis, design and performance assessment because recordings of these components are unavailable. A number of procedures have been proposed to extract rotational components of ground motion from translational time series recorded at multiple, closely spaced recording stations. In this paper, a new procedure that is capable of capturing higher frequency content in rotational time‐series is presented. The frequencies at which numerical errors are introduced in the solution, which are a function of apparent wave velocity and array dimension, are identified. Results are presented for the proposed procedure, the widely accepted geodetic method, and a single‐station procedure developed by the authors, all using data recorded at the Lotung array in Taiwan. Copyright © 2012 John Wiley & Sons, Ltd.     

Spatial reconstruction of rainfall fields from rain gauge and radar data

Rainfall is a phenomenon difficult to model and predict, for the strong spatial and temporal heterogeneity and the presence of many zero values. We deal with hourly rainfall data provided by rain gauges, sparsely distributed on the ground, and radar data available on a fine grid of pixels. Radar data overcome the problem of sparseness of the rain gauge network, but are not reliable for the assessment of rain amounts. In this work we investigate how to calibrate radar measurements via rain gauge data and make spatial predictions for hourly rainfall, by means of Monte Carlo Markov Chain algorithms in a Bayesian hierarchical framework. We use zero-inflated distributions for taking zero-measurements into account. Several models are compared both in terms of data fitting and predictive performances on a set of validation sites. Finally, rainfall fields are reconstructed and standard error estimates at each prediction site are shown via easy-to-read spatial maps.     

MIDAS-W: a workstation-based incoherent scatter radar data acquisition system

The Millstone Hill Incoherent Scatter Data Acquisition System (MIDAS) is based on an abstract model of an incoherent scatter radar. This model is implemented in a hierarchical software system, which serves to isolate hardware and low-level software implementation details from higher levels of the system. Inherent in this is the idea that implementation details can easily be changed in response to technological advances. MIDAS is an evolutionary system, and the MIDAS hardware has, in fact, evolved while the basic software model has remained unchanged. From the earliest days of MIDAS, it was realized that some functions implemented in specialized hardware might eventually be implemented by software in a general-purpose computer. MIDAS-W is the realization of this concept. The core component of MIDAS-W is a Sun Microsystems UltraSparc 10 workstation equipped with an Ultrarad 1280 PCI bus analog to digital (A/D) converter board. In the current implementation, a 2.25 MHz intermediate frequency (IF) is bandpass sampled at 1 s intervals and these samples are multicast over a high-speed Ethernet which serves as a raw data bus. A second workstation receives the samples, converts them to filtered, decimated, complex baseband samples and computes the lag-profile matrix of the decimated samples. Overall performance is approximately ten times better than the previous MIDAS system, which utilizes a custom digital filtering module and array processor based correlator. A major advantage of MIDAS-W is its flexibility. A portable, single-workstation data acquisition system can be implemented by moving the software receiver and correlator programs to the workstation with the A/D converter. When the data samples are multicast, additional data processing systems, for example for raw data recording, can be implemented simply by adding another workstation with suitable software to the high-speed network. Testing of new data processing software is also greatly simplified, because a workstation with the new software can be added to the network without impacting the production system. MIDAS-W has been operated in parallel with the existing MIDAS-1 system to verify that incoherent scatter measurements by the two systems agree. MIDAS-W has also been used in a high-bandwidth mode to collect data on the November, 1999, Leonid meteor shower.     

Mapping hydraulic biotopes using terrestrial laser scan data of water surface properties

For more than a decade, habitat mapping using biotopes (in‐channel hydraulically‐defined habitats) has underpinned aquatic conservation in the UK through (a) providing baseline information on system complexity and (b) allowing environmental and ecological change to be monitored and evaluated. The traditional method used is the subjective river habitat or corridor survey. This has recently been revised to include the floodplain via GeoRHS, but issues still exist concerning development of a national database due to the labour intensive nature of the data collection, subjectivity issues between samplers, temporal changes, the fuzzy nature of perceived habitats and habitat boundaries. This paper takes an innovative approach to biotope definition using high resolution spatial data to define water surface roughness for two representative reaches of the River South Tyne, Cumbria, and the River Rede, Northumberland, UK. Data was collected using a terrestrial laser scanner (TLS) and hydraulic variability simply expressed through assigning a local standard deviation value to a set of adjacent water surface values. Statistical linkage of these data with biotope locations defined visually in the field allowed complete mapping of the surveyed reach defining habitat and biotope areas to the fine scale resolution of the TLS data. Despite issues of data loss due to absorption and transmission through the water, the reflected signal generated an extremely detailed and objective map of the water surface roughness, which may be compared with known biotope locations as defined by visual identification in the field. The TLS accuracy achieved in the present study is comparable with those obtained using hyperspectral imagery: with 84% of the pool/glide/marginal deadwater amalgamated biotope, 88% of riffles, 57% of runs and 50% of the amalgamated cascade/rapid biotope successfully plotted. It is clear from this exercise that biotope distribution is more complex than previously mapped using subjective techniques, and based upon the water surface roughness delimiters presented in this study, the amalgamation of pools with glides and marginal deadwaters, riffles with unbroken standing waves, and cascades with rapids, is proposed. Copyright © 2010 John Wiley & Sons, Ltd     

Origin and compensation of Chagos-Laccadive ridge, Indian ocean, from admittance analysis of gravity and bathymetry data

The Chagos-Laccadive ridge (CLR) is a prominent aseismic, volcanic ridge in the northern Indian ocean. The ridge, together with the Southern Mascarene plateau (SMP), to which it is genetically related, is considered as a volcanic trace of the Reunion hotspot. We have examined the isostatic compensation of the CLR through transfer function analysis of gravity and bathymetry data along seven profiles. The analysis suggests that the CLR is compensated locally, with an Airy crustal thickness (T_c) of 20 km. The rather low elastic plate thickness (T_e) of about 4 km implies that the volcanism of the ridge took place very near a spreading centre. The proximity of the Chagos fracture zone indicates that the emplacement was probably near a spreading centre-transform junction.     

Estimation of effective source signatures from marine VSP data1

We present a simple method for estimating an effective source wavelet from the first arrival in marine vertical seismic profiling (VSP) data. The method, which utilizes the free-space Green's function of the Helmholz equation, is simple and very computer efficient. We show examples from synthetic and real offset and walkaway VSP data. In the synthetic examples, we show that data modelled with the estimated wavelet give small residuals when subtracted from the reference data. In the real data examples, we show that when modelling with the wavelet estimated from the real data, in a smooth macromodel, we obtain a good fit between the first arrivals in the real and modelled data.     

Evapotranspiration modelling from climatic data using a neural computing technique

Evapotranspiration is one of the basic components of the hydrologic cycle and essential for estimating irrigation water requirements. This paper investigates the modelling of evapotranspiration using the feed‐forward artificial neural network (ANN) technique with the Levenberg–Marquardt (LM) training algorithm. The LM algorithm has never been used in evapotranspiration estimation before. The LM is used for the optimization of network weights, since this algorithm is more powerful and faster than the conventional gradient descent. Various combinations of daily climatic data, i.e. wind speed, air temperature, relative humidity and solar radiation, from three stations in Los Angeles, USA, are used as inputs to the ANN so as to evaluate the degree of effect of each of these variables on evapotranspiration. A comparison is made between the estimates provided by the ANN and those of the following empirical models: Penman, Hargreaves, Turc. Mean square error, mean absolute error and determination coefficient statistics are used as comparing criteria for the evaluation of the models' performances. Based on the comparisons, it was found that the neural computing technique could be employed successfully in modelling evapotranspiration process from the available climatic data. The results also indicate that the Hargreaves method provides better performance than the Penman and Turc methods in estimation of the evapotranspiration. The accuracy of the ANN technique in evapotranspiration estimation using nearby station data was also investigated. Copyright © 2007 John Wiley & Sons, Ltd.     

Time-domain identification of dynamic properties of layered soil by using extended Kalman filter and recorded seismic data

A novel time-domain identification technique is developed for the seismic response analysis of soil-structure interaction. A two-degree-of-freedom (2DOF) model with eight lumped parameters is adopted to model the frequency-dependent behavior of soils. For layered soil, the equivalent eight parameters of the 2DOF model arc identified by the extended Kalman filter (EKF) method using recorded seismic data. The polynomial approximations for derivation of state estimators are applied in the EKF procedure. A realistic identification example is given for the layered-soil of a building site in Anchorage, Alaska in the United States. Results of the example demonstrate the feasibility and practicality of the proposed identification technique. The 2DOF soil model and the identification technique can be used for nonlinear response analysis of soil-structure interaction in the time-domain for layered of complex soil conditions. The identified parameters can be stored in a database for use in other similar soil conditions. If a universal database that covers information related to most soil conditions is developed in the future, engineers could conveniently perform time history analyses of soil-structural interaction.     

Hydrological modelling using data from monthly GCMs in a regional catchment

This study demonstrates the use of spatially downscaled, monthly general circulation model (GCM) rainfall and temperature data to drive the established HyMOD hydrological model to evaluate the prospective effects of climate change on the fluvial run‐off of the River Derwent basin in the UK. The evaluation results of this monthly hydrological model using readily available, monthly GCM data are consistent with studies on nearby catchments employing high‐temporal resolution data, indicating that useful hydro‐climatic planning studies may be possible using standard datasets and modest computational resources. HyMOD was calibrated against 5 km² gridded UK Climate Projections dataset data and then driven using monthly spatially interpolated (~5 km²) outputs from Hadley Centre Coupled Model, version 3 and the Canadian Centre for Climate Modelling and Analysis for Intergovernmental Panel on Climate Change Special Report on Emissions Scenarios (IPCC‐SRES) A2a and B2a covering the 2020s, 2050s and 2080s. Results for both GCMs project a decrease in annual run‐off in both GCM models and scenarios with higher values in the summer/autumn months, whereas an increase in the later winter months. Both Hadley Centre Coupled Model, version 3 and the Canadian Centre for Climate Modelling and Analysis show higher ranges of uncertainty during the winter season with higher values of run‐off associated with December in all three simulation periods and two scenarios. A seasonal comparison of run‐off simulations shows that both GCMs give similar results in summer and autumn, whereas disparities due to GCM uncertainties are more conspicuous in winter and spring. In this study, both the GCMs under A2a scenario have demonstrated the high possibility of time shift in monthly average peak run‐offs in the Derwent River by 2080s in comparison with the early 21st century. Copyright © 2013 John Wiley & Sons, Ltd.     

Long-term diagnostics of precipitation estimates and the development of radar hardware monitoring within a radar product data quality management system

Abstract

Quality is key to ensuring that the potential offered by weather radar networks is realized. To achieve optimum quality, a comprehensive radar data quality management system, designed to monitor the end-to-end radar data processing chain and evaluate product quality, is being developed at the UK Met Office. Three contrasting elements of this system are described: monitoring of key radar hardware performance indicators; generation of long-term integrations of radar products; and monitoring of radar reflectivity factor using synthesized observations from numerical weather prediction model fields. Examples of each component are presented and ways in which the different types of monitoring information have been used to both identify issues with the radar product data quality and help formulate solutions are given.

Editor Z.W. Kundzewicz; Guest editor R.J. Moore

Citation Harrison, D., Georgiou, S., Gaussiat, N., and Curtis, A., 2013. Long-term diagnostics of precipitation estimates and the development of radar hardware monitoring within a radar product data quality management system. Hydrological Sciences Journal, 59 (7), 1327–1342. http://dx.doi.org/10.1080/02626667.2013.841316     

基于VSP直达波资料的粘弹介质参数遗传算法反演

粘弹性参数变得越来越重要,其反演算法也逐渐成为众多研究者的研究热点。而遗传算法是一种随机、自适应、启发式的算法, 具有很好的鲁棒性和全局收敛性, 本文基于VSP直达波方程,引入了遗传算法来进行粘弹参数反演, 首先将频率域直达波方程表示为复速度的函数,然后通过遗传算法反演出复速度。而复速度和品质因子又是复速度的函数,从而便可很容易的得出。但若直接反演复速度, 反演参数太多, 不容易实现, 所以又将复速度表示成参数C0和C∞的函数,以减少反演参数数量。最后给出了理论模型实验,以证明该算法的有效性。     

Detection of wave disturbances in the mesosphere using a MF-HF radar

Theoretical backgrounds of the active and passive methods for detecting the set of acoustic gravity waves in the upper mesosphere are presented. The observation and processing techniques based on the registration, statistical and spectral analysis of the envelopes of backscatter signals and radio noise are presented. It has been indicated that wave disturbances with periods of 5–120 min were observed during the entire observation period. The average relative amplitude of these disturbances was 2–4%. The dependences of the wave disturbance parameters on the time of day, season, and magnetic activity level have been studied.