Minimum, maximum and mean annual ground temperatures in permafrost around Abisko ANS between 2008-2009

Monitoring of permafrost has been ongoing since 1978 in the Abisko area, northernmost Sweden, when measurements of active layer thickness started. In 1980, boreholes were drilled in three mires in the area to record permafrost temperatures. Recordings were made twice per year, and the last data were obtained in 2002. During the International Polar Year (2007-2008), new boreholes were drilled within the 'Back to the Future' (BTF) and 'Thermal State of Permafrost' (TSP) projects that enabled year-round temperature monitoring. Mean annual ground temperatures (MAGT) in the mires are close to 0°C, ranging from -0.16 to -0.47°C at 5 m depth. Data from the boreholes show increasing ground temperatures in the upper and lower part by 0.4 to 1°C between 1980 and 2002. At one mire, permafrost thickness has decreased from 15 m in 1980 to ca. 9 m in 2009, with an accelerating thawing trend during the last decade.

I-82 Prosser Landslide Investigation Benton City, WA

The scope of this technical report is to establish the mechanisms by which the eastbound lanes of Interstate 82 at mile post (MP) 91.9 near Benton City continue to deform. Within the Washington State Department of Transportation (WSDOT), the area is known as the Prosser Landslide and has been an ongoing concern since the 1980s. Results from previous technical investigations have been conflicted or inconclusive as to whether landslide movement persists beneath or through the shear key-buttress or that pavement distress is related to swelling of a clay-rich unit that underlies the slope and interstate. For this report, the following steps were taken. First, I conducted a desk review of archived reports, memos, data, and drill logs from the original construction of I-82 and previous geotechnical investigations commissioned by WSDOT. Findings of this desk review are reported in Part III. Second, WSDOT drillers drilled two new boreholes at the Prosser Landslide site above the buttress and instrumentation was installed within the boreholes. Borehole logs produced from the 2013 drilling can be found in Appendix A of this report. Material retrieved from the suspected failure zone during drilling was tested at the WSDOT Materials Lab by WSDOT personnel for its mechanical properties including Atterberg limits, grain-size analysis, and residual shear strength (Appendix B). Samples were also analyzed for mineral content using X -ray powder diffraction (XRD). These data and observations are reported in Part III and Appendix C. Finally, using drill logs produced by WSDOT from the latest drilling and from historic drilling campaigns, I constructed a 2-D geologic model of the landslide site. This model is the basis for slope stability analysis reported in Part IV and Appendix D. This study concludes that the deformation observed in the eastbound lanes of I-82 could be the result of continued landslide movement, despite previous remediation efforts.

Locating a strand of the Seattle Fault Zone in southeast Seattle, Washington through topographic analysis, borehole analysis and ground penetrating radar

Contributing to the evaluation of seismic hazards, a previously unmapped strand of the Seattle Fault Zone (SFZ), cutting across the southwest side of Lake Washington and southeast Seattle, is located and characterized on the basis of bathymetry, borehole logs, and ground penetrating radar (GPR). Previous geologic mapping and geophysical analysis of the Seattle area have generally mapped the locations of some strands of the SFZ, though a complete and accurate understanding of locations of all individual strands of the fault system is still incomplete. A bathymetric scarp-like feature and co-linear aeromagnetic anomaly lineament defined the extent of the study area. A 2-dimensional lithology cross-section was constructed using six boreholes, chosen from suitable boreholes in the study area. In addition, two GPR transects, oblique to the proposed fault trend, served to identify physical differences in subsurface materials. The proposed fault trace follows the previously mapped contact between the Oligocene Blakeley Formation and Quaternary deposits, and topographic changes in slope. GPR profiles in Seward Park and across the proposed fault location show the contact between the Blakeley Formation and unconsolidated glacial deposits, but it does not constrain an offset. However, north-dipping beds in the Blakely Formation are consistent with previous interpretations of P-wave seismic profiles on Mercer Island and Bellevue, Washington. The profiles show the mapped location of the aeromagnetic lineament in Lake Washington and the inferred location of the steeply-dipping, high-amplitude bedrock reflector, representing a fault strand. This north-dipping reflector is likely the same feature identified in my analysis. I characterize the strand as a splay fault, antithetic to the frontal fault of the SFZ. This new fault may pose a geologic hazard to the region.

An investigation of contaminant distribution in subsurface marine wood waste at Mill-A South Terminal, Port of Everett, Washington

This project investigates the correlation between contaminants and the wood waste present in marine sediments off the shore of the Port of Everett in the former Weyerhaeuser Mill-A pulp mill site. The investigation includes the results of two field studies, which tested contaminant levels in 22 boreholes as well as several surface samples. The contaminants include heavy metals and wood waste byproducts. These results, along with 14 other bore logs, provide the framework for a three-dimensional site model, interpolating the full extent of the depositional units and organic and inorganic chemicals found at Mill-A. The sediments of interest are divided into five depositional units defined by the percent wood content and type of wood: native material (<5% wood), intermediate (<30% wood), sawdust (<30% wood), woodchips (<30% wood), and poorly sorted sands with silt (SM-SP) (0% wood). The contaminants include arsenic, 2,4-dimethylphenol, and total organic carbon. Three-dimensional modeling software, RockWorks, interpolated the discrete borehole data of sediment and contaminants assuming horizontal continuity between sampling locations. The sediment distribution was calculated within concentration ranges for each contaminant of concern. The lowest detection limits, the screening levels, and the cleanup levels defined these ranges. Total organic carbon served as a proxy to estimate the quantity of wood waste in the sediment. As a known byproduct of wood decomposition, 2,4-dimethylphenol was expected to be more prevalent in the depositional units with more wood waste. Finally, arsenic was a proxy for other contaminants to determine if contaminants at Mill-A are dominant in sediments with high percentages of wood waste. The volumetric distribution established that high levels of total organic carbon are present in the sediment with higher percentages of wood waste. This correlation was stronger in the decomposing sawdust-rich sediment than the woodchip-rich sediment. The 2,4-dimethylphenol concentrations above cleanup standards were dominant in the sawdust-rich, intermediate and native sediments. Concentrations of 2,4-dimethylphenol below cleanup levels characterized the native sediment. The distribution of arsenic showed no statistically significant correlation to wood content in sediment. These results do not support the hypothesis of contaminant-rich wood waste, as many of the high concentrations of contaminants were not in the wood-rich sediments. This suggests that the contaminants are more distributed among all depositional units at Mill-A rather than focused within sediments with a high percent of wood waste. Understanding the distribution of potentially toxic compounds with wood waste is important for restoring the Puget Sound waterways to a more habitable environment. Future studies should include new data to validate these results and to limit the uncertainty of the extent of contaminants. Future studies may also find motive in looking for a correlation between contaminants and grain size based on previous studies linking these characteristics. These investigations will benefit the current cleanup effort as well as future cleanup efforts at similarly contaminated waterways.

Palynostratigraphy of Ordovician strata, Canning Basin, Western Australia. Part One: acritarchs and prasinophytes

This paper is the initial part of a comprehensive bipartite monograph of palynomorphs (viz., acritarchs, prasinophyte phycomata, and chitinozoans) that are represented profusely in marine lower Palaeozoic strata of the Canning Basin, Western Australia. The prime aim is to establish a palynologically based zonal scheme for the Ordovician sequence as represented in five cored boreholes drilled through the Lower to Middle Ordovician strata of the central-northeastern Canning Basin. These strata embrace the Oepikodus communis through Phragmodus-Plectodina conodont zonal interval and comprise (in ascending order) the Willara, Goldwyer, and Nita formations, of inferred early Arenig to Llanvirn age. All three formations contain moderately diverse and variably preserved palynomorphs. The palynomorph taxa, detailed systematically in the current Part One of this monograph, comprise 66 species of acritarchs and six of prasinophytes. Of these, two species of prasinophytes and 11 of acritarchs are newly established: Cymatiosphaera meandrica and Pterospermella franciniae; Aremoricanium hyalinum, A. solaris, Baltisphaeridium tenuicomatum, Gorgonisphaeridium crebrum, Lophosphaeridium aequalium, L. aspersum, Micrhystridium infrequens, Pylantios hadrus, Sertulidium amplexum, Striatotheca indistincta, and Tribulidium globosum. Pylantios (typified by P. hadrus), Sertulidium (typified by S. amplexum), and Tribulidium (typified by T globosum); are defined as new acritarch genera. Three new combinations are instituted: Baltisphaeridium pugiatum (PLAYFORD & MARTIN 1984), Polygonium canningianum (COMRAZ & PENIGUEL 1972), and Sacculidium furtivum (PLAYFORD & MARTIN 1984); and Ammonidium macilentum PLAYFORD & MARTIN 1984 and Sacculidium furtivum (PLAYFORD & MARTIN 1984) are emended. An appreciable number of palynomorph species are not formally named owing to lack of sufficient or adequately preserved specimens; others are compared but not positively identified with previously instituted species. The ensuing Part Two of this study will complete the systematic-descriptive documentation, i.e., chitinozoans, and evaluate the Canning Basin palynoflora in terms of its chronological and stratigraphic-correlative significance.

Indagine modellistica sulla tenuta idraulica delle difese arginali deputate alla protezione idraulica del territorio

La prima parte di questo lavoro di tesi tratta dell’interazione tra un bacino di laminazione e il sottostante acquifero: è in fase di progetto, infatti, la costruzione di una cassa di espansione sul torrente Baganza, a monte della città di Parma. L’obiettivo di tale intervento è di ridurre il rischio di esondazione immagazzinando temporaneamente, in un serbatoio artificiale, la parte più pericolosa del volume di piena che verrebbe rilasciata successivamente con portate che possono essere agevolmente contenute nel tratto cittadino del torrente. L’acquifero è stato preliminarmente indagato e monitorato permettendone la caratterizzazione litostratigrafica. La stratigrafia si può riassumere in una sequenza di strati ghiaioso-sabbiosi con successione di lenti d’argilla più o meno spesse e continue, distinguendo due acquiferi differenti (uno freatico ed uno confinato). Nel presente studio si fa riferimento al solo acquifero superficiale che è stato modellato numericamente, alle differenze finite, per mezzo del software MODFLOW_2005. L'obiettivo del presente lavoro è di rappresentare il sistema acquifero nelle condizioni attuali (in assenza di alcuna opera) e di progetto. La calibrazione è stata condotta in condizioni stazionarie utilizzando i livelli piezometrici raccolti nei punti d’osservazione durante la primavera del 2013. I valori di conducibilità idraulica sono stati stimati per mezzo di un approccio geostatistico Bayesiano. Il codice utilizzato per la stima è il bgaPEST, un software gratuito per la soluzione di problemi inversi fortemente parametrizzati, sviluppato sulla base dei protocolli del software PEST. La metodologia inversa stima il campo di conducibilità idraulica combinando osservazioni sullo stato del sistema (livelli piezometrici nel caso in esame) e informazioni a-priori sulla struttura dei parametri incogniti. La procedura inversa richiede il calcolo della sensitività di ciascuna osservazione a ciascuno dei parametri stimati; questa è stata valutata in maniera efficiente facendo ricorso ad una formulazione agli stati aggiunti del codice in avanti MODFLOW_2005_Adjoint. I risultati della metodologia sono coerenti con la natura alluvionale dell'acquifero indagato e con le informazioni raccolte nei punti di osservazione. Il modello calibrato può quindi essere utilizzato come supporto alla progettazione e gestione dell’opera di laminazione. La seconda parte di questa tesi tratta l'analisi delle sollecitazioni indotte dai percorsi di flusso preferenziali causati da fenomeni di piping all’interno dei rilevati arginali. Tali percorsi preferenziali possono essere dovuti alla presenza di gallerie scavate da animali selvatici. Questo studio è stato ispirato dal crollo del rilevato arginale del Fiume Secchia (Modena), che si è verificato in gennaio 2014 a seguito di un evento alluvionale, durante il quale il livello dell'acqua non ha mai raggiunto la sommità arginale. La commissione scientifica, la cui relazione finale fornisce i dati utilizzati per questo studio, ha attribuito, con molta probabilità, il crollo del rilevato alla presenza di tane di animali. Con lo scopo di analizzare il comportamento del rilevato in condizioni integre e in condizioni modificate dall'esistenza di un tunnel che attraversa il manufatto arginale, è stato realizzato un modello numerico 3D dell’argine mediante i noti software Femwater e Feflow. I modelli descrivono le infiltrazioni all'interno del rilevato considerando il terreno in entrambe le porzioni sature ed insature, adottando la tecnica agli elementi finiti. La tana è stata rappresentata da elementi con elevata permeabilità e porosità, i cui valori sono stati modificati al fine di valutare le diverse influenze sui flussi e sui contenuti idrici. Per valutare se le situazioni analizzate presentino o meno il verificarsi del fenomeno di erosione, sono stati calcolati i valori del fattore di sicurezza. Questo è stato valutato in differenti modi, tra cui quello recentemente proposto da Richards e Reddy (2014), che si riferisce al criterio di energia cinetica critica. In ultima analisi è stato utilizzato il modello di Bonelli (2007) per calcolare il tempo di erosione ed il tempo rimanente al collasso del rilevato.

The impact on agriculture of the drawdown of shallow watertables

The promoters of the large groundwater developments implemented in the 1970's paid little attention to the effects of pumping on soil moisture. A field study, conducted in 1979 in the Tern Area of the Shropshire Groundwater Scheme, revealed that significant quantities of the available moisture could be removed from the root zone of vegetation when drawdown of shallow watertables occurred. Arguments to this effect, supported by the field study evidence, were successfully presented at the Shropshire Groundwater Scheme public inquiry. The aim of this study has been to expand the work which was undertaken in connection with the Shropshire Groundwater Scheme, and to develop a method whereby the effects of groundwater pumping on vegetation can be assessed, and hence the impacts minimised. Two concepts, the critical height and the soil sensitivity depth, formulated during the initial work are at the core of the Environmental Impact Assessment method whose development is described. A programme of laboratory experiments on soil columns is described, as is the derivation of relationships for determining critical heights and field capacity moisture profiles. These relationships are subsequently employed in evaluating the effects of groundwater drawdown. In employing the environmental assessment technique, digitised maps of relevant features of the Tern Area are combined to produce composite maps delineating the extent of the areas which are potentially sensitive to groundwater drawdown. A series of crop yield/moisture loss functions are then employed to estimate the impact of simulated pumping events on the agricultural community of the Tern Area. Finally, guidelines, based on experience gained through evaluation of the Tern Area case study, are presented for use in the design of soil moisture monitoring systems and in the siting of boreholes. In addition recommendations are made for development of the EIA technique, and further research needs are identified.

Efficient Computation of Electromagnetic Waves in Hydrocarbon Exploration Using the Improved Numerical Mode Matching (NMM) Method

In this study, we developed and improved the numerical mode matching (NMM) method which has previously been shown to be a fast and robust semi-analytical solver to investigate the propagation of electromagnetic (EM) waves in an isotropic layered medium. The applicable models, such as cylindrical waveguide, optical fiber, and borehole with earth geological formation, are generally modeled as an axisymmetric structure which is an orthogonal-plano-cylindrically layered (OPCL) medium consisting of materials stratified planarly and layered concentrically in the orthogonal directions.

In this report, several important improvements have been made to extend applications of this efficient solver to the anisotropic OCPL medium. The formulas for anisotropic media with three different diagonal elements in the cylindrical coordinate system are deduced to expand its application to more general materials. The perfectly matched layer (PML) is incorporated along the radial direction as an absorbing boundary condition (ABC) to make the NMM method more accurate and efficient for wave diffusion problems in unbounded media and applicable to scattering problems with lossless media. We manipulate the weak form of Maxwell's equations and impose the correct boundary conditions at the cylindrical axis to solve the singularity problem which is ignored by all previous researchers. The spectral element method (SEM) is introduced to more efficiently compute the eigenmodes of higher accuracy with less unknowns, achieving a faster mode matching procedure between different horizontal layers. We also prove the relationship of the field between opposite mode indices for different types of excitations, which can reduce the computational time by half. The formulas for computing EM fields excited by an electric or magnetic dipole located at any position with an arbitrary orientation are deduced. And the excitation are generalized to line and surface current sources which can extend the application of NMM to the simulations of controlled source electromagnetic techniques. Numerical simulations have demonstrated the efficiency and accuracy of this method.

Finally, the improved numerical mode matching (NMM) method is introduced to efficiently compute the electromagnetic response of the induction tool from orthogonal transverse hydraulic fractures in open or cased boreholes in hydrocarbon exploration. The hydraulic fracture is modeled as a slim circular disk which is symmetric with respect to the borehole axis and filled with electrically conductive or magnetic proppant. The NMM solver is first validated by comparing the normalized secondary field with experimental measurements and a commercial software. Then we analyze quantitatively the induction response sensitivity of the fracture with different parameters, such as length, conductivity and permeability of the filled proppant, to evaluate the effectiveness of the induction logging tool for fracture detection and mapping. Casings with different thicknesses, conductivities and permeabilities are modeled together with the fractures in boreholes to investigate their effects for fracture detection. It reveals that the normalized secondary field will not be weakened at low frequencies, ensuring the induction tool is still applicable for fracture detection, though the attenuation of electromagnetic field through the casing is significant. A hybrid approach combining the NMM method and BCGS-FFT solver based integral equation has been proposed to efficiently simulate the open or cased borehole with tilted fractures which is a non-axisymmetric model.

(Table 3) Geochemistry of borehole fluids of ODP Hole 137-504B

Circulation of seawater through basaltic basement for several million years after crustal emplacement has been inferred from studies of surface heat flow, and may play a significant role in the exchange of elements between the oceanic crust and seawater. Without direct observation of the fluid chemistry, interpretations regarding the extent and timing of this exchange must be based on the integrated signal of alteration found in sampled basalts. Much interest has thus been expressed in obtaining and analyzing fluids directly from basaltic formations. It has been proposed that open oceanic boreholes can be used as oceanic groundwater wells to obtain fluids that are circulating within the formation. Water samples were collected from the open borehole in Hole 504B prior to drilling operations on Leg 137, with the original intention of collecting formation fluids from the surrounding basaltic rocks. Past results have yielded ambiguous conclusions as to the origin of the fluids recovered-specifically, whether or not the fluids were true formation fluids or merely the result of reaction of seawater in the borehole environment. The chemistry of eight borehole fluid samples collected during Leg 137 is discussed in this paper. Large changes in major, minor, and isotopic compositions relative to unaltered seawater were observed in the borehole fluids. Compositional changes increased with depth in the borehole. The samples exhibit the effect of simple mixing of seawater, throughout the borehole, with a single reacted fluid component. Analysis and interpretation of the results from Leg 137 in light of past results suggest that the chemical signals observed may originate predominantly from reaction with basaltic rubble residing at the bottom of the hole during the interim between drilling legs. Although this endeavor apparently did not recover formation waters, information on the nature of reaction between seawater and basalt at the prevalent temperatures in Hole 504B (>160°C) has been gained that can be related to reconstruction of the alteration history of the oceanic crust. Isotopic analyses allow calculation of element-specific water/rock mass ratios (Li and Sr) and are related to the extent of chemical exchange between the borehole fluids and basalt.

Analysis of Evidence from Core Sample Fractures in Stewardson Lake Uranium Prospect

This study was performed to characterize evidence of potential unconformity-type U mineralizing fluids in drill core fractures from the Stewardson Lake prospect, in the Athabasca Basin, located in Northern Saskatchewan and Alberta, Canada. Fractures were visually classified into eight varieties. This classification scheme was improved with the use of mineralogical characterization through SEM (Scanning Electron Microscope) and XRD analyses of the fracture fills and resulted in the identification of various oxides, hydroxides, sulfides, and clays or clay-sized minerals. Fractures were tallied to a total of ten categories with some commonalities in color. The oxidative, reductive or mixed nature of the fluids interacting with each fracture was determined based on its fill mineralogy. The measured Pb isotopic signature of samples was used to distinguish fractures affected solely by fluids emanating from a U mineralization source, from those affected by mixed fluids. Anomalies in U and U-pathfinder elements detected in fractures assisted with attributing them to the secondary dispersion halo of potential mineralization. Three types of fracture functions (chimney, composite and drain) were defined based on their interpreted flow vector and history. A secondary dispersion halo boundary with a zone of dominance of infiltrating fluids was suggested for two boreholes. The control of fill mineralogy on fracture color was investigated and the indicative and non-indicative colors and minerals, with respect to a secondary dispersion halo, were formally described. The fracture colors and fills indicative of proximity to the basement host of the potential mineralization were also identified. In addition, three zones of interest were delineated in the boreholes with respect to their geochemical dynamics and their relationship to the potential mineralization: a shallow barren overburden zone, a dispersion and alteration zone at intermediate depth, and a second deeper zone of dispersion and alteration.

Elfas und Ahlsburg - die Struktur der Schollengrenze zwischen Hils und Solling

The Elfas, Ahlsburg, and Salzderhelden overthrusts have been considered for a long time to be halotectonic structures caused by diapiric rise of Zechstein salt into thrust planes and other zones of weakness in the overlying rocks. Much of the Markoldendorf syncline between the overthrusts is covered by Quaternary deposits. The structure of the syncline and its western and southern boundaries were mapped with the help of many boreholes and micropalaeontological dating during remapping of the 1 : 25 000 geological sheet 4124 (Dassel). The Ahlsburg overthrust has now been shown to continue to the NW and to lead into the Lüthorst graben, which truncates the syncline in the west and passes into the crestal fault of the Elfas anticline. Figure 6 shows the NW-SE trending faults bounding three blocks which, according to seismic interpretation (Geotectonic Atlas 1996), have moved several times, in some areas in opposite senses. The NNE-SSW-trending grabens act as hinges, separating the NW-SE faults into sections that have undergone different movement.

Die Solling-Folge (Mittlerer Buntsandstein) im Grenzgebiet Niedersachsen - Thüringen - Hessen

This paper deals with the lithostratigraphic structure of the Solling sequence (Lower Triassic, Middle Buntsandstein) in the area between the Weser river region in the west and the Thuringian Eichsfeld region in the east. Lithologic profile mapping and the gamma-ray logs of several boreholes and 40 exposures have been used to define the lithostratigraphic Classification of the Solling sequence, to mark the facies zones and to find the connection between Sediments of the Thuringian basin in the east and the Weser fault trough via the crest of the Eichsfeld-Altmark Ridge. Tectonically controlled movements of synsedimentary character are the reason for the extreme convergence within the Solling sequence and the extreme Stratigraphie gap at its base (Hardegsen unconformity, Trusheim 1961) in the region of the swells. The discussion also demonstrates the importance of fault bundles active during Triassic and responsible for the thickness pattem of the Solling sequence between the Weser fault trough and the Eichsfeld-Altmark Ridge. The largest Stratigraphie gap is present at the line Brehme (Ohm Mountains) - Beuren - Treffurt where the Solling sequence covers Av/cn/a-bearing layers of the Volpriehausen sequence. In paiticular the Ridge sequences prove the existence of a further erosion unconformity within the Solling sequence (Solling unconformity, Kunz 1965) below the Thuringian Chirotheriensandstein as found by Rohling (1986) in the North German basin at the Stratigraphie level of the Karlshafen layers.

AMPLITUDE-VERSUS-ANGLE ANALYSIS AND WIDE-ANGLE-INVERSION OF CROSSWELL SEISMIC DATA IN A CARBONATE RESERVOIR

Crosswell data set contains a range of angles limited only by the geometry of the source and receiver configuration, the separation of the boreholes and the depth to the target. However, the wide angles reflections present in crosswell imaging result in amplitude-versus-angle (AVA) features not usually observed in surface data. These features include reflections from angles that are near critical and beyond critical for many of the interfaces; some of these reflections are visible only for a small range of angles, presumably near their critical angle. High-resolution crosswell seismic surveys were conducted over a Silurian (Niagaran) reef at two fields in northern Michigan, Springdale and Coldspring. The Springdale wells extended to much greater depths than the reef, and imaging was conducted from above and from beneath the reef. Combining the results from images obtained from above with those from beneath provides additional information, by exhibiting ranges of angles that are different for the two images, especially for reflectors at shallow depths, and second, by providing additional constraints on the solutions for Zoeppritz equations. Inversion of seismic data for impedance has become a standard part of the workflow for quantitative reservoir characterization. Inversion of crosswell data using either deterministic or geostatistical methods can lead to poor results with phase change beyond the critical angle, however, the simultaneous pre-stack inversion of partial angle stacks may be best conducted with restrictions to angles less than critical. Deterministic inversion is designed to yield only a single model of elastic properties (best-fit), while the geostatistical inversion produces multiple models (realizations) of elastic properties, lithology and reservoir properties. Geostatistical inversion produces results with far more detail than deterministic inversion. The magnitude of difference in details between both types of inversion becomes increasingly pronounced for thinner reservoirs, particularly those beyond the vertical resolution of the seismic. For any interface imaged from above and from beneath, the results AVA characters must result from identical contrasts in elastic properties in the two sets of images, albeit in reverse order. An inversion approach to handle both datasets simultaneously, at pre-critical angles, is demonstrated in this work. The main exploration problem for carbonate reefs is determining the porosity distribution. Images of elastic properties, obtained from deterministic and geostatistical simultaneous inversion of a high-resolution crosswell seismic survey were used to obtain the internal structure and reservoir properties (porosity) of Niagaran Michigan reef. The images obtained are the best of any Niagaran pinnacle reef to date.

Inversion of ambient seismic noise HVSR to evaluate velocity and structural models of the Lower Tagus Basin, Portugal

During its history, several significant earthquakes have shaken the Lower Tagus Valley (Portugal). These earthquakes were destructive; some strong earthquakes were produced by large ruptures in offshore structures located southwest of the Portuguese coastline, and other moderate earthquakes were produced by local faults. In recent years, several studies have successfully obtained strong-ground motion syntheses for the Lower Tagus Valley using the finite difference method. To confirm the velocity model of this sedimentary basin obtained from geophysical and geological data, we analysed the ambient seismic noise measurements by applying the horizontal to vertical spectral ratio (HVSR) method. This study reveals the dependence of the frequency and amplitude of the low-frequency (HVSR) peaks (0.2–2 Hz) on the sediment thickness. We have obtained the depth of the Cenozoic basement along a profile transversal to the basin by the inversion of these ratios, imposing constraints from seismic reflection, boreholes, seismic sounding and gravimetric and magnetic potentials. This technique enables us to improve the existing three-dimensional model of the Lower Tagus Valley structure. The improved model will be decisive for the improvement of strong motion predictions in the earthquake hazard analysis of this highly populated basin. The methodology discussed can be applied to any other sedimentary basin.