Desarrollo de métodos numérico-interpretativos para la realización de ensayos de bombeo

Los ensayos de bombeo son, sin lugar a dudas, una de las pruebas más fiables y de mayor interés que se hacen en el medio físico. No son pruebas estrictamente puntuales, dado que el bombeo atrae flujo desde distancias lejanas al pozo, la prueba tiene una excelente representatividad espacial. Los métodos de interpretación mediante ensayos de bombeo se empezaron a plantear en la primera mitad del pasado siglo. Con los ensayos de bombeo se puede calcular la transmisividad y coeficiente de almacenamiento de las formaciones acuíferas y suministran información sobre el tipo de acuífero, la calidad constructiva del pozo de extracción, la existencia de barreras impermeable o bordes de recarga próximos, e incluso en algunas circunstancias permiten el cálculo del área de embalse subterráneo. Desde mediados del siglo 20 existe una eficaz y abundante gama de métodos analítico-interpretativos de ensayos de bombeo, tanto en régimen permanente como transitorio. Estos métodos son ampliamente conocidos y están muy experimentados a lo largo de muchos países, sin embargo, hoy día, podrían utilizarse modelos de flujo para la interpretación, logrando la misma fiabilidad e incluso mejores posibilidades de análisis. Muchos ensayos que no pueden interpretarse porque las configuraciones del medio son demasiado complejas y no están disponibles, o no es posible, el desarrollo de métodos analíticos, tienen buena adaptación y en ocasiones muy fácil solución haciendo uso de los métodos numéricos de simulación del flujo. En esta tesis se ha buscado una vía de interpretar ensayos de bombeo haciendo uso de modelos de simulación del flujo. Se utiliza el modelo universal MODFLOW del United States Geological Survey, en el cual se configura una celda de simulación y mallado particularmente adecuados para el problema a tratar, se valida con los métodos analíticos existentes. Con la célula convenientemente validada se simulan otros casos en los que no existen métodos analíticos desarrollados dada la complejidad del medio físico a tratar y se sacan las oportunas conclusiones. Por último se desarrolla un modelo específico y la correspondiente aplicación de uso general para la interpretación numérica de ensayos de bombeo tanto con las configuraciones normales como con configuraciones complejas del medio físico. ABSTRACT Pumping tests are, without doubt, one of the most reliable and most interesting tests done in the physical environment. They are not strictly anecdotal evidence, since pumping flow attracts from far distances to the well, the test has excellent spatial representation. Methods of interpretation by pumping tests began to arise in the first half of last century. With pumping tests, can be calculated transmissivity and storage coefficient of the aquifer formations, and provide information on the type of aquifer, the construction quality of the well, the existence of waterproof barriers or borders next recharge, and even in some circumstances allow calculating the area of underground reservoir. Since the mid-20th century there is effective and abundant range of analytical interpretative pumping tests, both in steady state and transient methods. These methods are very widely known and experienced over many countries, however, nowadays, may flow models used for interpretation, obtaining equally reliable or even better possibilities for analysis. Many trials cannot be interpreted as environmental settings are too complex and are not available, or not possible, the development of analytical methods, have good adaptation and sometimes very easily solved using numerical flow simulation methods. This thesis has sought a way to interpret pumping tests using flow simulation models. MODFLOW universal model of United States Geological Survey, in which a simulation cell and meshing particularly suitable for the problem to be treated, is validated with existing analytical methods used is set. With suitably validated cell other cases where there are no analytical methods developed given the complexity of the physical environment to try and draw appropriate conclusions are simulated. Finally, a specific model and the corresponding application commonly used for numerical interpretation of pumping tests both with normal settings as complex configurations of the physical environment is developed.

Amelioration of a degraded acid soil from SW Spain by no-tillage and Ca-amendment

En los suelos, el exceso de acidez lleva asociado deficiencias en ciertos nutrientes y una alta disponibilidad de aluminio, tóxico para los cultivos propios del ambiente mediterráneo. Su laboreo, provoca la pérdida de materia orgánica (MO), deteriora su estructura y reduce la actividad biológica, provocando en última instancia una menor calidad del suelo. Es de esperar pues que cuando se labran suelos ácidos, sus problemáticas particulares tiendan a agravarse. En nuestra zona de estudio, la “raña” de Cañamero (Extremadura, España), predominan los suelos muy ácidos y degradados por un laboreo inadecuado. Las rañas constituyen amplias plataformas casi horizontales, con unos suelos muy viejos (Palexerults), que se caracterizan por tener el complejo de cambio dominado por el aluminio, y un pH ácido que decrece en profundidad. Poseen un potente horizonte Bt rico en arcillas caoliníticas, que propicia que en periodos con exceso de lluvia, se generen capas colgadas de agua cercanas a la superficie. En torno a los años 1940’s estos suelos, que previamente sostenían un alcornocal, o su matorral de sustitución, se pusieron en cultivo. El laboreo aceleró la mineralización de la materia orgánica, agravó los problemas derivados del exceso de acidez y condujo al abandono de los campos cultivados por falta de productividad. Para recuperar la calidad de estos suelos degradados y obtener unos rendimientos compatibles con su uso agrícola es necesario, por un lado, aplicar enmiendas que eleven el pH y reduzcan la toxicidad del aluminio y, por otro, favorecer el incremento en el contenido en MO. En 2005 se implantó en esta raña un ensayo de campo para estudiar la influencia del no laboreo y de la utilización de una enmienda cálcica en parámetros relacionados con la calidad del suelo en un cultivo forrajero. El diseño experimental fue en parcelas divididas con cuatro repeticiones donde el factor principal fue el tipo de laboreo, no laboreo (NL) frente a laboreo convencional (LC), y el factor secundario el uso o no de una enmienda cálcica. La enmienda consistió básicamente en una mezcla de espuma de azucarería y yeso rojo y se incorporó al comienzo del ensayo hasta los 7 cm de profundidad. Desde el comienzo del ensayo el NL influyó positivamente en el contenido de carbono orgánico total (COT) y particulado (COP), mientras que la enmienda tuvo una ligera influencia al principio del ensayo en ambos pero su efecto positivo se desvaneció con el paso del tiempo. Los mayores contenidos en COT y POC se observaron cuando se combinó el NL con la enmienda. La enmienda incrementó con rapidez el pH, y el Ca, y disminuyó el contenido en aluminio hasta una profundidad de 50 cm, incluso en NL, y mejoró ligeramente la agregación del suelo. El NL por sí solo, gracias al aumento en POC, TOC y las proteínas del suelo relacionadas con la glomalina (PSRG), que son capaces de formar compuestos estables no tóxicos con el aluminio, también contribuyó a la reducción de la toxicidad de aluminio en la capa más superficial. Cuando en las campañas con exceso de precipitaciones se generaron capas colgadas de agua próximas a la superficie, el NL generó unas condiciones más favorables para la germinación y desarrollo del cultivo, resultando en una producción más alta que el LC. A ello contribuyó la mayor capacidad de almacenamiento de agua y la mayor transmisividad de esta hacia abajo, en la capa más superficial (0-5 cm) que propició una menor saturación por agua que el LC. Respecto a los parámetros relacionados con la agregación, el NL aumentó los macroagregados hasta los 10 cm de profundidad y favoreció la acumulación de CO y N en todas las fracciones de tamaño de agregados. Sin embargo, la recuperación del grado de macroagregación tras el cese del laboreo resulta lenta en comparación con otros suelos, posiblemente debido al bajo contenido en arcilla en el horizonte Ap. En comparación con el NL, la enmienda mostró también un efecto positivo, aunque muy ligero, en la agregación del suelo. En contradicción con otros estudios en suelos ácidos, nuestros resultados indican la existencia de una jerarquía de agregados, y destacan el papel importante de la MO en la mejora de la agregación. Tanto el NL como la enmienda favorecieron por separado varias propiedades químicas, físicas y biológicas del suelo, pero, en general, encontramos los mayores beneficios con su uso combinado. Además, a largo plazo el efecto positivo de NL en las propiedades del suelo fue en aumento, mientras que el efecto beneficioso de la enmienda se limitó básicamente a las propiedades químicas y se desvaneció en pocos años. Destacamos que las condiciones meteorológicas a lo largo del ensayo beneficiaron la producción de biomasa en NL, y en consecuencia las propiedades relacionadas con la materia orgánica, por lo que son un factor a tener en cuenta a la hora de evaluar los efectos de la enmienda y el laboreo sobre las propiedades del suelo, especialmente en zonas donde esas condiciones son muy variables entre una campaña y otra. Los resultados de este estudio han puesto de manifiesto que el NL no ha mermado la eficacia de la enmienda caliza, posiblemente gracias a la alta solubilidad de la enmienda aplicada, es más, el manejo con NL y enmienda es el que ha favorecido en mayor medida ciertos parámetros de calidad del suelo. Por el contrario el LC sí parece anular los beneficios de la enmienda en relación con las propiedades relacionadas con la MO. Por tanto, cabe concluir que la combinación de NL y la enmienda es una práctica adecuada para mejorar las propiedades químicas y físicas de suelos ácidos degradados por el laboreo. ABSTRACT Excessive acidity in soils is associated with deficiencies in certain nutrients and high concentrations of available aluminum, which is toxic for most Mediterranean crops. Tilling these soils results in the loss of soil organic matter (SOM), damages soil structure and reduces biological activity, ultimately degrading soil quality. It is expected, therefore, that when acid soils are tilled, their particular problems will tend to get worse. In our study area, the "Cañamero’s Raña” (Extremadura, Spain), acid soils degraded by an inappropriate tillage prevail. Rañas are large and flat platforms with very old soils (Palexerults), which are characterized by an exchange complex dominated by aluminum and an acid pH which decreases with depth. These soils have a strong Bt horizon rich in kaolinite clays, which encourages the formation of perched water-tables near the soil surface during periods of excessive rain. During the first third of the 20th century, these soils, that previously supported cork oak or its scrub replacement, were cultivated. Tillage accelerated the mineralization of the SOM, aggravating the problems of excessive acidity, which finally led to the abandonment of the land due to low productivity. To recover the quality of these degraded soils and to obtain consistent yields it is necessary, first, to apply amendments to raise the pH and reduce aluminum toxicity, and second to encourage the accumulation of SOM. In 2005 a field trial was established in the Raña to study the influence of no-tillage and the use of a Ca-amendment on soil quality related parameters in a forage crop agrosystem. The experimental design was a split-plot with four replicates where the main factor was tillage type, no-tillage (NT) versus traditional tillage (TT) and the secondary factor was the use or not of a Ca-amendment. The Ca-amendment was a mixture of sugar foam and red gypsum that was incorporated into the top 7 cm of the soil. Since the beginning of the experiment, NT had a positive influence on total and particulate organic carbon (TOC and POC, respectively), while the Ca-amendment had a small positive influence at the beginning of the study but its effect diminished with time. The highest TOC and POC contents were observed when NT and the Ca-amendment were combined. The Ca-amendment, even under NT, rapidly increased pH and Ca, and decreased the aluminum content to a depth of 50 cm, as well as improving soil aggregation slightly. NT, due to the increased POC, TOC and Glomalin-related soil proteins (GRSP), which can form stable non-toxic compounds with aluminum, also contributed to the reduction of aluminum toxicity in the upper layer. When perched water-tables near the soil surface were formed in campaigns with excessive rainfall, NT provided more favorable conditions for germination and crop development, resulting in higher yields compared with TT. This was directly related to the higher water storage capacity and the greater transmissivity of the water downwards from the upper layers, which led to lower water saturation under NT compared with TT. With regards to the aggregation-related parameters, NT increased macroaggregation to a depth of 10 cm and favored the accumulation of OC and N in all aggregate size fractions. However, the degree of recovery of macroaggregation after tillage ceased was slow compared with other soils, possibly due to the low clay content in the Ap horizon. Compared with NT, the Ca-amendment had a slight positive effect on soil aggregation. In contrast to other studies in acid soils, our results indicate the existence of an aggregate hierarchy, and highlight the important role of SOM in improving aggregation. Both NT and the Ca-amendment separately favored various chemical, physical and biological soil properties, but in general we found the greatest benefits when the two treatments were combined. In addition, the positive effect of NT on soil properties increased with time, while the beneficial effect of the Ca-amendment, which was limited to the chemical properties, vanished after a few years. It is important to note that the meteorological conditions throughout the experiment benefited biomass production under NT and, as a consequence, organic matter related properties. This suggests that meteorological conditions are a factor to consider when evaluating the effects of Ca-amendments and tillage on soil properties, especially in areas where such conditions vary significantly from one campaign to another. The results of this study show that NT did not diminish the effectiveness of the Ca-amendment, possibly due to the high solubility of the selected amendment. Moreover, the combination of NT and the Ca-amendment was actually the management that favored certain soil quality parameters the most. By contrast, TT seemed to nullify the benefits of the Ca-amendment with regards to the OM related properties. In conclusion, the combination of NT and the application of a Ca-amendment is an advisable practice for improving the chemical and physical properties of acid soils degraded by tillage.

Concepção de um receptor de cavidade para concentração de energia solar para aplicação em reatores químicos.

Este trabalho dimensionou um receptor de cavidade para uso como reator químico de um ciclo de conversão de energia solar para energia química. O vetor energético proposto é o hidrogênio. Isso implica que a energia solar é concentrada em um dispositivo que absorve a radiação térmica e a transforma em energia térmica para ativar uma reação química endotérmica. Essa reação transforma o calor útil em gás hidrogênio, que por sua vez pode ser utilizado posteriormente para geração de outras formas de energia. O primeiro passo foi levantar os pares metal/óxido estudados na literatura, cuja finalidade é ativar um ciclo termoquímico que possibilite produção de hidrogênio. Esses pares foram comparados com base em quatro parâmetros, cuja importância determina o dimensionamento de um receptor de cavidade. São eles: temperatura da reação; estado físico de reagentes e produtos; desgaste do material em ciclos; taxa de reação de hidrólise e outros aspectos. O par escolhido com a melhor avaliação no conjunto dos parâmetros foi o tungstênio e o trióxido de tungstênio (W/WO3). Com base na literatura, foi determinado um reator padrão, cujas características foram analisadas e suas consequências no funcionamento do receptor de cavidade. Com essa análise, determinaram-se os principais parâmetros de projeto, ou seja, a abertura da cavidade, a transmissividade da janela, e as dimensões da cavidade. Com base nos resultados anteriores, estabeleceu-se um modelo de dimensionamento do sistema de conversão de energia solar em energia útil para um processo químico. Ao se analisar um perfil de concentração de energia solar, calculou-se as eficiências de absorção e de perdas do receptor, em função da área de abertura de um campo de coleta de energia solar e da radiação solar disponível. Esse método pode ser empregado em conjunto com metodologias consagradas e dados de previsão de disponibilidade solar para estudos de concentradores de sistemas de produção de hidrogênio a partir de ciclos termoquímicos.

Metodologias utilizadas no estudo do escoamento em aquíferos cársicos e o caso prático da captação do Olho de Mira (Maciço Calcário Estremenho)

Tese de mestrado em Geologia Aplicada, apresentada à Universidade de Lisboa, através da Faculdade de Ciências, 2016

Determining Spatial Distribution and Physical Properties of the Vashon Advance Outwash near Mountlake Terrace, Washington

This study is aimed at determining the spatial distribution, physical properties, and groundwater conditions of the Vashon advance outwash (Qva) in the Mountlake Terrace, WA area. The Qva is correlative with the Esperance Sand, as defined at its type section; however, local variations in the Qva are not well-characterized (Mullineaux, 1965). While the Qva is a dense glacial unit with low compressibility and high frictional shear strength (Gurtowski and Boirum, 1989), the strength of this unit can be reduced when it becomes saturated (Tubbs, 1974). This can lead to caving or flowing in excavations, and on a larger scale, can lead to slope failures and mass-wasting when intersected by steep slopes. By studying the Qva, we can better predict how it will behave under certain conditions, which will be beneficial to geologists, hydrogeologists, engineers, and environmental scientists during site assessments and early phases of project planning. In this study, I use data from 27 geotechnical borings from previous field investigations and C-Tech Corporation’s EnterVol software to create three-dimensional models of the subsurface geology in the study area. These models made it possible to visualize the spatial distribution of the Qva in relation to other geologic units. I also conducted a comparative study between data from the borings and generalized published data on the spatial distribution, relative density, soil classification, grain-size distribution, moisture content, groundwater conditions, and aquifer properties of the Qva. I found that the elevation of the top of the Qva ranges from 247 to 477 ft. I found that the Qva is thickest where the modern topography is high, and is thinnest where the topography is low. The thickness of the Qva ranges from absent to 242 ft. Along the northern, east-west trending transect, the Qva thins to the east as it rises above a ridge composed of Pre- Vashon glacial deposits. Along the southern, east-west trending transect, the Qva pinches out against a ridge composed of pre-Vashon interglacial deposits. Two plausible explanations for this ridge are paleotopography and active faulting associated with the Southern Whidbey Fault Zone. Further investigations should be done using geophysical methods and the modeling methods described in this study to determine the nature of this ridge. The relative density of the Qva in the study area ranges from loose to very dense, with the loose end of the spectrum probably relating to heave in saturated sands. I found subtle correlations between density and depth. Volumetric analysis of the soil groups listed in the boring logs indicate that the Qva in the study area is composed of approximately 9.5% gravel, 89.3% sand, and 1.2% silt and clay. The natural moisture content ranges from 3.0 to 35.4% in select samples from the Qva. The moisture content appears to increase with depth and fines content. The water table in the study area ranges in elevation from 231.9 to 458 ft, based on observations and measurements recorded in the boring logs. The results from rising-head and falling-head slug tests done at a single well in the study area indicate that the geometric mean of hydraulic conductivity is 15.93 ft/d (5.62 x 10-03 cm/s), the storativity is 3.28x10-03, and the estimated transmissivity is 738.58 ft2/d in the vicinity of this observation well. At this location, there was 1.73 ft of seasonal variation in groundwater elevation between August 2014 and March 2015.

Development of a Conceptual Site Model for Evaluating Seawater Intrusion under Sea Level Rise Scenarios Using Analytical Methods

Island County is located in the Puget Sound of Washington State and includes several islands, the largest of which is Whidbey Island. Central Whidbey Island was chosen as the project site, as residents use groundwater for their water supply and seawater intrusion near the coast is known to contaminate this resource. In 1989, Island County adopted a Saltwater Intrusion Policy and used chloride concentrations in existing wells in order to define and map “risk zones.” In 2005, this method of defining vulnerability was updated with the use of water level elevations in conjunction with chloride concentrations. The result of this work was a revised map of seawater intrusion vulnerability that is currently in use by Island County. This groundwater management strategy is defined as trigger-level management and is largely a reactive tool. In order to evaluate trends in the hydrogeologic processes at the site, including seawater intrusion under sea level rise scenarios, this report presents a workflow where groundwater flow and discharge to the sea are quantified using a revised conceptual site model. The revised conceptual site model used several simplifying assumptions that allow for first-order quantitative predictions of seawater intrusion using analytical methods. Data from water well reports included lithologic and well construction information, static water levels, and aquifer tests for specific capacity. Results from specific capacity tests define the relationship between discharge and drawdown and were input for a modified Theis equation to solve for transmissivity (Arihood, 2009). Components of the conceptual site model were created in ArcGIS and included interpolation of water level elevation, creation of groundwater basins, and the calculation of net recharge and groundwater discharge for each basin. The revised conceptual site model was then used to hypothesize regarding hydrogeologic processes based on observed trends in groundwater flow. Hypotheses used to explain a reduction in aquifer thickness and hydraulic gradient were: (1) A large increase in transmissivity occurring near the coast. (2) The reduced aquifer thickness and hydraulic gradient were the result of seawater intrusion. (3) Data used to create the conceptual site model were insufficient to resolve trends in groundwater flow. For Hypothesis 2, analytical solutions for groundwater flow under Dupuit assumptions were applied in order to evaluate seawater intrusion under projected sea level rise scenarios. Results indicated that a rise in sea level has little impact on the position of a saltwater wedge; however, a reduction in recharge has significant consequences. Future work should evaluate groundwater flow using an expanded monitoring well network and aquifer recharge should be promoted by reducing surface water runoff.

Analysis of the Hydrogeological Conditions and the Effects of Development on the Hydraulic Head at Redmond Ridge East

Redmond Ridge East (RRE) is a large-scale master plan community in East King County, WA. In this report, I evaluate the spatial variability of the Quaternary Advance Outwash (Qva) at RRE and the time-series data for 16 water wells with the intent to better understand groundwater below the RRE area. I investigate changes between pre- and post-development conditions through the determination of temporal changes in annual water level, annual water level fluctuations, hydraulic head response to precipitation, and ambient drainage of the aquifer. I also perform a basic analysis of the annual aquifer recharge and a determination for the storage through the implementation of the water table fluctuation (WTF) method. Associated Earth Sciences (AESI) was tasked with monitoring the geological and environmental impacts during the development of RRE and collected the data I use in this report. AESI involvement in monitoring began in 1998 and extends to the present. Sixteen wells were identified in the RRE area with adequate temporal data to conduct the analysis. A comparison of the well logs and aquifer testing data allowed local variations in the Qva to be mapped. The WTF was used to determine a range of reasonable specific yield values for locations where the Qva was unconfined. Yearly average of the seasonal water level high and lows, and the fluctuations were quantified. Temporal relationships were established through linear regression. The average water level was found to be increasing in some locations, and the corresponding fluctuations were found to decrease. However, no clear change between pre- and post-development was observed. The response of hydraulic head to precipitation was investigated through an analysis of hydrographs for ten wells. Periods of consistent response and the corresponding precipitation during each period were delineated. A linear relationship between precipitation and water level change was determined. The threshold precipitation under which there is a positive response in the hydraulic head was established. No observable changes were apparent between pre- and post-development conditions. The ambient drainage for the Qva was calculated using recessional periods on the hydrograph. The transmissivity of Qva varies with thickness of the overlying lodgment till and thickness of the Qva, itself. Water level fluctuations observed in the Qva are consistent with regional observations. Localized areas in the Qva display the large 10 foot fluctuations and these anomalies are likely due to a combination of the local variability in the storativity as well as the concentration and channeling of water due to geographical variations in the Qva and the overlying topography. All trends seen in the RRE area remained relatively constant through time. There was no evidence showing an effect of development on the hydraulic head at RRE. This implies that the style and distribution of infiltration has not changed as a result of development, and that any measures in place are properly mitigating the effects of development on the RRE region.

Ground water model calibration using pilot points and regularization

Use of nonlinear parameter estimation techniques is now commonplace in ground water model calibration. However, there is still ample room for further development of these techniques in order to enable them to extract more information from calibration datasets, to more thoroughly explore the uncertainty associated with model predictions, and to make them easier to implement in various modeling contexts. This paper describes the use of pilot points as a methodology for spatial hydraulic property characterization. When used in conjunction with nonlinear parameter estimation software that incorporates advanced regularization functionality (such as PEST), use of pilot points can add a great deal of flexibility to the calibration process at the same time as it makes this process easier to implement. Pilot points can be used either as a substitute for zones of piecewise parameter uniformity, or in conjunction with such zones. In either case, they allow the disposition of areas of high and low hydraulic property value to be inferred through the calibration process, without the need for the modeler to guess the geometry of such areas prior to estimating the parameters that pertain to them. Pilot points and regularization can also be used as an adjunct to geostatistically based stochastic parameterization methods. Using the techniques described herein, a series of hydraulic property fields can be generated, all of which recognize the stochastic characterization of an area at the same time that they satisfy the constraints imposed on hydraulic property values by the need to ensure that model outputs match field measurements. Model predictions can then be made using all of these fields as a mechanism for exploring predictive uncertainty.

The cost of uniqueness in groundwater model calibration

Calibration of a groundwater model requires that hydraulic properties be estimated throughout a model domain. This generally constitutes an underdetermined inverse problem, for which a Solution can only be found when some kind of regularization device is included in the inversion process. Inclusion of regularization in the calibration process can be implicit, for example through the use of zones of constant parameter value, or explicit, for example through solution of a constrained minimization problem in which parameters are made to respect preferred values, or preferred relationships, to the degree necessary for a unique solution to be obtained. The cost of uniqueness is this: no matter which regularization methodology is employed, the inevitable consequence of its use is a loss of detail in the calibrated field. This, ill turn, can lead to erroneous predictions made by a model that is ostensibly well calibrated. Information made available as a by-product of the regularized inversion process allows the reasons for this loss of detail to be better understood. In particular, it is easily demonstrated that the estimated value for an hydraulic property at any point within a model domain is, in fact, a weighted average of the true hydraulic property over a much larger area. This averaging process causes loss of resolution in the estimated field. Where hydraulic conductivity is the hydraulic property being estimated, high averaging weights exist in areas that are strategically disposed with respect to measurement wells, while other areas may contribute very little to the estimated hydraulic conductivity at any point within the model domain, this possibly making the detection of hydraulic conductivity anomalies in these latter areas almost impossible. A study of the post-calibration parameter field covariance matrix allows further insights into the loss of system detail incurred through the calibration process to be gained. A comparison of pre- and post-calibration parameter covariance matrices shows that the latter often possess a much smaller spectral bandwidth than the former. It is also demonstrated that, as all inevitable consequence of the fact that a calibrated model cannot replicate every detail of the true system, model-to-measurement residuals can show a high degree of spatial correlation, a fact which must be taken into account when assessing these residuals either qualitatively, or quantitatively in the exploration of model predictive uncertainty. These principles are demonstrated using a synthetic case in which spatial parameter definition is based oil pilot points, and calibration is Implemented using both zones of piecewise constancy and constrained minimization regularization. (C) 2005 Elsevier Ltd. All rights reserved.

Facilitating soil-and plant-water research through the use of TDR

The development of TDR for measurement of soil water content and electrical conductivity has resulted in a large shift in measurement methods for a breadth of soil and hydrological characterization efforts. TDR has also opened new possibilities for soil and plant research. Five examples show how TDR has enhanced our ability to conduct our soil- and plant-water research. (i) Oxygen is necessary for healthy root growth and plant development but quantitative evaluation of the factors controlling oxygen supply in soil depends on knowledge of the soil water content by TDR. With water content information we have modeled successfully some impact of tillage methods on oxygen supply to roots and their growth response. (ii) For field assessment of soil mechanical properties influencing crop growth, water content capability was added to two portable soil strength measuring devices; (a) A TDT (Time Domain Transmittivity)-equipped soil cone penetrometer was used to evaluate seasonal soil strengthwater content relationships. In conventional tillage systems the relationships are dynamic and achieve the more stable no-tillage relationships only relatively late in each growing season; (b) A small TDR transmission line was added to a modified sheargraph that allowed shear strength and water content to be measured simultaneously on the same sample. In addition, the conventional graphing procedure for data acquisition was converted to datalogging using strain gauges. Data acquisition rate was improved by more than a factor of three with improved data quality. (iii) How do drought tolerant plants maintain leaf water content? Non-destructive measurement of TDR water content using a flat serpentine triple wire transmission line replaces more lengthy procedures of measuring relative water content. Two challenges remain: drought-stressed leaves alter salt content, changing electrical conductivity, and drought induced changes in leaf morphology affect TDR measurements. (iv) Remote radar signals are reflected from within the first 2 cm of soil. Appropriate calibration of radar imaging for soil water content can be achieved by a parallel pair of blades separated by 8 cm, reaching 1.7 cm into soil and forming a 20 cm TDR transmission line. The correlation between apparent relative permittivity from TDR and synthetic aperture radar (SAR) backscatter coefficient was 0.57 from an airborne flyover. These five examples highlight the diversity in the application of TDR in soil and plant research.