Database of multi-parametric geophysical data from the TOMO-DEC experiment on Deception Island, Antarctica

Deception Island volcano (Antarctica) is one of the most closely monitored and studied volcanoes on the region. In January 2005, a multi-parametric international experiment was conducted that encompassed both Deception Island and its surrounding waters. We performed this experiment from aboard the Spanish oceanographic vessel ‘Hespérides’, and from five land-based locations on Deception Island (the Spanish scientific Antarctic base ‘Gabriel de Castilla’ and four temporary camps). This experiment allowed us to record active seismic signals using a large network of seismic stations that were deployed both on land and on the seafloor. In addition, other geophysical data were acquired, including bathymetric high precision multi-beam data, and gravimetric and magnetic profiles. To date, the seismic and bathymetric data have been analysed but the magnetic and gravimetric data have not. We provide P-wave arrival-time picks and seismic tomography results in velocity and attenuation. In this manuscript, we describe the main characteristics of the experiment, the instruments, the data, and the repositories from which data and information can be obtained.


Background & Summary
Multi-parametric geophysical experiments using marine and terrestrial techniques on active volcanoes are rare 1 . The TOMO-DEC experiment is a pioneering study that has combined several techniques to provide multi-disciplinary information that is more complete than previous studies at Vesuvius volcano 2 and Campi Flegrei 3 , Italy. We performed this experiment at Deception Island volcano (Antarctica; Fig. 1) using the Spanish oceanographic research vessel 'Hespérides' as a platform for generating active seismic sources and as a recording system for other geophysical data, including magnetic and gravimetric fields. We recorded seismic signals using a temporary network of 95 land-based seismometers deployed on Deception Island and 14 ocean-bottom seismometers (OBS) deployed on the floor of a flooded caldera (Port Foster) and in the Bransfield Strait around the island. These included 7 Lennartz Marslite seismic stations, 4 M24 instruments, and 11 seismic arrays, all working in continuous recording mode. We generated active seismic sources using an air gun array with a maximum capacity of 3,520 cubic inches (c.i.). Together with the seismic waveforms, we collected a final dataset consisting of travel times for more than 70,000 crustal P-wave first arrivals. We recorded bathymetric data using both multi-beam and single-beam sounders. We obtained magnetic profiles using a Marine Magnetics SEASPY marine magnetometer on the base of an Overhauser sensor. Finally, we acquired gravimetric data using the marine gravimeter BELL AEROSPACE-TEXTRON BGM-3. In total, the various data types accounted for more than 1,000 km of profiles.
Deception Island (Antarctica) is an active volcano of considerable scientific interest. Its origin and evolution remains controversial; however, most recently it has been interpreted as a collapse caldera resulting from the explosive eruption of basaltic-to-andesitic magmas 4 . A reconstruction of paleovolcanic geometries 5 reveals an initial stratovolcano island that reached 640 m above sea level. The present caldera hosts a very active geothermal system, which is responsible for most of the present-day seismic activity [6][7][8][9][10] . The importance of research at Deception Island is widely accepted and a recent review of multidisciplinary research performed on the island 11 showed that between 1964 and 2012 at least 173 impact papers have been published, with the number of new publications growing exponentially. This scientific interest is related to the continuous presence of multiple research teams, including those stationed at permanent bases (e.g., the Spanish Antarctic base 'Gabriel de Castilla' and the Argentinean base 'Deception'), or in temporary camps. The large number of researchers working in the volcano is a clear increasing of the volcanic risk 12 . However, this increase in volcanic risk is primarily the result of increasing visitor numbers, with thousands of tourists accessing the island each year 13 .
In this manuscript, we present a dataset that has been analysed using different approaches and techniques and which has provided new insights into the nature and structure of the island. Tomographic and structural images of seismic velocity 14,15 and attenuation 16,17 have been used to help in the interpretation of the inner structure of the volcano, sometimes combined with bathymetric analysis 18 . The interpretation of these results indicates the presence of strong lateral structural variations that affect seismic wave propagation 19,20 . The presence of many volcano-clastic deposits also introduces nearsurface low velocity effects 21 . This velocity heterogeneity can produce changes in the seismic-waveforms of the classic volcanic signals used to assess volcanic risk (e.g., LP events 22 ). Finally, when analysing the data, researchers considered the impact of unique local weather conditions, in which warm water increases the presence of biological activity 23 .
Although considerable information has already been extracted from the dataset, there remains the potential for additional research given the high quality of the data, the volume of unanalysed (magnetic, gravity) or partially-analysed data (bathymetry), and the recent advances in the field of signal processing, from automatic signal recognition techniques [24][25][26][27] to new methods for the determination of first arrivals [28][29][30] . We expect that our data is of broad interest and may yield new and high-impact results that will add to the existing scientific impact of the experiment. The database includes not only the raw data, but also seismological models, such as 3-D velocity and attenuation structures.

Methods
Deception Island volcano is a remote location where complex experiments require a large number of people, expensive logistical complexities, and complex material support including transport and accommodation. The success of the TOMO-DEC experiment extends beyond the above-cited scientific and social implications. Moreover, the results have provided additional support to other multidisciplinary approaches, including those in deformation and geodesy 31,32 , magnetic and gravimetric models 33,34 , and other tomographic and structural studies [35][36][37] . On this basis, we believe that this multidisciplinary database must be opened to the scientific community. Finally, it is important to remark that we provide derived models associated with the structure of Deception Island, mainly 3-D velocity structure and attenuation models, which will be important for scholars wishing to apply new advanced relocation techniques, such as has been done in several volcanic environments 38 .
Despite the large number of geological and geophysical studies performed at Deception Island, existing volcanological models lack complete information on the inner structure of the island. This shortcoming was acutely felt during a seismo-volcanic crisis that occurred during the 1998-1999 Austral summerautumn 8 . In the absence of realistic velocity and volcanological models, it was not possible to distinguish the real triggering process of this crisis and the nature (volcanic or tectonic) of the observed seismic stress migration. We conceived the TOMO-DEC experiment to improve our knowledge, and by implementing   a robust interpretation of the results 9 , we were better able to inform the emergency and contingency volcanic eruption plan 12 . This experiment was organised in collaboration with several research institutions. For fieldwork and data collection, the following institutions participated (each providing seismic instruments and/or researchers): Instituto Andaluz De Geofísica, Spain (

TOMO-DEC experiment
The TOMO-DEC experiment was organised into three main stages: (1) a search for sites to deploy seismic stations and establish temporary camps to accommodate researchers during the austral summer of 2003-2004; (2) field-based data collection during the austral summer of 2004-2005; and finally (3) laboratory-based data gathering, organisation, and analysis. Stage 3 remains on going and papers continue to be published. The preliminary analysis of data involved an initial phase of data gathering from instruments of different types, followed by their organisation in a joint database, their conversion to a common format, and finally the determination of the first P-wave arrival times.
Stage 1 was performed during two summer field expeditions (2002-2003 and 2003-2004). The whole island was surveyed (including glacial areas) in order to identify potential sites for seismic stations and camps. We assembled a handbook for use during the 2004-2005 summer deployment that included precise locations (by hand-held GPS) and cross-island access routes. We identified four sites for temporary camps (Whalers Bay, Pendulum Cove, Sea Wolf Cove, and Obsidian Valley; Fig. 1), each to house 6-8 researchers and include an adapted shelter tent, beds, food, and other necessities (Fig. 2). We identified 170 potential sites for seismic stations, including flat areas where seismic antennas could be deployed. This number was larger than the final number of sites used in the experiment.

Geometry of the experiment and observations
In December 2004, we deployed four radio-telemetered seismic stations at different sites on the island (one 3-component station and three vertical-component stations) to monitor the background level of seismic activity. In January 2005, 11 individual seismic stations and 10 arrays with a maximum of 12 channels each were deployed. This work included the resurvey of pre-established station sites, their installation, and a status check. For the single stations, this check consisted of a true plug-and-play and a real-time check in the field (Fig. 2). The array installation comprised individual sensors connected to the central hub of the acquisition system by 50-to 200-m-long cables. Data could be checked in a quasi-real time after conversion to a readable format (Fig. 2). Each seismometer was buried as deeply as possible in what were typically pyroclastic deposits over permafrost. In most cases the burial depth was~20 cm, and the deepest hole was 50 cm (Fig. 2). In glaciated areas, holes were dug to depths sufficient to cover the sensor and minimise wind noise. The central hubs of the acquisition systems of many stations were protected inside tents (Fig. 2). The position of each seismic station was initially determined using handheld GPSs; however, during the experiment, improved positions were obtained using the differential GPS system of the 'Hespérides' for both individual seismic stations and for the central hub of every seismic array. A laser theodolite was then used to measure the height of each array node with respect to the differential GPS solution for its central hub.
At the same time, a total of 14 OBSs were deployed. Of these, 11 were placed offshore of the island at depths between 200 and 686 m below sea level (b.s.l.), while 3 were placed inside Port Foster Bay at a maximum depth of 150 m. This procedure included on-deck instrument preparation, deployment on the seafloor, and acoustic verification of the operating status. The entire deployment of instruments took 3 days. During days 8 to 11, the first pattern of air gun shots was fired, both inside and outside the bay of Post Foster. After the first shooting phase, 10 of the autonomous seismic stations were relocated in order to increase the number of recording sites to 21. At the same time, 10 of the 14 OBSs were also recovered and redeployed at new locations. A second phase of air gun shooting was performed between 16 January and 18 January. The positions of all seismic stations deployed during the experiment are shown in Fig. 3. Following the second phase of air gun shooting, most of the seismic stations were recovered. Four M24type land stations and 12 OBS were left in their final positions until late February 2005 in order to record natural seismicity in and around Deception Island volcano.
The active sources for the seismic experiment were generated using air guns towed behind the R/V 'Hespérides' of the Spanish Research Council, and conducted by the Spanish Army. This process was conducted by the Spanish Army. The air gun shot lines were designed to provide dense coverage inside and around the island. In the interior of Port Foster, the tracks followed a dense grid of perpendicular lines with a shot spacing of 120 m on a 0.5 km grid and with shots fired every 60 s (Fig. 4). A safety distance from the coast of~500 m was respected. Offshore of the island, shots were fired every 120 s at intervals of 170-340 m along three main trajectories ( Fig. 4): (1) two straight tracks for refraction profiles, one 92 km long and oriented NNW-SSE and the other 55 km long and oriented WSW-ENE; (2) concentric lines around the island (out to 20 km from the centre of the island); and (3) lines in a 'star' configuration. The theoretical tracks were planned on the basis of station distribution, with the aim of sampling known tectonic structures with an optimum source and receiver configuration. For example, the seismic line oriented NNW-SSE was designed to obtain a 2-D section of the expected velocity contrasts NW of Deception Island 15 . For the inner part of the island (Port Foster), two rounds of shooting were performed. In the first round, only four guns were fired with a total capacity of 1520 c.i. (24.91 l or 0.024908 m 3 ). In the second round inside the bay, the total capacity was increased from 1520 to 2020 c.i. (33.1 l or 0.033101 m 3 ). Approximately 380 nautical miles (more than 700 km) of air gun shooting were collected. This configuration, along with the operational capability of the air guns allowed the recording of high-energy P-wave signals at distances as great as 40 km. The dataset consists of 6,630 useful shots.
While the ship was shooting, deploying, or recovering, additional geophysical and marine data were also recorded. In particular marine gravimetric data, magnetic measurements, and multi-beam ecosounder signals (using two types of instruments) were collected. Active seismic source measurements Seismic sources. Active seismic sources were generated using six BOLT 1,500LL air guns in an array configuration 12-m long ( Fig. 5) that was chosen to optimise the directivity of the source and minimise the water bubble. The bubble was generated using four Hamworthy 4TH190W70 compressors producing a maximum pressure of 140 bars. The air guns had capacities of 500, 1,000, 500, 255, 265, and 1,000 c.i., giving a total volume of 3,520 c.i. (57.7 l), and had inter-gun spacings of 0.8-2.5 m. The power of the system was decreased to 1,520 and 2,020 c.i. (as indicated above) for shallow-water shooting inside Port Foster. Each shot was recorded by a near-field hydrophone. The average of depth of the air gun array was 9 m b.s.l. In the link https://data.aad.gov.au/metadata/records/TOMODEC_2005_UTM-SPAIN we provide the whole information of the temporal and spectral characteristic of the air gun source, with the technical features of the air gun system detailed in Table 1. It is noteworthy that during the whole experiment no problems were detected and all of the planned shots were fired.
The air guns were controlled by a Hydra Minipulse system, which synchronised the guns and triggered the shots. Using information from a sensor located in each gun, the system can correct the moment of fire in order to synchronise shots to within 1 ms. The shot command is initially supplied by the GPS Seapath clock to the KONMAP navigation system (fix point) with a delay of 1 ms. This information is then sent to the gun control system within 60 ms. Summing up these time intervals, the real shot time has a delay of 61 ms with respect to the time written in the telegram supplied for every shot. This telegram contains information about the fix point (e.g., latitude, longitude, date, direction of vessel, speed, depth, name of the line of shots, and number of the fix point). The location information applies to the position acquired 60 ms earlier (after correction for the distance between the GPS and gun). After the shots, another message is appended to these data by the Minipulse system, detailing the number of working air guns, and the total capacity and delay of each gun with respect to the aiming point. All of these data are supplied for every shot and are recorded on a PC in text format. The whole telegram is made up of two parts: (1) data supplied by the Minipulse system (starting with the $HYDRA sequence), and (2) data from the navigation system ($HESSIS sequence). Owing to the huge amount of data, the output files were filtered in order to make them more manageable. Filtered files contain 6 initial fields (shot number, latitude, longitude, date, hour, fix), of which the first field was extracted from the shot telegram, while the others were extracted from the navigation telegram. In addition, a log file was supplied detailing system failures and errors. As a result of this organisation, the data appear as reported in Table 2.
The gun positions were assumed to be 125 m astern and 5 m to port of the GPS antenna, with a constant depth of 9 m. In the UTM conversion to geographical coordinates, the UTM region was 20E.  Seismic instruments. Lennartz Marslite autonomous stations: The 7 Marslite seismic stations used in the experiment were Lennartz LE-3Dlite MkII short-period sensors that record 125 samples per second using a 3-channel, 20-bit data acquisition system. This system includes a preamplifier with four 12 dB steps. Since digitisation always takes place at a fixed sampling rate of 4 kHz, the digital data stream needed to be decimated to the final sampling interval, after the signal had been digitally low-pass filtered by a FIR (Finite Impulse Response) filter. The system was completed with a GPSlite time signal receiver. The data were recorded in binary format on rewriteable 540 Mb magneto-optical disks. The acquisition systems were programmed to record in continuous mode. The power supply was provided from an external 12 V DC supply, later stabilised and regulated to the voltages required by the system.
Lennartz M24 seismic stations: Four of the stations used Lennartz M24 Compact/LP portable seismometers with 24 bits of dynamic range working with 3 built-in channels. There is one separate over sampled A/D converter per input channel, which digitises the signal coming from the preamplifier at very high frequency (20 kHz), until such point that the data stream is decimated down to the user's desired sampling rate: in our case, a sampling frequency of 100 Hz was configured. A single chip works as both a general-purpose CPU (Central Processing Unit) and as a DSP (Digital Signal Processor). The latter functionality is required to achieve the high-speed, high-precision digital filtering commonly associated with over sampled ADCs. GPS timing is performed thanks to a built-in receiver combined with an exterior antenna connected through a cable. The selected data format was of binary type, and recording was on a removable 2.5'' IDE 20 GB hard disk. For power supply, external 12 V DC power is accepted (in our case, a 100 A/h battery in combination with a 47 W solar panel) and converted to internally required voltages. The acquisition system is equipped with a Lennartz LE-3D (20 s) seismometer, which has a corner frequency of 0.05 Hz, an output sensitivity of 1,000 V/m/s, and damping of 0.707.
Seismic arrays: The seismic array modules were high resolution (24 nominal bits) with an acquisition system that can work with a maximum of 12 sensors, in array configuration. In this experiment, the acquisition modules were used with 3-12 channels (see below). Sampling frequency was configured to 100 samples/s and the gain was set to 1. Data acquisition was performed by four SEISAD18 plates, each managing three analog channels. The plates were synchronised thanks to a PLL plate, whose time pulse was supplied by a GPS (Garmin 35-HVS) connected with a serial port. Digital data from four plates were passed to a commercial PC (Lippert Cool Roadrunner II) with low energy consumption, and were finally stored on a 30 Gb USB hard disk. Moreover, a 10 GB hard disc was used for the operative system and for temporary data storage (moved every 6 h to the external hard disk). Central control of the process was achieved by means of a low-power industrial PC, working under MS-DOS. The acquisition software generated raw files in a binary.ARR format. The array box was powered using 12 V batteries. The array modules used vertical L-28B Mark Products sensors as seismometers, with preamplifier and a natural frequency of 4.5 Hz. Electronic extension permitted all of the sensors to achieve a flat response curve in the 1-50 Hz frequency interval. In this way, their response was almost identical to other 1 Hz commercial sensors. Some of the systems were additionally equipped with 3-component L-4C Mark Products sensors, without preamplifiers.
Ocean bottom seismometers (OBSs): The OBSs used in this experiment were provided by the Lamont-Doherty Laboratory of Columbia University. These instruments have low power consumption, recording capacities exceeding 18 GB, employ 24 bit digitizers, and are equipped with broadband 3-component seismometers with noise levels below 10-16 (m/s 2 )/Hz in the band from 0.01 to 50 Hz. These 'broadband' seismometers are short period sensors (1 Hz geophones) coupled to efficient and very low noise amplifiers. A broadband hydrophone (differential pressure gauge) provides an additional seismo-acoustic channel. The pressure sensor is used primarily to provide a redundant channel (looking very much like vertical velocity) in case the seismic sensor fails. The pressure sensor is considered to provide good Rayleigh and body wave records, as well as good active source and micro-earthquake data at shorter periods. The instruments sample continuously at 125 Hz, and are easily converted to standard PASSCAL sampling rates. During the experiment two of the OBSs were not recovered.
Radio-telemeters seismic stations: In order to control background seismic activity, four analog radiotelemetered seismic stations were deployed at strategic sites on the island, with deployment designed to cover as wide a range as possible (The Spanish base with a 3-component seismometer, Fumarole Bay, Obsidians Bay, and Pendulum Cove). These stations were equipped with Mark Products L4C sensors that have a natural frequency of 1 Hz. The telemetry system consisted of a radio-frequency transmitter with a directional antenna that sent the sensor signal to the Spanish base, where it was received and recorded digitally using a de-modulator and a 16-bit A/D converter 10 .

P-waves first arrivals
To determine P-wave first arrivals, we used the SAC routine named APK. With this algorithm, the detection of a pick is based upon abrupt changes in the ratio of a short term and long term running average of the signal. Once detected, the pick is subjected to an optional validation phase that attempts to distinguish a true event from seismic noise. Once validated, the pick is further evaluated to determine other characteristics of the event. To assess the performance of the APK routine, we compared automatically obtained picks with the real traces. Usually, the results for detected traces differed by a mean value of 0.02 s with respect to the handpicked arrivals. This can be considered an acceptable performance of the automatic routine. However, to be sure we did not use 'falsely triggered' arrivals, and to take in account as many phase arrivals as possible, we visually checked the automatic algorithm results. For this manual task, we used the UPICKER software, which permitted trace gather plotting, starting 2 s before shot time and with duration of 20 s. This choice depended on the known interval between shots, which was 120 s. We checked the error introduced by filtering the traces in SAC and found out that it is usually below 1 sample. The final database includes delay times from 70,411 crustal P-waves first arrivals automatically picked and manually checked. On average, each station recorded almost 1,500 shots.

Nomenclature of seismic stations
In order to provide a clear nomenclature of the different seismic stations, each type of sensor received a different name. Seismic stations are summarise in Table 3 and in the excel file of the repository data with all positions.

Multi-beam sounding
Multi-beam instruments. EM120 multi-bean sounder: To estimate the thickness of the water column below each shot, bathymetric data were collected throughout the period of the experiment using a multibeam sounding EM120, which was able to reach up to 11.000 m b.s.l. This information was essential for producing accurate seismic tomography of the region. The main features of the multi-beam sounding EM120 are detailed in Table 4. Bathymetric profiles were collected outside of Port Foster, allowing us to obtain a high definition map of the sea floor near Deception Island volcano. The bathymetric data also revealed potential new submarine volcanic edifices and their relationships with important geodynamic characteristics 18 . EM-1002 multi-beam sounder: An EM-1002 multi-beam sounder was used to collect very high precision sea floor maps. Designed to work in shallow waters (2-700 m depth), the system was operated simultaneously with the EM102 when the 'Hesperides' was operating inside Port Foster. Data collection stopped automatically when the depth of the seawater column was greater than 700 m. The main features of the multi-beam sounding EM-1002 are detailed in Table 5.   Potential field measurements SeaSPY marine magnetometer. The SeaSPY magnetometer (Marine Magnetics) is a highly sensitive/ precise acquisition system and omnidirectional Overhauser sensor. This sensor, which allows measurement of the total magnetic field, remained stable for the duration of the survey period. The system has a high-precision clock (1 ppm) that can be synchronised with the onboard GPS clock. It is thermally stabilised, permitting it to operate in all type of seawater, including Antarctic waters, and owing to its design there are no heading errors. The technical features of the SeaSPY magnetometer are summarised in Table 6. The magnetometer was used during the first phase of shooting outside of Post Foster; however, its 200-m long tow cable precluded its use within the Deception Island inner bay.
BGM-3 marine gravimeter. The gravimeter BELL AEROSPACE-TEXTRON BGM-3 (Lockheed Martin Federal Systems) is an acquisition system that can be used with both aerial and marine vehicles. This system is composed of a gravimeter sensor mounted on a gyro-stabilised platform, and an acquisition system. Raw data are recorded and processed using a HP-486/50 computer running the BGM software, which provides two types of data, raw and previously analysed.
Gravimeter calibration. The gravimeter is regularly calibrated by the manufacturer. In addition, before every field survey, it is calibrated in the starting harbour to avoid potential shifts. The mean accuracy of the system is ±2.0 mGal. All calibration controls are performed using a WORDEN mod. MASTER portable gravimeter. In this experiment, the control was performed at the Argentinean AeroNaval (Navy   and Air) base of Ushuaia on 29 December 2004, three days before the start of the field experiment. Figure 6 shows the reported calibration folder, written in Spanish, as reported by the onboard operator.

Water measurements
In order to obtain accurate sound velocities, we also collected water-temperature profiles. The SIPPICAN MK-21 bathy-thermo-graphic system was connected to several sensors, including temperature data (XBT sensor), sound velocity (XSV), conductivity and salinity (XCTD). This system acquires information in quasi-real time. The nominal resolution of the system is about 2% of water depth or ±0.15°C.

Data Records
The present database is a multidisciplinary package that contains raw data and derived results. The total size of the database is approximately 106 GB. The full data set is available from the Australian Antarctic Data Centre (AADC), which is a Public Scientific Data Repository. The AADC hosts parent metadata descriptor records that describe the whole experiment. These records are entitled 'TOMODEC_2005_-PROJECT-SPAIN' and can be accessed from the URL 'https://data.aad.gov.au/metadata/records/ TOMODEC_2005_PROJECT-SPAIN'. The main database is connected to the parent metadata records as child metadata records with their own DOI numbers. The child data records are divided according to discipline and instrument (see table Table 7 for a complete description of all child links). Additional information on the organisation of folders and data is provided in the 'Usage Notes' section of this manuscript. Continuous data from OBSs are also available in raw format from IRIS DMC (http://ds.iris.edu/gmap/ XU?timewindow = 2005-2005). Additional information of bathymetric data can be obtained from the Lamont-Doherty Earth Observatory IEDA Marine Geoscience Data System (http://www.marine-geo.org/). A summary of downloads that have been made of these data and documents is available from http://www. marine-geo.org/about/downloadreport/person/Ibanez_Jesus/2016A.

Seismic data formats
The four seismic station types (M24, Marslite, Seismic arrays, and OBSs) produced different data formats; therefore, we used different programs to handle the data. The original format of the M24 and Marslite seismic stations was binary type, but data were stored as SAC files following recovery from each station.
Since the digitizer used with the M24 stations introduces a delay in the conversion from analog to digital signals, the start time of these data is always −0.2288 s before the nominal start time of the file. This was taken in account during conversion of the data to SAC format. Data from seismic array instruments were in.ARR format and were easily converted and managed using the common seismic software SEISAN, which was also used to visualise records and plot some spectral analyses. The OBS systems stored raw data in binary format organised into a file of 8,192 records. The data in each of these 1 Mb records follow a 28-byte header that contains information on the instrument, timing, and sample rate. The data were easily converted to SEGY format after retrieval. In summary, all seismic data in the database were converted to SAC and SEGY formats. The SEGY format was selected because it permits the gathering of traces in time windows after each shot and in the display of the record section. This facilitates the comparison (and identification) of waveforms. Moreover, because only 120 s of data are stored in each file, this format results in efficient manipulation and display of the data (Fig. 7).

Seismic data quality
The quality of data registered during the seismic experiment was very good, with low noise recorded by both on-land and ocean bottom seismometers. Owing to the nature of the emplacement sites, cultural noise was completely absent. Record quality decreased for on-land stations during high winds and for OBS records when sea waves were high 23 . The clearest signals in seismic traces were those produced by shots (Fig. 8), although the records also included a variety of seismic signals that could be classified according to their shape, magnitude, and frequency content. For example, possible natural signals included ice-quakes, local and regional earthquakes, and tremor/long-period events. The air gun signals have P-wave arrivals from within a 40 km range. We clearly recognised crustal phases (Pg), except when they were masked by water-wave phases. These waves are easily identified owing to their travel-times (assuming a velocity of 1.5 km/s and that the shot-receiver distance is known) and to their larger amplitudes (Fig. 9). At greater distances, depending on crustal thickness, first arrivals associated with mantle refraction phases (Pn) were also identified. Interestingly, the waveforms of crustal phases with similar ray paths showed strong variability in their shapes, which was probably due to closely spaced attenuating heterogeneities in the medium (Fig. 9).

Technical Validation
As discussed, many of the recorded data have been used to produce derived products that have been described in published manuscripts. The database contains both raw and derived data. All of the derived data are associated with seismic analyses, including P-wave first arrivals, the velocity tomographic model, 2-D scattering and an intrinsic attenuation model, and a 3-D coda-wave normalisation attenuation model. The following provides a brief description of each derived data product.

Seismic velocity tomography
Seismic velocity tomography was the first significant output of the TOMO-DEC experiment 14 and was performed using a well-tested method for inverting P-wave travel times 39,40 . This method applies the LSQR inversion algorithm and uses the shortest ray tracing time 40 . To obtain the 3-D tomography of Deception Island, we followed a two-step approach using two grid configurations. Firstly, we applied the method to the study of a larger region, encompassing Deception Island and surroundings. The volume, which was geometrically represented in an x, y, z Cartesian system, was 53 × 52 × 12 km wide and centred in the middle of Port Foster (62°58'S, 60°40'W). It was parameterised using 0.25-km grid-node spacing for the ray tracing and 0.5-km grid-node spacing for the perturbational grid. Secondly, we focused on Port Foster and reduced the studied area to a sub-region of 12 × 14 × 7 km, centred in the middle of the bay (62°57.2'S, 60°37.2'W). We increased the density of nodes in the parameterisation grids by using a 0.1-km grid for the ray path tracing and a 0.2 km grid for the velocity perturbation.

Intrinsic attenuation model and 2-D scattering
We obtained 2-D regional maps of inverse intrinsic (Qi − 1 ), scattering (Qs − 1), and total (Qt − 1 ) quality factors for the volcanic environment of Deception Island 16 on the basis of diffusion approximation. These maps were obtained using the above described active seismic database. Attenuation parameters were estimated by fitting observed energy envelopes to the diffusion model (Fig. 10). Then, the spatial distribution of the attenuation parameters was obtained by averaging all single source-receiver couples in 1 × 1 km cells. Finally, using a back-projection method based on a Gaussian-type weighting function, 2-D intrinsic-and scattering-Q images were produced. The depth for which information is obtained can be approximately estimated assuming depth equal to the minor axis 16 ; in this case, 6 km on average.

Normalisation attenuation model using 3-D coda-waves
We obtained a high-resolution 3-D P-wave attenuation tomography model using the coda normalisation method 17 , where the energy ratios of 20,293 high-quality waveforms were analysed in a single-step inversion 41 . We applied the Thurber-modified ray-bending approach 42 in the 3-D velocity model described above. Owing to observational data associated with incoherent estimates of ray paths, the final model is restricted to depths of 1-4 km. The inverted area is a region of 20 × 20 km centred on Deception Island. The seismic tomography was obtained using cells of 2 × 2 km horizontal resolution and 1 km depth resolution with a maximum depth of 4 km b.s.l.

Usage Notes
The TOMO-DEC experiment employed a variety of instruments and produced a variety of data types. To efficiently manage experimental information and data, the folder structure for the different metadata records deposited in the AADC is divided first according to the nature of data, and then within each folder according to instrument or derived product. The root folder, which is named TOMODEC-KNOWAVES, contains six sub-folders: SEISMIC, MAGNETIC, GRAVIMETRIC, BATHYMETRY, UTM, and MODELS (see Table 7 for complete information). The UTM folder contains the original report provided by the science officer of the O/V 'Hesperides', written in Spanish and containing a description of all onboard instruments, the details of different tests, calibration, and other useful notes associated with general navigation actions.
An excel file is provided with full information about seismometers, including type of system (arrays are named MalIAG), name of the unit, channels used, component of the seismometers (3-C or V), latitude, longitude, UTM position East and North, elevation above sea level, start and end of the recording period, and shooting phases recorded. This excel file is named stations_data_05022006.xls.
For the derived data (MODELS folder), seismic velocity tomography results for 'Dense grid' and 'Small grid' models are located in the folder 'V models'. The associated P-wave onset database is located in the folder 'P travel time'. Seismic attenuation results, including the bi-dimensional model separating intrinsic and scattering Q ('Qi-Qs') and the 3-D total Q model ('TotalQ-3D') are located in the folder 'Q models'.

Author Contributions
J.M.I. was the coordinator of the TOMODEC experiment. He was the principal researcher for all of the Spanish projects involved in the experiment, and was a main researcher for the remaining projects. He was the advisor of A.D.-M., J.P., D.Z., J.A., and A.G.-Y., whose doctoral theses have contributed to the present database and to improving knowledge of Deception Island. He contributed to the structure and text of this manuscript, and prepared the figures. He coordinated data analysis and actively participated in the derived models and in the data storage. A.D.-M. actively participated in writing this manuscript, in the attenuation derived models, in data organisation, and in the storage of the data. J.P. actively participated in writing this manuscript, and in editing and improving the figures. Through her PhD thesis, she was responsible for all works associated with the attenuation models and provided these derived products. D.Z. actively participated in the success of the TOMODEC experiment, which formed the main part of her PhD thesis. She actively participated in all phases of the experiment, including the preliminary search for sites, the planning of the experiment, the field campaign, and the derivation of velocity seismic tomography of Deception Island. W.W. was responsible for the NSF project that provided OBSs to the experiment, and actively participated in the TOMODEC field campaign. He advised the velocity-derived models, including the 3-D seismic tomography. His support during the preparation phase of the experiment, during the experiment itself, and after the experiment was crucial to the creation of this high quality database. A.B. was, together with W.W., responsible for the successful use of OBSs, participating in their deployment and data quality control. He participated in on-land activities during the experiment. He was actively involved in the derivation of structural velocity models of Deception Island. He was responsible for bathymetric data and prepared the high-resolution sea floor map of the region, which is necessary for any analysis using this database. He has performed a deep revision of the English version of this manuscript. J.A. was the co-advisor of the doctoral thesis of D.Z. He actively participated in the planning and development of the TOMODEC experiment and in the velocity derived models, including the determination of the first arrival, quality control, and the publications of