Source wavelet estimation for waveform inversion of crosshole georadar data

AbstractFor a wide range of environmental, hydrological, and engineering applications there is a fast growing need for high-resolution imaging. In this context, waveform tomographic imaging of crosshole georadar data is a powerful method able to provide images of pertinent electrical properties in near-surface environments with unprecedented spatial resolution. In contrast, conventional ray-based tomographic methods, which consider only a very limited part of the recorded signal (first-arrival traveltimes and maximum first-cycle amplitudes), suffer from inherent limitations in resolution and may prove to be inadequate in complex environments. For a typical crosshole georadar survey the potential improvement in resolution when using waveform-based approaches instead of ray-based approaches is in the range of one order-of- magnitude. Moreover, the spatial resolution of waveform-based inversions is comparable to that of common logging methods. While in exploration seismology waveform tomographic imaging has become well established over the past two decades, it is comparably still underdeveloped in the georadar domain despite corresponding needs. Recently, different groups have presented finite-difference time-domain waveform inversion schemes for crosshole georadar data, which are adaptations and extensions of Tarantola's seminal nonlinear generalized least-squares approach developed for the seismic case. First applications of these new crosshole georadar waveform inversion schemes on synthetic and field data have shown promising results. However, there is little known about the limits and performance of such schemes in complex environments. To this end, the general motivation of my thesis is the evaluation of the robustness and limitations of waveform inversion algorithms for crosshole georadar data in order to apply such schemes to a wide range of real world problems.One crucial issue to making applicable and effective any waveform scheme to real-world crosshole georadar problems is the accurate estimation of the source wavelet, which is unknown in reality. Waveform inversion schemes for crosshole georadar data require forward simulations of the wavefield in order to iteratively solve the inverse problem. Therefore, accurate knowledge of the source wavelet is critically important for successful application of such schemes. Relatively small differences in the estimated source wavelet shape can lead to large differences in the resulting tomograms. In the first part of my thesis, I explore the viability and robustness of a relatively simple iterative deconvolution technique that incorporates the estimation of the source wavelet into the waveform inversion procedure rather than adding additional model parameters into the inversion problem. Extensive tests indicate that this source wavelet estimation technique is simple yet effective, and is able to provide remarkably accurate and robust estimates of the source wavelet in the presence of strong heterogeneity in both the dielectric permittivity and electrical conductivity as well as significant ambient noise in the recorded data. Furthermore, our tests also indicate that the approach is insensitive to the phase characteristics of the starting wavelet, which is not the case when directly incorporating the wavelet estimation into the inverse problem.Another critical issue with crosshole georadar waveform inversion schemes which clearly needs to be investigated is the consequence of the common assumption of frequency- independent electromagnetic constitutive parameters. This is crucial since in reality, these parameters are known to be frequency-dependent and complex and thus recorded georadar data may show significant dispersive behaviour. In particular, in the presence of water, there is a wide body of evidence showing that the dielectric permittivity can be significantly frequency dependent over the GPR frequency range, due to a variety of relaxation processes. The second part of my thesis is therefore dedicated to the evaluation of the reconstruction limits of a non-dispersive crosshole georadar waveform inversion scheme in the presence of varying degrees of dielectric dispersion. I show that the inversion algorithm, combined with the iterative deconvolution-based source wavelet estimation procedure that is partially able to account for the frequency-dependent effects through an "effective" wavelet, performs remarkably well in weakly to moderately dispersive environments and has the ability to provide adequate tomographic reconstructions.

Parameter estimation for the distribution of single cell lag times.

In Quantitative Microbial Risk Assessment, it is vital to understand how lag times of individual cells are distributed over a bacterial population. Such identified distributions can be used to predict the time by which, in a growth-supporting environment, a few pathogenic cells can multiply to a poisoning concentration level. We model the lag time of a single cell, inoculated into a new environment, by the delay of the growth function characterizing the generated subpopulation. We introduce an easy-to-implement procedure, based on the method of moments, to estimate the parameters of the distribution of single cell lag times. The advantage of the method is especially apparent for cases where the initial number of cells is small and random, and the culture is detectable only in the exponential growth phase.

Panel index VAR models: Specification, Estimation, Testing and Leading Indicators

This paper proposes a method to conduct inference in panel VAR models with cross unit interdependencies and time variations in the coefficients. The approach can be used to obtain multi-unit forecasts and leading indicators and to conduct policy analysis in a multiunit setups. The framework of analysis is Bayesian and MCMC methods are used to estimate the posterior distribution of the features of interest. The model is reparametrized to resemble an observable index model and specification searches are discussed. As an example, we construct leading indicators for inflation and GDP growth in the Euro area using G-7 information.

Body surface area and body weight predict total liver volume in Western adults.

Computed tomography (CT) is used increasingly to measure liver volume in patients undergoing evaluation for transplantation or resection. This study is designed to determine a formula predicting total liver volume (TLV) based on body surface area (BSA) or body weight in Western adults. TLV was measured in 292 patients from four Western centers. Liver volumes were calculated from helical computed tomographic scans obtained for conditions unrelated to the hepatobiliary system. BSA was calculated based on height and weight. Each center used a different established method of three-dimensional volume reconstruction. Using regression analysis, measurements were compared, and formulas correlating BSA or body weight to TLV were established. A linear regression formula to estimate TLV based on BSA was obtained: TLV = -794.41 + 1,267.28 x BSA (square meters; r(2) = 0.46; P <.0001). A formula based on patient weight also was derived: TLV = 191.80 + 18.51 x weight (kilograms; r(2) = 0.49; P <.0001). The newly derived TLV formula based on BSA was compared with previously reported formulas. The application of a formula obtained from healthy Japanese individuals underestimated TLV. Two formulas derived from autopsy data for Western populations were similar to the newly derived BSA formula, with a slight overestimation of TLV. In conclusion, hepatic three-dimensional volume reconstruction based on helical CT predicts TLV based on BSA or body weight. The new formulas derived from this correlation should contribute to the estimation of TLV before liver transplantation or major hepatic resection.

Primary Health Care from the perception of women living in a rural area

AbstractOBJECTIVEUnderstanding the perception of women living in a rural area about the actions and services of Primary Health Care (PHC) in a municipality of southern Brazil, which is the only one regarded as predominantly rural.METHODA descriptive study of qualitative approach, carried out with women who lived in the countryside and required health services in the 15 days prior to collection.RESULTSThe results registered low fidelity to PHC attributes, focusing its functional axis on sickness, transforming the unit into small points of emergency care and a bureaucratic place where patients are referred to other types of services. The quality of service offered is compromised to offering quick, fragmented and unequal treatment in the rural context.CONCLUSIONThe findings of this study highlight the need for greater efforts in order to adequate the new care model in the development of appropriate actions as designated by PHC in the rural context studied.

Eleven-year physical activity trends in a Swiss urban area.

OBJECTIVE: Regular physical activity is a major health determinant. Little is known about physical activity trends. We evaluated whether adult physical activity levels are changing in a Swiss urban state (Geneva). METHOD: We analyzed 11-year trends of physical activity indicators, including 3+MET-minutes per week and physical activity outside working hours, in population representative adults (n=9320, aged 35-74years, 50% women), relating declared physical activity to socioeconomic status, lifestyle, and clinical and blood markers. RESULTS: Combining yearly cohorts from 1999 to 2009, we found a significant trend for increased physical activity levels. Weekly age and sex adjusted 3+MET-minutes per week increased from 3023 to 3752, between 1999 and 2009 (P=0.02). The increase also concerned physical activity outside working hours (+18kcal/day/year). There was a shift from low levels of physical activity levels towards higher activities. Physical activity indicators were associated with socioeconomic status, comorbidities, and biological and anthropometric measures. The trend for increased physical activity was more prominent over the latter 5years. CONCLUSION: We found that physical activity levels have increased in an urban Swiss state. The increase is significant but small, and further efforts to promote physical activity are therefore warranted.

Mandatory accounting disclosure by small private companies

This article analyzes how mandatory accounting disclosure is grounded on differentrationales for private and public companies. It also explores technological changes, such ascomputerised databases and the Internet, which have recently made disclosure of companyaccounts by small companies potentially less costly and more valuable, thanks to electronicfiling and universal online access to credit information systems. These recent developmentsfavour policies that would expand the scope of mandatory publication for small companies incountries where it is voluntary. They also encourage policies to reduce the costs and enhancethe value of disclosure through administrative reforms of filing, archive and retrieval systems.Survey and registry evidence on how the information in the accounts is valued and used bycompanies is consistent with these claims about the evolution of the tradeoff of costs andbenefits that should guide policy in this area.

Late Proterozoic thrust tectonics, high-pressure metamorphism and Uranium mineralization in the domes area, Lufilian arc, northwestern Zambia

The following main lithostratigraphic units have been distinguished in the Domes Area. The Kibaran basement complex composed of gneisses, migmatites with amphibolite bands and metagranites is exposed in dome structures; metamorphic features of Kibaran age have been almost completely obliterated by extensive Lufilian reactivation. The post-Kibaran cover sequence is subdivided into the Lower Roan Group consisting of well-preserved quartzites with high Mg content, talc-bearing, extremely foliated schists intercalated with pseudo-conglomerates of tectonic origin and the Upper Roan Group including dolomitic marbles with rare stromatolites, metapelites and a sequence of detrital metasediments, with local volcano-sedimentary components and interlayered banded ironstones. The sediments of the Lower Roan Group are interpreted as continental to lagoonal-evaporitic deposits partly converted into the talc-kyanite + garnet assemblage characteristic of ``white schists''. The dolomites and metapelites of the Upper Roan Group are attributed to a carbonate platform sequence progressively subsiding under terrigenous deposits, whilst the detrital metasediments and BIF may be interpreted as a basinal sequence, probably deposited on oceanic crust grading laterally into marbles. Metagabbros and metabasalts are considered as remnants of an ocean-floor-type crustal unit probably related to small basins. Alkaline stocks of Silurian age intruded the post-Kibaran cover. Significant ancestral tectonic discontinuities promoted the development of a nappe pile that underwent high-pressure metamorphism during the Lufilian orogeny and all lithostratigraphic units. Rb-Sr and K-Ar and U-Pb data indicate an age of 700 Ma for the highest grade metamorphism and 500 Ma for blocking of the K-Ar and Rb-Sr system in micas, corresponding to the time when the temperature dropped below 350-degrees-400-degrees-C and to an age of about 400 Ma for the emplacement of hypabyssal syenitic bodies. A first phase of crustal shortening by decoupling of basement and cover slices along shallow shear zones has been recognized. Fluid-rich tectonic slabs of cover sediments were thus able to transport fluids into the anhydrous metamorphic basement or mafic units. During the subsequent metamorphic re-equilibration stage of high pressure, pre-existing thrusts horizons were converted into recrystallized mylonites. Due to uplift, rocks were re-equilibrated into assemblages compatible with lower pressures and slightly lower temperatures. This stage occurs under a decompressional (nearly adiabatic) regime, with P(fluid) almost-equal-to P(lithostatic). It is accompanied by metasomatic development of minerals, activated by injection of hot fluids. New or reactivated shear zones and mylonitic belts were the preferred conduits of fluids. The most evident regional-scale effect of these processes is the intense metasomatic scapolitization of formerly plagioclase-rich lithologies. Uraninite mineralization can probably be assigned to the beginning of the decompressional stage. A third regional deformation phase characterized by open folds and local foliation is not accompanied by significant growth of new minerals. However, pitchblende mineralization can be ascribed to this phase as late-stage, short-range remobilization of previously existing deposits. Finally, shallow alkaline massifs were emplaced when the level of the Domes Area now exposed was already subjected to exchange with meteoric circuits, activated by residual geothermal gradients generally related to intrusions or rifting. Most of the superficial U-showings with U-oxidation products were probably generated during this relatively recent phase.

Drosophilids (Diptera) from an Atlantic Forest Area in Santa Catarina, Southern Brazil

The present work aims at knowing the faunal composition of drosophilids in forest areas of southern Brazil. Besides, estimation of species richness for this fauna is briefly discussed. The sampling were carried out in three well-preserved areas of the Atlantic Rain Forest in the State of Santa Catarina. In this study, 136,931 specimens were captured and 96.6% of them were identified in the specific level. The observed species richness (153 species) is the largest that has been registered in faunal inventories conducted in Brazil. Sixty-three of the captured species did not fit to the available descriptions, and we believe that most of them are non-described species. The incidence-based estimators tended to give rise to the largest richness estimates while the abundance based give rise to the smallest ones. Such estimators suggest the presence from 172.28 to 220.65 species in the studied area. Based on these values, from 69.35 to 88.81% of the expected species richness were sampled. We suggest that the large richness recorded in this study is a consequence of the large sampling effort, the capture method, recent advances in the taxonomy of drosophilids, the high preservation level and the large extension of the sampled fragment and the high complexity of the Atlantic Rain forest. Finally, our data set suggest that the employment of estimators of richness for drosophilid assemblages is useful but it requires caution.

Firm investment in imperfect capital markets: A structural estimation

We set up a dynamic model of firm investment in which liquidity constraintsenter explicity into the firm's maximization problem. The optimal policyrules are incorporated into a maximum likelihood procedure which estimatesthe structural parameters of the model. Investment is positively related tothe firm's internal financial position when the firm is relatively poor. This relationship disappears for wealthy firms, which can reach theirdesired level of investment. Borrowing is an increasing function of financial position for poor firms. This relationship is reversed as a firm's financial position improves, and large firms hold little debt.Liquidity constrained firms may be unused credits lines and the capacity toinvest further if they desire. However the fear that liquidity constraintswill become binding in the future induces them to invest only when internalresources increase.We estimate the structural parameters of the model and use them to quantifythe importance of liquidity constraints on firms' investment. We find thatliquidity constraints matter significantly for the investment decisions of firms. If firms can finance investment by issuing fresh equity, rather than with internal funds or debt, average capital stock is almost 35% higher overa period of 20 years. Transitory shocks to internal funds have a sustained effect on the capital stock. This effect lasts for several periods and ismore persistent for small firms than for large firms. A 10% negative shock to firm fundamentals reduces the capital stock of firms which face liquidityconstraints by almost 8% over a period as opposed to only 3.5% for firms which do not face these constraints.

SwissParam: a fast force field generation tool for small organic molecules.

The drug discovery process has been deeply transformed recently by the use of computational ligand-based or structure-based methods, helping the lead compounds identification and optimization, and finally the delivery of new drug candidates more quickly and at lower cost. Structure-based computational methods for drug discovery mainly involve ligand-protein docking and rapid binding free energy estimation, both of which require force field parameterization for many drug candidates. Here, we present a fast force field generation tool, called SwissParam, able to generate, for arbitrary small organic molecule, topologies, and parameters based on the Merck molecular force field, but in a functional form that is compatible with the CHARMM force field. Output files can be used with CHARMM or GROMACS. The topologies and parameters generated by SwissParam are used by the docking software EADock2 and EADock DSS to describe the small molecules to be docked, whereas the protein is described by the CHARMM force field, and allow them to reach success rates ranging from 56 to 78%. We have also developed a rapid binding free energy estimation approach, using SwissParam for ligands and CHARMM22/27 for proteins, which requires only a short minimization to reproduce the experimental binding free energy of 214 ligand-protein complexes involving 62 different proteins, with a standard error of 2.0 kcal mol(-1), and a correlation coefficient of 0.74. Together, these results demonstrate the relevance of using SwissParam topologies and parameters to describe small organic molecules in computer-aided drug design applications, together with a CHARMM22/27 description of the target protein. SwissParam is available free of charge for academic users at www.swissparam.ch.

Phlebotomines (Diptera, Psychodidae) in the Speleological Province of the Ribeira Valley: 3. Serra district - area of hostels for tourists who visit the Parque Estadual do Alto Ribeira (PETAR), state of São Paulo, Brazil

Phlebotomines (Diptera, Psychodidae) in the Speleological Province of the Ribeira Valley: 3. Area of hostels for tourists who visit the Parque Estadual do Alto Ribeira (PETAR), state of São Paulo, Brazil. The study characterizes some ecological aspects of the phlebotomine fauna in an endemic area of cutaneous leishmaniasis (CL) situated in the Serra district, Iporanga municipality where the hostels for tourists visiting the PETAR are located. Captures were undertaken on a smallholding and a small farm situated near the hostels, monthly between January/2001 and December/2003 with automatic light traps (ALT) in pigsty, hen-house and veranda of a domicile at the two sites, and in peridomicile of the small farm also with black/white Shannon traps. With the ALT a total of 87,224 phlebotomines representing 19 species and also two hybrids of Nyssomyia intermedia (Lutz & Neiva) and Nyssomyia neivai (Pinto) and two anomalous specimens were captured. The standardized index species abundance was for Ny. intermedia = 1.0 and Ny. neivai = 0.935. The highest frequencies of the smallholding occurred in the pigsty, the Williams' mean/capture for Ny. intermedia being 63.7 specimens and for Ny. neivai 29.2, and on the small farm, in the hen-house, Ny. intermedia 402.6 and Ny. neivai 116.2. A total of 863 phlebotomines (Ny. intermedia: 75.4%; Ny. neivai: 24.3%) were captured with black/white Shannon traps; females of both species being predominant in the white trap. The high frequencies of Ny. intermedia and Ny. neivai, both implicated in CL transmission, indicate the areas presenting risk of the disease.

Convergence in panel data: Evidence from the skipping estimation

This paper demonstrates that, unlike what the conventional wisdom says, measurement error biases in panel data estimation of convergence using OLS with fixed effects are huge, not trivial. It does so by way of the "skipping estimation"': taking data from every m years of the sample (where m is an integer greater than or equal to 2), as opposed to every single year. It is shown that the estimated speed of convergence from the OLS with fixed effects is biased upwards by as much as 7 to 15%.

Quantitative analysis of crystal/grain sizes and their distributions in 2D and 3D

We review methods to estimate the average crystal (grain) size and the crystal (grain) size distribution in solid rocks. Average grain sizes often provide the base for stress estimates or rheological calculations requiring the quantification of grain sizes in a rock's microstructure. The primary data for grain size data are either 1D (i.e. line intercept methods), 2D (area analysis) or 3D (e.g., computed tomography, serial sectioning). These data have been used for different data treatments over the years, whereas several studies assume a certain probability function (e.g., logarithm, square root) to calculate statistical parameters as the mean, median, mode or the skewness of a crystal size distribution. The finally calculated average grain sizes have to be compatible between the different grain size estimation approaches in order to be properly applied, for example, in paleo-piezometers or grain size sensitive flow laws. Such compatibility is tested for different data treatments using one- and two-dimensional measurements. We propose an empirical conversion matrix for different datasets. These conversion factors provide the option to make different datasets compatible with each other, although the primary calculations were obtained in different ways. In order to present an average grain size, we propose to use the area-weighted and volume-weighted mean in the case of unimodal grain size distributions, respectively, for 2D and 3D measurements. The shape of the crystal size distribution is important for studies of nucleation and growth of minerals. The shape of the crystal size distribution of garnet populations is compared between different 2D and 3D measurements, which are serial sectioning and computed tomography. The comparison of different direct measured 3D data; stereological data and direct presented 20 data show the problems of the quality of the smallest grain sizes and the overestimation of small grain sizes in stereological tools, depending on the type of CSD. (C) 2011 Published by Elsevier Ltd.

Estimation of the extreme flash flood evolution in Barcelona county from 1351 to 2005

Every year, flash floods cause economic losses and major problems for undertaking daily activity in the Catalonia region (NE Spain). Sometimes catastrophic damage and casualties occur. When a long term analysis of floods is undertaken, a question arises regarding the changing role of the vulnerability and the hazard in risk evolution. This paper sets out to give some information to deal with this question, on the basis of analysis of all the floods that have occurred in Barcelona county (Catalonia) since the 14th century, as well as the flooded area, urban evolution, impacts and the weather conditions for any of most severe events. With this objective, the identification and classification of historical floods, and characterisation of flash-floods among these, have been undertaken. Besides this, the main meteorological factors associated with recent flash floods in this city and neighbouring regions are well-known. On the other hand, the identification of rainfall trends that could explain the historical evolution of flood hazard occurrence in this city has been analysed. Finally, identification of the influence of urban development on the vulnerability to floods has been carried out. Barcelona city has been selected thanks to its long continuous data series (daily rainfall data series, since 1854; one of the longest rainfall rate series of Europe, since 1921) and for the accurate historical archive information that is available (since the Roman Empire for the urban evolution). The evolution of flood occurrence shows the existence of oscillations in the earlier and later modern-age periods that can be attributed to climatic variability, evolution of the perception threshold and changes in vulnerability. A great increase of vulnerability can be assumed for the period 1850¿1900. The analysis of the time evolution for the Barcelona rainfall series (1854¿2000) shows that no trend exists, although, due to changes in urban planning, flash-floods impact has altered over this time. The number of catastrophic flash floods has diminished, although the extraordinary ones have increased.