Corrigendum: Topography and Tilt at Volcanoes

In the original article, there was an oversight, resulting in the previous work of Beauducel and Carbone (2015) and Albino et al. (2011) not being cited. A correction has been made to Section 2. Influence of the Original Slope Angle, paragraph 3. Inflation or deflation of a magma reservoir is commonly modeled analytically using a point source (Mogi, 1958), where the orientation of the displacement field at the surface varies with horizontal distance from the source. The point source approximation is valid for a spherical source where the source radius is small with respect to the distance from that source. Addressing the influence of topography, Cayol and Cornet (1998) showed that displacement produced by reservoir pressure is greater at lower elevation, where the surface is closer to the source. Delving further, Ronchin et al. (2015) found that the slope angle θ is a secondary control on the displacement. Through numerical modeling, they showed that the maximum displacement is produced where the surface is perpendicular to a line between the surface and source, where the exposure to the pressure source is greatest. We refer to the slope angle with the maximum exposure to the source at any point as θmaxexp, and the deviation from this angle as θ − θmaxexp. Hence, tilt generated by a point or spherical source is influenced by the slope angle. We vary this over a 90 range around a point centered at x = 3 km, z = 3 km (Figure 1C). We calculate u and w at each point along the slope using it’s distance from the centre of a Mogi source (x = 0 km, z =−1 km)

In the original article, there was an oversight, resulting in the previous work of Beauducel and Carbone (2015) and Albino et al. (2011) not being cited.
A correction has been made to Section 2. Influence of the Original Slope Angle, paragraph 3. Inflation or deflation of a magma reservoir is commonly modeled analytically using a point source (Mogi, 1958), where the orientation of the displacement field at the surface varies with horizontal distance from the source. The point source approximation is valid for a spherical source where the source radius is small with respect to the distance from that source. Addressing the influence of topography, Cayol and Cornet (1998) showed that displacement produced by reservoir pressure is greater at lower elevation, where the surface is closer to the source. Delving further, Ronchin et al. (2015) found that the slope angle θ is a secondary control on the displacement. Through numerical modeling, they showed that the maximum displacement is produced where the surface is perpendicular to a line between the surface and source, where the exposure to the pressure source is greatest. We refer to the slope angle with the maximum exposure to the source at any point as θ maxexp , and the deviation from this angle as θ − θ maxexp . Hence, tilt generated by a point or spherical source is influenced by the slope angle. We vary this over a 90 • range around a point centered at x = 3 km, z = 3 km ( Figure 1C). We calculate u and w at each point along the slope using it's distance from the centre of a Mogi source ( where P is the pressure, α is the source radius, d is the vertical distance between each point and the source center, ν is the Poisson's ratio and E m is the Young's modulus. v is the displacement in the y direction, which is not used in this 2D example, but is considered in the 3D modeling later in section 5. The new slope angle can be computed as the derivative of the new coordinates for each point along the slope, x + u and z + w. The tilt θ is then calculated as the change in slope angle. Beauducel and Carbone (2015) suggested that tilt induced by a spherical pressure source is greater where the slope of the edifice is steeper. However, we find that the tilt θ/ θ max is greatest where the surface is 45 • from the angle of maximum exposure, i.e., where θ − θ maxexp = 45 • ( Figure 1C). In other words, no tilt is produced where the surface is parallel or perpendicular to the displacement field, and the highest tilt is generated where the surface is 45 • from this. This, unsurprisingly, is as previously shown for a purely horizontal displacement field.
A correction has been made to Section 5. Real Topography, paragraphs 1 and 2.
The 2D axisymmetric models in sections 3, 4 are useful in demonstrating how exaggerated features of simplified topography can affect the tilt. However, it is important to consider whether the effect of real topography is significant or negligible. To do this, we present results of 3D deformation modeling, using 10 m × 10 m digital elevation models (DEMs) of Tungurahua volcano, Ecuador and Soufrire Hills volcano, Montserrat (SHV), and compare how tilt varies spatially across both DEMs where topography is the only variable. We again asses the topographic effect on tilt induced by each of the sources described in section 1. The mesh must be suitably fine such that the resolution of the DEM is matched. Beauducel and Carbone (2015) investigated the influence of real topography on tilt produced by a spherical pressure source. However, since they use a spatially variant mesh with a minimum element size of 200 m, they thus do not consider the influence of local topography, which we show can significantly influence the tilt (section 4). Here, we use spatially variant triangular mesh across the entire DEM, to a distance of almost 8 km in x and y from the conduit, with a minimum element size of ∼1 m where the topography is most complex. The model is extrapolated to distance of 40 km in x and y and a depth of 50 km, to ensure that the deformation within the region that the DEM covers is not affected by the boundary constraints, described in section 3. A much coarser mesh (maximum element size ∼7 km) is used for this extrapolation. We consider the maximum tilt for θ , which is always positive and not necessarily radial to the source.
In defining the reference tilt, θ r , we followed the common practice where topography is not considered. Therefore, θ AMP is the factor by which the modeled tilt is scaled when topography is considered relative to the common practice when it is not. For tilt generated by reservoir pressure, an analytical solution was used (Mogi, 1958) (Equation 7), taking the x, y, z coordinates of each individual point across the mesh to calculate θ r . For tilt produced by conduit pressure or shear stress, in the absence of an available analytical solution, we obtain θ r from 2D axisymmetric models using a constant slope, based on the average dip angle radially away from the summit for each edifice. This mimics the approach of Albino et al. (2011), who represent topography using only this first-order approximation. The computation of this average dip angle, and a full description of these 2D axisymmetric models, is included in the Supplementary Material, Section S1.