Influence of end-expiratory level and tidal volume on gravitational ventilation distribution during tidal breathing in healthy adults

Our understanding of regional filling of the lung and regional ventilation distribution is based on studies using stepwise inhalation of radiolabelled tracer gases, magnetic resonance imaging and positron emission tomography. We aimed to investigate whether these differences in ventilation distribution at different end-expiratory levels (EELs) and tidal volumes (V (T)s) held also true during tidal breathing. Electrical impedance tomography (EIT) measurements were performed in ten healthy adults in the right lateral position. Five different EELs with four different V (T)s at each EEL were tested in random order, resulting in 19 combinations. There were no measurements for the combination of the highest EEL/highest V (T). EEL and V (T) were controlled by visual feedback based on airflow. The fraction of ventilation directed to different slices of the lung (VENT(RL1)-VENT(RL8)) and the rate of the regional filling of each slice versus the total lung were analysed. With increasing EEL but normal tidal volume, ventilation was preferentially distributed to the dependent lung and the filling of the right and left lung was more homogeneous. With increasing V (T) and maintained normal EEL (FRC), ventilation was preferentially distributed to the dependent lung and regional filling became more inhomogeneous (p < 0.05). We could demonstrate that regional and temporal ventilation distribution during tidal breathing was highly influenced by EEL and V (T).

Regional and overall ventilation inhomogeneities in preterm and term-born infants

OBJECTIVES: We compared ventilation inhomogeneity assessed by electrical impedance tomography (EIT) and multiple breath washout (MBW) in preterm and term-born infants. We hypothesised that EIT measurements in spontaneously breathing infants are repeatable and that differences in regional ventilation distribution measured by EIT can distinguish between preterm and term-born infants. DESIGN: Cross-sectional group comparison study. SETTING: Lung function laboratory at a University Children's Hospital. PARTICIPANTS: Seventeen healthy term-born and 15 preterm infants at a matched postmenstrual age of 44 weeks. MEASUREMENTS AND RESULTS: We concurrently measured ventilation inhomogeneity by EIT, ventilation inhomogeneity (LCI) and functional residual capacity (FRC) by MBW and tidal breathing variables during unsedated quiet sleep. EIT measurements were highly repeatable (coefficient of variation 3.6%). Preterm infants showed significantly more ventilation of the independent parts of the lungs compared to healthy term-born infants assessed by EIT (mean difference 5.0, 95 CI 1.3-8%). Whereas the two groups showed no differences in lung volumes or ventilation inhomogeneities assessed by MBW, EIT discriminated better between term and preterm infants. (FRC/kg: mean difference 1.1 mL, 95% CI -1.4-3.8 mL; LCI: mean difference 0.03, 95% CI -0.32-0.25). CONCLUSIONS: EIT shows distinct differences in ventilation distribution between preterm and term-born infants, which cannot be detected by MBW. Although preterm infants are capable of dynamically maintaining overall functional residual volume and ventilation distribution, they show some spatial differences from fullterm infants.

Effect of different breathing AIDS on ventilation distribution in adults with cystic fibrosis.

BACKGROUND AND OBJECTIVES We investigated the effect of different breathing aids on ventilation distribution in healthy adults and subjects with cystic fibrosis (CF). METHODS In 11 healthy adults and 9 adults with CF electrical impedance tomography measurements were performed during spontaneous breathing, continuous positive airway pressure (CPAP) and positive expiratory pressure (PEP) therapy randomly applied in upright and lateral position. Spatial and temporal ventilation distribution was assessed. RESULTS The proportion of ventilation directed to the dependent lung significantly increased in lateral position compared to upright in healthy and CF. This effect was enhanced with CPAP but neutralised with PEP, whereas the effect of PEP was larger in the healthy group. Temporal ventilation distribution showed exactly the opposite with homogenisation during CPAP and increased inhomogeneity with PEP. CONCLUSIONS PEP shows distinct differences to CPAP with respect to its impact on ventilation distribution in healthy adults and CF subjects EIT might be used to individualise respiratory physiotherapy.

Influence of inspiration to expiration ratio on cyclic recruitment and derecruitment of atelectasis in a saline lavage model of acute respiratory distress syndrome*

OBJECTIVE Cyclic recruitment and derecruitment of atelectasis can occur during mechanical ventilation, especially in injured lungs. Experimentally, cyclic recruitment and derecruitment can be quantified by respiration-dependent changes in PaO2 (ΔPaO2), reflecting the varying intrapulmonary shunt fraction within the respiratory cycle. This study investigated the effect of inspiration to expiration ratio upon ΔPaO2 and Horowitz index. DESIGN Prospective randomized study. SETTING Laboratory investigation. SUBJECTS Piglets, average weight 30 ± 2 kg. INTERVENTIONS At respiratory rate 6 breaths/min, end-inspiratory pressure (Pendinsp) 40 cm H2O, positive end-expiratory pressure 5 cm H2O, and FIO2 1.0, measurements were performed at randomly set inspiration to expiration ratios during baseline healthy and mild surfactant depletion injury. Lung damage was titrated by repetitive surfactant washout to induce maximal cyclic recruitment and derecruitment as measured by multifrequency phase fluorimetry. Regional ventilation distribution was evaluated by electrical impedance tomography. Step changes in airway pressure from 5 to 40 cm H2O and vice versa were performed after lavage to calculate PO2-based recruitment and derecruitment time constants (TAU). MEASUREMENTS AND MAIN RESULTS In baseline healthy, cyclic recruitment and derecruitment could not be provoked, whereas in model acute respiratory distress syndrome, the highest ΔPaO2 were routinely detected at an inspiration to expiration ratio of 1:4 (range, 52-277 torr [6.9-36.9 kPa]). Shorter expiration time reduced cyclic recruitment and derecruitment significantly (158 ± 85 torr [21.1 ± 11.3 kPa] [inspiration to expiration ratio, 1:4]; 25 ± 12 torr [3.3 ± 1.6 kPa] [inspiration to expiration ratio, 4:1]; p < 0.0001), whereas the PaO2/FIO2 ratio increased (267 ± 50 [inspiration to expiration ratio, 1:4]; 424 ± 53 [inspiration to expiration ratio, 4:1]; p < 0.0001). Correspondingly, regional ventilation redistributed toward dependent lung regions (p < 0.0001). Recruitment was much faster (TAU: fast 1.6 s [78%]; slow 9.2 s) than derecruitment (TAU: fast 3.1 s [87%]; slow 17.7 s) (p = 0.0078). CONCLUSIONS Inverse ratio ventilation minimizes cyclic recruitment and derecruitment of atelectasis in an experimental model of surfactant-depleted pigs. Time constants for recruitment and derecruitment, and regional ventilation distribution, reflect these findings and highlight the time dependency of cyclic recruitment and derecruitment.

Effect of smoke inhalation on viscoelastic properties and ventilation distribution in sheep

Smoke inhalation injuries are the leading cause of mortality from burn injury. Airway obstruction due to mucus plugging and bronchoconstriction can cause severe ventilation inhomogeneity and worsen hypoxia. Studies describing changes of viscoelastic characteristics of the lung after smoke inhalation are missing. We present results of a new smoke inhalation device in sheep and describe pathophysiological changes after smoke exposure. Fifteen female Merino ewes were anesthetized and intubated. Baseline data using electrical impedance tomography and multiple-breath inert-gas washout were obtained by measuring ventilation distribution, functional residual capacity, lung clearance index, dynamic compliance, and stress index. Ten sheep were exposed to standardized cotton smoke insufflations and five sheep to sham smoke insufflations. Measured carboxyhemoglobin before inhalation was 3.87 +/- 0.28% and 5 min after smoke was 61.5 +/- 2.1%, range 50-69.4% ( P < 0.001). Two hours after smoke functional residual capacity decreased from 1,773 +/- 226 to 1,006 +/- 129 ml and lung clearance index increased from 10.4 +/- 0.4 to 14.2 +/- 0.9. Dynamic compliance decreased from 56.6 +/- 5.5 to 32.8 +/- 3.2 ml/ cmH(2)O. Stress index increased from 0.994 +/- 0.009 to 1.081 +/- 0.011 ( P < 0.01) ( all means +/- SE, P < 0.05). Electrical impedance tomography showed a shift of ventilation from the dependent to the independent lung after smoke exposure. No significant change was seen in the sham group. Smoke inhalation caused immediate onset in pulmonary dysfunction and significant ventilation inhomogeneity. The smoke inhalation device as presented may be useful for interventional studies.

A wavenumber independent boundary element method for an acoustic scattering problem

In this paper we consider the impedance boundary value problem for the Helmholtz equation in a half-plane with piecewise constant boundary data, a problem which models, for example, outdoor sound propagation over inhomogeneous. at terrain. To achieve good approximation at high frequencies with a relatively low number of degrees of freedom, we propose a novel Galerkin boundary element method, using a graded mesh with smaller elements adjacent to discontinuities in impedance and a special set of basis functions so that, on each element, the approximation space contains polynomials ( of degree.) multiplied by traces of plane waves on the boundary. We prove stability and convergence and show that the error in computing the total acoustic field is O( N-(v+1) log(1/2) N), where the number of degrees of freedom is proportional to N logN. This error estimate is independent of the wavenumber, and thus the number of degrees of freedom required to achieve a prescribed level of accuracy does not increase as the wavenumber tends to infinity.

A high-wavenumber boundary-element method for an acoustic scattering problem

In this paper we show stability and convergence for a novel Galerkin boundary element method approach to the impedance boundary value problem for the Helmholtz equation in a half-plane with piecewise constant boundary data. This problem models, for example, outdoor sound propagation over inhomogeneous flat terrain. To achieve a good approximation with a relatively low number of degrees of freedom we employ a graded mesh with smaller elements adjacent to discontinuities in impedance, and a special set of basis functions for the Galerkin method so that, on each element, the approximation space consists of polynomials (of degree $\nu$) multiplied by traces of plane waves on the boundary. In the case where the impedance is constant outside an interval $[a,b]$, which only requires the discretization of $[a,b]$, we show theoretically and experimentally that the $L_2$ error in computing the acoustic field on $[a,b]$ is ${\cal O}(\log^{\nu+3/2}|k(b-a)| M^{-(\nu+1)})$, where $M$ is the number of degrees of freedom and $k$ is the wavenumber. This indicates that the proposed method is especially commendable for large intervals or a high wavenumber. In a final section we sketch how the same methodology extends to more general scattering problems.

Debonding detection of FRP strengthened concrete beams by using impedance measurements and an ensemble PSO adaptive spectral model

An impedance-based midspan debonding identification method for RC beams strengthened with FRP strips is presented in this paper using piezoelectric ceramic (PZT) sensor?actuators. To reach this purpose, firstly, a two-dimensional electromechanical impedance model is proposed to predict the electrical admittance of the PZT transducer bonded to the FRP strips of an RC beam. Considering the impedance is measured in high frequencies, a spectral element model of the bonded-PZT?FRP strengthened beam is developed. This model, in conjunction with experimental measurements of PZT transducers, is used to present an updating methodology to quantitatively detect interfacial debonding of these kinds of structures. To improve the performance and accuracy of the detection algorithm in a challenging problem such as ours, the structural health monitoring approach is solved with an ensemble process based on particle of swarm. An adaptive mesh scheme has also been developed to increase the reliability in locating the area in which debonding initiates. Predictions carried out with experimental results have showed the effectiveness and potential of the proposed method to detect prematurely at its earliest stages a critical failure mode such as that due to midspan debonding of the FRP strip.

A simplified spectral approachfor impedance-based damage identification of frp-strengthened rc beams

Hoy en día, el refuerzo y reparación de estructuras de hormigón armado mediante el pegado de bandas de polímeros reforzados con fibras (FRP) se emplea cada vez con más frecuencia a causa de sus numerosas ventajas. Sin embargo, las vigas reforzadas con esta técnica pueden experimentar un modo de fallo frágil a causa del despegue repentino de la banda de FRP a partir de una fisura intermedia. A pesar de su importancia, el número de trabajos que abordan el estudio de este mecanismo de fallo y su monitorización es muy limitado. Por ello, el desarrollo de metodologías capaces de monitorizar a largo plazo la adherencia de este refuerzo a las estructuras de hormigón e identificar cuándo se inicia el despegue de la banda constituyen un importante desafío a abordar. El principal objetivo de esta tesis es la implementación de una metodología fiable y efectiva, capaz de detectar el despegue de una banda de FRP en una viga de hormigón armado a partir de una fisura intermedia. Para alcanzar este objetivo se ha implementado un procedimiento de calibración numérica a partir de ensayos experimentales. Para ello, en primer lugar, se ha desarrollado un modelo numérico unidimensional simple y no costoso representativo del comportamiento de este tipo vigas de hormigón reforzadas con FRP, basado en un modelo de fisura discreta para el hormigón y el método de elementos espectrales. La formación progresiva de fisuras a flexion y el consiguiente despegue en la interface entre el hormigón y el FRP se formulan mediante la introducción de un nuevo elemento capaz de representar ambos fenómenos simultáneamente sin afectar al procedimiento numérico. Además, con el modelo propuesto, se puede obtener de una forma sencilla la respuesta dinámica en altas frecuencias de este tipo de estructuras, lo cual puede hacer muy útil su uso como herramienta de diagnosis y detección del despegue en su fase inicial mediante una monitorización de la variación de las características dinámicas locales de la estructura. Un método de evaluación no destructivo muy prometedor para la monitorización local de las estructuras es el método de la impedancia usando sensores-actuadores piezoeléctricos (PZT). La impedancia eléctrica de los sensores PZT se puede relacionar con la impedancia mecánica de las estructuras donde se encuentran adheridos Ya que la impedancia mecánica de una estructura se verá afectada por su deterioro, se pueden implementar indicadores de daño mediante una comparación del espectro de admitancia (inversa de la impedancia) a lo largo de distintas etapas durante el periodo de servicio de una estructura. Cualquier cambio en el espectro se podría interpretar como una variación en la integridad de la estructura. La impedancia eléctrica se mide a altas frecuencias con lo cual esta metodología debería ser muy sensible a la detección de estados de daño incipiente local, tal como se desea en la aplicación de este trabajo. Se ha implementado un elemento espectral PZT-FRP como extensión del modelo previamente desarrollado, con el objetivo de poder calcular numéricamente la impedancia eléctrica de sensores PZT adheridos a bandas de FRP sobre una viga de hormigón armado. El modelo, combinado con medidas experimentales captadas mediante sensores PZT, se implementa en el marco de una metodología de calibración de modelos para detectar cuantitativamente el despegue en la interfase entre una banda de FRP y una viga de hormigón. El procedimiento de optimización se resuelve empleando el método del enjambre cooperativo con un algoritmo bagging. Los resultados muestran una gran aproximación en la estimación del daño para el problema propuesto. Adicionalmente, se ha desarrollado también un método adaptativo para el mallado de elementos espectrales con el objetivo de localizar las zonas dañadas a partir de los resultados experimentales, el cual contribuye a aumentar la robustez y efectividad del método propuesto a la hora de identificar daños incipientes en su aparición inicial. Finalmente, se ha llevado a cabo un procedimiento de optimización multi-objetivo para detectar el despegue inicial en una viga de hormigón a escala real reforzada con FRP a partir de las impedancias captadas con una red de sensores PZT instrumentada a lo largo de la longitud de la viga. Cada sensor aporta los datos para definir cada una de las funciones objetivo que definen el procedimiento. Combinando el modelo previo de elementos espectrales con un algoritmo PSO multi-objetivo el procedimiento de detección de daño resultante proporciona resultados satisfactorios considerando la escala de la estructura y todas las incertidumbres características ligadas a este proceso. Los resultados obtenidos prueban la viabilidad y capacidad de los métodos antes mencionados y también su potencial en aplicaciones reales. Abstract Nowadays, the external bonding of fibre reinforced polymer (FRP) plates or sheets is increasingly used for the strengthening and retrofitting of reinforced concrete (RC) structures due to its numerous advantages. However, this kind of strengthening often leads to brittle failure modes being the most dominant failure mode the debonding induced by an intermediate crack (IC). In spite of its importance, the number of studies regarding the IC debonding mechanism and bond health monitoring is very limited. Methodologies able to monitor the long-term efficiency of bonding and successfully identify the initiation of FRP debonding constitute a challenge to be met. The main purpose of this thesisis the implementation of a reliable and effective methodology of damage identification able to detect intermediate crack debonding in FRP-strengthened RC beams. To achieve this goal, a model updating procedure based on numerical simulations and experimental tests has been implemented. For it, firstly, a simple and non-expensive one-dimensional model based on the discrete crack approach for concrete and the spectral element method has been developed. The progressive formation of flexural cracks and subsequent concrete-FRP interfacial debonding is formulated by the introduction of a new element able to represent both phenomena simultaneously without perturbing the numerical procedure. Furthermore, with the proposed model, high frequency dynamic response for these kinds of structures can also be obtained in a very simple and non-expensive way, which makes this procedure very useful as a tool for diagnoses and detection of debonding in its initial stage by monitoring the change in local dynamic characteristics. One very promising active non-destructive evaluation method for local monitoring is impedance-based structural health monitoring(SHM)using piezoelectric ceramic (PZT) sensor-actuators. The electrical impedance of the PZT can be directly related to the mechanical impedance of the host structural component where the PZT transducers are attached. Since the structural mechanical impedance will be affected by the presence of structural damage, comparisons of admittance (inverse of impedance) spectra at various times during the service period of the structure can be used as damage indicator. Any change in the spectra might be an indication of a change in the structural integrity. The electrical impedance is measured at high frequencies with which this methodology appears to be very sensitive to incipient damage in structural systems as desired for our application. Abonded-PZT-FRP spectral beam element approach based on an extension of the previous discrete crack approach is implemented in the calculation of the electrical impedance of the PZT transducer bonded to the FRP plates of a RC beam. This approach in conjunction with the experimental measurements of PZT actuator-sensors mounted on the structure is used to present an updating methodology to quantitatively detect interfacial debonding between a FRP strip and the host RC structure. The updating procedure is solved by using an ensemble particle swarm optimization approach with abagging algorithm, and the results demonstrate a big improvement for the performance and accuracy of the damage detection in the proposed problem. Additionally, an adaptive strategy of spectral element mesh has been also developed to detect damage location with experimental results, which shows the robustness and effectiveness of the proposed method to identify initial and incipient damages at its early stage. Lastly, multi-objective optimization has been carried out to detect debonding damage in a real scale FRP-strengthened RC beam by using impedance signatures. A net of PZT sensors is distributed along the beam to construct impedance-based multiple objectives under gradually induced damage scenario. By combining the spectral element model presented previously and an ensemble multi-objective PSO algorithm, the implemented damage detection process yields satisfactory predictions considering the scale and uncertainties of the structure. The obtained results prove the feasibility and capability of the aforementioned methods and also their potentials in real engineering applications.

Data analysis in multiple-frequency bioelectrical impedance analysis

The performance of three analytical methods for multiple-frequency bioelectrical impedance analysis (MFBIA) data was assessed. The methods were the established method of Cole and Cole, the newly proposed method of Siconolfi and co-workers and a modification of this procedure. Method performance was assessed from the adequacy of the curve fitting techniques, as judged by the correlation coefficient and standard error of the estimate, and the accuracy of the different methods in determining the theoretical values of impedance parameters describing a set of model electrical circuits. The experimental data were well fitted by all curve-fitting procedures (r = 0.9 with SEE 0.3 to 3.5% or better for most circuit-procedure combinations). Cole-Cole modelling provided the most accurate estimates of circuit impedance values, generally within 1-2% of the theoretical values, followed by the Siconolfi procedure using a sixth-order polynomial regression (1-6% variation). None of the methods, however, accurately estimated circuit parameters when the measured impedances were low (<20 Omega) reflecting the electronic limits of the impedance meter used. These data suggest that Cole-Cole modelling remains the preferred method for the analysis of MFBIA data.

Factorization of elliptic boundary value problems by invariant embedding and application to overdetermined problems

Dissertação para obtenção do Grau de Doutor em Matemática

A systematic and feasible method for computing nuclear contributions to electrical properties of polyatomic molecules

An analytic method to evaluate nuclear contributions to electrical properties of polyatomic molecules is presented. Such contributions control changes induced by an electric field on equilibrium geometry (nuclear relaxation contribution) and vibrational motion (vibrational contribution) of a molecular system. Expressions to compute the nuclear contributions have been derived from a power series expansion of the potential energy. These contributions to the electrical properties are given in terms of energy derivatives with respect to normal coordinates, electric field intensity or both. Only one calculation of such derivatives at the field-free equilibrium geometry is required. To show the useful efficiency of the analytical evaluation of electrical properties (the so-called AEEP method), results for calculations on water and pyridine at the SCF/TZ2P and the MP2/TZ2P levels of theory are reported. The results obtained are compared with previous theoretical calculations and with experimental values

Rate and processes of river network rearrangement during incipient faulting: The case of the Cahabon river, Guatemala

Deeply incised river networks are generally regarded as robust features that are not easily modified by erosion or tectonics. Although the reorganization of deeply incised drainage systems has been documented, the corresponding importance with regard to the overall landscape evolution of mountain ranges and the factors that permit such reorganizations are poorly understood. To address this problem, we have explored the rapid drainage reorganization that affected the Cahabon River in Guatemala during the Quaternary. Sediment-provenance analysis, field mapping, and electrical resistivity tomography (ERT) imaging are used to reconstruct the geometry of the valley before the river was captured. Dating of the abandoned valley sediments by the Be-10-Al-26 burial method and geomagnetic polarity analysis allow us to determine the age of the capture events and then to quantify several processes, such as the rate of tectonic deformation of the paleovalley, the rate of propagation of post-capture drainage reversal, and the rate at which canyons that formed at the capture sites have propagated along the paleovalley. Transtensional faulting started 1 to 3 million years ago, produced ground tilting and ground faulting along the Cahabon River, and thus generated differential uplift rate of 0.3 +/- 0.1 up to 0.7 +/- 0.4 mm . y(-1) along the river's course. The river responded to faulting by incising the areas of relative uplift and depositing a few tens of meters of sediment above the areas of relative subsidence. Then, the river experienced two captures and one avulsion between 700 ky and 100 ky. The captures breached high-standing ridges that separate the Cahabon River from its captors. Captures occurred at specific points where ridges are made permeable by fault damage zones and/or soluble rocks. Groundwater flow from the Cahabon River down to its captors likely increased the erosive power of the captors thus promoting focused erosion of the ridges. Valley-fill formation and capture occurred in close temporal succession, suggesting a genetic link between the two. We suggest that the aquifers accumulated within the valley-fills, increased the head along the subterraneous system connecting the Cahabon River to its captors, and promoted their development. Upon capture, the breached valley experienced widespread drainage reversal toward the capture sites. We attribute the generalized reversal to combined effects of groundwater sapping in the valley-fill, axial drainage obstruction by lateral fans, and tectonic tilting. Drainage reversal increased the size of the captured areas by a factor of 4 to 6. At the capture sites, 500 m deep canyons have been incised into the bedrock and are propagating upstream at a rate of 3 to 11 mm . y(-1) deepening at a rate of 0.7 to 1 5 mm . y(-1). At this rate, 1 to 2 million years will be necessary for headward erosion to completely erase the topographic expression of the paleovalley. It is concluded that the rapid reorganization of this drainage system was made possible by the way the river adjusted to the new tectonic strain field, which involved transient sedimentation along the river's course. If the river had escaped its early reorganization and had been given the time necessary to reach a new dynamic equilibrium, then the transient conditions that promoted capture would have vanished and its vulnerability to capture would have been strongly reduced.

A filtering method to correct time-lapse 3D ERT data and improve imaging of natural aquifer dynamics

We have developed a processing methodology that allows crosshole ERT (electrical resistivity tomography) monitoring data to be used to derive temporal fluctuations of groundwater electrical resistivity and thereby characterize the dynamics of groundwater in a gravel aquifer as it is infiltrated by river water. Temporal variations of the raw ERT apparent-resistivity data were mainly sensitive to the resistivity (salinity), temperature and height of the groundwater, with the relative contributions of these effects depending on the time and the electrode configuration. To resolve the changes in groundwater resistivity, we first expressed fluctuations of temperature-detrended apparent-resistivity data as linear superpositions of (i) time series of riverwater-resistivity variations convolved with suitable filter functions and (ii) linear and quadratic representations of river-water-height variations multiplied by appropriate sensitivity factors; river-water height was determined to be a reliable proxy for groundwater height. Individual filter functions and sensitivity factors were obtained for each electrode configuration via deconvolution using a one month calibration period and then the predicted contributions related to changes in water height were removed prior to inversion of the temperature-detrended apparent-resistivity data. Applications of the filter functions and sensitivity factors accurately predicted the apparent-resistivity variations (the correlation coefficient was 0.98). Furthermore, the filtered ERT monitoring data and resultant time-lapse resistivity models correlated closely with independently measured groundwater electrical resistivity monitoring data and only weakly with the groundwater-height fluctuations. The inversion results based on the filtered ERT data also showed significantly less inversion artefacts than the raw data inversions. We observed resistivity increases of up to 10% and the arrival time peaks in the time-lapse resistivity models matched those in the groundwater resistivity monitoring data.

Integration of local-scale hydrological and regional-scale geophysical based on a nonlinear Bayesian sequential simulation approach

Significant progress has been made with regard to the quantitative integration of geophysical and hydrological data at the local scale. However, extending the corresponding approaches to the regional scale represents a major, and as-of-yet largely unresolved, challenge. To address this problem, we have developed a downscaling procedure based on a non-linear Bayesian sequential simulation approach. The basic objective of this algorithm is to estimate the value of the sparsely sampled hydraulic conductivity at non-sampled locations based on its relation to the electrical conductivity, which is available throughout the model space. The in situ relationship between the hydraulic and electrical conductivities is described through a non-parametric multivariate kernel density function. This method is then applied to the stochastic integration of low-resolution, re- gional-scale electrical resistivity tomography (ERT) data in combination with high-resolution, local-scale downhole measurements of the hydraulic and electrical conductivities. Finally, the overall viability of this downscaling approach is tested and verified by performing and comparing flow and transport simulation through the original and the downscaled hydraulic conductivity fields. Our results indicate that the proposed procedure does indeed allow for obtaining remarkably faithful estimates of the regional-scale hydraulic conductivity structure and correspondingly reliable predictions of the transport characteristics over relatively long distances.