14.11.2014 Views

Geophysical data acquisition - OGS

Geophysical data acquisition - OGS

Geophysical data acquisition - OGS

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

Istituto Nazionale<br />

di Oceanografia e di Geofisica Sperimentale - <strong>OGS</strong><br />

GEOPHYSICS<br />

OF THE<br />

LITHOSPHERE<br />

DEPARTMENT<br />

Director: Giuliano BRANCOLINI<br />

2 0 0 0<br />

A N N U A L R E P O R T<br />

3


About <strong>OGS</strong><br />

The Istituto Nazionale di Oceanografia e di Geofisica Sperimentale -<br />

<strong>OGS</strong> (formerly Osservatorio Geofisico Sperimentale di Trieste)<br />

is a research institute financed by the Italian Ministry of Universities<br />

and Research. Its function is research in geology, geophysics and<br />

oceanography. More specific tasks are: crustal studies; the search for<br />

oil, gas and minerals; earthquake seismology; environmental<br />

geophysics; hydrogeology; hydrodynamics and ecology of the seas<br />

and oceans.<br />

These activities are carried out by the three Departments of<br />

Geophysics of the Lithosphere, Oceanography, and Seismology, which<br />

employ 56 researchers, 19 senior technicians and 63 technicians.<br />

Although established in 1949, the origin of <strong>OGS</strong> can be dated back<br />

to 1841, when the Osservatorio Meteorologico was founded at Trieste.<br />

The institution publishes the results of its studies, exploration, and<br />

investigations, and preserves for study and reference the geophysical<br />

<strong>data</strong> collected by the r/v <strong>OGS</strong>-Explora during the seven Antarctic<br />

campaigns since 1988, together with the historical archives of the<br />

Trieste seismographic station.<br />

<strong>OGS</strong> is concerned with transferring the results of its research<br />

activities to industry, and it is open to cooperation with scientists<br />

from academic and research institutions, as well as to partnership<br />

with industrial research centers.<br />

Istituto Nazionale<br />

di Oceanografia e di Geofisica Sperimentale - <strong>OGS</strong><br />

GEOPHYSICS OF THE LITHOSPHERE DEPARTMENT<br />

Borgo Grotta Gigante 42/C<br />

34010 Sgonico<br />

Trieste, Italy<br />

Tel: 0039-040 21401<br />

Fax: 0039-040 327521<br />

E-mail: gbrancolini@ogs.trieste.it<br />

Web: http://www.ogs.trieste.it<br />

4


Contents<br />

Introduction . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 1<br />

Measurements while drilling (<strong>data</strong> <strong>acquisition</strong>) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 5<br />

<strong>Geophysical</strong> <strong>data</strong> <strong>acquisition</strong> . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 10<br />

Seismic <strong>data</strong> processing . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 17<br />

Measurements while drilling (R & D) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 23<br />

Wave modeling . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 25<br />

• Seismic modeling for exploration geophysics . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 25<br />

• Ground penetrating radar for environmental problems . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 28<br />

• Simulation of the ground motion caused by earthquakes and site response analysis . . . . 28<br />

Seismic inversion . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 33<br />

• The Cat-3D tomographic software . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 33<br />

• Joint 3D inversion of P, S and converted waves . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 34<br />

• Time-lapse 3D tomography . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 34<br />

• Seismic tomography for environmental studies . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 36<br />

<strong>Geophysical</strong> interpretation . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 38<br />

• Joint Italian/Australian Marine Geoscience Expedition to the George V<br />

Land Region of East Antarctica (Wilkes Land Glacial History, WEGA project) . . . . . . . . . . . . 39<br />

• Physical properties and seismic stratigraphy of ODP Leg 178 well sites,<br />

Antarctic Peninsula Pacific margin . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 40<br />

• Orbitally-Controlled rhythmic sedimentation in the Wild Drift,<br />

Antarctica (ODP Leg 188, Site 1165) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 42<br />

• Subsidence at the Cape Roberts drill sites (Ross Sea, Antarctica)<br />

from backstripping techniques . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 42<br />

• Cenozoic Evolution of the South Orkney Microcontinent . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 45<br />

• Structure and Cenozoic evolution of the South America - Scotia plate<br />

boundary in the Tierra del Fuego region . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 46<br />

• Mapping the BSR on the South Shetland Margin (Antarctica)<br />

and assessing gas hydrate and free gas quantities . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 46<br />

• Gas hydrate physical properties imaging by multi-attribute analysis - Blake<br />

Ridge BSR Case History . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 48<br />

• Physical properties of sediment cores from the Antarctic continental margins . . . . . . . . . . 51<br />

• Backstripping modelling in the frame of the Stratigraphical Development<br />

of the Glaciated European Margin (STRATAGEM) - EU project . . . . . . . . . . . . . . . . . . . . . . . . . . . 52<br />

• Earth gravity field: measurements, <strong>data</strong> processing and interpretation . . . . . . . . . . . . . . . . . 53<br />

• Synthetic Aperture Radar (SAR) remote sensing . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 55<br />

Administration . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 56<br />

Publications . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 57<br />

Presentations at meetings and conventions . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 62<br />

Book reviews . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 66<br />

Educational video . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 66<br />

Visitors . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 67<br />

International seminars in solid earth geophysics . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 67<br />

5


Introduction<br />

The Department Geophysics of the Lithosphere (GdL) carries out studies and<br />

researches in a wide range of applied and theoretical geophysics. It is organised<br />

into four research groups, three operative groups and one support group. It<br />

employees 20 researchers, 16 senior technician and 27 technician. During 2000<br />

the Department gave hospitality to 7 grants.<br />

Main activity of the GdL is the application of geophysical methods to the<br />

knowledge of the underground. This knowledge has always been important in<br />

human history, for finding out minerals and fresh water, but only in the last 50<br />

years the need for accurate reconstructions of the underground geology drastically<br />

increased, due to the world wide utilisation of the hydrocarbon as a primary energy<br />

source.<br />

Only in the last years however, grew the consciousness that the underground is a<br />

masterpiece in the global environment and that its knowledge is of importance not<br />

only for the exploitation of natural resources, but also for a sustainable<br />

management of the global environment.<br />

GdL activities during the 2000 were based on a deep consciousness of the central<br />

role that the correct and sustainable managment of the underground resources<br />

plays on the human development and impaction on the natural systems. The main<br />

fields of our researches can be grouped in four main fields:<br />

1) to recognise the presence of natural resources, with particular attention to<br />

hydrocarbon,<br />

2) to reconstruct the dynamic of the natural systems through the study of the<br />

sediments,<br />

3) to evaluate the impact of human activities to the underground system,<br />

4) to produce high tecnology services for oil industry.<br />

In all these fields, GdL takes advantage of its long experience in multichannel<br />

seismic <strong>data</strong> <strong>acquisition</strong> and processing for oil industry.<br />

In the hydrocarbon detection, a major improvement was reached by the first truly<br />

3D Seisbit survey during the drilling of the Vallazza well. The Seisbit techniques is<br />

a trade mark <strong>OGS</strong>-Agip and it allows the reconstruction of a direct and reverse<br />

Vertical Seismic Profile during the drilling by listening to the noise produced by<br />

the drilling bit. The 3D multi-offset in the Vallazza well produced, in a-quasi-real<br />

time, an accurate and clear reconstruction of the geological strata below the<br />

drilling bit.<br />

The 2000 budget for natural resources’ study was 795,000 EURO, which derives<br />

from research contracts with Agip and European Community<br />

The study of the natural systems through the sediments was focussed in the high<br />

latitude environments (Artic and Antartic). High resolution seismic, correlated to<br />

sediment cores and drilling <strong>data</strong>, were used for detailed reconstruction of the last<br />

glacial-interglacial cycles and to infer type and characteristics and dynamics of the<br />

ice caps. Study have been carried out on two marine cruises, one in the Weddel<br />

1


Sea, close to the Antarctic Peninsula, in co-operation with the US Antarctic<br />

Program, the other off the Wilkes Land, in co-operation with the Australian<br />

Antarctic Program.<br />

The 200 budget for these studies was 450,000 EURO, mainly supported by the<br />

Italian Antarctic Program.<br />

The main problem we are facing in the application of the traditional onshore<br />

seismic to environmental and hydrogeological studies, is the need to drastically<br />

increase the resolution of the method: our approach follows three convergent<br />

directions: i) by a theoretical approach based on the study and modelling of the<br />

propagation of P and S waves in heterogenous-visco-elastic media; ii) studying the<br />

source, by the <strong>acquisition</strong> of an high frequencies (up to 400 Hz) vibroseis system;<br />

iii) the processing, by the application of evolved statics computation techniques.<br />

The 2000 budget for these researches was about 500,000 EURO, mainly supported<br />

by the Ministry for University and Research and by Fondo Trieste.<br />

During the 2000, the GdL was also involved in significant activities for the oil<br />

industry, particularly Agip. These were based on the Seisbit tecnology that was<br />

applied for monitoring the bit position and carring direct and reverse WSP while<br />

drilling.<br />

The total budget for these activities was 1,230,000 EURO.<br />

Giuliano BRANCOLINI<br />

Department Director<br />

2


PRESIDENCY<br />

IGINIO MARSON<br />

ADMINISTRATION<br />

AND<br />

SERVICESI<br />

DEPARTMENT<br />

OF<br />

GEOPHYSICS<br />

DEPARTMENT<br />

OF<br />

OCEANOGRAPHY<br />

DEPARTMENT<br />

OF<br />

SEISMOLOGY<br />

GIULIANO BRANCOLINI<br />

RENZO MOSETTI<br />

ALBERTO MICHELINI<br />

ADMINISTRATION<br />

OFFICE<br />

DARIO COLONNELLO<br />

WAVE<br />

MODELLING<br />

GÉZA SERIANI<br />

GEOPHYSICAL<br />

DATA ACQUISITION<br />

DANIEL NIETO YABAR<br />

SEISMIC<br />

INVERSION<br />

ALDO VESNAVER<br />

MEASUREMENTS<br />

WHILE DRILLING<br />

(DATA ACQUISITION)<br />

GIULIANO DORDOLO<br />

MEASUREMENTS<br />

WHILE DRILLING<br />

(R & D)<br />

FLAVIO POLETTO<br />

SIGNAL<br />

PROCESSING<br />

NIGEL WARDELL<br />

GEOPHYSICAL<br />

INTERPRETATION<br />

ANGELO CAMERLENGHI<br />

3


Measurements while drilling (<strong>data</strong> <strong>acquisition</strong>)<br />

Coordinator: Giuliano DORDOLO<br />

G. CAPELLI<br />

B. CATTANI<br />

P. COMELLI<br />

G. CRISTOFANO<br />

P. GHIDINI<br />

M. GIORGI<br />

B. MOIMAS<br />

V. PASCIULLO<br />

A. SCHLEIFER<br />

G. VASCOTTO<br />

F. ZGAUC<br />

In the year 2000 the ASTI group, keeping on the industrial application of the<br />

Seisbit ® System 1, begins the field tests and research applications of the new<br />

version of Seisbit ® System 2.<br />

Starting the application of the new system means a complete revision and a<br />

hardware and software redesign of the older version. It was necessary indeed to<br />

increase the operational capabilities like interfacing with other systems. Most of<br />

the hardware problems due to the synchronization between different systems has<br />

been solved developing new hardware specifically for our purpose. Using these new<br />

products with completely new software solutions, the new system is able to<br />

manage the timing problems merging <strong>acquisition</strong> <strong>data</strong> form different systems with<br />

different sample rates and different sampling characteristics.<br />

Redesigning the system makes also possible to manage a huge amount of different<br />

types of <strong>acquisition</strong> units such as Sercel SN348, SN368 and SN368E mixed almost<br />

without limits.<br />

The analogic section of the system has been modified and improved too, so in the<br />

<strong>data</strong> transmission, as well as in the safety and in the assembling features.<br />

Seisbit ® System 2 scheme.<br />

5


Field operations showed a normal course, so it has been possible to get all the<br />

targets:<br />

Cerro F.- service well: seismic while drilling <strong>acquisition</strong> with ENI-AGIP. 55 channel<br />

seismic survey. Main target top of cretaceous platform, deviation geometry control,<br />

refraction survey for static corrections, indications on well casing operations.<br />

Monte A.- service well: seismic while drilling <strong>acquisition</strong> with ENI AGIP. 60<br />

channel seismic survey. Main target top of carbonate platform, deviation geometry<br />

control, refraction survey for static corrections, indications on well casing<br />

operations, receivers pattern test.<br />

Cerro F.- research well: seismic while drilling <strong>acquisition</strong> with ENI-AGIP in the<br />

frame of a research project. Down-hole instrumentation testing.<br />

Geothermic service well: seismic while drilling <strong>acquisition</strong> service in Larderello<br />

area for ENEL. Well geometry tracing. 164 channel pseudo-3D seismic survey.<br />

Geothermic well<br />

recording plan.<br />

Rig: 12 pilot channels<br />

L1 :22 channels, L2: 26 channels, L3: 26 channels<br />

L4: 26 channels, L5: 26 channels, L6: 26 channels<br />

First channel offset from centre of well: 325 m,<br />

distance between two station units: 75 m.<br />

6


Vallazza research well: 3D-RVSP seismic survey. – CEE Research Project. This<br />

Project required a large amount of technical and human resource, to realize an<br />

experimental survey, with a 500 station units spread configuration. The main goal<br />

was the validation of the 3D Seismic While Drilling. The spread configuration is<br />

surely unusual: it consists in 2 circles centred on the well and 2 branches crossed<br />

on the well too. The units pattern follow a saw-teeth shape, each segment<br />

orthogonal to the radius of the circles, of 1000 and 2000 m length respectively. The<br />

operation scheduling, people, instrumentation, <strong>acquisition</strong>s and <strong>data</strong> management<br />

followed the expected flow without problems. There was some difficulty solving the<br />

<strong>data</strong> management problems, due to the huge amount of them: over 130 Gbyte.<br />

3D-RVSP recording<br />

pattern.<br />

7


After these research and testing phases, the new <strong>acquisition</strong> and pre-processing<br />

techniques are validated and have been applied in the service and industrial<br />

surveys.<br />

Cerro Fal. Service well: 54 channels seismic survey, main target top oil trap,<br />

refraction survey for static corrections, using stationary noise too, indications on<br />

well casing operations.<br />

Before <strong>data</strong> <strong>acquisition</strong>, it has to perform some preparatory activities, such as<br />

scouting, topography units location, positioning of the recording materials,<br />

mounting of the sensor pilots on the rig and the short refraction survey to have<br />

the static corrections.<br />

Elaboration and application of an accurate static correction is a very important<br />

step in the processing of VSP and VSP CDP mapping. Data correctly elaborated are<br />

useful in drilling too, i.e. to schedule the casing operations.<br />

Example of a while drilling seismic<br />

survey area.<br />

8


Therefore, having beside a high elevation spread and a datum plane quite deeper,<br />

a lot of effort has gone into studying and designing the survey, first of all using the<br />

already existing informations (seismic lines, up-hole survey, stacking chart,<br />

dromochrones and thematic maps) and then the foreseen stratigraphy. After the<br />

ray-tracing, this set of elements allowed to plan the <strong>data</strong> recording geometries,<br />

fitting the basis extension with high channel number (48), using the Summit<br />

telemetric system and a Hydrapulse as source.<br />

Moreover, together with this new approach using short refraction for a while<br />

drilling survey, <strong>data</strong> are dynamically correct with Seisbit ® information until the<br />

datum plane depth.<br />

The ASTI group worked in collaboration with more crews up to 24 people, most of<br />

them <strong>OGS</strong> senior technicians and engineers and some people from the contractors<br />

companies.<br />

After the testing phases of the Seisbit ® 2, innovation studying and designing are<br />

pursued, to improve the automation capabilities and the <strong>data</strong> storage and<br />

elaboration features in field.<br />

SHORT REFRACTION SURVEY RECORDING PLAN<br />

Spread configuration<br />

SEISBIT<br />

channel<br />

SEISBIT<br />

channel<br />

SEISBIT<br />

channel<br />

SEISBIT<br />

channel<br />

SEISBIT<br />

channel<br />

SEISBIT<br />

channel<br />

P.s.: shot point<br />

G1 - G48: geophones of the seismic short refraction profile<br />

Short refraction spread example.<br />

9


<strong>Geophysical</strong> <strong>data</strong> <strong>acquisition</strong><br />

Coordinator: Daniel NIETO YABAR<br />

S. BARBAGALLO<br />

L. BARADELLO<br />

R. BOLIS<br />

A. BRATUS<br />

G. COVA<br />

C. D’AMICANTONIO<br />

E. DEL NEGRO<br />

F. FANZUTTI<br />

M. GROSSI<br />

B. MARINO<br />

P. PAGANINI<br />

M. POROPAT<br />

G. VISNOVIC<br />

The wide and profitable group activities, supporting the other research groups, are<br />

listed below. Acquiring appliances have been implemented in all the operating<br />

fields, such as GPR, geoelectric, onshore and offshore seismic investigations.<br />

A methodological research project has been presented and approved at the Fondo<br />

Trieste, with the results of a further technological development of the <strong>acquisition</strong><br />

systems and an improvement of scientific technicians.<br />

Acquiring appliances are:<br />

SEISMIC ON-SHORE:<br />

Positioning<br />

• Total Station SOKKIA<br />

• DGPS<br />

• Infrared level LEIKA 3003<br />

Energization<br />

• MiniBang ISOTTA<br />

• Seismic Sources Power Weight Drop PWD-80<br />

• IDROBANG for well energization<br />

• Vibroseis-MINIVIB (to be completed)<br />

• Hammer<br />

Geophones<br />

• Single 100 Hz<br />

• Array of 6x20 Hz<br />

• Array of 12x10 Hz<br />

Recording<br />

• SUMMIT telemetric 140 channels<br />

• OYO DAS-1 96 channels<br />

SEISMIC OFFSHORE:<br />

Positioning<br />

• DGPS<br />

• Communication Technology Navigation System<br />

Energization<br />

• PULSAR 2000 Power Unit for Uniboom<br />

• IDROBANG for water energization<br />

• GI-Guns<br />

• Water-Gun<br />

• Sure Shot gun Controller<br />

• Bauer I 28.0 – 75, high-pressure air compressor<br />

10


• Standard 20’ container modified to host the compressor, the high-pressure air<br />

control and the airgun workshop, complete with all the required electric<br />

wiring.<br />

Streamer<br />

• 1200 m ITI solid state<br />

• Monochannel 3.5 m Geometrics<br />

• Single hydrophones and array hydrophones (10 single in array)<br />

• Streamer Tester A-2000<br />

• Digicourse streamer position control system<br />

• ARDEA streamer winch<br />

Recording<br />

• OYO DAS-1 recording system, expanded to 96 channels, with three IMB<br />

3480/3490 recording modules.<br />

• DELPH II System (version n° 2 seismic channels)<br />

DATASONICS - CHIRP-II (sub-bottom profiler)<br />

ERT (Earth Resistivity Tomography)<br />

• Resistivimeter Syscal-R2, IRIS Instruments<br />

• Output : 800V- 2,5A<br />

• Ground energization power supply: external<br />

- 250W using 12V input<br />

- 1200W using 220V input<br />

• Multinode System of IRIS, to control from 32 to 256 intelligent nodes<br />

GPR (Ground Penetrating Radar)<br />

• SIR-2000 GSSI<br />

• Antennas:<br />

• SUBECHO 35 MHz<br />

- SUBECHO 70 MHz<br />

- GSSI 100 MHz bistatic<br />

- GSSI 200 MHz<br />

MAGNETOMETRIC ACQUISITION:<br />

• Cesium Gradiometer G-858, Geometrics<br />

• Proton Magnetometer G-856, Geometrics<br />

The group’s activities are:<br />

• Project Tasman Sea. Geophysics research project in the Tasman Sea, with the<br />

R/V research ship Polar Duke, in cooperation with BGR of Hannover.<br />

Multichannel seismic lines and magnetometry were displayed.<br />

• High resolution seismic in Mica. The high-resolution 3D seismic program of<br />

the Mica Project has been acquired using the new Summit telemetric<br />

<strong>acquisition</strong> system and comprehends the sites of S. Pier d’Isonzo and Iamiano<br />

11


(GO). The energizing systems were an Isotta gun and PWD. In the Mica project<br />

127 transepts of 20-channels have been acquired, for a total amount of 994<br />

points, while in Iamiano the acquirement grid was provided by 600 channels<br />

and 530 energization points.<br />

• Integrated metodology in the Mica project. Combined to the high resolution<br />

seismic, the two sites of S. Pier d’Isonzo and Iamiano were characterised using<br />

GPR, ERT (earth resistivity tomography) and magnetometry. The 2D<br />

resistivity profiles had length from 31 to 315 m.<br />

• TRUCK Project. In February 2000, for the AGIRE srl society a magnetic survey<br />

was carried out on an area near PERPIGNAN (F). The aim of the survey was<br />

the determination of magnetic anomalies connected with the supposed<br />

presence of a buried vehicle. The area was delimited with a total station<br />

SOKIA. The vertexes of the area were joined with the local topography. The<br />

high-resolution magnetic <strong>acquisition</strong> was carried out with a portable cesium<br />

gradiometer GEOMETRICS mod. G-858 along parallel profiles with a 2-m<br />

range. Magnetic <strong>data</strong> were reduced by the field reference model IGRF 2000.<br />

• KRSKO High Resolution. In February 2000 three high-resolution seismic<br />

profiles, with a total length of 4 km were recorded in the surroundings of the<br />

KRSKO nuclear-power plant (SLO). This <strong>acquisition</strong> was performed in the<br />

framework of the EU-programme PHARE to complete <strong>data</strong> acquired in 1999.<br />

• Project LARSEN. The NSF project “Paleohistory of the Larsen Ice Shelf:<br />

Evidence from the Marine Record”, with the aim of reconstructing the<br />

history and collecting sedimentological, biostratigraphical, biological and<br />

oceanographical information on the environmental change of the sea floor in<br />

the area that, until five years ago was still covered by the Larsen ice platform<br />

(Eastern part of the Anctartic peninsula). Among the various applied surveys<br />

the group took part in <strong>data</strong> <strong>acquisition</strong> with multibeam echosounder, Side<br />

Scan Sonar and monochannel reflection seismic. The expedition has been<br />

carried out in May 2000, on the icebreaker N.B. Palmer (USA).<br />

• Extended Program <strong>Geophysical</strong> Research in the surroundings of the Krsko<br />

NPP. Three new seismic lines were made with a double purpose: to clarify the<br />

position of some faults zones and correlate the seismic <strong>data</strong> with logs of the<br />

drill hole Drnovo 1.<br />

• Project STRATAGEM. Research project with the R/V research ship Dana, in<br />

three different investigation areas. In one of these, the Faeroe-Shetland<br />

margin, the main aim was the construction of a mid –to late Cenozoic<br />

stratigraphic framework for the Faeroe-Shetland and the setting up an<br />

evolution model for the Faeroe-Shetland margin, with particular emphasis on<br />

the development of shelf-margin progradational wedges. During the fieldwork<br />

4 monochannel high-resolution reflection profiles have been recorded, for a<br />

total amount of 300 km. The <strong>acquisition</strong> system was composed by a GI-Gun<br />

(90 in 3 ) checked by a Real Time System Sure Shot, one streamer with an array<br />

of 10 Hydrophones 1.6 m space and one recording system Elics Delph-2x. The<br />

positioning was supplied by SHIPMATE GPS equipped with differential GPS<br />

corrector.<br />

12


• Offshore gravimetric survey “CERVETERI”. The gravimetric offshore program<br />

in the area covered by the IGM 1:50000 “CERVETERI” sheet was executed for<br />

the “Servizio Geologico Nazionale” with a floor gravimeter Lacoste &<br />

Romberg on the ship N/R VEGA owned by Sopromar. 300 points, with a 1-km<br />

interval between them, have been measured.<br />

• Integrated methodologies at the Doria Cave. This test was carried on to define<br />

the applicability of the ERT, GPR and magnetometry in detecting underground<br />

caves in Karst area.<br />

• High-resolution mono-channel seismic Acque Profonde. The investigated area<br />

comprehended sea/littoral, lagoon and continental/fluvial multiple<br />

environments in the Marano Lagoon. The survey was carried out by the<br />

employment of two high-resolution <strong>acquisition</strong> systems linked to a satellite<br />

positioning system. The first was composed by a Delph monochannel<br />

acquiring system synchronized with a Uniboom high frequency<br />

electrodynamical impulse source. The second by an integrated acquiringenergizing<br />

Chirp-<strong>data</strong>sonic system using a non-impulsive high frequency<br />

source, allowing centimetric resolutions.<br />

• High Resolution Survey in the Barcis Lake. Research project for the “Regione<br />

Friuli Venezia Giulia” with monochannel seismic. A Uniboom energizing<br />

source, a Delph acquiring system and a DGPS positioning system were<br />

mounted on a boat provided by the Barcis municipality.<br />

During this year Salvatore Barbagallo and Giorgio Cova have retired: we would<br />

like to thank them for their work.<br />

High resolution<br />

seismic in Mica. 3D<br />

<strong>data</strong> <strong>acquisition</strong> in<br />

the Iamiano site.<br />

13


High-resolution mono-channel<br />

seismic Acque Profonde.<br />

The off-shore workstation.<br />

Integrated methodologies at the Doria Cave.<br />

The ERT (Earth Resistivity Tomography) image.<br />

14


The GPR instruments.<br />

Integrated methodologies in the Mica project.<br />

Geoelectrical <strong>data</strong> <strong>acquisition</strong>.<br />

Extended Program<br />

<strong>Geophysical</strong> Research<br />

in the surroundings of<br />

the Krsko NPP.<br />

Acquisition near the<br />

nuclear power plant.<br />

15


Extended Program<br />

<strong>Geophysical</strong> Research<br />

in the surroundings<br />

of the Krsko NPP.<br />

Seismic line.<br />

Integrated<br />

methodologies at the<br />

Doria Cave. GPR<br />

(Ground Penetrating<br />

Radar) section.<br />

TRUCK project. High resolution<br />

gradiometric <strong>acquisition</strong>.<br />

Mono-channel<br />

seismic offshore.<br />

Data <strong>acquisition</strong>.<br />

16


Seismic <strong>data</strong> processing<br />

Coordinator: Nigel WARDELL<br />

G. CENTONZE<br />

L. CERNOBORI<br />

P. DIVIACCO<br />

M. MARCHI<br />

R. OLIVOTTI<br />

C. PELOS<br />

M. ROMANELLI<br />

R. SINCERI<br />

L. SORMANI<br />

F. ZGUR<br />

During the year 2000, the Processing group was involved in a variety of projects,<br />

from deep crustal studies to very high resolution land and marine surveys in both<br />

2-D and 3-D. However, the year was marred by the tragic loss of Licio Cernobori<br />

who died suddenly after a short and terrible illness. His experience, self-motivation<br />

and enthusiasm had made him a key member of the group. He will be greatly<br />

missed, as a colleague and also as a friend.<br />

One of the projects in which Licio was very much involved, was “<strong>Geophysical</strong><br />

research in the surroundings of the Krsko Nuclear Power plant”, to study the<br />

geological stability of the area around the nuclear power plant at Krsko in<br />

Slovenia. The three high resolution lines that had had to be postponed, in 1999,<br />

due to inclement weather conditions, were acquired in the spring. Although field<br />

tests had been conducted on the first line, during the winter, further testing was<br />

performed to take into consideration the different weather conditions for the three<br />

remaining lines. The group utilised the Vista processing package installed on a<br />

laptop computer to analyse these field tests. The <strong>data</strong> were subsequently finalised<br />

in the processing centre in Trieste. An elaborate processing sequence was used to<br />

attempt to overcome problems introduced by ploughed fields and shallow peat<br />

layers. These near surface conditions caused static problems and absorption of the<br />

high frequency component of the frequency spectrum.<br />

To complete the interpretation of the area and to better define the geological<br />

structures, three additional lines were planned in the vicinity of the power plant.<br />

These lines were semi-regional in character; the <strong>acquisition</strong> parameters were midway<br />

between those of the regional and the high resolution lines. Since the Isotta<br />

rifle, which was used for the high resolution, was not considered to have the<br />

penetration necessary for the target structures, a new source Hydrapulse, a<br />

powered weight drop, was used for these lines. This choice entailed more testing<br />

in the field to optimise the recording parameters for the new source and objectives.<br />

The processing group provided the technical expertise during these tests. This was<br />

followed by quality control (QC) and processing in the field to produce preliminary<br />

stacks. Since these surveys were very close to the nuclear power plant, 50 Hz noise<br />

from the power cables emanating from the power station, was very evident on the<br />

field records. The QC and processing in the field was important to monitor that the<br />

level of this noise was not saturating the signal and that it could be attenuated in<br />

the processing phase.<br />

The final processing of these <strong>data</strong>, which was performed at the processing centre<br />

in Trieste, was tailored to detail the particular objectives and to match the regional<br />

lines recorded the previous year. To assist the interpretation, the final stack<br />

sections of both <strong>data</strong>sets were migrated and converted to depth. All the processing<br />

17


sections were included, in digital form on a CD-ROM, with the interpreted results<br />

and conclusions in the final report.<br />

Another important project in which the group was involved was the CROP or Deep<br />

Crustal project. During the year, the processing of first part the seismic line CROP-<br />

11 (From Lazio to Abruzzo) was finalised. A non-standard processing sequence<br />

involving noise reduction, coherency enhancement, array simulation and<br />

refraction statics had been derived to improve continuity and enhance the signal<br />

to noise ratio, especially at deeper depths. This sequence is also being applied to<br />

the second part of this line (from Abruzzo to the Adriatic Coast) which the<br />

processing group had been involved in the field QC and processing in 1999.<br />

The group was also involved in QC and field processing in the marine environment<br />

in the WEGA project, offshore Wilkes Land in Antarctica. As part of this project to<br />

study the glacial history, analysts from the processing group participated in this<br />

cruise. In the spring of 2000, 1500 kms of multi-channel seismic were acquired by<br />

an Australian research vessel. Following initial processing on board, the <strong>data</strong> were<br />

transported back to Trieste for further processing and finalisation. The final <strong>data</strong><br />

were presented at an International Conference in Tasmania in December.<br />

The research activities of the group continued with the EU project Very high<br />

resolution marine 3D seismic method for detailed site investigation (VHR3D)<br />

which is aimed at a cost-effective detailed 3-D reconnaissance of the seabed<br />

sediment properties for geological, geotechnical and environmental site<br />

investigation purposes. The group plays a major role in this three-year project<br />

coordinating <strong>OGS</strong>’s contributions in processing, tomography, geotechnical<br />

studies, and seismic modelling. In this third year of the project, limitations in the<br />

<strong>acquisition</strong> technology, restricting the available navigation information and the<br />

range of offsets, reduced the involvement of the other <strong>OGS</strong> contributors. However,<br />

additional research had to be undertaken by the processing group to include preprocessing<br />

corrections in 3-D for wave-motion and tidal effects.<br />

Initial work had already been completed by the group, defining a methodology in<br />

2-D that used static corrections to correct for wave motion. This had been<br />

extended in the second year to include geometry regularisation to take into<br />

account variations in cable positions, and hence offset, which were not able to be<br />

recorded by the navigation system. This methodology was presented at the EAGE<br />

meeting in Glasgow in June 2000. The 3-D corrections were based on a similar<br />

statistical analysis of the first break arrivals after common offset spatial averaging<br />

that had been used in 2-D. In the 3-D case, the spatial averaging had to be<br />

performed in an areal sense to include not only the component of the wave motion<br />

but also the component due to tidal differences between lines. First results on a<br />

<strong>data</strong>set recorded in the Dover strait were encouraging.<br />

Members of the group participated in the next VHR3D cruise in St. Austell Bay in<br />

Cornwall where the two different multi-cable <strong>acquisition</strong> systems, proposed in the<br />

project, were used. After seeing these <strong>acquisition</strong> systems in operation a number<br />

of modifications were made to improve the technique. Since the cables were seen<br />

to act independently at times, a separate component was derived for each cable<br />

rather than for each shot. Also, different 3-D spatial filters had to be introduced to<br />

18


allow for the different spatial sampling in the in-line and cross-line directions.<br />

Results from the St. Austell Bay survey showed that the methodology worked well.<br />

3-D stack cubes using a one metre bin size were produced on <strong>data</strong>sets from both<br />

the <strong>acquisition</strong> systems. The complex channeling in the area was well delineated<br />

with a vertical resolution of less than half a metre. A presentation on this 3D<br />

methodology has been accepted for the 63rd Annual Meeting of the EAGE in<br />

Amsterdam in June 2001.<br />

The group was also involved in a research project, funded by ENI (Agip Division)<br />

through the University of Parma to study crustal movements and the processes of<br />

convergence and sedimentation in the Eastern Mediterranean. Within the project,<br />

1000 kms of seismic <strong>data</strong>, recorded by <strong>OGS</strong> in the seventies, was reprocessed and<br />

interpreted in order to design scaled physical models for sediment deformation<br />

experiments in the laboratory. This physical modelling was performed by a team<br />

from the University of Parma (Department of Earth Sciences) and IGN-CNR<br />

Bologna. The reprocessed <strong>data</strong> crossed the eastern Mediterranean ridge<br />

accretionary complex and extended to the Herodotus foredeep and African foreland<br />

(Nile river deep sea fan). The reprocessing involved re-generation of stack sections<br />

from the original tapes, migration and depth conversion. Pre-stack depth<br />

migration was also performed on selected parts of the profiles with Geodepth ©<br />

software. The migration velocity focusing analysis and grid tomography in<br />

Geodepth © was also used to study the velocity distribution and define possible<br />

facies changes in the salt sequence.<br />

The development of the group’s capabilities in 3-D processing continued with the<br />

MICA project. The aim of this project, funded by the ‘Fondo di Trieste’, was to<br />

define the flow of underground water in the zone of the Carso. Two high resolution<br />

land 3-D surveys were planned in the project; one was acquired in the early part of<br />

the year, whilst the second towards the end of the year. The processing group was<br />

involved both in the preparatory phase and in the processing of the recorded <strong>data</strong>.<br />

Notwithstanding the problem of lack of high frequencies in <strong>data</strong> recorded on land,<br />

the 3-D stack cube of the first survey produced good results in the area of interest<br />

after the application of residual static correction routines.<br />

The group has continued its collaboration with OCSA (Orellana Consultores S.A.<br />

Madrid), processing a number of high resolution land lines to plan the best routing<br />

for rail and road tunnels under mountainous regions. These <strong>data</strong> tend to be<br />

inherently noisy and lacking in continuous reflections due to the deformed, often<br />

metamorphic, areas in which they are recorded. Careful noise reduction and<br />

continuity enhancement routines have been used in the final sections. Continued<br />

close ties with other academic institutions produced a number of small processing<br />

projects. Two 3-D surveys were processed with the University of Trieste (DINMA),<br />

one line in the Acqua Profonda project for the local Friuli Venezia Giulia region<br />

(also with DINMA), and a line acquired in Sardegna by the University of Cagliari.<br />

The recovery and archiving of old multi-channel seismic <strong>data</strong> has continued to be<br />

an important part of the group’s activities. All the Mediterranean <strong>data</strong> recorded by<br />

<strong>OGS</strong> in the seventies, and about 70% of the <strong>data</strong> recorded by <strong>OGS</strong> in the Antarctic<br />

in recent years, have been transcribed from their original field format to SEG-Y on<br />

19


3480/90 cartridges. The group also performed a transcription service for<br />

ENEL/ERGA, copying their entire seismic field <strong>data</strong> library onto cartridges and<br />

CD-ROM.<br />

Processed stack <strong>data</strong> from the Antarctic is also being transcribed onto CD-ROM as<br />

part of the Seismic Data Library System (SDLS) project. The Antarctic Seismic<br />

Data Library System (SDLS) provides open access to multi-channel seismic<br />

reflection <strong>data</strong> collected by all countries in the Antarctica, to facilitate large-scale<br />

cooperative research projects. The SDLS has 11 library branches that are located<br />

in 10 countries world-wide. Researchers may go to any library branch to inspect<br />

Antarctic multi-channel seismic reflection <strong>data</strong>. The processing group performs<br />

any necessary reformatting, filtering and scaling to the stack <strong>data</strong> before<br />

transcribing them, in SEG-Y format, onto CD-ROM together with a visualisation<br />

program.<br />

VHR3D - Common offset surfaces (Dover)<br />

Results of the 3D<br />

methodology for<br />

determining static<br />

corrections for tide and<br />

wave motion effects.<br />

A common offset<br />

surface of the water<br />

bottom is shown<br />

without<br />

corrections (top),<br />

with tidal corrections<br />

(middle) and with tidal<br />

and shot corrections<br />

(bottom).<br />

20


Results of application<br />

of the static<br />

corrections to a<br />

stacked 3D <strong>data</strong><br />

cube. The original<br />

<strong>data</strong> is on the left<br />

and the static<br />

corrected cube on<br />

the right.<br />

VHR3D - Dover 3D cube<br />

St. Austell Bay 3D Cube<br />

Time slice at 33 ms.<br />

Time slices from the<br />

processed<br />

3D cube from<br />

St. Austell Bay<br />

after application of the<br />

static corrections for<br />

tide and wave motion.<br />

The high resolution<br />

definition of the<br />

channels is clearly<br />

evident.<br />

WAVE MODELING<br />

St. Austell Bay 3D Cube<br />

Time slice at 35 ms.<br />

21


S<br />

MS-53 SP 1175-1475 Velocity Analysis “semblance”<br />

N<br />

An example of the migration velocity focusing analysis<br />

using the GeoDepth Pre-stack depth imaging package.<br />

A line extracted<br />

from the Mica<br />

project 3D cube<br />

showing the<br />

improvement<br />

obtained by the<br />

application of 3D<br />

residual statics<br />

(bottom) over the<br />

original stack<br />

<strong>data</strong> (top).<br />

22


Measurements while drilling (R & D)<br />

Coordinator: Flavio POLETTO<br />

C. BELLEZZA<br />

P. CORUBOLO<br />

A. CRAGLIETTO<br />

M. LOVO<br />

M. MALUSA<br />

L. PETRONIO<br />

G. PINNA<br />

U. TINIVELLA<br />

S. TINONIN<br />

The research was primarily aimed at studying in deep and extending the applicability<br />

of the seismic while drilling technology as well as increasing the industrial<br />

potential by raising the number of geological information given by the method.<br />

The projects developed by the SERE group in collaboration with ENI/AGIP are the<br />

following:<br />

A) Acquisition of a 3D reverse VSP using the drill-bit source during the drilling of<br />

an ENI/AGIP well in Sicily (Italy).<br />

Such a method, which is the only capable of acquiring an onshore 3D VSP on<br />

a large number of <strong>acquisition</strong> levels, was applied in cooperation with the “while<br />

drilling” <strong>data</strong> <strong>acquisition</strong> group (ASTI).<br />

The experiment, funded with the contribution of the 3D-RVSP European Union<br />

programme (Contract Thermie OG 278/98 IT/UK, partners <strong>OGS</strong>, ENI/AGIP and<br />

Prosol Technology), allowed us to acquire a good quality 3D Seisbit <strong>data</strong>set.<br />

The receivers were placed on a 15 square Km area and the investigation was<br />

carried out on a 3 km long well section.<br />

We used radial and circular seismic lines in a geometry suited to discriminate<br />

the signal and noise arrivals (first Figure). During the survey setup we<br />

computed some 3D elastic models and after the <strong>acquisition</strong> phase <strong>data</strong> have<br />

been processed as multioffset vsp sections along the cross shaped radial seismic<br />

lines. We are still processing the circular lines with the aim of obtainig a 3D<br />

reverse vsp imaging.<br />

B) Theoretical analysis and experimental study of signals measured in the drill<br />

pipes and of reflection coefficients at the drill-bit / rock interface. This allowed<br />

us to measure the acoustic impedance of the drilled formations and to gather<br />

information about the signal amplitude in relation to changes in the drillstring<br />

and drilled rock (formation evaluation while drilling).<br />

C) Feasibility study for the Geosteering project (ENI/AGIP). This project is aimed<br />

at steering the drilling operations basing on “while drilling” information.<br />

Particularly, the research is aimed at the tuning of a method and an acoustic<br />

system to be used downhole for monitoring the lithology in the vicinity of the<br />

drill bit for the purposes of drillers and geopysicysts.<br />

The feasibility study includes the use of downhole instruments and the<br />

syncronization with surface measurements. Field and laboratory tests have<br />

been carried out with the prototype instrument called Instrumented Sub<br />

ENI/AGIP. Such measurements allowed us to detect the drill bit signal also in<br />

highly unfavorable conditions, even for the surface measurements (second<br />

Figure). Numerical modelling methods were developed to analyze the signals<br />

propagation throughout the drill string. Moreover, the research was extended<br />

23


to study methods suited to predict overpressure zones using SWD and to the<br />

analysis of electromagnetic signals in a well.<br />

D) Feasibility study and preparation of the technology for the extension of the<br />

method to the deep sea and downhole environments, in collaboration with<br />

ENI/AGIP and Tecnomare.<br />

E) In deep study and analysis of <strong>data</strong>sets collected during a SWD experiment<br />

carried out in a tunnel. Study of the applicative contexts, of the tunnel boring<br />

machine signal resolution and of the signal and noise components.<br />

F) Study and multioffset processing of a <strong>data</strong>set acquired in a geothermal area<br />

along seismic lines laid out in a multiradial geometry with respect to the well<br />

head (in collaboration with ENEL). Analysis and preparation of a <strong>data</strong>set for<br />

tomographic inversion, in collaboration with REDS.<br />

24


Wave modeling<br />

Coordinator: Géza SERIANI<br />

Seismic modeling for exploration geophysics<br />

J.M. CARCIONE<br />

F. CAVALLINI<br />

G. SERIANI<br />

The numerical modeling of seismic waves plays a key role in exploration<br />

geophysics, reservoir engineering, and environmental studies. In the year 2000,<br />

<strong>OGS</strong> has continued this kind of activity under the aegis of the European<br />

Community (research program “Detection of Overpressure Zones from Seismic<br />

and Well Data - ODS”), with emphasis on the following topics.<br />

The first one is the estimation of gas-hydrate concentration and free-gas<br />

saturation. When no direct measurements are available, a detailed knowledge of<br />

the compressional and shear velocity is essential for the quantitative estimation of<br />

gas hydrate and free gas in bottom-simulating reflectors (BSR). Discrepancies<br />

between experimental velocity profiles and the velocity for water-filled sediments<br />

reveal the presence of gas hydrate (positive anomaly) or free gas (negative<br />

anomaly). A three-phase Biot-type theory yields wave velocities of sediments<br />

saturated with water and gas hydrate, while the Hill average is used to model the<br />

patchy free-gas saturation below the BSR. The model has been applied to field <strong>data</strong><br />

acquired by the r/v <strong>OGS</strong> Explora in Antarctica.<br />

Velocity field across a section parallel<br />

to the South-Shetland Margin in Antarctica.<br />

The BSR is evident at the left where<br />

a low-velocity layer, caused by the presence<br />

of free gas, is embedded in a higher<br />

velocity background.<br />

A line extracted<br />

from the Mica<br />

project 3D cube<br />

showing the<br />

improvement<br />

obtained by the<br />

application of 3D<br />

residual statics<br />

(bottom) over the<br />

original stack<br />

<strong>data</strong> (top).<br />

Concentration map of gas hydrate<br />

(positive values) and free gas (negative values)<br />

corresponding to the BSR offshore the<br />

South-Shetland Islands. The figure shows<br />

the hydrate concentration (and free gas<br />

saturation) multiplied by the porosity<br />

(i.e., the volume concentration). In this way,<br />

we can compare the content of hydrate and<br />

free gas between zones of different porosity.<br />

25


Poisson’s ratio as an indicator of overpressure has been investigated. Poisson’s<br />

ratio values of dry samples are significantly smaller than those of fluid-saturated<br />

samples. The values are anomalously high for high pore pressure, with the<br />

possibility of differentiating between gas-saturated, brine-saturated and oilsaturated<br />

porous rocks. Two overpressure models, based on oil/gas conversion and<br />

disequilibrium compaction, have been developed to obtain Poisson’s ratio versus<br />

differential pressure. Poisson’s ratio is approximately constant at high differential<br />

pressures and increases (decreases) for saturated (dry) rocks at low differential<br />

pressures. Fluid type can be determined at all differential pressures from Poisson’s<br />

ratio. Moreover, the analysis is extended to the transversely isotropic case by<br />

computing the three Poisson’s ratios. Experiments performed on cores, under<br />

different pressure conditions, and calibration of the models with these <strong>data</strong>,<br />

provide a tool for inverting pore pressure from seismic <strong>data</strong>.<br />

1<br />

0.2<br />

0.15<br />

0.1<br />

gas<br />

water<br />

oil<br />

(a)<br />

2<br />

0.3<br />

0.25<br />

0.2<br />

0.15<br />

0.1<br />

gas<br />

water<br />

oil<br />

(b)<br />

3<br />

0.25<br />

0.2<br />

0.15<br />

0.1<br />

gas<br />

water<br />

oil<br />

(c)<br />

0.05<br />

0 20 40 60 80 100<br />

Effective pressure (MPa)<br />

0.05<br />

0 20 40 60 80 100<br />

Effective pressure (MPa)<br />

0.05<br />

0 20 40 60 80 100<br />

Effective pressure (MPa)<br />

Anisotropic Poisson’s ratioes σ 1<br />

(a), σ 2<br />

(b) and σ 3<br />

(c), for brine-, oil- and gas-saturated Berea<br />

sandstone versus effective pressure, compared with experimental dry-rock Poisson’s ratios.<br />

The dashed line is the best-fit curve to the dry-rock <strong>data</strong>.<br />

well 1B<br />

NW<br />

SE<br />

formations<br />

vel. P (km/s)<br />

Snapshot of the<br />

seismic wave field<br />

propagating in a<br />

complex geological<br />

structure during a 3D<br />

simulation of a<br />

seismic-while drilling<br />

experiment.<br />

depth (km)<br />

distance (km)<br />

26


Phase velocities<br />

of the five wave<br />

modes<br />

propagating in<br />

partially frozen<br />

Berea sandstone<br />

versus water<br />

proportion.<br />

Phase velocity (km/s)<br />

5<br />

4<br />

3<br />

2<br />

1<br />

0<br />

P1<br />

S1<br />

P2<br />

S2<br />

P3<br />

(a)<br />

0.05 0.1 0.15 0.2<br />

Water proportion<br />

Snapshot of the rock-frame<br />

vertical particle velocity<br />

corresponding to the wave<br />

modes illustrated<br />

in the previous figure.<br />

The compressional waves<br />

are labeled P1, P2 and P3,<br />

and the shear waves<br />

are labeled S1 and S2.<br />

Finally investigation on dynamics of frozen porous media has been conducted. The<br />

knowledge of the physical properties of frozen soils is essential, in polar areas, for<br />

the exploitation of mineral resources and for the construction of highways and<br />

pipelines. A three-phase Biot-type theory has been developed to describe a seismic<br />

wave propagating in a porous rock filled with ice and water. The model predicts<br />

three compressional waves and two shear waves, and takes into account energy<br />

dissipation and wave dispersion as observed in rocks. Attenuation is introduced<br />

with exponential relaxation functions, which allow a differential fomulation based<br />

on memory variables. The wavefield is computed with a grid method based on the<br />

Fourier differential operator in space and a Runge-Kutta time-integration<br />

algorithm. The presence of slow quasi static modes makes the differential<br />

equations stiff and hence numerically untractable as such. But a splitting timeintegration<br />

algorithm has allowed to solve the stiff part analytically. The algorithm<br />

is second-order accurate with respect to time and has spectral accuracy in the<br />

computation of the spatial derivatives.<br />

27


Ground penetrating radar for environmental problems<br />

J.M. CARCIONE<br />

F. CAVALLINI<br />

G. PADOAN<br />

G. SERIANI<br />

The ground-penetrating radar (GPR) has unique capabilities as a subsurface<br />

exploration device; as such, its applications (since 1992) include archaeology,<br />

mining, glaciology, hydrogeology, and environmental remediation. But its<br />

performance is highly sensitive on the knowledge of the physics of the medium and<br />

of the instrument itself. Therefore, basic research in this field is of crucial<br />

importance for pracical applications as well. Many concepts and techniques that<br />

were devised for seismic waves can be applied to the processing of GPR <strong>data</strong>, in<br />

virtue of the mathematical analogy between electromagnetic and elastodynamic<br />

equations. This research, funded by the Regional Government of Friuli - Venezia<br />

Giulia (northeastern Italy), aims at exploiting these methods for the monitoring of<br />

a polluted area.<br />

The exploding-reflector method, originally developed for seismic waves, has been<br />

adapted to produce zero-offset synthetic radargrams without need to compute<br />

common-shot records. The basic underlying idea consists in assuming that each<br />

reflecting point acts at the initial time as an instantaneous source of motion with<br />

a magnitude proportional to the normal-incidence reflection coefficient. The<br />

magnetic permeability is used as a free parameter to obtain a constant-impedance<br />

model and, so, to avoid multiple reflections. Moreover, the condition that the phase<br />

velocity remain unchanged also requires the scaling of the permittivity and of the<br />

conductivity. Thus, the method generates normal-incidence reflections, i.e., those<br />

having identical downgoing and upgoing wave paths. In an example with<br />

transverse-magnetic equations, this method has shown more accuracy, but less<br />

efficiency, than the plane-wave technique.<br />

Analytical solution for lossy anisotropic media has been developed. Microstructural<br />

features, fine layering, and fluid-filled cracks give rise to electric and magnetic<br />

anisotropy, while mineralized water in the fractures makes conductivity<br />

anisotropic. Water relaxation and ferromagnetism are responsible of energy<br />

dissipation. This motivates an orthotropic constitutive law in which attenuation<br />

stems from the Debye model. The corresponding solution is obtained in the<br />

frequency domain from a change of coordinates that turns Maxwell’s equations<br />

into Helmholtz equations, which are solved analytically. Then, the solution in time<br />

domain is recovered by numerically computing the inverse Fourier transform. The<br />

results are in agreement with a plane-wave analysis of the slowness, the<br />

attenuation and the energy velocity.<br />

Simulation of the ground motion caused by earthquakes<br />

and site response analysis<br />

E. PRIOLO<br />

G. LAURENZANO<br />

Synthetic seismograms computed as a solution of the full-wave propagation<br />

through a realistic geological structure can reproduce accurately much of the<br />

effects of medium heterogeneity and local soil conditions on the ground motion<br />

caused by earthquakes. They are a powerful tool for building seismic hazard<br />

scenarios as well as for performing site response analyses. For this kind of studies,<br />

we currently use two methods, i.e. the 2-D Chebyshev spectral element method<br />

28


(SPEM) and the Wavenumber Integration Method (WIM), which have very<br />

different features.<br />

The 2-D Chebyshev spectral element method (SPEM) is a high-order finite element<br />

technique, which solves the variational formulation of the equation. It has been<br />

entirely developed at <strong>OGS</strong> by E. Priolo and G. Seriani. The use of irregular meshes,<br />

as well as the high computational accuracy, which derives from the use of highorder<br />

Chebyshev polynomials, make the SPEM particularly suitable to solve<br />

numerically the seismic wave propagation through complex geological structures.<br />

The WIM solves the 3-D elastic full-wave equation in a plane layer medium. It has<br />

been developed by R. Herrmann at St. Louis University. Synthetic seismograms<br />

contain all wavefield phases, in both the near- and the far-fields. Furthermore,<br />

simulations can be set up very simply, since the method needs only few parameters to<br />

be defined. These properties make WIM suitable for performing predictions at a regional<br />

scale, as well as for kinematic modelling of the rupture process along the fault.<br />

The methodological development has been focused on improving the capability of<br />

the SPEM for building complex models. In the standard implementation, the<br />

spectral element method decomposes the whole model into a patch of sub-regions<br />

having constant physical properties. This approach is inadequate to represent geomodels,<br />

since they usually feature small scale heterogeneities and medium<br />

properties that change continuously in the space. This fact motivated the<br />

development of a different approach, in which 1) the geological structure is builtup<br />

through an interpolator which handles discontinuities, and using a minimum<br />

number of control points and lines, and 2) the mesh size adapts continuously to<br />

the medium properties and is controlled by few geometric constraints.<br />

In general, the activity has been developed mainly within a number of national<br />

research projects funded or co-ordinated by the GNDT1 and GNV2. During year<br />

2000, one three-years project – the Catania Project – finished its activity, while five<br />

new research projects started. In particular, the latter include a one-year project –<br />

the Marche Microzonation Project – and four three-years projects, namely:<br />

1) Development and Comparison between Methodologies for the Evaluation of<br />

Seismic Hazard in Seismogenic Areas: Application to the Central and Southern<br />

Apennines; 2) Damage Scenarios in the Veneto-Friuli Area; 3) Detailed Scenarios<br />

and Actions for Seismic Prevention of Damage in the Urban Area of Catania; and<br />

4) Integrated Seismic Methods Applied to the Investigation of the Active Volcano<br />

Structure. An Application to the Phlegrean Fields Caldera. We will shortly present<br />

some of those projects in the following.<br />

The Catania Project aimed at evaluating the seismic risk of a highly urbanised area<br />

typical of the Mediterranean region. The <strong>OGS</strong> contribution of the last year was<br />

twofold. In a first part, we simulated a recent event, the M 5.8 earthquake, which<br />

struck Eastern Sicily on December 13, 1990, and was recorded by the Catania<br />

ENEA-ENEL accelerometric station. Seismograms computed using the SPEM<br />

agree very well with those recorded by the accelerometric station. On the contrary,<br />

any plane layer representation, which simplifies the complex structure used by<br />

SPEM, does not provide a comparable agreement. In this case, the synthetic<br />

seismograms are computed by the WIM. This study has several outcomes. Firstly,<br />

29


it demonstrates that the whole approach, based on the SPEM, previously used to<br />

simulate the ground shaking for a destructive scenario earthquake, provides<br />

reliable results. Secondly, it shows that the model used to represent the crustal<br />

structure beneath this area is realistic. Indeed, simplified models may not be<br />

adequate for predicting the site response. Moreover, the high amplitude displayed<br />

by the Catania station during the 1990 earthquake can be explained as a combined<br />

effect of site and structure-path.<br />

The second contribution to the project was the analysis of the environmental<br />

seismic noise <strong>data</strong> (microtremors) acquired within the Catania municipal area in<br />

May 1999, with the aim of improving the prediction of the seismic ground motion<br />

locally. To this end, we followed Nakamura’s approach which, as proven, provides<br />

the main features of the dynamic ground response through the calculation of the<br />

spectral ratio between the horizontal and the vertical components (i. e., H/V ratio)<br />

of background microtremors. We found (first Figure) that several sites exhibit<br />

irrelevant or weak amplification, i. e., sites which are located either on lava (for<br />

example, the sites in the Catania centre or in the northern part of the municipal<br />

area) or well-consolidated sedimentary<br />

soils (Western districts of the city).<br />

The only sites which bear evidence of<br />

some amplification are located on<br />

either fillings soils lying over lava, or<br />

on the fine alluvial deposits of the<br />

Catania Plain.<br />

Catania Project. Map of the Catania<br />

municipal area showing the pseudoamplification<br />

estimated from background<br />

microtremor horizontal to vertical spectral<br />

ratio (HVSR).<br />

The color and size of the circles indicate<br />

peak frequency and amplitude of the<br />

HVSR, respectively. Blue and white circles<br />

indicate sites with very low pseudoamplification<br />

(< 2) and flat response in<br />

the low frequency range, respectively.<br />

Contour lines indicate the depth of the<br />

top of the light-blue clays formation,<br />

which has been assumed as a possible<br />

seismic bedrock. The background map<br />

shows the simplified geotechnical<br />

zonation.<br />

30


The Marche Microzonation Project is a two-year research project, whose general<br />

objective is the detailed microzonation of some selected cities of the Marche<br />

Region. The <strong>OGS</strong> contribution develops within the second year of the project. The<br />

aim is to perform a detailed simulation of the ground motion and site response<br />

analysis for the cities of Treia (MC) and Cagli (PS). Numerical simulations are<br />

performed using the SPEM. The reference event is a M=5.7 earthquake, and it is<br />

associated to two normal faults located beneath the cities, respectively. The second<br />

Figure shows a wavefield snapshot for the simulations performed for Treia.<br />

Marche<br />

Microzonation<br />

Project. Wavefield<br />

snapshot (P-SV<br />

acceleration) of the<br />

SPEM simulation<br />

performed for the<br />

town of Treia<br />

(Macerata, Marche,<br />

Italy). The locations<br />

of the town and<br />

source are<br />

indicated at the top<br />

and by the couple<br />

of arrows,<br />

respectively.<br />

Within the project entitled Damage Scenarios in the Veneto-Friuli Area, the<br />

activity of the first year aims at building-up a ground shaking scenario at a regional<br />

scale for the area surrounding the town of Vittorio Veneto. The reference<br />

earthquake is the M=5.8 event, which occurred in “Cansiglio” on October 14, 1936.<br />

The associated fault mechanism is oblique, with a strong character of reverse fault.<br />

Both point source models and rupture propagation along an extended fault are<br />

considered. The dominant feature of the geological structure is the thrust<br />

corresponding to the Alpago-Cansiglio highlands, the southern edge of which<br />

features a very steep change of elevation. The reference earthquake is associated to<br />

this thrust. The particular location of the source, which is rather deep and right in<br />

the middle of a structural discontinuity, makes it possible to tackle this study with<br />

a non classical approach, that is differentiating the surface structure and top<br />

elevation in different zones, i.e., the plane area, the foot-hill zone, the alpine<br />

valleys, and the mountain area, respectively. The areas of maximum ground<br />

shaking predicted in this way agree very well with those observed in the<br />

macroseismic field (third Figure) even just by using a point source.<br />

31


Collaborations: National Group for the Defence Against Earthquakes (GNDT) and<br />

National Group of Volcanology (GNV), belonging formerly to the National Council<br />

of Research (CNR) and currently to the National Institute of Geophysics and<br />

Volcanology (INGeV); Politecnico di Milano; Università della Basilicata (Potenza);<br />

Università di Catania; Università di Camerino; Università La Sapienza (Roma);<br />

Università di Napoli; Università di Udine; Università di Trieste; National Institute of<br />

Geophysics and Volcanology.<br />

A<br />

B<br />

Damage Scenarios in the Veneto-Friuli Area.<br />

The October 14, 1936 (M=5.8) “Cansiglio” earthquake.<br />

(A) Observed macroseismic field, and (B) peak ground<br />

accelerations predicted numerically by the WIM.<br />

32


Seismic inversion<br />

Coordinator: Aldo VESNAVER<br />

The Cat3D tomographic software<br />

G. BÖHM<br />

M. PERONIO<br />

The software package Cat3D allows the 3D traveltime inversion by adaptive<br />

irregular grids, for arbitrary recording geometry and combinations of different<br />

wave types: direct, reflected, refracted and diffracted arrivals. Since it allows, for<br />

example, the joint inversion of surface and VSP <strong>data</strong>, it is a valuable tool for<br />

calibrating the seismic surveys with the well information.<br />

During the year 2000, the software commercialization was supported by the <strong>OGS</strong><br />

partners: Paneura (Trieste, Italy) and Fact (Houston, USA). Paneura had a booth at<br />

the EAGE International Conference and Exhibition in Glasgow. Various copies of<br />

Cat3D were installed in Europe.<br />

A major improvement of the software efficiency was an aggressive use of the<br />

dynamic memory allocation, which allowed increasing the computational speed,<br />

the size of processed <strong>data</strong> and the model complexity. Several graphic features were<br />

improved during the year, and new Import/Export model formats Surfer ® and<br />

GeoQuest ® , which were added to the existing ones (GoCad ® and Jason<br />

Geosystems ® ). The software installation was made easier and its documentation<br />

significantly extended. A particular effort was spent for the software porting for all<br />

major Linux distributions (Slackware, Red Hat, Mandrake, Caldera, Debian) in<br />

low-cost personal computers; of course, the Cat3D package runs also in different<br />

hardware platforms (as IBM, Sun and SGI) under the Unix operating system.<br />

Furthermore, we added a user-friendly menu for defining grids with a circular<br />

symmetry, which we prepared and tested for 3D VSP’s.<br />

Cover image of<br />

the user manual for<br />

the tomographic<br />

Cat3D package.<br />

33


Joint 3D inversion of P, S and converted waves<br />

A. VESNAVER<br />

G. ROSSI<br />

During the last few years, ocean-bottom cables (OBC) were used at a large scale for<br />

3D marine surveys for recording the three components of the particle velocity and<br />

the pressure waves. Also on land surface surveys and in VSP’s this technology is<br />

expanding, since it allows detect and separate P, S and converted waves. P and S<br />

velocity and, in particular, their ratio, are important lithological parameters to<br />

reconstruct porosity and pore pressure in the rocks. This information is important<br />

to delineate boundaries and preferential migration paths both in hydrocarbon<br />

reservoirs and in shallow water-saturated rocks.<br />

In most cases, rock interfaces reflect both P and S waves, and convert part of their<br />

energy from one wave type into the other one. When these events are observable<br />

and can be reliably picked, we can enhance a lot the 3D lithological description of<br />

the Earth, because of the <strong>data</strong> redundancy with respect to the estimated<br />

parameters. In fact, the unknowns describing the interfaces do not increase, and<br />

those ones for the velocities just double; viceversa, if P, S, P-S and S-P waves are<br />

available, including both reflected and head waves, the traveltime <strong>data</strong> increase by<br />

a factor much larger than two. Thus, the Earth model is much better constrained<br />

by the experimental <strong>data</strong>, and more robust with respect to errors in the<br />

traveltimes.<br />

The joint inversion of P, S and converted waves has an additional advantage. The<br />

sought velocity fields of P and S waves are almost decoupled, when considering<br />

pure P and S arrivals: their only connection are the possible common reflecting<br />

interfaces in the Earth. Converted waves provide new equations in the tomography<br />

inversion, which directly relates the two velocity fields.<br />

We set up a software prototype for the elastic inversion by generalizing our<br />

previous acoustic code. The ray tracing algorithm is based on the Fermat’s<br />

principle, and determines the ray path for a given ray signature by a bending type<br />

approach. For example, in a P-S converted wave, we use the P velocities up to the<br />

conversion point, where we switch to the S field, for computing and minimizing<br />

the traveltime along the ray path. The P and S velocities are estimated by an ART,<br />

SIRT or weighted SIRT algorithm.<br />

Time-lapse 3D tomography<br />

A. VESNAVER<br />

F. ACCAINO<br />

G. DAL MORO<br />

G. MADRUSSANI<br />

G. ROSSI<br />

The response in a seismic survey over a producing hydrocarbon reservoir changes<br />

during the years, due to the different pressure conditions of gas and the<br />

movements of the gas/oil/water interfaces. This information is crucial for<br />

optimizing the hydrocarbon production, locating the new wells in the unswept<br />

areas, and where the hydraulic conductivity is expected to be adequate, as in<br />

fractured or permeable formations. The 4D-TAIL Project is aimed at detecting<br />

these variations of lithologic parameters by Amplitude Versus Offset (AVO) analysis<br />

and 3D tomographic imaging. The <strong>OGS</strong> partners are two oil companies, i.e.<br />

TotalFinaElf UK and Norsk Hydro, and the University of Milan. This project is<br />

supported by the European Union in the Thermie Programme.<br />

34


The joint inversion of<br />

direct, reflected and<br />

refracted P waves<br />

allows reconstructing<br />

the velocity and depth<br />

in a model (above)<br />

without any<br />

redundancy, unlike in<br />

the elastic case (below).<br />

For a proper comparison of the differences in the seismic response, one has to<br />

compensate all changes that do not depend on the reservoir itself, as the recording<br />

geometry and equipment. A factor usually neglected that we studied is the seasonal<br />

variation of the seawater velocity, due to changing currents and the temperature<br />

of the mixed layer. These variations can be of the same order of magnitude (or even<br />

larger) than those expected at the reservoir; furthermore, since they are spatially<br />

organized, they can be attributed to the reservoir, so totally distorting the timelapse<br />

analysis. In the Figure, we see that two 3D surveys acquired at the North sea<br />

in 1989 and 1992 provide very different estimates for the sound speed in the sea<br />

water, obtained by the joint inversion of reflected and head waves. (We remark that<br />

these plots can be quite interesting for oceanographic studies too).<br />

Outside the reservoir, we do not expect variations of the seismic velocities and the<br />

reflector structure. Thus, we can impose that the Earth model in depth, obtained<br />

35


from different <strong>data</strong> vintages, is unique. At the reservoir, the P and S velocities<br />

should not be constrained at all, because their variations (and the resulting ratio)<br />

is the primary information that we are looking for. Thus, we will have a set of<br />

coupled models, one of each vintage year, that are mostly identical, except in the<br />

upper layer and at the reservoir. When a proper amplitude-preserving surfaceconsistent<br />

processing is carried out, including a vintage cross-calibration, we can<br />

carry out a pre-stack depth migration of each vintage, and compare the reflectivity<br />

changes at the reservoir. AVO can add further details, and guide the design of the<br />

tomographic grid.<br />

Sound speed in the seawater estimated by the joint inversion of<br />

reflected and head waves in the year 1989 (above) and 1992 (below)<br />

in the same area at the North Sea.<br />

Seismic tomography for environmental studies<br />

G. ROSSI<br />

F. ACCAINO<br />

G. BÖHM<br />

G. DAL MORO<br />

G. MADRUSSANI<br />

M. PERONIO<br />

A. VESNAVER<br />

The oil and gas industry pushed seismic technology at advanced levels by<br />

significant investments, which rarely are available for environmental studies.<br />

However, except for a scale factor, many practical problems encountered in<br />

hydrology within the shallowest Earth layers are very similar to those ones<br />

considered in the hydrocarbon reservoirs. Seismic tomography is a possible<br />

example: the inversion of velocity anomalies can be generally related to lateral<br />

facies variations of the geological formations and, sometimes, to variation of other<br />

properties of hydraulic interest – as permeability and porosity. If both P and S<br />

waves can be jointly inverted, various analytic expressions exists which related<br />

these velocity fields to the fluids’ pressure and saturation.<br />

36


Transferring technology from the oil and gas industry to environmental<br />

application is a goal we pursued within the MICA Project, funded by the “Fondo<br />

Trieste”. We acquired a high-resolution 3D survey close to the Trieste airport, in an<br />

area where the city aqueduct has several catchment wells. The tomographic<br />

inversion and imaging of this <strong>data</strong> showed not only a good correlation with the<br />

available stratigraphy from 3 wells, but also remarked the significant vertical and<br />

lateral variations of the layers overlying the water-saturated formations. This fact<br />

proves the viability of 3D surveys for the protection of groundwater from pollution,<br />

and also their need: sparse 2D profiles would miss the area complexity;<br />

furthermore, drilling sparse exploration wells would be not only even less<br />

significant (and probably more expensive), but could break the impermeable<br />

covers over the groundwater pool, allowing it to be polluted by industrial or<br />

agricultural activities.<br />

Complementary to the seismic surveys, we carried out also a Georadar and a geoelectric<br />

survey, for integrating the elastic and electro-magnetic parameters. Of<br />

course, the penetration scale of the different techniques is quite different (but<br />

complementary), and their integration into a consistent physical model is<br />

ongoing. Partner of <strong>OGS</strong> within the MICA Project is GeoKarst, which is a<br />

geochemical company located in the Trieste Area Science Park. This company is<br />

complementing the <strong>OGS</strong> activity by analyzing the water origin and fluxes by<br />

measuring the isotopes it contains.<br />

The stratigraphy<br />

of 3 water well (above)<br />

matches the pre-stack<br />

depth migrated section<br />

(below) obtained by<br />

tomographic velocities<br />

in a high-resolution<br />

3D survey.<br />

37


Geophisical interpretation<br />

Coordinator: Angelo CAMERLENGHI<br />

The activity of the group in 2000 was characterised by intensive <strong>data</strong> <strong>acquisition</strong> at<br />

sea and on land. At the same time <strong>data</strong> processing and interpretation have been<br />

carried out in the <strong>OGS</strong> headquarters. Data <strong>acquisition</strong>, focussed to the<br />

understanding of the geology of continental margins, was mainly in Antarctica,<br />

within the projects WEGA (Wilkes Land Glacial History), LARSEN (Deglacial<br />

History of the Larsen Ice Shelf - Weddell Sea), ODP Leg 188 in Prydz Bay<br />

(M. Rebesco was shipboard sedimentologist), and TESAC (Tectonic and Cenozoic<br />

Evolution of the South America-Scotia Plate Boundaries) sponsored by PNRA. In<br />

addition, we participated in a cruise of the EU project STRATAGEM<br />

(Stratigraphical Development of the Glaciated European Margin) on the Faroe-<br />

Shetland margin. The <strong>data</strong> we collected span from single and multi channel<br />

seismic reflection, chirp sonar, core logging <strong>data</strong>, gravity at sea, to structural<br />

geology <strong>data</strong> and lake bathymetry on land. Of particular relevance for the group<br />

was the continuing activity with the recently acquired Geotek multi-sensor core<br />

logger, and the utilisation of the Datasonic chirp sonar in Antarctica.<br />

Processing and interpretation were on seismic, gravity, magnetic <strong>data</strong>, core<br />

samples, core logging, downhole logging, and oceanographic <strong>data</strong> collected in the<br />

previous years within the research programs sponsored by PNRA on the Pacific<br />

Margin of the Antarctic Peninsula (Projects ODP Leg 178, SEDANO, Sediment<br />

Drifts of the Antarctic Offshore, and BSR, Bottom Simulating Reflectors), in the<br />

South Scotia Sea (Crustal Structure and Evolution of the Powel Basin), and in the<br />

Ross Sea (Cape Roberts Drilling Project and Evolution of the West Antarctic Ice<br />

Sheet). Research was conducted with the highest possible level of international<br />

cooperation and by supporting research fellowships at other research institutions<br />

such as INGV, Rome, and the University of Trieste.<br />

The research activity involving the processing of Synthetic Aperture Radar (SAR)<br />

satellite images has continued with applications to areas such as Antarctica, Tierra<br />

del Fuego and the surroundings of city of Trieste. In parallel, we started a new<br />

project as a joint venture with SOPROMAR in a contract with the Italian<br />

Geological Survey to produce test sheets of the marine gravity maps of the coastal<br />

zones of Italy.<br />

In the field of environmental geology and geophysics, we participated in research,<br />

within broader <strong>OGS</strong> projects, on the aquifers identification and mapping in the<br />

Friuli Venezia Giulia region. Another project has involved the identification of<br />

karstic cavities with micro-gravimetric techniques.<br />

Research applied to petroleum exploration included one project funded by AGIP on<br />

the kinematics of salt deformation in the Eastern Mediterranean and distal Nile<br />

cone (in cooperation with the University of Parma and IGM, Bologna), and the<br />

participation in the Ormen Lange Verification project sponsored by NorskHydro<br />

through SINTEFF on the theme of gas hydrates.<br />

38


At the end of the year, we started a doctoral Program in Polar Sciences (applied to<br />

polar continental margin evolution) in cooperation with the University of Siena.<br />

Michele Rebesco hosted the Workshop “Seismic expression of contourites<br />

and related deposits” in the framework of the IUGS-UNESCO International<br />

Geological Correlation Programme n. 432 (Contourites, Bottom Currents and<br />

Palaeocirculation).<br />

We illustrate below the major scientific results achieved in the year 2000.<br />

Joint Italian/Australian Marine Geoscience Expedition<br />

to the George V Land Region of East Antarctica<br />

(Wilkes Land Glacial History, WEGA project)<br />

G. BRANCOLINI<br />

M. BUSETTI<br />

C. PELOS<br />

L. SORMANI<br />

R. VIDMAR<br />

A collaborative Italian PNRA/Australian AGSO-Antarctic CRC marine geoscience<br />

research voyage to the George V Land sector of the East Antarctic continental<br />

margin was carried out in February-March, 2000, on board the of the RV Tangaroa.<br />

A total of 1827 km of multi-channel seismic <strong>data</strong> (2 x 150 cu.in. GI airguns, 600 m<br />

streamer length), 562 km of Chirper (2.5-7kHz and 8-21 kHz transducers, 1000m<br />

armoured cable) sonar <strong>data</strong>, 11 gravity cores, 28 piston cores, 18 surface grabs and<br />

11 short trigger cores were collected on the voyage. Water profile (CTD)<br />

measurements and water samples were collected at nine stations and seabed<br />

bottom photographs were made at 11 stations.<br />

The expedition discovered and mapped a shelf current-derived, sediment drift<br />

deposit called the “Mertz Drift”, covering about 400 km 2 lying in an >800m deep<br />

section of the George V basin west of the Mertz Glacier.<br />

On the continental rise multi-channel seismic <strong>data</strong> were taken across contourite<br />

drift deposits and a submarine canyon system in 2500 to 3500 m water depth.<br />

Piston cores were collected along the profile of one drift deposit which gave a<br />

preliminary Mid-Pliocene age to truncated strata that crop out on the drift’s<br />

steeper lee side. These <strong>data</strong> will provide crucial information about the Antarctic<br />

late Cenozoic glaciations and useful site-survey support of a proposal sent to the<br />

Ocean Drilling Program under the auspices of the SCAR-ANTOSTRAT project for<br />

drilling key sites along the Antarctic margin.<br />

Deep Tow Chirper profile W-16 in the George V basin (continental shelf) and<br />

core stations across the “Mertz Drift”. The drift is over 35 m thick and it is composed<br />

of laminated, anoxic, gelatinous olive green, siliceous mud and diatom ooze (SMO).<br />

39


Physical properties and seismic stratigraphy of ODP<br />

Leg 178 well sites, Antarctic Peninsula Pacific margin<br />

V. VOLPI<br />

A. CAMERLENGHI<br />

M. REBESCO<br />

P. CORUBOLO<br />

U. TINIVELLA<br />

C. DE CILLIA<br />

In this work, sponsored by PNRA as participation to the scientific activity of the<br />

Ocean Drilling Program, we have re-analysed the porosity, bulk density and<br />

seismic velocity information collected from three bore holes on the continental<br />

rise of the pacific Margin of the Antarctic Peninsula. The purpose is to provide a<br />

comprehensive, composite digital <strong>data</strong> set of <strong>data</strong> readily available for future<br />

studies aimed at well-seismic correlation. The work originates from the<br />

occurrence of overlapping sets of physical parameters and acoustic velocity<br />

collected with different methods (downhole logging, core logging, laboratory<br />

determination, derivation from seismic <strong>data</strong>), and in different holes of the same<br />

site. These <strong>data</strong> not always provide the same information.<br />

Composite vertical profiles of velocity (sonic logs, core logs, and measurements on<br />

samples), density (RHOM - LDS corrected bulk density) and porosity (APLC - APS<br />

Near/Array limestone porosity corrected) have been obtained by combining <strong>data</strong><br />

from different instruments and different holes. The comparison between core and<br />

log porosity for site 1095 shows that, except for the lower part of the section,<br />

downhole porosity is systematically larger than core sample porosity due to<br />

the poor open hole conditions encountered in a fine grained, generally<br />

underconsolidated formation.<br />

At site 1095 additional information on acoustic velocity comes from the velocity<br />

check-shots, obtained from a vertical seismic profile (VSP), while at site 1096<br />

additional information comes from acoustic tomographic inversion of travel times.<br />

These interval velocities can be compared with the nearest available stacking<br />

velocity and the velocities obtained from core samples, at site 1095. The in situ<br />

velocity check-shots provide a more reliable velocity information than the stacking<br />

velocity throughout the section, while the vertical seismic profile (VSP) provides a<br />

reliable tie to the site survey MCS.<br />

40


ODP Leg 178 - Site Location<br />

Site 1095 Porosity<br />

(g/cm 3 )<br />

A<br />

B<br />

Site 1095 Line 195-135A Site 1095<br />

NW<br />

SE<br />

mbsf<br />

TWT time (s)<br />

interval velocity (m/s)<br />

C<br />

D<br />

A) Location map of ODP Sites 1095, 1096, and 1101 (highlighted), together with all<br />

other sites drilled during ODP Leg 178. B) Comparison between the downhole<br />

logging (APLC, APS Near/Array limestone porosity corrected) porosity (red curve )<br />

and porosity from the index properties (IP) measurements in the laboratory.<br />

C) Comparison among interval velocities obtained with the in situ velocity check<br />

shots (VSP), stacking velocities and core log velocity. D) Vertical Seismic Profile and tie<br />

in two way travel times between lithostratigraphic and seismostratigraphic units.<br />

41


Orbitally-Controlled rhythmic sedimentation<br />

in the Wild Drift, Antarctica (ODP Leg 188, Site 1165)<br />

M. REBESCO<br />

This research, conducted in collaboration with J. Gruetzner (Bremen University,<br />

Germany) is an outcome of the participation of an <strong>OGS</strong> researcher (Michele<br />

Rebesco) to the ODP Leg 188 in Prydz Bay – Cooperation Sea (Antarctica). Leg 188<br />

cruise began in Fremantle (Australia) on 10 January 2000 and ended in Hobart<br />

(Australia) on 11 March 2000. Three sites were drilled on continental shelf, slope<br />

and rise to document onset and fluctuations of East-Antarctic glaciations.<br />

Site 1165 is situated in a water depth of 3357m on the continental rise in front of<br />

the outlet of the Lambert glacier-Amery Ice Shelf system that today drains 22% of<br />

East Antarctica. The site cored a 999-m lower Miocene-Holocene section into an<br />

elongate sediment body (Wild Drift) formed by the interaction of westward-flowing<br />

currents with the sediment supplied from the shelf. Alternations with wavelengths<br />

ranging from cm to m size between a greenish grey diatom bearing clay facies and<br />

dark grey clay facies with silt laminations are apparent throughout the hole back<br />

to early Miocene time. Furthermore the greenish intervals are characterised by<br />

lower density, susceptibility and iron content. The dark grey intervals are<br />

interpreted as contouritic facies deposited during maximum ice advances whereas<br />

the greenish sediments indicate hemipelagic sedimentation under warmer climate<br />

conditions.<br />

Analysis of high resolution colour photo-spectrometer <strong>data</strong> reveals that the colour<br />

cycles are best described by the ratio of the reflectivity in the green colour band<br />

and the average reflectivity (grey).<br />

Spectral analyses on depth and time series of the investigated parameters over<br />

selected intervals demonstrate that variance is dominated by orbital frequencies as<br />

predicted by the Milankovitch theory. The detected obliquity and precession cycles<br />

allow a refined evaluation of sedimentation rates.<br />

Subsidence at the Cape Roberts drill sites<br />

(Ross Sea, Antarctica) from backstripping techniques<br />

L. DE SANTIS<br />

G. BRANCOLINI<br />

The tectonic subsidence of the western margin of the Victoria Land basin has been<br />

estimated from the physical properties and ages of the sediments in the Cape<br />

Roberts Project drill cores 2/2A and 3, using backstripping techniques, assuming a<br />

local isostatic compensation. The sediment load effects is removed from the total<br />

subsidence and the tectonic contribution through time at each location is<br />

calculated.<br />

The analysis indicates a total tectonic subsidence of about 660 m at this location<br />

between 34 Ma and the present time. Two main trends are defined, i) about 230<br />

m/m.y. from 34 Ma to 32.5 Ma, and ii) about 23 m/m.y. from 32.5 Ma to 21 Ma.<br />

Since 21 Ma, the subsidence is not well constrained. Extrapolation indicates a very<br />

low subsidence rate, but uplift within this period may have greatly affected the<br />

estimate.<br />

42


Core 1165B-14H<br />

Composite diagram showing the cyclicity at different scales nested together in the<br />

colour record from Core 1165-14H. The cm-scale cycles are evidenced by colour<br />

banding and lamination (see the interpreted black and white photo in the left<br />

bottom). Such cycles are included within dm-scale cycles evidenced by lighteningup<br />

intervals. In turn, these cycles are included within dark grey facies (see the black<br />

and white photo in the centre and the synthetic interpretation on its right). Finally<br />

m-scale cycles are constituted by a couplet of dark grey and greenish grey facies.<br />

The dm- and m-scale cycles are precisely recorded by colour photo-spectrometer<br />

<strong>data</strong> (see the ratio of the reflectivity in the green colour band versus the average<br />

reflectivity on the right side of the diagram). Moreover (not shown here) larger<br />

scale (tens to hundreds of m) cycles are produced by the variation in ratio between<br />

the thickness of the two facies of the m-scale couplets.<br />

43


The age and the calculated fast rates of the tectonic subsidence affecting the<br />

western Ross sea until 21 Ma is consistent with previous studies made in the<br />

Central and Eastern Ross Sea that suggest an Oligocene age of the basin opening<br />

phase in those regions. Since 20 Ma, extrapolation of the tectonic subsidence<br />

curves indicates a period of very low subsidence. Seismic reflection <strong>data</strong> across the<br />

VLB indicates that extensional tectonics was diachronous within the VLB and it<br />

was progressively younger toward east and possibly toward south. We believe that<br />

the apparent major slow down of the overall subsidence rate after about 20 Ma may<br />

be the result of subsequent uplift of the region and migration of the extensional<br />

tectonics toward east and south.<br />

A<br />

TECTONIC SUBSIDENCE<br />

A. Tectonic subsidence curves obtained<br />

using backstripping technique at the CRP-2<br />

(squares) and the CRP-3 (circles) drill sites.<br />

The dashed curve represents the<br />

extrapolated tectonic subsidence of the<br />

basement at CRP-3 site between 31 Ma<br />

and present. This curve is inferred from the<br />

trend at CRP-2 (from 30 to 21 Ma) and the<br />

depth of the basement caused by the<br />

tectonics at present time at the CRP-3 site.<br />

The error bars indicate the range of the<br />

uncertainties in the paleo-depth<br />

information and age.<br />

B. line drawing of a composite seismic<br />

section made by IT69, US 403 and US 404<br />

across the VLB about 5 km north of the<br />

CRP drill sites.<br />

WEST<br />

EAST<br />

B<br />

44


Cenozoic Evolution of the South Orkney Microcontinent<br />

M. BUSETTI<br />

A. MARCHETTI<br />

The South Orkney Microcontinent is one of the fragments of the South Scotia<br />

Ridge between the Antarctic and Scotia Plate. It is a key area for both Cenozoic<br />

paleoclimatic and geodynamic studies as it is a remnant of the separation between<br />

Antarctica and South America about 26 Ma. The separation permitted the onset of<br />

the Circum-Antarctic current and the following climate isolation of the Antarctic<br />

continent.<br />

Due to the complexity of the geodynamic setting of the area and the different<br />

tectonic regime (transform, convergence and divergence) in the surroundings, the<br />

small microcontinent has complicate geological evolution. Basins formation on<br />

the microcontinent is related to the South Scotia Ridge fragmentation in the<br />

?Eocene/Oligocene time. The northern margin of the South Orkney<br />

Microcontinent is an obliquely convergent plate boundary dominated by<br />

transcurrent condition exhibiting strain partitioning of convergent motion<br />

accommodating by a thrust zone in the oceanic area and vertical strike-slip zone<br />

at the border of the steep escarpment of the SOM. Between the microcontinent and<br />

the Bruce Bank the plate boundary changes, exhibiting transcurrent and<br />

extensional regime.<br />

Structural scheme of the South Orkney Microcontinent and bathymetric<br />

contours of the area. The continental shelf of the microcontinent is delimited<br />

by the 1000-meter contour. Extensional tectonic in the Eocene/Oligocene time<br />

produced several north-south elongated basins. The northern margin of the<br />

microcontinent is the plate boundary between Antarctic and Scotia Plates.<br />

45


Structure and Cenozoic evolution of the South America –<br />

Scotia plate boundary in the Tierra del Fuego region<br />

E. LODOLO<br />

R. GELETTI<br />

Onshore field geological and geophysical studies, and a multichannel seismic<br />

survey have been conducted in the last two years in the Tierra del Fuego region, in<br />

the frame of an Argentinean-Italian scientific research. The main aim of this<br />

project was to analyse the regional geological setting of the Island and reconstruct<br />

the Cenozoic geodynamic evolution of an important segment of the South<br />

America-Scotia plate boundary, called Magallanes-Fagnano fault system. This is a<br />

mainly wrench lineament which cuts across the Island and runs from the Pacific<br />

entrance of the Magallanes Strait to the Atlantic coast of the Island. The Lago<br />

Fagnano, located in the central part of the Tierra del Fuego, is an E-W-trending<br />

major depression which hides part of the fault, as revealed by the bathymetric map,<br />

which shows the presence of significant and steep scarps in correspondence of the<br />

onshore prosecution of the lineament. In cross-section, this tectonic lineament is<br />

represented by sub-vertical faults and associated asymmetric basins, generated by<br />

simultaneous strike-slip motion and transform-normal extension, as imaged by<br />

the seismic profiles acquired off the Atlantic coast of the Island and in the central<br />

and western Magallanes Strait.<br />

Data analyses support the interpretation that the Magallanes-Fagnano fault system<br />

is remarkably transtensive in nature, and is structurally and temporally<br />

superposed on the older tectonic framework of the Tierra del Fuego (i.e., the<br />

contractional system of the Magallanes fold and thrust belt), even if the<br />

displacement history of this fault system is unclear. The near parallelism among<br />

the younger and older lineaments suggests that the development of the<br />

transtensional structures may have reactivated pre-existing weakened zones<br />

formed by the Cretaceous-Tertiary shortening.<br />

Mapping the BSR on the South Shetland Margin<br />

(Antarctica) and assessing gas hydrate and free gas<br />

quantities<br />

E. LODOLO<br />

A. CAMERLENGHI<br />

G. MADRUSSANI<br />

G. ROSSI<br />

U. TINIVELLA<br />

Bottom Simulating Reflectors (BSRs) along the South Shetland continental<br />

margin have been first identified by <strong>OGS</strong> researchers on two multichannel seismic<br />

profiles acquired by the R/V <strong>OGS</strong>-Explora (1990 Antarctic Campaign). A dedicated<br />

survey was conducted on 1997 to purposely map the extent of the BSR on this<br />

margin, and study the relationships between geological structure and gas hydrate<br />

and free gas distribution. Processing and interpretation of the collected grid of<br />

high-resolution multichannel seismic reflection profiles have been completed<br />

(about 700 km of <strong>data</strong>), and have allowed us to map the lateral extent of the BSRs.<br />

The South Shetland continental margin, an accretionary wedge located off the<br />

northern tip of the Antarctic Peninsula, consists of two distinct and superimposed<br />

tectonic regimes: an older regime is related to Mesozoic - Middle Cenozoic<br />

subduction-related tectonism; a younger one is associated with a mainly<br />

46


Magallanes-Fagnano<br />

master fault<br />

TWT (s)<br />

Tectonic sketch of the Magallanes-Fagnano fault system across the<br />

Tierra del Fuego Island, as revealed from the interpretation of seismic<br />

profiles acquired at the Pacific entrance<br />

of the Magellan Strait (section A), in the central part of the Strait<br />

(section B), and off the Atlantic coast of the Island (section C).<br />

The profile on the right (section C) shows the presence of a principal<br />

sub-vertical fault (the master fault), with associated asymmetric<br />

basin, generated by simultaneous strike-slip motion and transformnormal<br />

extension along the strike of the fault system.<br />

47


extensional tectonic phase, and related to the Oligocene development of the<br />

western Scotia Sea. The occurrence of the BSR appears to be controlled by the<br />

geological structure of the margin. The BSR lacks continuity near basement<br />

structures, main geological discontinuities and faults. On the contrary, the<br />

amplitude and continuity of the BSR are not affected by the presence of folded<br />

structures and undeformed sedimentary layering. We found that the BSR is mostly<br />

confined to the N-E sector of the South Shetland Margin, where propagation of<br />

faulting associated to the Shackleton Fracture Zone may have driven migration of<br />

natural gas towards the surface and created the conditions for a BSR to appear. The<br />

application of reflection tomography techniques allowed us to reconstruct the<br />

averaged seismic velocity field between the seafloor and BSR in order to map the<br />

depth of BSR. By averaging the observed velocity structure above and below the<br />

BSR, and applying a theoretical model of elastic wave propagation in porous media,<br />

we attempted as rigorously as possible a quantitative assessment of the natural gas<br />

present as gas hydrate above the BSR and as free gas between the BSR and the Base<br />

of Gas Reflector (BGR).<br />

We found that in this reservoir, where free gas layers are widespread and of<br />

significant thickness, the amount of free gas trapped beneath the gas hydrate layer<br />

is about two orders of magnitude less than the natural gas stored as hydrate phase.<br />

The obtained values, compared with global estimates of total natural gas trapped<br />

in and beneath gas hydrate reservoirs proposed by several authors, represent<br />

1/10000 of the global estimate and are about twice that estimated in the Prudhoe-<br />

Kuparuk oil field (Alaska), and one order of magnitude less than the Bering Sea<br />

reservoir, which extends however over an area extremely larger than the South<br />

Shetland margin reservoir.<br />

Gas hydrate physical properties imaging by multi-attribute<br />

analysis - Blake Ridge BSR Case History<br />

F. COREN<br />

V. VOLPI<br />

U. TINIVELLA<br />

In this project, partially supported by the European Union under the project<br />

‘Application of Seismic Attribute Analysis for Reservoir Characterisation’ (CEE/NIS<br />

(INCO/COPERNICUS 1995 – 1996), we analysed the Bottom Simulating Reflector<br />

(BSR) of the Blake Ridge (Figure A) through the “multi-attribute” analysis, by<br />

using the EMERGE (Hampson – Russel) software. This technique allows<br />

prediction of petrophysical/geological parameters along seismic lines (velocity,<br />

porosity, density, resistivity etc.) starting from log <strong>data</strong>. Its application consists of<br />

three steps: i) tie the well log to seismic <strong>data</strong>; ii) train the seismic to predict the<br />

reservoir parameter of interest at the tie locations and iii) apply the result of the<br />

training to the seismic <strong>data</strong>set. The algorithm that allows the correlation between<br />

logs and seismics is the generalised multiple linear regression. The resulting<br />

function is then applied to the seismic profile, generating a target log-predicted<br />

section. We obtained five sections (VSP and P-wave velocity, density, porosity and<br />

resistivity), each of them representing the distribution of the corresponding<br />

property along the profile. The <strong>data</strong> set is constituted by two seismic lines and<br />

48


(Top): Structural map<br />

of the BSRs occurrence<br />

along the South<br />

Shetland margin,<br />

where the differences<br />

in strength of the BSRs<br />

have been highlighted.<br />

(Bottom): Part of a<br />

seismic profile where<br />

the BSR crosses an<br />

anticline fold, without<br />

loss in amplitude and<br />

strength.<br />

W<br />

E<br />

49


three wells in the Blake Ridge, offshore southern United States. Here, a strong<br />

‘BSR’ marks a transition from a hydrate rich sediments zone above, to a free gas<br />

bearing sediments below. This transition is reflected in the velocity profile with a<br />

boundary at 4150 ms between a high velocity region (1670 m/s above) and a low<br />

velocity one (1500/1600 m/s below) (Figure B).<br />

A velocity structure obtained by prediction can be translated in gas hydrate and<br />

free gas concentration structure. We used the method proposed by Tinivella<br />

(1999). The concentration is estimated by fitting the theoretical velocity to the<br />

experimental P/wave velocity (derived by the prediction of the VSP in our case).<br />

The discrepancies between the inverted velocity profile and the velocity for water<br />

filled marine sediments are interpreted as due to the presence of the gas hydrate<br />

(where positive anomalies are present) and free gas (where negative anomalies are<br />

present). Figure C shows the distribution of the two phases: positive values (red<br />

colours) are the gas hydrate concentration, while negative concentrations (blue<br />

colours) are related to the percentage of volume occupied by free gas.<br />

A<br />

VSP velocity Section<br />

B<br />

A) Position map of the<br />

study area. Wells<br />

location and seismic<br />

line position are shown<br />

in the close-up<br />

window.<br />

B) VSP velocity panel.<br />

C) Distribution of the<br />

two gas phases;<br />

positive values (red<br />

colours) are the gas<br />

hydrate concentration,<br />

while the negative<br />

concentrations (blue<br />

colour).<br />

C<br />

50


Physical properties of sediment cores<br />

from the Antarctic continental margins<br />

M. BUSETTI<br />

Antarctic paleoclimate and paleoenvironmental studies have a key role in the<br />

understanding of the global climate changing. Waxing and waning dynamic of the<br />

Antarctic ice coverage is the most important effect at high latitudes as response to<br />

climatic variations. As sediment deposited on the continental shelf and rise record<br />

the ice fluctuations, they are an essential task in the paleoclimate study. Analysing<br />

the acoustical and physical properties of sediment cores it is possible to investigate<br />

glacial and interglacial stages. P-wave velocity, bulk density and magnetic<br />

susceptibility are closely related to sediment composition and may reflect changes<br />

in grain-size distribution or in the ratio of terrigenous (quartz and clay) and<br />

biogenic components. Generally terrigenous sediment characterise glacial<br />

deposits, while biogenic material is present in the interglacial deposits.<br />

Acoustical and physical property logs of a marine piston core collected on the<br />

continental rise of the Antarctic margin. P-wave velocity and bulk density reflect<br />

the grain size distribution: the higher values in the upper part of the core are<br />

related to silt, and even greater to sand, the lower values in the deeper part<br />

of the core are related to the clay fraction. Ice rafted debris (black spots in the<br />

visual log) produce clear spikes in the magnetic susceptibility measurements.<br />

51


Backstripping modelling in the frame of the<br />

Stratigraphical Development of the Glaciated European<br />

Margin (STRATAGEM) - EU project<br />

L. DE SANTIS<br />

S. CERAMICOLA<br />

A. CAMERLENGHI<br />

The Project is funded by the European Community within the 5 th framework and<br />

it is part of the Ocean Margin Deep Water Research Consortium (OMARC) cluster.<br />

It is Co-ordinated by Dan Evans (British Geological Survey, UK).<br />

The principal objectives of STRATAGEM (web site: www.stratagem-europe.org) are<br />

to address problems related to the mid-Cenozoic to Recent development of the<br />

glaciated north European margin, that extends from northern Norway to central<br />

Ireland. The end product of the project will be the development of a Margin<br />

Evolution Model, to assess the interactions of sedimentary processes that have<br />

produced the character of the present day margin and the forcing factors (such as<br />

climate, tectonics and palaeoceanography).<br />

The seismic <strong>data</strong> to be interpreted will comprise existing <strong>data</strong> held by the partners,<br />

new <strong>data</strong> to be acquired in STRATAGEM, and <strong>data</strong> from industry, and from the<br />

IMAGES project.<br />

The <strong>OGS</strong> INTE group contributes to the project in:<br />

1) 300 km of Single channel (90 in 3 GI gun, 16 m streamer length) seismic <strong>data</strong><br />

collection (and processing) in the Faroe Island continental margin (R/V DANA<br />

cruise);<br />

2) 2D flexural, post-rift backstripping modelling on several transects along the<br />

length of the continental margin, in particular in the Vøring, in Shetlands-<br />

Faroe Islands and in the Hebrides margins. The models will constrain the<br />

more-qualitative Margin Evolution Model, and provide insights into factors<br />

such as uplift and subsidence that have been among the crucial controls on<br />

margin development. The observed stratigraphy is modelled using flexural<br />

backstripping combined with decompaction and reverse thermal subsidence<br />

calculations by mean of a commercial software package (Flex-Decomp by<br />

Badleys Ltd.-UK). Reliable direct palaeobathymetry information from drill sites<br />

and indirect palaeobathymetry markers, such us erosion surfaces, are used to<br />

constrain the subsidence history. During the first year of the project several<br />

tests have been carried out to verify the input parameters (e.g. stretching factor<br />

and the lithosphere elastic thickness).<br />

3) the comparison with the Antarctic glaciated margins, that will place the results<br />

of the STRATAGEM study in a global perspective.<br />

52


Bathymetric map of the North Atlantic continental margin studied in the STRATAGEM<br />

frame. WP1 is the Work Package 1 study area (the Norwegian margin); WP2 is the Work<br />

Package 2 study area (the South Shetland-Faeroe margins); WP3 is the Work Package 3<br />

study area (the Rockall Trough and Porcupine Basin).<br />

Earth gravity field: measurements, <strong>data</strong> processing<br />

and interpretation<br />

C. ZANOLLA<br />

F. PALMIERI<br />

F. COREN<br />

C. DE CILLIA<br />

<strong>OGS</strong> has a long tradition in gravimetry (i.e. the international gravity station<br />

network IGSN71, Italian gravity map, marine gravity surveys in the Mediterranean<br />

sea). Gravity surveys (land or marine) can be managed from planning, to field<br />

<strong>acquisition</strong>, processing and interpretation. The INTE group has the responsibility<br />

to manage and operate three LaCoste-Romberg gravity meters (one model D and<br />

two model G), one underwater gravity meter LCR model H, one marine surface<br />

gravity meter Bodenseewerk KSS31. These allow us to cover all the possible<br />

applications of gravimetry:<br />

At present our group is involved in the following projects:<br />

1) A joint venture group <strong>OGS</strong> and SO.PRO.MAR has been entrusted by Servizio<br />

Geologico Nazionale to carry out two marine gravity surveys in areas located<br />

near Rome.<br />

53


The purpose is to merge the land and marine gravity <strong>data</strong> to obtain a higher<br />

spatial definition of the gravity anomalies near the coastline and a gravity map<br />

that is overlapped with the geological maps at scale 1:50.000. In particular we<br />

have been required an underwater gravity survey and a surface marine gravity<br />

survey. The first survey has been completed while the second one will be<br />

carried out in february 2001.<br />

2) Two micro-gravimetry nets have been established in order to study the<br />

temporal gravity variations associated with: a) water table fluctuations to study<br />

the effective porosity of the geological formations involved (Fagagna area); b)<br />

the subsidence problems that affect some districts of the city of Trieste; this<br />

survey is jointly carried out with SAR techniques.<br />

In the near future we will repeat measurements in a network located in<br />

seismically active areas of Friuli-Venezia Giulia.<br />

3) Several micro-gravity surveys have been planned in order to detect the<br />

presence of cavities in the Carso area surrounding the city of Trieste: One<br />

survey is ongoing in the “Grotta Doria”, with the aim to model the gravity<br />

anomalies and to compare the results obtained with several geophysical<br />

techniques applied to the same geological feature.<br />

4) Our group is active also in <strong>data</strong> Antarctica were we acquired a gravity transect<br />

(see Figure) crossing the Wilkes Basin (East Antarctica).<br />

AHVRR image of<br />

Antarctica with<br />

superimpose the ITASE<br />

traverse along which<br />

the gravity profile has<br />

been acquired (above).<br />

Free air gravity profile<br />

acquired across the<br />

Wilkes Basin, East<br />

Antarctica (below).<br />

Free Air gravity anomaly (mGal)<br />

54


Synthetic Aperture Radar (SAR) remote sensing<br />

F. COREN<br />

R. VIDMAR<br />

P. STERZAI<br />

Horizontal velocity<br />

field of the<br />

Ligosullo landslide<br />

(northeast Italy)<br />

derived by SAR<br />

interferometry<br />

(yellow contours)<br />

superimposed<br />

to a aerial<br />

image of the area.<br />

Since 1996 <strong>OGS</strong> developed within the INTE group its own research line in<br />

Synthetic Aperture Radar (SAR) remote sensing. The activity is now mainly<br />

addressed to interferometry and analytical products generation using ERS-1 and<br />

ERS-2 satellites.<br />

Research activity have been focused on two main tasks:<br />

1) Monitoring terrain deformation (landslide) using differential interferometry<br />

techniques and analysis of the hydrological setting using SAR remote sensed<br />

<strong>data</strong>. In this field TS-SAR is the major project where SAR interferometry has<br />

been used to monitor ground deformation in the area of Trieste (Italy) in the<br />

period between 1996-2000 for the Civil Protection Department of the<br />

Municipality. A specific investigation has been carried out in industrialised and<br />

urban area to assess the possibility of a continuous monitoring of inferred<br />

subsiding phenomena. A set of complex interferograms has been computed in<br />

order to identify sliding and subsiding areas and to verify the general urban<br />

stability. A validation procedure on the interferogram has been applied<br />

considering only pixels characterised by high coherence in all the<br />

combinations. These pixels represent permanent backscatter and are mainly<br />

associated to civil structures and buildings. The displacement history of three<br />

warehouses has been computed on the temporal basis of 23 month. At the<br />

present, a validation has been carried out only by ground observation without<br />

elevation measurements. Other minor project has been also carried out to<br />

monitor and detect landslides in the Friuli area (see figure).<br />

2) Study of the glacial setting of East Antarctica (STARGLASS project financed by<br />

PNRA Italian National Antarctic Program) and Northwestern sector of<br />

Antarctic Peninsula (a self sustained research activity). Main goal of both<br />

projects is the computation of the ice velocity fields and digital elevation model<br />

of outcrop areas. Polar regions play an important role in the global<br />

environment. The potentiality of<br />

SAR interferometry for the<br />

monitoring of high latitude areas is<br />

today a well understood geophysical<br />

tool and represents a flexible and<br />

powerful method to study large polar<br />

sectors at low specific cost. In<br />

STARGLASS project we have used<br />

pairs of tandem images of satellites<br />

ERS-1 and ERS-2 synthetic aperture<br />

radar (SAR) of the David glacier, East<br />

Antarctica, for the purpose to outline<br />

the grounding line. All the SAR radar<br />

<strong>data</strong> have been provided under the<br />

VECTRA project cover (European<br />

Space Agency Announcement<br />

Opportunity 3.108) in raw format.<br />

55


Administration<br />

Coordinator: Dario Colonnello<br />

F. PETRONIO<br />

The activities of the administration office were, basically, to manage the purchase of<br />

raw material, scientific instruments, services and consumables necessary to carry out<br />

the various projects; personnel expenses for traveling to working sites and scientific<br />

meetings and conventions; invoices and statements of accounts and the legal aspects<br />

of the projects; and the Department’s inventory.<br />

The total turnover of the Department in 2000 was approximately 3.7 milion Euro’s of<br />

which we effectively spent 1.7 milion Euros. The real income related to specific<br />

programs amounted to more than 3.2 milion Euros.<br />

56


Publications<br />

Wave modeling<br />

ARNTSEN, B., AND CARCIONE, J. M., 2000, A new insigth into the reciprocity principle,<br />

Geophysics, 65, 1604-1612.<br />

ARNTSEN, B., AND CARCIONE, J. M., 2000, Numerical simulation of the Biot slow wave<br />

in water-saturated Nivelsteiner sandstone, submitted to Geophysics.<br />

CARCIONE, J. M., ARNTSEN, B., AND CAVALLINI, F., 2000, Simulation of ultrasonic waves<br />

in a natural sandstone. In Bermudez, A. et al. (Eds.), Mathematical and Numerical Aspects<br />

of Wave propagation, SIAM, 128-132.<br />

CARCIONE, J. M., AND SERIANI, G., 2000, Numerical simulation of wave propagation in<br />

frozen porous media. In Bermudez, A. et al. (Eds.), Mathematical and Numerical Aspects<br />

of Wave propagation, SIAM, 771-775.<br />

CARCIONE, J. M., GUREVICH, B. AND CAVALLINI, F., 2000, A generalized Biot-<br />

Gassmann model for the acoustic properties of clayley sandstones, Geophys. Prosp., 48,<br />

539-557.<br />

CARCIONE, J. M., AND GANGI, A., 2000, Non-equilibrium compaction and abnormal porefluid<br />

pressures: effects on seismic attributes, Geophys. Prosp., 48, 521-537.<br />

CARCIONE, J. M., AND POLETTO, F., 2000, Sound velocity of drilling mud saturated with<br />

reservoir gas, Geophysics, 65, 646-651.<br />

CARCIONE, J. M., AND GANGI, A., 2000, Gas generation and overpressure: effects on<br />

seismic attributes, Geophysics, 65, 1769-1769.<br />

CARCIONE, J. M., 2000, A model for seismic velocity and attenuation in petroleum source<br />

rocks, Geophysics, 66, 1080-1092.<br />

CARCIONE, J. M., AND SCHOENBERG, M., 2000, 3-D ground-penetrating radar<br />

simulation and plane wave theory, Geophysics, bf 65, 1527-1541.<br />

CARCIONE, J. M., 2000, AVO effects of a hydrocarbon source-rock layer, submitted to<br />

Geophysics.<br />

CARCIONE, J. M., AND CAVALLINI, F., 2000, Abnormal pore pressure and Poisson’s ratio,<br />

submitted to Geophys. Prosp.<br />

CARCIONE, J. M., 2000, Amplitude variations with offset of pressure-seal reflections,<br />

Geophysics, in print.<br />

CARCIONE, J. M., 2000, Energy balance and fundamental relations in dynamic anisotropic<br />

poro-viscoelasticity, Proc. Roy. Soc. London A, 457, 331-348.<br />

CARCIONE, J. M., AND SERIANI, G., 2000, Wave simulation in frozen sediments,<br />

J. Comput. Phys., in print.<br />

CARCIONE, J. M., FELICIANGELI, L., AND ZAMPARO. M., 2000, The exploding-reflector<br />

concept for ground penetrating radar modeling, submitted to Geophysics.<br />

CARCIONE, J. M., AND TINIVELLA, U., 2000, The seismic response to overpressure: a<br />

modeling methodology based on laboratory, well and seismic <strong>data</strong>, submitted to<br />

Geophys. Prosp.<br />

CARCIONE, J. M., CAVALLINI, F., AND MAINARDI, F., AND HANYGA, A., 2000, Timedomain<br />

seismic modeling of constant Q-wave propagation using fractional<br />

derivatives, Pure and Applied Geophysics, in print.<br />

CARCIONE, J. M., AND CAVALLINI, F., 2000, A semi-analytical solution for the propagation<br />

of electro-magnetic waves in 3-D lossy orthotropic media, Geophysics, in print.<br />

CARCIONE, J. M., AND POLETTO, F., 2000, Simulation of stress waves in attenuating drill<br />

strings, including piezoelectric sources and sensors, J. Acoust. Soc. Am., 108(1), 53-64.<br />

CARCIONE, J. M., PADOAN, G., AND CAVALLINI. F., 2000, Synthetic seismograms of the<br />

sea-bottom under different streamers conditions, Boll. Geof. Teor. Appl., 41, 21-29.<br />

57


CARCIONE, J. M., AND HERMAN, G., AND TEN KROODE, F. P. E., 2000, Seismic<br />

modeling, A review for Geophysics, submitted.<br />

CARCIONE, J. M., AND GEI, D., 2001, A seismic modeling study of Vostok lake, submitted<br />

to Journal of Glaciology.<br />

CARCIONE, J. M., AND CAVALLINI, F., 2000, Poisson’s ratio at high pore pressure, Norsk<br />

Hydro, E&P research centre, Bergen, NH-report R-089643.<br />

CARCIONE, J. M., MARCAK, H., SERIANI, G., AND PADOAN, G., 2000, GPR modeling<br />

study in a contaminated area of Krzywa airbase, Geophysics, 65, 521-525.<br />

CARCIONE, J. M., AND TINIVELLA, U., 2000, A modeling study based on laboratory, well<br />

and seismic <strong>data</strong>, Norsk Hydro, E&P research centre, Bergen, NH-report R-089737.<br />

CARNEVALE, G.F., CAVALLINI, F., and F. Crisciani, 2000, Dynamic boundary conditions<br />

revisited. J. Phys. Oceanogr., In press.<br />

CAVALLINI, F., Reply to comment by K. Helbig on “The best isotropic approximation of an<br />

anisotropic Hooke’s law” by F. Cavallini, Boll. Geof. Teor. Appl. 41, 1, 2000, 89-90.<br />

CAVALLINI, F., AND CRISCIANI, F., A generalized 2-D Poincare inequality. J. of Inequal. &<br />

Appl., 5:343-349, 2000.<br />

GUREVICH, B. AND CARCIONE, J. M., 2000, Gassmann modeling of acoustic properties of<br />

sand/clay mixtures, Pure and Applied Geophysics, 157, 811-827.<br />

HANYGA, A., AND CARCIONE, J.M., 2000, Numerical solutions of a poro-acoustic wave<br />

equation with generalized fractional integral operators. In Bermudez, A. et al. (Eds.),<br />

Mathematical and Numerical Aspects of Wave propagation, SIAM, 163-167.<br />

PHAM, N. H., CARCIONE, J. M., Helle, H. B., 2001, Poro-viscoelastic representation of<br />

shaley sandstones, 63th Ann. Internat. Mtg. Europ. Assoc. Expl. Geophys., Expanded<br />

Abstracts.<br />

POLETTO, F., CARCIONE, J. M., LOVO, M., AND MIRANDA, F., 2000, Acoustic velocity of<br />

SWD bore-hole guided waves, submitted to Geophysics.<br />

POLETTO, F., AND CARCIONE, J. M., 2000, On the group velocity of guided waves in drill<br />

strings, Submitted to J. Acoust. Soc. Am..<br />

PRIOLO, E., 2000, Deterministic computation of the reference ground motion in Fabriano<br />

(Marche, Italy). Ital. Geotech. J., in print.<br />

PRIOLO, E., 2000, Earthquake ground motion simulation through the 2-D spectral<br />

element method. J. Comp. Acoustics, in print.<br />

PRIOLO, E., 2000, 2-D spectral element simulation of the ground motion for a<br />

catastrophic earthquake. In: E. Faccioli and V. Pessina (Eds.), The Catania Project:<br />

Earthquake Damage Scenarios for High Risk Area in the Mediterranean. CNR-GNDT,<br />

Rome (Italy), 67-73.<br />

PRIOLO, E., AND MICHELINI, A., 2000, Measurments of environmental seismic noise for<br />

site response prediction. In: E. Faccioli and V. Pessina (Eds.), The Catania Project:<br />

Earthquake Damage Scenarios for High Risk Area in the Mediterranean, CNR-GNDT,<br />

Rome (Italy), 73-75.<br />

TINIVELLA, U., AND CARCIONE, J. M., 2000, Estimation of gas-hydrate concentration and<br />

free-gas saturation from log and sesimic <strong>data</strong>, The Leading Edge, in print.<br />

SERIANI, G., 2000, An iterative time-stepping method for solving first-order time<br />

dependent problems and its application to the wave equation, J. of Comp. Acoustics, 8(1),<br />

241-255.<br />

SERIANI, G., AND PRIOLO, E., 2000, Heterogeneous Chebyshev spectral elements for<br />

acoustic wave modelling. In: Onate, E. et al (Eds.), ECCOMAS 2000: Finite Element<br />

Schemes for Waves Problems. CIMNE, Barcelona (Spain), 13 pp., CD-ROM.<br />

VALLE, S., AND CARCIONE, J. M., 2000, Detection of liquid contaminants in the subsoil<br />

using the GPR technique, submitted to J. Appl. Geophys.<br />

58


Seismic inversion<br />

ROSSI, G., AND VESNAVER, A., 2000, Joint 3D traveltime inversion of P, S and converted<br />

waves, Journal of Computational Acoustics, 8, (in stampa).<br />

ROSSI, G., VESNAVER, A., AND PETERSEN, S., 2000, Anisotropy detection by tomography<br />

and polarization analysis in a 3D three-component VSP, First Break, (in stampa).<br />

BÖHM, G., GALUPPO, P., AND VESNAVER, A., 2000, 3D adaptive tomography by Delaunay<br />

triangles and Voronoi polygons, <strong>Geophysical</strong> Prospecting, 48, 723-744.<br />

BÖHM, G., MADRUSSANI, G., ROSSI, G., AND VESNAVER, A., 2001, Ray footprint and<br />

redundancy in seismic tomography, Journal of Seismic Exploration, 10, (in stampa).<br />

ROSSI, G., MADRUSSANI, G., AND VESNAVER, A., 2000, Adaptive 3D joint inversion of<br />

direct, reflected and refracted arrivals, in: Caiti, A., Hermand, J. P., Jesus, S. M., and Porter,<br />

M. B., Eds., Experimental Acoustic Inversion Methods for Exploration of the Shallow Water<br />

Environment, Kluwer, Dordrecht, 235-248.<br />

VESNAVER, A., AND BÖHM, G., 2000, Staggered or adapted grids for seismic tomography?,<br />

The Leading Edge, 19, 944-950.<br />

VESNAVER, A., BÖHM, G., MADRUSSANI, G., ROSSI, G., AND GRANSER, H., 2000, Depth<br />

imaging and velocity calibration by 3D adaptive tomography, First Break, 18, 303-312.<br />

<strong>Geophysical</strong> interpretation<br />

BARKER, P.F., CAMERLENGHI, A., and the ODP Leg 178 Shipboard Scientific Party, in<br />

press. Antarctic Glacial history, step 1: the continental margin drilled by ODP Leg 178. In:<br />

J. Gamble, D. Skinner, & S. Henrys (Ed), Proceedings of the VIII° International<br />

Symposium on Antarctic Earth Sciences, New Zealand Journal of Geology and Geophysics,<br />

Royal Society of New Zealand.<br />

BONACCORSI, R., BRAMBATI, A., BUSETTI, M., FANZUTTI, G.P., in press. Relationship<br />

among X-Ray Lithofacies, Magnetic Susceptibility, P-wave Velocity and Bulk Density in<br />

Core ANTA95-89C (Ross Sea, Antarctica): First Results. Proceedings of the Workshop<br />

“Ricostruzioni paleo-climatiche dai sedimenti marini del MarediRoss (Antartide) e<br />

dell’Oceano Meridionale”,Trieste, 26-27 novembre 1998. Terra Antartica.<br />

BUSETTI, M., MARCHETTI, A., ZANOLLA, C., DE CILLIA, C. and BELYAEV, V., in press:<br />

Seismic Structure and Stratigraphy of the South Orkney Microcontinent. In: J. Gamble, D.<br />

Skinner, & S. Henrys (Ed), Proceedings of the VIII° International Symposium on Antarctic<br />

Earth Sciences, New Zealand Journal of Geology and Geophysics, Royal Society of New<br />

Zealand.<br />

BUSETTI, M., ZANOLLA, C. and MARCHETTI, A. in press. Geological Structure of the<br />

South Orkney Microcontinenvt. Proceedings of the workshop: “Broad Band Observations<br />

and the Geodynamics of the Scotia Sea Region, Antarctica”, 25-26 October, 1999, Trieste<br />

(Italy), Terra Antarctica.<br />

CAMERLENGHI, A., REBESCO, M., DE SANTIS, L., VOLPI, V., in press. The Antarctic<br />

Peninsula Pacific Margin: modelling flexure and decompaction with constraints from ODP<br />

Leg 178 initial results. New Zealand Journal of Geology and Geophysics.<br />

COREN, F., LODOLO, E., CECCONE, G. submitted. Age Constraints for the Evolution of<br />

the Northern Powell Basin (Antarctica). Bollettino di Geofisica Teorica ed Applicata.<br />

DE SANTIS, L., DAVEY, F., PRATO, S., and BRANCOLINI, G,. submitted. Subsidence at the<br />

CRP drillsites from backstripping techniques. Terra Antartica, Scientific Report on CRP-3.<br />

DI VINCENZO, G., CABURLOTTO, A., and CAMERLENGHI, A., submitted. An 40 Ar- 39 Ar<br />

investigation of volcanic clasts in glaciogenic sediments at Sites 1097 and 1103 (ODP Leg<br />

178, Antarctic Peninsula). In Barker, P.F., Camerlenghi A., Acton, G.A. and Ramsay, T.(Eds.)<br />

Proc. ODP, Sci. Results, 178.<br />

59


FERRACIOLI, F., COREN, F., BOZZO, E., FREZZOTTI, M., ZANOLLA, C., GANDOLFI, S.<br />

AND TABACCO, I. submitted. Rifted(?) crust at the East Antarctic Craton margin: gravity<br />

and magnetic interpretation along a traverse across the Wilkes Basin. Earth and Planetary<br />

Sciences Letters.<br />

HARRIS, P.T., BRANCOLINI, G., ARMAND, L., BUSETTI, M., BEAMAN, R.J., GIORGETTI,<br />

G., PRESTI, M. and TRINCARDI, F., submitted. Continental shelf drift deposits indicate<br />

non-steady state Antarctic bottom water production in the Holocene. Nature.<br />

LA MACCHIA, C. and DE SANTIS, L., in press.Seismostratigraphic sequences analysis in<br />

the Prydz Bay area (East Antarctica). In the Proceeding volume of the Italian workshop on<br />

Antarctic Paleoclimate, Trieste Nov. 1998. Terra Antartica.<br />

LODOLO, E. and CAMERLENGHI, A., 2000. The occurrence of BSRs on the Antarctic<br />

Margin. In: M.D. Max (Ed.): Natural Gas Hydrate in Oceanic and Permafrost<br />

Environments, Kluwer Ac. Pub., 199-212.<br />

LODOLO, E., TASSONE, A., MENICHETTI, M., STERZAI, P. AND COREN, F., submitted.<br />

Superposed tectonic styles in the Tierra del Fuego region (southernmost South America).<br />

Terra Nova.<br />

LODOLO, E., CAMERLENGHI, A., MADRUSSANI, G., TINIVELLA, U. AND ROSSI, G., in<br />

press. Assessment of gas hydrate and free gas distribution on the South Shetland margin<br />

(Antarctica), based on multichannel seismic reflection <strong>data</strong>. Geophys. Journ. Intl.<br />

LUCCHI, R.G., REBESCO, M., BUSETTI, M., CABURLOTTO, A., COLIZZA, E., AND<br />

FONTOLAN, G., in press, Sedimentary Processes and Glacial Cycles on the Sediment Drifts<br />

of the Antarctic Peninsula Pacific Margin: Preliminary Results of SEDANO-II Project, In:<br />

J. Gamble, D. Skinner, & S. Henrys (Ed), Proceedings of the VIII° International<br />

Symposium on Antarctic Earth Sciences, New Zealand Journal of Geology and<br />

Geophysics, Royal Society of New Zealand.<br />

M. BRAUN, F. RAU, F. COREN AND H. SAURER. Submitted. Delimiting glacier drainage<br />

basins using remote sensing <strong>data</strong> of various sensor types and digital elevation models of<br />

different accuracies. Journal of Glaciology.<br />

PROTOPSALTI, I., IMMORDINO, F., DE SANTIS, L. FANZUTTI, G. P., in press. Sediment<br />

grain size and quartz grain morphology from Cape Roberts 1 core sample (Ross Sea):<br />

proxies for transport and depositional processes. In: J. Gamble, D. Skinner, & S. Henrys<br />

(Ed), Proceedings of the VIII° International Symposium on Antarctic Earth Sciences, New<br />

Zealand Journal of Geology and Geophysics, Royal Society of New Zealand.<br />

REBESCO, M., COOPER, A.K., O’BRIEN, P.E., and the shipboard Scientific Party, 2000.<br />

Southern Ocean Contourites - Preliminary Results from ODP Leg 188 in Prydz Bay,<br />

Antarctica. Comtourite Watch, issue 3, IGCP 432 newsletter, Southhampton<br />

Oceanography Centre, U.K.<br />

REBESCO, M., DELLA VEDOVA, B., CERNOBORI, L., AND ALOISI, G., 2000. Acoustic<br />

Facies of Holocene Megaturbidites in the Eastern Mediterranean. In: Shiki T., Cita M.,<br />

Gorsline D. (Ed), Sedimentary Features of Seismities, Seismo-turbidites and Tsunamites,<br />

Sedimentary Geology 135, 1/4 (Special Issue), 65-74.<br />

REBESCO, M., PUDSEY, C., CANALS, M., CAMERLENGHI, A., BARKER, P., ESTRADA, F.,<br />

GIORGETTI, A., in press, Sediment Drift and Deep-Sea Channel Systems, Antarctic<br />

Peninsula Pacific Margin. In: Stow D.A.V., Pudsey C.J., Howe J., & Faugeres J.C. (Ed), Atlas<br />

of Deep-Water Contourite Systems. Memoir of the Geological Society, Special publication.<br />

SAGNOTTI, L., MACRÌ, P., CAMERLENGHI, A., AND REBESCO, M., submitted.<br />

Environmental magnetism of Antarctic Pleistocene sediments and interhemispheric<br />

correlation of climatic events. Earth Planet. Sci. Lett<br />

SHIPBOARD SCIENTIFIC PARTY, 2000. Leg 188 Preliminary Report: Prydz Bay -<br />

Cooperation Sea, Antarctica: glacial history and paleoceanography. ODP Preliminary<br />

Report, 188 [online] Available from:<br />

<br />

60


TINIVELLA, U. AND LODOLO, E., 2000. The Blake Ridge BSR transect: tomographic<br />

velocity field and theoretical model to estimate methane hydrate and free gas quantities.<br />

Proceedings of the Ocean Drilling Program, Scientific Results, vol. 164. College Station<br />

(TX): 273-281.<br />

TINIVELLA, U., CAMERLENGHI, A., AND REBESCO M., submitted. Sesimic velocity<br />

analysis on the continental shelf transect, ODP Leg 178, Antarctic Peninsula. In Barker,<br />

P.F., Camerlenghi A., Acton, G.A. and Ramsay, T. (Eds.) Proc. ODP, Sci. Results, 178.<br />

VOLPI, V. CAMERLENGHI, A., MOERZ, T., CORUBOLO, P., REBESCO, M., AND<br />

TINIVELLA, U., submitted. Physical properties and seismic stratigraphy, continental rise<br />

sites 1095, 1096, and 1101, ODP Leg 178, Antarctic Peninsula. In Barker, P.F., Camerlenghi<br />

A., Acton, G.A. and Ramsay, T. (Eds.) Proc. ODP, Sci. Results, 178.<br />

61


Presentations at meetings and conventions<br />

Seismic <strong>data</strong> processing<br />

WARDELL N., DIVIACCO P., SINCERI R., 2000, Pre-Processing Corrections on Very High<br />

Resolution 3D Marine Seismic Data, 62 nd Annual International Meeting EAGE, Glasgow<br />

2000, Expanded Abstracts, P-130.<br />

DIVIACCO P., SINCERI R., WARDELL N., 2000, Estensione 3D per tecniche di<br />

preprocessing di dati sismici marini ad alta risoluzione, GNGTS, Rome 2000.<br />

Wave Modeling<br />

CARCIONE, J. M., GUREVICH, B., CAVALLINI, F., AND SERIANI. G., 2000, A generalized<br />

Biot-Gassmann model for the acoustic properties of shaley sandstones, 62th Ann. Internat.<br />

Mtg. Europ. Assoc. Expl. Geophys., Glasgow (UK), Expanded Abstracts, D35.<br />

CARCIONE, J. M., ARNTSEN, B., AND CAVALLINI, F., 2000, Simulation of ultrasonic waves<br />

in a natural sandstone, Fifth International Conference on Mathematical and Numerical<br />

Aspects of Wave propagation, Santiago de Compostela (Spain), July 10-14<br />

CARCIONE, J. M., AND SERIANI, G., 2000, Numerical simulation of wave propagation in<br />

frozen porous media, Fifth International Conference on Mathematical and Numerical<br />

Aspects of Wave propagation, Santiago de Compostela (Spain), July 10-14.<br />

CARCIONE, J. M., CAVALLINI, F., AND MAINARDI, F., 2000, Modeling constant-Q wave<br />

propagation with fractional derivatives, 70th Ann. Internat. Mtg., Soc. Expl.<br />

Geophys.,Calgary (Canada), Expanded Abstracts, 2345-2348.<br />

CARCIONE, J.M., CAVALLINI, F., MAINARDI, F., AND HANYGA, A., 2000, Time-domain<br />

seismic modeling of constant Q-wave propagation using fractional derivatives, Workshop<br />

Meeting on Seismic Wave in Laterally Inhomogeneous Media V, Zahradky, Czech Republic,<br />

June 5-9.<br />

CARCIONE, J. M., AND GANGI, A., 2000, Gas generation, overpressure and seismic<br />

properties, 70th Ann. Internat. Mtg., Soc. Expl. Geophys., Expanded Abstracts, 2405-2408.<br />

CARCIONE, J. M., GANGI, A., AND H. B. HELLE, 2000, Gas generation, overpressure and<br />

seismic properties (CD-Rom), III Conferencia Latinoamericana de Geofisica, Villahermosa,<br />

Mexico.<br />

CARCIONE, J. M., CAVALLINI, F., GUREVICH, B., AND SERIANI, G., 2000, A generalized<br />

Biot-Gassmann model for the acoustic properties of frozen porous media, Workshop on<br />

Seismic Signatures of Fluid Transport, Berlin, Germany.<br />

CARCIONE, J. M., CAVALLINI, F., AND SERIANI, G., 2000, Biot-type three-phase modeling<br />

of seismic wave propagation, EGS XXV General Assembly, Nice, France.<br />

CARCIONE, J. M., AND SERIANI, G., 2000, Electromagnetic properties of fluid<br />

contaminated soils using composite models, EGS XXV General Assembly, Nice, France.<br />

CAVALLINI, F., BOBBIO, M., PETTENATI, F., AND SIROVICH, L., 2000, ConVor, A newgeneration<br />

methodology for tracing objective and reproducible iso-seismals: the case of<br />

Feb. 28, 1925 Charlevoix earthquake in Canada. In EOS, Proceedings of AGU Spring<br />

Meeting, May 30 - June 3 2000, Washington, DC.<br />

FOSTER, M., LODOLO, E., TASSONE, A., GELLETTI, R., CARCIONE, J. M., 2000, Seismic<br />

structure and sedimentary setting of the souther Magallanes Basin off the Tierra del Fuego<br />

Island, Workshop Contiental Shelf, Buenos Aires, Argentina.<br />

HANYGA, A., AND CARCIONE, 2000, Numerical solutions of a poro-acoustic wave equation<br />

with generalized fractional integral operators, Fifth International Conference on<br />

Mathematical and Numerical Aspects of Wave propagation, Santiago de Compostela<br />

(Spain), July 10-14.<br />

HANYGA, A., AND CARCIONE, J. M., 2000, Numerical study of pulse delay effects in a poroacoustic<br />

wave equation, 70th Ann. Internat. Mtg., Soc. Expl. Geophys., Expanded<br />

Abstracts, 2337-2340.<br />

62


POLETTO, F., LOVO, M., AND CARCIONE, J. M., 2000, Acoustic velocity of drilling mud<br />

and SWD borehole guided waves, 70th Ann. Internat. Mtg., Soc. Expl. Geophys., Expanded<br />

Abstracts, 1763-1766.<br />

PRIOLO, E. 2000, Numerical simulation of the reference ground motion in Fabriano<br />

(Marche, Italy). Proc. 12th World Conf. on Earthq. Eng. (12WCEE), 30 January - 4<br />

February, 2000, Auckland, New Zealand. 8 pp., CD-ROM.<br />

PRIOLO, E., MICHELINI, A., AND LAURENZANO, G. 2000, Rapporti spettrali H/V di<br />

rumore sismico ambientale nel Comune di Catania. XIX Convegno GNGTS, Roma 7-9<br />

novembre.<br />

PRIOLO, E. 2000, Modellazioni numeriche del moto sismico del suolo a Catania e misure<br />

di rumore. Corso CISM-APT: “La riduzione del rischio sismico nella pianificazione del<br />

territorio (l’input sismico, le modellazioni, i valori di amplificazione), Lucca, 15-17<br />

novembre.<br />

TINIVELLA, U., AND CARCIONE, J. M., 2000, Estimation of gas hydrate concentration and<br />

free gas saturation from log and seismic <strong>data</strong>, 70th Ann. Internat. Mtg.,<br />

Soc. Expl. Geophys., Expanded Abstracts, 1568-1571.<br />

SERIANI, G., PRIOLO, E. 2000, Heterogeneous Chebyshev spectral elements for acoustic<br />

wave modelling. ECCOMAS 2000 - European Cong. on Comp. Meth. in Applied Sciences &<br />

Engng., Barcellona (Spain), 11-14 September.<br />

SIROVICH, L., CAVALLINI, F., PETTENATI, F., AND BOBBIO, M., 2000, ConVor: Un codice<br />

grafico per tracciare isosiste obiettive e riproducibili. In Convegno Nazionale GNGTS,<br />

Roma, 7 - 9 novembre.<br />

Seismic inversion<br />

BÖHM, G., GALUPPO, P., AND VESNAVER, A., 2000, Multiresolution in 3D seismic<br />

tomography within physical limits, Proceeding of ICTCA ’99 Conference, (in stampa).<br />

BÖHM, G., AND VESNAVER, A., 2000, The ray footprint in the joint 3D inversion of surface<br />

and well <strong>data</strong>: 62th Mtg. Eur. Assoc. Expl. Geophys., Extended Abstracts, Glasgow, P-160.<br />

VESNAVER, A., AND BÖHM, G., 2000, OBC versus conventional seismic <strong>data</strong> in 3D<br />

adaptive tomography: 70th Annual Internat. Mtg., Soc. Expl. Geophys., Expanded<br />

Abstracts, Calgary, 2269-2272.<br />

ROSSI, G., MADRUSSANI, G., AND VESNAVER, A., 2000, Tomographic inversion of the<br />

water layer in the 4D analysis: 70th Annual Internat. Mtg., Soc. Expl. Geophys., Expanded<br />

Abstracts, Calgary, 1287-1290.<br />

TINIVELLA, U., ACCAINO, F., CAMERLENGHI, A., 2000, Gas hydrate and free gas<br />

distribution from inversion of seismic <strong>data</strong> on the South Shetland margin (Antarctica).<br />

Sottomesso a <strong>Geophysical</strong> Prospecting.<br />

ACCAINO, F., BATINI, F., CORUBOLO, P., LOVO, M., PETRONIO, L., POLETTO, F., ROSSI,<br />

G., AND VESNAVER, A., 2000, Tomografia SWD con griglie sfalsate a simmetria radiale:<br />

Atti 19° Convegno Nazionale GNGTS, Roma.<br />

ACCAINO, F., BÖHM, G., MADRUSSANI, G., ROSSI, G., AND VESNAVER, A., 2000, Il<br />

problema “acqua” nella tomografia 4D: Atti 19° Convegno Nazionale GNGTS, Roma.<br />

DAL MORO, G., ACCAINO, F., BÖHM, G., MADRUSSANI, G., ROSSI, G., AND VESNAVER,<br />

A., 2000, Tomografia sismica 4D nel Mare del Nord: Atti 19° Convegno Nazionale GNGTS,<br />

Roma.<br />

DELLA MORETTA, D., MAZZOTTI, A., AND VESNAVER, A., 2000, Individuazione di<br />

anomalie di velocità tramite tomografia a riflessione: Atti 19° Convegno Nazionale GNGTS,<br />

Roma.<br />

ROBEIN, E., LAFOND, C., MAZZOTTI, A., AND VESNAVER, A., 2000, Time-lapse analysis<br />

by AVO and tomographic inversion at a producing field in the North Sea: Atti 19°<br />

Convegno Nazionale GNGTS.<br />

63


ROSSI, G., BÖHM, G., MADRUSSANI, G., AND VESNAVER, A., 2000, Seguendo le impronte<br />

dei raggi …..: Atti 19° Convegno Nazionale GNGTS, Roma.<br />

VESNAVER, A., AND BÖHM, G., 2000, Il modello iniziale nella tomografia sismica: Atti 19°<br />

Convegno Nazionale GNGTS, Roma.<br />

ROSSI, G., BUSETTI, M., BALLARIN, L., PIPAN M. E GRUPPO DI LAVORO MICA, 2000,<br />

Integrazione dei metodi geochimici e geofisici per lo studio idrogeologico: esempio di<br />

applicazione nella piana alluvionale dell’Isonzo: Riassunti dell’80 Riunione estiva della<br />

Società Geologica Italiana, Trieste, 411-413.<br />

ROSSI, G., ZADRO M. e EBBLIN, C., 2000. Processi geodinamici nell’Italia Nord-orientale:<br />

osservazioni e modellazione: Riassunti dell’80 Riunione estiva della Società Geologica<br />

Italiana, Trieste, 414-415.<br />

<strong>Geophysical</strong> interpretation<br />

BUSETTI, M., 2000. WEGA (Wilkes Land) Site Survey for ODP proposal 482 AGU 2000<br />

Spring Meeting, Washington, DC 30/05-03/06/2000, Supplement to Eos, May 9, 2000,<br />

p. S267.<br />

BUSETTI, M, 2000. Physical Properties from cores on the continental rise. WEGA Post-<br />

Cruise Workshop, Hobart (Tasmania, Australia), 6-11 December, 2000.<br />

CAMERLENGHI A., REBESCO M., DE SANTIS L., VOLPI V., DE ROSSI A. (2000),<br />

Modelling Flexure and Decompaction on the Antarctic Peninsula Pacific Margin with<br />

Constraints from ODP Leg 178, AGU 2000 Spring Meeting, Washington, DC 30/05-<br />

03/06/2000, Supplement to Eos, May 9, 2000, p. S267-268.<br />

CAMERLENGHI, A., COSTA, E., POLONIA, A., COOPER, C., FABRETTI, P., MOSCONI, A.,<br />

MURELLI, P., ROMANELLI, M., SORMANI, L., AND WARDELL, N., 2000. New insights on<br />

the mechanisms of deformation of the Eastern Mediterranean Ridge. EAGE Conference on<br />

Geology and Petroleum Geology of the Mediterranean and circum-Mediterranean Basins,<br />

Malta 1-4 October 2000. Extended Abstract Book.<br />

GRUETZNER, J., FORSBERG, C., REBESCO, M., 2000. Orbitally Controlled Sedimentation<br />

at the East Antarctic Continental Rise: Evidence from ODP Site 1165 (Leg 188, Prydz Bay)<br />

AGU 2000 Fall Meeting, December 15-19, 2000, San Francisco, California, Supplement to<br />

Eos, p. OS22A-05.<br />

LODOLO, E. AND TASSONE, A., 2000. The South America-Scotia Plate Boundary in the<br />

Tierra del Fuego Island: A <strong>Geophysical</strong> and Geological Study. 31th International<br />

Geological Congress, Rio de Janeiro, August 2000.<br />

LODOLO, E. TASSONE, A. MENICHETTI, M. COREN, F. STERZAI, P., 2000. Deciphering<br />

the morphostructure of the Tierra del Fuego region from remote-sensing and geophysical<br />

<strong>data</strong>. European <strong>Geophysical</strong> Society, XXV General Assembly, Nice, April 2000.<br />

MACRÌ, P., L. SAGNOTTI, A. CAMERLENGHI, M. REBESCO, F. FLORINDO, A.P.<br />

ROBERTS, AND A. WINKLER (2000), Environmental Magnetism and Paleomagnetism of<br />

Sediment Drifts from the Western Continental Rise of the Antarctic Peninsula, 25° EGS<br />

Assembly (Nice, 25-29/04/00, Abstracts).<br />

ODP LEG 188 SHIPBOARD SCIENCE PARTY (2000), Lithostratigraphy of Continental<br />

Shelf, Trough-Mouth Fan and Sediment Drift Deposits, ODP Leg 188, Prydz Bay, East<br />

Antarctica, AGU 2000 Spring Meeting, Washington, DC 30/05-03/06/2000, Supplement to<br />

Eos, May 9, 2000, p. S273-274.<br />

ODP LEG 188 SHIPBOARD SCIENCE PARTY (2000), Physical Property Changes as a Proxy<br />

for East Antarctic Sedimentation: First Results From ODP Leg 188 (Prydz Bay), AGU 2000<br />

Spring Meeting, Washington, DC 30/05-03/06/2000, Supplement to Eos, May 9, 2000,<br />

p. S272.<br />

REBESCO, M., CITA, M.B., HIEKE, W., DELLA VEDOVA, B., ALOISI, G., WERNER, F.,<br />

CERNOBORI, L., 2000. Deep-water Megaturbidites in the Eastern Mediterranenan, . EAGE<br />

Conference on Geology and Petroleum Geology of the Mediterranean and circum-<br />

Mediterranean Basins, Malta 1-4 October 2000. Extended Abstract Book.<br />

64


REBESCO, M., CITA, M.B., HIEKE, W., DELLA VEDOVA, B., ALOISI, G., WERNER, F.,<br />

CERNOBORI, L., 2000. Megatorbiditi Abissali Oloceniche Prodotte da Onda di Tsunami nel<br />

Mare Mediterraneo Orientale, Riassunti delle comunicazioni orali e dei poster, 80°<br />

Riunione Estiva della Società Geologica Italiana (Trieste, 6-8/9/2000), 401-402.<br />

REBESCO, M., COOPER, A.K., O’BRIEN, P.E., AND THE SHIPBOARD SCIENTIFIC PARTY<br />

(2000) Southern Ocean Contourites - Preliminary Results from ODP Leg 188 in Prydz Bay,<br />

Antarctica. Contourite Watch, issue 3, IGCP 432 newsletter, Southhampton Oceanography<br />

Centre, U.K.<br />

COREN, F., VIDMAR, R., STERZAI, P., 2000. Utilizzo di dati SAR per applicazioni di<br />

protezione civile nel comune di Trieste: il progetto TS-SAR – Atti della 3 Conferenza<br />

Nazionale ASITA – Napoli – Vol 1 pp. 627 – 632<br />

CAPRA, A., COREN, F., FREZZOTTI, M., MANCINI, F., STERZAI, P., VIDMAR, R., 2000.<br />

Verso Il Monitoraggio Ambientale dell’Antartide A Scala Globale –Il Progetto Vectra – Atti<br />

della 3 Conferenza Nazionale ASITA – Napoli – Vol 1 pp. 489 – 496<br />

COREN, F., STERZAI, P. VIDMAR, R., 2000. Interferometric Analysis of David Glacier (East<br />

Antarctica) – ERS ENVISAT Symposium 2000 – Goteborg – ESA.<br />

65


Book reviews<br />

VESNAVER, A., 2000, Review of the book “Numerical methods for wave equations in<br />

geophysical fluid dynamics” by Dale R. Durran, The Leading Edge, (in stampa).<br />

VESNAVER, A., 2000, Review of the book “Processing near-surface seismic-reflection <strong>data</strong>:<br />

a primer” by Gregory S. Baker, The Leading Edge, (in stampa).<br />

Educational video<br />

<strong>OGS</strong> co-produced an educational video entitled “PERCHÉ L’ANTARTIDE”, in which the<br />

main scientific activities in the field of earth sciences carried out in Antarctica are<br />

described with the aid of original video material. The video was produced also by the<br />

National Antarctic Museum (MNA), and CNR. The english version of the video has been<br />

presented at the Italian stand at the XXXI International Geological Congress held in Rio<br />

de Janeiro in August 2000. It is now available for sale at the Antarctic Museum<br />

(http://www.mna/it). Copies can be obtained also at <strong>OGS</strong> (spersoglia@ogs.trieste.it).<br />

Summary of the video:<br />

Title of the Italian version: Perché l’Antartide<br />

Title of the English version: Why Antarctica<br />

Produced by: CNR-IRPI-RCS/<strong>OGS</strong>/MNA<br />

Production: CNR-IRPI Reparto di Cinamatogrtafia Sceintifica<br />

Directed by: Teodoro Mercuri<br />

Science advisors: A. Camerlenghi, M. Parotto, F. Salvini, F. Talarico<br />

Editing: E. Valente<br />

Original Music: V. Ricca<br />

Archive Images: MURST - PNRA, ENEA, <strong>OGS</strong><br />

Computer Graphics: IMMAGINE SaS - Cosenza<br />

VHS PAL - Color - 37 Minutes<br />

Printed in the year 2000<br />

66


Visitors<br />

Peter F. BARKER, British Antarctic Survey and Xavier F. Molina, Cadiz University,<br />

visiting A. Camerlenghi within the project ODP Leg 178, Antarctic Peninsula.<br />

M.Y. MOSKALEVSKY, Institute of Geography, Russian Academy of Sciences,<br />

Moscow, Russia, visiting F. Coren within the VECTRA Project.<br />

Fred J. DAVEY, Inst. of Geological and Nuclear Science, Wellington, NZ, visiting<br />

Giuliano Brancolini within the Cape Roberts Drilling Project.<br />

Belinda BROWN, School of Earth Sciences, University of Sydney, Sydney,<br />

AUSTRALIA, visiting Laura DeSantis within the WEGA project<br />

Atle NYGÅRD, Department of Geology, University of Bergen, Bergen (Norway),<br />

visiting Laura De Santis within the STRATAGEM project.<br />

International seminars<br />

in solid earth geophysics<br />

Chairman: Fabio CAVALLINI<br />

Giovanni SANTARATO (University of Ferrara, Italy)<br />

Electrical tomography for environmental geophysics and cultural heritage: some<br />

examples<br />

Michele REBESCO (<strong>OGS</strong>, Trieste, Italy)<br />

Glacial history and paleoceanography: preliminary results of the cruise “ODP Leg<br />

188” and “WEGA” in Antarctica<br />

Jurgen MIENERT (University of Tromso, Norway)<br />

Storegga slide gas hydrate drilling on the mid-Norwegian margin<br />

Jeno GAZDAG (<strong>OGS</strong>, Trieste, Italy)<br />

The effects of regularization on 3-D pre-stack migration<br />

Steven R. PRIDE (University of Rennes, France)<br />

Electroseismic wave phenomena<br />

Steven R. PRIDE (University of Rennes, France)<br />

The theory of poroelasticity applied to exploration seismology<br />

Giovanni P. GREGORI, IFA (CNR, Roma, Italy)<br />

The origin of the magnetic field and of the endogenous energy of the Earth and<br />

of celestial bodies<br />

Gabriele PAPARO, IDAC (CNR, Roma, Italy)<br />

Acoustic emission as a diagnostic tool in geophysics<br />

Alfredo MAZZOTTI (University of Milan, Italy)<br />

Principles of AVO exploration<br />

67


Fred DAVEY (Institute of Geological & Nuclear Sciences, Wellington, New<br />

Zealand)<br />

<strong>Geophysical</strong> Investigation of a Modern Continental Transpressional Orogen: the<br />

Southern Alps, New Zealand<br />

Professor Sierd CLOETINGH (Vrije Universiteit, Amsterdam, The Netherlands)<br />

Intraplate tectonics and continental lithosphere evolution: models and constraint<br />

Professor Sierd CLOETINGH (Vrije Universiteit, Amsterdam, The Netherlands)<br />

Sedimentary basins and continental topography: from the Mediterranean to the<br />

Carpathian region<br />

Professor Sierd CLOETINGH (Vrije Universiteit, Amsterdam, The Netherlands)<br />

Continental rifts and rifted continental margins<br />

Dr. Maurizio BATTAGLIA (Stanford University, USA)<br />

GPS applications in the Earth Sciences<br />

Dr. Maurizio BATTAGLIA (Stanford University, USA)<br />

Temporal gravity investigations at Long Valley caldera<br />

Dr. Maxim Yu. MOSKALEWSKY (Institute of Geography, Moscow, Russia)<br />

Ice formation zones: probably the most sensitive indicators of short-term global<br />

changes<br />

Dr. Maxim Yu. MOSKALEWSKY (Institute of Geography, Moscow, Russia)<br />

Active Antarctic coastal zones as objects for remote sensing monitoring of<br />

environmental changes<br />

Peter F. BARKER (British Antarctic Survey, Cambridge, UK)<br />

The Antarctic Circumpolar Current and Antarctic glaciation<br />

Klaus HELBIG (Hannover, Germany)<br />

Seismic anisotropy for the rest of us<br />

Klaus HELBIG (Hannover, Germany)<br />

Singularities of the phase velocity of anisotropic media: specific examples for<br />

orthorhombic media.<br />

68


About Trieste, Italy<br />

Overwiew<br />

Trieste is geographically at the center of Europe and is a crossroads between<br />

the Central-European and Mediterranean cultures. Located on a strip of land<br />

where the white Karst cliffs plunge abruptly into the sea, Trieste has been strongly<br />

affected by its history both in terms of architecture and life style.<br />

Evidence of the contact between different ethnic groups can be found in the local<br />

dialect, in the family names of the inhabitants, and in the local cuisine that cannot<br />

be found elsewhere in Italy. In its welcoming restaurants, the savory flavors of mid-<br />

European dishes, or the simple cooking typical of the Northern Adriatic, the<br />

healthy and colourful cuisine of the south, or delicate seafood dishes can be<br />

experienced.<br />

The spirit of Trieste can be found in its cafes, pubs and buffets, which are the<br />

traditional meeting places in the life of the city. Cafes in fact are parlours where<br />

you can rest, read national or international magazines, and meet people.<br />

They have been the cultural and political meeting points for writers and artist,<br />

such as Italo Svevo and James Joyce.<br />

Museums and theaters<br />

There are many historical buildings, the most famous is probably the romantic<br />

Miramare castle, the residence of the luckless Maximilian of Habsburg<br />

and Charlotte of Belgium. Moreover, Trieste is rich in prestigious museums<br />

boasting numerous collections of great artistic value, and in theaters,<br />

such as Giuseppe Verdi Opera theater, opened by 1801, which is one of<br />

the centers of cultural life with its winter opera and concert season, and<br />

the Operetta Festival during the summer.<br />

69


Trieste and science<br />

Trieste is also a city of science, with a prestigious University and many worldrenowned<br />

laboratories, like the International Center for Theoretical Physics, the<br />

International Center for Genetic Engineering and Biothechnology, the Research<br />

Area Park with its Synchrotron, as well the School of Modern Languages for<br />

Interpreters and Translators, and the Osservatorio Astronomico. The <strong>OGS</strong> National<br />

Institute is located in a beautiful green<br />

area on the Karst plateau near the Giant<br />

Cave, visited by thousands of tourist<br />

each year.<br />

A bit of history<br />

The Roman origins of Trieste are still<br />

visible in the well preserved remains<br />

near San Giusto Castle, which is the<br />

symbol of the city, together with the<br />

Romanesque Cathedral, and the<br />

remains of the Roman Theater erected<br />

by Quinto Petronio Modesto in the II<br />

Century A.D.<br />

Because of its geographical position, as<br />

a crossroads between important trade<br />

routes, possession of the city was contended by Venice, France and then Austria.<br />

After realizing the commercial importance of Trieste, and as a consequence of the<br />

decline of the Venetian Republic, the Austrian Empire proclaimed it a free port,<br />

thereby laying the foundation of modern Trieste and making it one of the most<br />

important ports in Europe.<br />

70


Near Trieste<br />

One cannot miss a visit to Venice, 150 km from Trieste. The splendid beaches of<br />

Sistiana, Grado and Lignano are 20, 50 and 100 km away, respectively. Near Grado<br />

is the ancient Roman town of Aquileia, with a museum, well preserved open-air<br />

ruins and the largest Byzantine mosaic floor, dated 314 A.D.<br />

The magnificent Dolomite mountains are only 200 km away. The town of Udine,<br />

with numerous frescoes by the 17th century painter Giambattista Tiepolo and its<br />

International Center for mechanical Sciences, is 80 km away.<br />

Transportation and travel to Trieste<br />

By air: Trieste International Airport is 33 km from the city. There are direct flights<br />

to/from Rome, Milan, Genoa, Munich. A shuttle bus is available to town. In<br />

addition, the airport offers rental car facilities and 24-hour taxi service.<br />

By train: Trieste is connected to the national and international networks.<br />

The Railway Station is located at walking distance from the city center and main<br />

hotels.<br />

By car: Trieste is connected to the national and international motorways, and the<br />

A4 Venice-Trieste and Trieste-Ljubliana highways.<br />

By sea: there is a seasonal ferry service to/from Greece.<br />

Photo by Gabriele Crozzoli - Trieste<br />

71


72<br />

TRIESTE

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!