Hydrothermal vent waters at 13°N on the East Pacific Rise: isotopic composition and gas concentration

Gas concentrations and isotopic compositions of water have been measured in hydrothermal waters from 13°N on the East Pacific Rise. In the most Mg-depleted samples ( 5 × 10⁻³ moles/kg) the gas concentrations are: 3–4.5 × 10⁻⁵ cm³ STP/kg helium, 0.62–1.24 cm³ STP/kg CH₄, 10.80–16.71 × 10⁻³ moles/kg CO₂. The samples contain large quantities (95–126 cm³/kg) of H₂ and some carbon monoxide (0.26–0.36 cm³/kg) which result from reaction with the titanium sampling bottles. δ¹³C in methane and CO₂ (−16.6 to −19.5 and −4.1 to −5.5 respectively) indicate temperatures between 475 and 550°C, whereas δ¹³C_CO is compatible with formation by reduction of CO₂ on Ti at 350°C close to the sampling temperature.³He/⁴He are very homogeneous at (7.5 ± 0.1)R_A(³He/⁴He = 1.0 × 10⁻⁵) and very similar to already published data as well as CH₄/³He ratios between 1.4 and 2.1 × 10⁶.¹⁸O and D in water show enrichments from 0.39 to 0.69‰ and from 0.62 to 1.49‰ respectively. These values correspond to W/R ratios of 0.4–7. The distinct¹⁸O enrichments indicate that the isotopic composition of the oceans is not completely buffered by the hydrothermal circulations. The³He-enthalpy relationship is discussed in terms of both hydrothermal heat flux and³He mantle flux.     

Dissolved carbon dioxide in basaltic glasses: concentrations and speciation

Carbon dioxide dissolved in both synthetic Ca±Mg-bearing silicate glasses and natural basaltic glasses has been characterized using infrared spectroscopy. CO₂ is inferred to be dissolved in these glasses as distorted Ca or Mg carbonate ionic complexes that result in unique infrared absorption bands at 1515 cm⁻¹ and 1435 cm⁻¹. This speciation contrasts with the case of CO₂-bearing sodium aluminosilicate glasses, which contain both dissolved molecular CO₂ and dissolved Na-carbonate ionic-complexes. The difference in speciation in Ca±Mg-bearing melts may result in part from a higher activity of oxygens that react with CO₂ molecules to produce carbonate.Dissolved CO₂ contents of natural basaltic glasses can be determined from the intensities of the carbonate absorption bands at 1515 cm⁻¹ and 1435 cm⁻¹. The uncertainty of the method is estimated to be ± 15% of the amount present. The infrared technique is a powerful tool for the measurement of dissolved CO₂ contents in natural basaltic glasses since it is non-destructive, can be aimed at regions of glass a few tens of microns in size, and can discriminate between dissolved carbonate and carbon present as carbonate alteration, contained in fluid inclusions, or adsorbed on the glass.A set of submarine basaltic glasses dredged from a variety of locations contain 0–400 ppm dissolved CO₂, measured using the infrared technique. These concentrations are lower than most previous reports for similar basaltic glasses. No general relationship is observed between dissolved CO₂ content and depth of magmatic eruption, although some correlation might be present in restricted geographic locales.     

Resonance in the gas discharge of the bocca nuova, etna (Italy), 1968–1969

During 1968–1969 the “Bocca Nuova”, a small vent, 7 m across, located in the northern part of Mt. Etna's summit crater showed a resonating-pipe phenomenon triggered by explosions at its bottom. The emitted gases were mixed with air. Temperatures, velocities, dynamic pressures and CO₂ and N₂ contents, were recorded.     

Quantification of the gas mass emitted during single explosions on Stromboli with the SO2 imaging camera

We performed measurements using an SO₂ imaging camera of the SO₂ gas mass emitted during five discrete explosive events on Stromboli volcano on 3 October 2006. The SO₂ gas mass released during discrete explosions was 15–40 kg per explosion, producing 3–8% of the total daily SO₂ gas emission, demonstrating that in terms of gas flux Strombolian explosions are a second-order phenomenon compared with quiescent degassing. Using the typical gas composition measured with OP-FTIR allows us to determine the total gas mass released during an explosion as 360–960 kg with a volume of 1500–4100 m³ at 1 bar. At the probable source pressure of gas slug formation of 75 MPa this gas amount would occupy a volume equivalent to a sphere with a radius of 0.8–1 m, comparable with estimates of Stromboli's conduit geometry.     

The relationship between vertical eddy diffusion and buoyancy gradient in the deep sea

Vertical eddy diffusivities (K_v's) have been estimated at fourteen widely separated locations from fourteen²²²Rn profiles and two²²⁸Ra profiles measured near the ocean floor as part of the Atlantic and Pacific GEOSECS programs. They show an inverse proportionality to the local buoyancy gradient [(g/?)(??_pot/?z)] calculated from hydrographic measurements. The negative of the constant of proportionality is the buoyancy flux [?K_v(g/?)(??_pot/?z)] which has a mean of ?4 × 10^?6 cm²/sec³. Our results suggest that the buoyancy flux varies very little near the ocean floor. K_v's for the interior of the deep Pacific calculated from the relationship K_v = (4 × 10^?6cm²/sec³)/[(g/?)(??_pot/?z)] agree well with published estimates. K_v's calculated for the pycnocline are one to two orders of magnitude smaller than upper limits estimated from tritium and⁷Be distributions.Heat fluxes calculated with the model K_v's obtained from the²²²Rn profiles average 31 μcal cm^?2 sec^?1 in the Atlantic Ocean and 8 μcal cm^?2 sec^?1 in the Pacific Ocean.     

On the heat flux related to stretching in the NW-Mediterranean continental margins

Summary The surface thermal flux of the continental margins of the northwestern Mediterranean Sea is interpreted on the basis of a 1-D instantaneous pure shear stretching model of the lithosphere in terms of three components: the background heat flowing out from the asthenosphere (38 mW m^–2), the transient contribution depending on the rift age and extension amount (35 mW m^–2 at the most), and the contribution due to the radiogenic elements of the lithosphere. The radiogenic component is estimated at the continental margins of the Ligurian-Provençal basin and Valencia trough, and in the surrounding mainland areas by means of available data of surface heat generation from Variscan Corsica, Maures-Estérel and the Central Massif along with a geophysical-petrological relationship between heat production and seismic velocity. The lithosphere radiogenic heat contribution q_l decreases with the thinning factor according to the exponential law: q_l() = a exp(-b), in which factor b is greater for that part of the lithosphere below the uppermost 10 km. Considering also the heat generated by radioactive isotopes in sediments, the stable Variscan lithosphere produces an average thermal flux of 30 mW m^–2 which decreases by about one half where the lithosphere is thinned by one third. Although the surface heat generation is 2·1 – 3·3 µW m^–3 in the Maures-Estérel massif — excepting small outcrops of dioritic rocks with lower heat production — and 1·8 µW m^–3 for most of Corsica, the radiogenic heating within the lithosphere for such areas is nearly the same and does not explain the higher heat flux of the Corsica margin. This asymmetric thermal pattern with surface heat flux which is 10 – 15 mW m^–2 higher than predictions is probably of upper mantle origin, or can be ascribed to penetrative magmatism.     

A crust-mantle model for the NORSAR area

Summary Elastic waves from explosions were recorded at NORSAR and at a number of field stations, and the data were used for determining a crust-mantle model under the array. The number of explosions was eleven distributed on seven shot points. The total number of recording points was fifty-one, and the interpretation was based on 350 individual records.The velocities obtained for the crustal phases were 6.2, 6.6 and 8.2 km/sec for theP _g ,P _g andP _n waves respectively. A deep crustal phase with a velocity of about 7.4 km/sec was observed. The mean depths to the discontinuities within the crust were determined to be 17 and 26 km. The depth to Moho varied greatly across the array from 31.5 km in the central part to 38 km under the C-ring. The maximum dip observed for the Moho was 12^o.Contribution No. 57 to Norwegian Geotraverse Project.     

Verification of land surface evapotranspiration estimation from remote sensing spatial contextual information

Estimation of evapotranspiration (ET) is of great significance in modeling the water and energy interactions between land and atmosphere. Negative correlation of surface temperature (T_s) versus vegetation index (VI) from remote sensing data provides diagnosis on the spatial pattern of surface soil moisture and ET. This study further examined the applicability of T_s–VI triangle method with a newly developed edges determination technique in estimating regional evaporative fraction (EF) and ET at MODIS pixel scale through comparison with large aperture scintillometer (LAS) and high‐level eddy covariance measurements collected at Changwu agro‐ecological experiment station from late June to late October, 2009. An algorithm with merely land and atmosphere products from MODIS onboard Terra satellite was used to estimate the surface net radiation (R_n) and soil heat flux. In most cases, the estimated instantaneous R_n was in good agreement with surface measurement with slight overestimation by 12 W/m². Validation results from LAS measurement showed that the root mean square error is 0.097 for instantaneous EF, 48 W/m² for instantaneous sensible heat flux, and 30 W/m² for daily latent heat flux. This paper successfully presents a miniature of the overall capability of T_s–VI triangle in estimating regional EF and ET from limited number of data. For a thorough interpretation, further comprehensive investigation needs to be done with more integration of remote sensing data and in‐situ surface measurements. Copyright © 2011 John Wiley & Sons, Ltd.     

Gas hazard assessment in a densely inhabited area of Colli Albani Volcano (Cava dei Selci, Roma)

The northwestern flank of the Colli Albani, a Quaternary volcanic complex near Rome, is characterised by high pCO₂ values and Rn activities in the groundwater and by the presence of zones with strong emission of gas from the soil. The most significant of these zones is Cava dei Selci where many houses are located very near to the gas emission site. The emitted gas consists mainly of CO₂ (up to 98 vol%) with an appreciable content of H₂S (0.8–2%). The He and C isotopic composition indicates, as for all fluids associated with the Quaternary Roman and Tuscany volcanic provinces, the presence of an upper mantle component contaminated by crustal fluids associated with subducted sediments and carbonates. An advective CO₂ flux of 37 tons/day has been estimated from the gas bubbles rising to the surface in a small drainage ditch and through a stagnant water pool, present in the rainy season in a topographically low central part of the area. A CO₂ soil flux survey with an accumulation chamber, carried out in February–March 2000 over a 12 000 m² surface with 242 measurement points, gave a total (mostly conductive) flux of 61 tons/day. CO₂ soil flux values vary by four orders of magnitude over a 160-m distance and by one order of magnitude over several metres. A fixed network of 114 points over 6350 m² has been installed in order to investigate temporal flux variations. Six surveys carried out from May 2000 to June 2001 have shown large variations of the total CO₂ soil flux (8–25 tons/day). The strong emission of CO₂ and H₂S, which are gases denser than air, produces dangerous accumulations in low areas which have caused a series of lethal accidents to animals and one to a man. The gas hazard near the houses has been assessed by continuously monitoring the CO₂ and H₂S concentration in the air at 75 cm from the ground by means of two automatic stations. Certain environmental parameters (wind direction and speed; atm P, T, humidity and rainfall) were also continuously recorded. At both stations, H₂S and CO₂ exceeded by several times the recommended concentration thresholds. The highest CO₂ and H₂S values were recorded always with wind speeds less than 1.5 m/s, mostly in the night hours. Our results indicate that there is a severe gas hazard for people living near the gas emission site of Cava dei Selci, and appropriate precautionary and prevention measures have been recommended both to residents and local authorities.     

Surface energy balance at a tropical station

Summary Studies of various fluxes, namely net radiation, soil heat, sensible heat and latent heat observed at a tropical station are presented in this paper. The time variation of these fluxes are examined in relation to various meteorological parameters and atmospheric conditions. The turbulent transfer coefficients have been evaluated to examine the applicability of the classical theory or the non-equivalence theory for eddy transport in the lower layers of the atmosphere. The energy balance at a tropical station is evaluated. It is found over year there is a net surplus of 94,000 ly. A detailed discussion of the disposal of this energy by various consuming processes is given.Symbols and notation All the quantities represented by symbols in the text of the paper are defined below C _p specific heat at constant pressure in cal. g^–1 °K - E evaporation in g cm^–2 hr^–1 - E ^* evaporative heat flux cal. cm^–2 sec^–1 (in Eq. 10) - e vapour pressure in millibars - e _z ,e _2z vapour pressure at heightsz and 2z - g acceleration due to gravity - H sensible heat flux cal cm^–2 sec^–1 (in Eq. 12) - K _M ,K _H ,K _W coefficients of eddy diffusivities of momentum, heat, and water vapour respectively in cm² sec^–1 - k von Karman' constant=0.4 - L Monin-Obukhov length (according to Monin and Obukhov [53] the structure of the turbulent boundary layer is determined by the non-dimensional variableZ/L whereL is defined byL=–(u _* ³ C _p T)/(kgH) - ly langleys - Q _c Q—sensible heat flux in langleys (in Eqs. 3 and 4) - Q _e E—latent heat flux in langleys (in Eqs. 3 and 4) - Q _s S—soil heat flux in langleys (in Eqs. 3 and 4) - Q _i Q _c +Q _e +Q _s whenK _M K _H K _W , (in Eq. 6) - Q' _i Q' _c +Q' _e +Q _s whenK _M =K _H =K _W (in Eq. 7) - qq mean specific humidity g kg^–1     

Correlated seasonal and climatic variations of trace constituents in the stratosphere

Summary Spectrometric experiments performed, in November 1976, within the framework of the Latitude Survey Mission on board the NASA Convair 990 from Ames Research Center are briefly deseribed. The results presented concern odd nitrogen molecules, HCl and water vapor. In terms of vertical column density, HNO₃ is predominant over NO+NO₂ at all latitudes higher than 40 degrees. A seasonal variation of NO₂ abundance is observed, with larger values in the summer hemisphere at high latitude. The mean zenith column density of HCl above 11 km is 1.5×10¹⁵ mol.cm^–2, with no evidence for any seasonal or climatic variation. Local number densities as high as 1.4×10¹⁰ mol.cm^–3 for HNO₃ and 5.4×10¹⁴ mol.cm^–3 for water vapor have been measured during the same flight near 11 km.     

Evaluation of CF2Cl2 spectral parameters for atmospheric sensing

Laboratory infrared absorption spectral measurements of CF₂Cl₂ in the 10.8 m band have been performed. Doppler limited resolution has been achieved using semiconductor tunable lasers. The line strength and the self and pressure broadening rates have been measured for the 918.562 and 918.621 cm^–1 lines. The results are used for considerations on the application of high resolution spectroscopy to the measurement of atmospheric CF₂Cl₂.     

A detailed energy budget analysis of river supercooling and the importance of accurately quantifying net radiation to predict ice formation

River supercooling and ice formation is a regular occurrence throughout the winter in northern countries. The resulting frazil ice production can obstruct the flow through intakes along the river, causing major problems for hydropower and water treatment facilities, among others. Therefore, river ice modellers attempt to calculate the river energy budget and predict when supercooling will occur in order to anticipate and mitigate the effects of potential intake blockages. Despite this, very few energy budget studies have taken place during freeze-up, and none have specifically analysed individual supercooling events. To improve our understanding of the freeze-up energy budget detailed measurements of air temperature, relative humidity, barometric pressure, wind speed and direction, short- and longwave radiation, and water temperature were made on the Dauphin River in Manitoba. During the river freeze-up period of late October to early November 2019, a total of six supercooling events were recorded. Analysis of the energy budget throughout the supercooling period revealed that the most significant heat source was net shortwave radiation, reaching up to 298 W/m², while the most significant heat loss was net longwave radiation, accounting for losses of up to 135 W/m². Longwave radiation was also the most significant heat flux overall during the individual supercooling events, accounting for up to 84% of the total heat flux irrespective of flux direction, highlighting the importance of properly quantifying this flux during energy budget calculations. Five different sensible (Q_h) and latent (Q_e) heat flux calculations were also compared, using the bulk aerodynamic method as the baseline. It was found that the Priestley and Taylor method most-closely matched the bulk aerodynamic method on a daily timescale with an average offset of 8.5 W/m² for Q_h and 10.1 W/m² for Q_e, while a Dalton-type equation provided by Webb and Zhang was the most similar on a sub-daily timescale with average offsets of 20.0 and 14.7 W/m² for Q_h and Q_e, respectively.     

The influence of plants on atmospheric carbon monoxide and dinitrogen oxide

It is shown by laboratory experiments and extensivein-situ measurements that higher plants (Vicia faba, Platanus acerifolia, Fagus silvatica, andPinus silvestris) produce carbon monoxide. The measurements were carried out under natural conditions with respect to the concentrations of O₂ and H₂O, and temperature. The CO₂- and CO-mixing ratios were varied in the ranges 350 to 530 ppm and 3–270 ppb, respectively. The CO-production rates were found to be light dependent with an average value per cm² of leaf area of 3×10^–13 g/sec for a radiation intensity of 5×10⁴ erg/cm² sec. The production rates are independent of the CO₂- and CO-mixing ratios employed in the test atmosphere. Considering the production rate of 3×10^–13 g/cm² sec to be representative for global conditions the total CO-production by plants is estimated to be 0.5–1.0×10¹⁴ g/year. In contrast to carbon monoxide atmospheric dinitrogen oxide is not influenced by plants in the same manner.     

Excess222Rn and the benthic boundary layer in the western and southern Indian Ocean

We present 9 bottom²²²Rn profiles measured from the western and southern Indian Ocean during the 1977–1978 GEOSECS expedition. These profiles can be grouped into three cypes: one-layer, two-layer, and irregular types. The one-layer profiles with quasi-exponential distributions allow one to estimate the apparent vertical eddy diffusivity,K_v, with a simple model. The two-layer profiles show that there is a benthic boundary layer of the order of 50–100 m in which the excess²²²Rn distribution shows a vertical gradient much smaller than that of the layer immediately above. Within the boundary layer, the STD potential temperature (θ) and density(σ₄) profiles are practically constant, and theK_v values are of the order of 1000 cm²/s. The STD profiles for the water column above the boundary layer show gradients of increasing stability, and theK_v values are of the order of 100 cm²/s. Modeling of the Rn data in the water column above the boundary layer indicates that there is a transition layer which effectively reduces the penetration of excess Rn from the benthic boundary layer into the upper layer.Sarmiento et al. [10] have shown that the buoyancy gradient or stability is inversely correlated with the apparent vertical eddy diffusivity, and the resulting buoyancy flux is fairly uniform, ranging from 1 to 14 × 10^?6 cm²/s³ in the Atlantic and Pacific Oceans. However, Sarmiento et al. [11] show that a much higher buoyancy flux is associated with an intensified flow of the bottom water through a passage. In the Indian Ocean basins, we have found that the buoyancy flux has a comparable range (3–14 × 10^?6 cm²/s³), except for a couple of stations where both stability and apparent vertical diffusivity are higher, resulting in a much higher buoyancy flux, probably indicative of rapid bottom water flow.     

A remote sensing model of CO2 flux for wheat and studying of regional distribution

A model for calculating CO₂ flux in the wheat field and an algorithm for estimating CO₂ flux in the mejonal scale were presented using the remote sensing data and supplementary micpo-met~orological data. First of all a-longertenn measurement wae carried out during winter wheat growing period in Yucheng Experimental Station udng the spectmradiometer system, the thermal infrared radiometer system, the Bowen-ratio device as well as the eddy-correlation device. Two kinds of issues concerning remote sensing and CO₂ flux can be obtained. Based on the obeervations a remote sensing model was estabilished. Then when the NOAA-AVHRR passed over the experimental area simultaneous measurements were carried out with the satellites. A regional distribution image for CO₂ flux over wheat canopy in North China (500×500 km²) was made using the supplementary ground data and NOAA-AVHRR remote sensing data which was calibrated by the synchronous observation. The sources and sinks for CO₂ fluxes in the region can be seen obviously. Project supported by the National Natural Science Foundation of China (Grant Nos. 49671058, 49890330)     

基于同步观测信息利用的高寒湖泊湍流通量数据插补方法

源区划分和质量过滤提高湖面涡动相关通量数据可靠性的同时,却降低了通量时间序列的连续性.为此,本文基于TensorFlow机器学习框架构建了一种超宽人工神经网络(ANN)模型.在选择输入ANN模型的特征变量信息时,我们采取了尽可能获取湍流输送过程中热力、动力学同步观测背景强迫信息的原则.通过ANN模型模拟通量的插补,本文实现了通量时间序列连续性的优化,插补后的羊卓雍错湖面通量数据的时间覆盖率从不足0.40提升至超过0.98.基于10次折叠交叉验证的ANN模型通量模拟性能检验则表明,各个检验组之间ANN模型的模拟误差波动较小,这显示出了较好的稳健性.具体地讲,感热通量、潜热通量和水汽通量原始观测平均值分别约为18.8 W/m~2、81.5 W/m~2和1.84 mmol/(s·m~2),10组交叉验证的插补感热通量、潜热通量和水汽通量平均绝对误差分别为5.4 W/m~2、15.7 W/m~2和0.35 mmol/(s·m~2).这表明本文所探索的ANN建模结构和同步观测变量筛选原则可更充分地利用观测点局地同步观测信息估算通量强度,有效地优化湍流通量数据的时间连续性,从而提升通量数据的可分析性.     

Solar activity dependence of ionospheric electron content and slab thickness using different solar indices

Ionospheric electron content (IEC) and slab thickness () data for the period 1977 to 1980 from Lunping (23.03°N; 121.90°E subionospheric) have been examined for their solar activity dependence. Local noontime monthly means as well as values for the 5 QQ days in a month have been examined separately with different solar indices, namely: solar EUV flux (170–190 Å),S _{10.7 cm} flux and sun spot number (SSN) on a seasonal basis. Both IEC and parameters exhibit better correlation with solar EUV andS _{10.7 cm} fluxes than with SSN for all seasons. IEC increases linearly with both EUV andS _{10.7 cm} flux whereas with SSN it shows a distinct nonlinear relationship during all seasons in both monthly mean and 5 QQ days' values. This study indicates that for correlating and predicting the variations (especially the medium term) in the ionospheric parameters, both EUV andS _{10.7 cm} fluxes have an advantage over SSN.     

Effects of soil compaction,rain exposure and their interaction on soil carbon dioxide emission

Soils release more carbon, primarily as carbon dioxide (CO₂), per annum than current global anthropogenic emissions. Soils emit CO₂ through mineralization and decomposition of organic matter and respiration of roots and soil organisms. Given this, the evaluation of the effects of abiotic factors on microbial activity is of major importance when considering the mitigation of greenhouse gases emissions. Previous studies demonstrate that soil CO₂ emission is significantly affected by temperature and soil water content. A limited number of studies have illustrated the importance of bulk density and soil surface characteristics as a result of exposure to rain on CO₂ emission, however, none examine their relative importance. Therefore, this study investigated the effects of soil compaction and exposure of the soil surface to rainfall and their interaction on CO₂ release. We conducted a factorial laboratory experiment with three soil types after sieving (clay, silt and sand soil), three different bulk densities (1·1 g cm^–3, 1·3 g cm^–3, 1·5 g cm^–3) and three different exposures to rainfall (no rain, 30 minutes and 90 minutes of rainfall). The results demonstrated CO₂ release varied significantly with bulk density, exposure to rain and time. The relationship between rain exposure and CO₂ is positive: CO₂ emission was 53% and 42% greater for the 90 minutes and 30 minutes rainfall exposure, respectively, compared to those not exposed to rain. Bulk density exhibited a negative relationship with CO₂ emission: soil compacted to a bulk density of 1·1 g cm^–3 emitted 32% more CO₂ than soil compacted to 1·5 g cm^–3. Furthermore we found that the magnitude of CO₂ effluxes depended on the interaction of these two abiotic factors. Given these results, understanding the influence of soil compaction and raindrop impact on CO₂ emission could lead to modified soil management practices which promote carbon sequestration. Copyright © 2012 John Wiley & Sons, Ltd.     

Hydrogen emissions from Erebus volcano, Antarctica

The continuous measurement of molecular hydrogen (H₂) emissions from passively degassing volcanoes has recently been made possible using a new generation of low-cost electrochemical sensors. We have used such sensors to measure H₂, along with SO₂, H₂O and CO₂, in the gas and aerosol plume emitted from the phonolite lava lake at Erebus volcano, Antarctica. The measurements were made at the crater rim between December 2010 and January 2011. Combined with measurements of the long-term SO₂ emission rate for Erebus, they indicate a characteristic H₂ flux of 0.03?kg s^–1 (2.8?Mg? day^–1). The observed H₂ content in the plume is consistent with previous estimates of redox conditions in the lava lake inferred from mineral compositions and the observed CO₂/CO ratio in the gas plume (～0.9 log units below the quartz–fayalite–magnetite buffer). These measurements suggest that H₂ does not combust at the surface of the lake, and that H₂ is kinetically inert in the gas/aerosol plume, retaining the signature of the high-temperature chemical equilibrium reached in the lava lake. We also observe a cyclical variation in the H₂/SO₂ ratio with a period of ～10?min. These cycles correspond to oscillatory patterns of surface motion of the lava lake that have been interpreted as signs of a pulsatory magma supply at the top of the magmatic conduit.