Geophysical data acquisition - OGS
Geophysical data acquisition - OGS
Geophysical data acquisition - OGS
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