Vertical to Horizontal UHS Ratios for Low to Medium Seismicity Regions with Deep Soil atop Deep Geological Sediments—An Example of the City of Osijek, Croatia

: The severity of vertical seismic ground motions is often factored into design regulations as a component of their horizontal counterparts. Furthermore, most design codes, including Eurocode 8, ignore the impact of local soil on vertical spectra. This paper investigates vertical pseudo-absolute acceleration spectral estimates, as well as the ratios of spectral estimates for strong motion in vertical and horizontal directions, for low to medium seismicity regions with deep local soil and deep geological sediments beneath the local soil. The case study region encompasses the city of Osijek in Croatia. New regional frequency-dependent empirical scaling equations are derived for the vertical spectra. According to these equations, for a 0.3 s spectral amplitude at deep soils atop deep geological sediments compared to the rock sites, the maximum ampliﬁcation is 1.48 times. The spectra of vertical components of various real strong motions recorded in the surrounding region are compared to the empirical vertical response spectra. The new empirical equations are used to construct a Uniform Hazard Spectra for Osijek. The ratios of vertical to horizontal Uniform Hazard Spectra are generated, examined, and compared to Eurocode 8 recommendations. All the results show that local soil and deep geology conditions have a signiﬁcant impact on vertical ground motions. The results also show that for deep soils atop deep geological strata, Eurocode 8 can underestimate the vertical to horizontal spectral ratios by a factor of three for Type 2 spectra while overestimating them by a factor of two for Type 1 spectra.


Introduction
In Croatia, as in many other countries, vertical design ground motions are calculated based on the horizontal design ground motions and do not depend on the soil conditions. In Eurocode 8 [1], which is in use in Croatia [2], vertical design Peak Ground Acceleration (hereinafter, PGA) values are determined by multiplying the horizontal PGA values spectra defined for the "rock" ground type (type A) by a constant. This scaling constant changes according to the predominant magnitudes; it is equal to 0.9 for Eurocode 8 Type 1 spectra (magnitudes larger than 5.5) and 0.45 for Type 2 spectra (magnitudes equal to or smaller than 5.5). In both cases, vertical pseudo-absolute acceleration spectra (hereinafter, PSA vert ) are the same for all ground types that are considered in Eurocode 8 [1]. In other words, the Eurocode 8 approach for vertical PGA assumes that seismic wave amplification happens only horizontally and not vertically. This concept underpins the so-called H/V methodology for calculating the natural period of soil vibrations based on microtremors, as suggested by Nogoshi and Igarashi [3] and later popularized by Nakamura [4][5][6]. However, many studies have demonstrated that this methodology does not always provide realistic popularized by Nakamura [4][5][6]. However, many studies have demonstrated tha methodology does not always provide realistic results for local soil vibration period 11]. The initial purpose of this research is to determine the amplification (o amplification) impacts of deep soil atop deep geological sediments on vertical PSA tra. The term "deep soil" will refer to the geotechnical site description on the sc several tens of meters, in this case to sites with more than 100 m of softer soil l above the first layer with an average shear wave velocity of more than 800 m/s [1 The term "deep geological sediments" will refer to geological settings on the scal few kilometers or at least hundreds of meters [14].
Vertical to horizontal PSA ratios obtained in this study will be compared to t tios defined in two different sets of empirical Ground Motion Prediction Equa (hereinafter, GMPEs). The first set was created for California [15], and these were th ly GMPEs we could find that consider the simultaneous impacts of deep geological ments and deep local soil on the vertical spectra. The second set [16] was created fo former Yugoslavia and is used to compare the vertical and horizontal spectral rat deep geological sediments only, as these GMPEs do not consider the effects of deep so In this paper, we first develop regional frequency-dependent empirical attenu equations that consider deep soil sites, then we generate probabilistic seismic h PSAvert estimates and compare them to Eurocode 8 spectra, and finally compare the ve and the horizontal PSA ratios to the ratios suggested by other authors and Eurocode In addition, seismic hazard maps for the various vibration periods and exceedance p bilities are created. Figure 1 depicts a research framework diagram for this study. This research builds on our previous work [17,18] on the severity of surface gr motion in low to medium seismicity regions with deep soils atop deeper geologica ta. We chose the city of Osijek, Croatia for the case study region since it is in the l medium seismicity region of the Pannonian basin (see Figure 2 in our accompa paper on the horizontal spectra [18]). Osijek is also an example of an area where geological deposits are found beneath deep soil sites. In our study in [17], we found vertical PGA is de-amplified by 10% at deep soil sites compared to rock soil site ~15% at deep geological sediments compared to geological rocks. This study will focus on PSAvert. This research builds on our previous work [17,18] on the severity of surface ground motion in low to medium seismicity regions with deep soils atop deeper geological strata. We chose the city of Osijek, Croatia for the case study region since it is in the low to medium seismicity region of the Pannonian basin (see Figure 2 in our accompanying paper on the horizontal spectra [18]). Osijek is also an example of an area where deep geological deposits are found beneath deep soil sites. In our study in [17], we found that vertical PGA is de-amplified by 10% at deep soil sites compared to rock soil sites and~15% at deep geological sediments compared to geological rocks. This study will now focus on PSA vert . Appl. Sci. 2021, 11,  Furthermore, in Eurocode 8 [1], the maximum spectral amplitudes are equal to the product of the vertical PGA and a constant value of 3.0. Another purpose of this study is to examine if this spectral acceleration factor provides realistic PSAvert estimates. We will also examine the ratios of the vertical to horizontal spectra. Table 1 presents the formulas for computing vertical to horizontal PSA ratios for Ground Type C using the spectra of Types 1 and 2 of Eurocode 8 [1]. As can be seen, the vertical to horizontal PSA ratio changes non-linearly with the vibration period, T, with a general tendency to become lower with T. For short periods, the so-called damping correction factor, η, also influences this ratio. Furthermore, in Eurocode 8 [1], the maximum spectral amplitudes are equal to the product of the vertical PGA and a constant value of 3.0. Another purpose of this study is to examine if this spectral acceleration factor provides realistic PSA vert estimates. We will also examine the ratios of the vertical to horizontal spectra. Table 1 presents the formulas for computing vertical to horizontal PSA ratios for Ground Type C using the spectra of Types 1 and 2 of Eurocode 8 [1]. As can be seen, the vertical to horizontal PSA ratio changes non-linearly with the vibration period, T, with a general tendency to become lower with T. For short periods, the so-called damping correction factor, η, also influences this ratio.  The GMPEs that will be generated in this study will account for the effects of deep soils and deep geological sediments on the vertical spectra, which is a unique feature of this research. Most other GMPEs for vertical spectra only address local soil conditions via the V S,30 parameter [19][20][21][22][23]. However, strong motion predictions based on scaling equations that account for local soil conditions (up to depths of 100 m and more) as well as deep geology, have been proven to be very accurate when compared to recorded ground motions and observed intensities in the north-western Balkans [24,25]. In addition, a series of recent regional seismic microzonation studies [26][27][28][29][30][31] found that changes in deep geological site conditions can considerably increase seismic hazard estimates when compared to hazard estimates based simply on local soil effects.

Empirical Frequency-Dependent Scaling Equations for Vertical Ground Motion
We start with the development of new empirical formulae for estimating vertical PSA values that can be used for deep soil sites atop deep geological deposits. The following mathematical form will be used to express the attenuation equations: where PSA vert (T) denotes the amplitudes of the vertical pseudo-absolute acceleration spectrum for a 5 percent viscous damping ratio, T is the vibration period of an undamped single-degree-of-freedom oscillator, M is the earthquake magnitude, and S L1 and S L2 (the local soil categorical variables) are 0 for the rock soil sites (further denoted as S L = 0), S L1 = 1 and S L2 = 0 for stiff soil (S L = 1), and S L1 = 0 and S L2 = 1 for deep soil (S L = 2). It should be noted that in this equation, local soil (i.e., shallow geology) is classified by Seed et al. [12,13], while deep geology is separated into three categories using Trifunac and Brady's [14] classification: (i) geological rock (further denoted as S = 2) with the corresponding categorical variables S G1 and S G2 equal to 0; (ii) deep geological sediments (S = 0) with S G1 = 0 and S G2 = 1; and (iii) intermediate (or complex) geological site surroundings (S = 1) for which S G1 = 1 and S G2 = 0. The unknown parameters resulting from the regression analysis are the coefficients c 1 through c 6 . Table 2 displays the coefficients c 1 to c 6 that can be used to predict the PSA vert values for various vibration periods. The regression analyses were conducted on the assumption that the studied data had a lognormal distribution. The standard deviation of the logarithm of PSA vert (T) is denoted by σ log (T), where p = 0 for median estimations and p = 1 for median plus one standard deviation estimations. Finally, R stands for hypocentral or epicentral distance, and R 0 (T) was iteratively calculated for each PSA(T) to minimize the standard deviation in the related prediction equation. The employed strong motion database comprises 218 vertical components of accelerograms obtained in the north-western Balkans from 112 earthquakes with magnitudes between 3 and 6.8. The majority (209) occurred between 1976 and 1987, with the remainder (9) occurring between 2010 and 2011. 203 are from the EQINFOS database [32] and 6 are from the ISESD database for the period 1976-1987 [33,34]. The Seismological Survey of Serbia used an accelerograph network to collect remaining 9 data in 2010 [35]. Because the majority of the accelerograms had a low signal-to-noise ratio for T > 2 s [16,36], most of the records were high-pass filtered with a relatively large corner frequency, and the results for T > 2 s are unreliable. Hence, we calculated scaling coefficients up to a 2 s vibration period. Bulajić et al. [17,24] provide more information on the used data.
The MATLAB ® function "regress" was used to perform two-phase multiple linear regression analyses, the first of which consisted of fitting the Equation (1) model to a dataset that only included 203 vertical acceleration components from the EQINFOS database [32], for which Trifunac et al. [37] established the soil categorization of the recording locations. Because there were no deep soil sites among the 203 components, we fitted Equation (1) without the coefficient c 5 . To lower the root mean squared error, the R 0 values were 5 of 16 adjusted iteratively. We ran a secondary analysis using only data gathered at epicentral distances of less than 30 km because the major portion of the data was obtained across short distances. We calibrated the resultant prediction model in the second phase by calculating the coefficient c 5 based on the additional 15 vertical acceleration components only observed at deep soil sites (again via linear regression).  Table 2 shows the calculated scaling coefficients for Equation (1) with R as the epicentral distance, based on the data recorded at all distances. Table 3 shows the calculated coefficients for Equation (1) with R as the hypocentral distance.    Figure 2 shows the empirical predictions at deep soil sites (sL = 2) for three different deep geology types. The solid lines reflect median empirical predictions based on the data recorded at all epicentral distances, whereas the dashed lines reflect predictions based on distances less than 30 km.
To calculate the exact differences between the attenuation curves for different site conditions, we can apply the coefficients corresponding to the categorical variables SL and SG. If the coefficients from Table 2 are utilized, the spectral amplitude for T = 0.05 s will be 1/10 −0.021 = 1.05 times larger at deep geological rock (s = 2) than at deep geological sediments (s = 0), as seen in the top left-hand plot in Figure 2. This is most likely because short-period waves travel through more compact rocks more easily. The opposite occurs for larger periods; the spectral amplitude for T = 0.50 s will be 1/10 0.147 = 1.40 times larger at deep geological sediments than at rock, as seen in the bottom right-hand plot in the same figure.
The exact differences between the PSA attenuation curves for different site conditions can be also calculated by using the coefficients c4, c5, c6, or c7. For example, the spectral amplitude for T = 0.05 s will be 1/10 0.052 = 1.13 times larger at deep soil sites (sL = 2) than at rock soil sites (sL = 0), while the amplitude for T = 0.30 s will be 1/10 0.088 = 1.23 times larger. The coefficient c5 will have positive values for all vibration periods up to 0.5 s. This indicates that deep soil sites produce amplification of vertical seismic waves for all vibration periods smaller or equal to 0.5. This is somewhat surprising because our investigation of the horizontal spectra revealed that short-period waves at the deep soil sites are de-amplified [18]. Coefficients c5 and c7 from Tables 2 and 3 can be used together to examine the combined influence of deep soils and deep geological sediments. When the coefficients from Table 2 Figure 4 compares the spectra of 15 vertical accelerograms recorded at 10 different deep soil sites to the empirical predictions obtained by Equation (1) and the scaling coefficients shown in Table 3. As can be observed, the actual and anticipated spectra have a high degree of agreement.  Figure 2 shows the empirical predictions at deep soil sites (S L = 2) for three different deep geology types. The solid lines reflect median empirical predictions based on the data recorded at all epicentral distances, whereas the dashed lines reflect predictions based on distances less than 30 km.
To calculate the exact differences between the attenuation curves for different site conditions, we can apply the coefficients corresponding to the categorical variables S L and S G . If the coefficients from Table 2 are utilized, the spectral amplitude for T = 0.05 s will be 1/10 −0.021 = 1.05 times larger at deep geological rock (S = 2) than at deep geological sediments (S = 0), as seen in the top left-hand plot in Figure 2. This is most likely because short-period waves travel through more compact rocks more easily. The opposite occurs for larger periods; the spectral amplitude for T = 0.50 s will be 1/10 0.147 = 1.40 times larger at deep geological sediments than at rock, as seen in the bottom right-hand plot in the same figure.
The exact differences between the PSA attenuation curves for different site conditions can be also calculated by using the coefficients c 4 , c 5 , c 6 , or c 7 . For example, the spectral amplitude for T = 0.05 s will be 1/10 0.052 = 1.13 times larger at deep soil sites (S L = 2) than at rock soil sites (S L = 0), while the amplitude for T = 0.30 s will be 1/10 0.088 = 1.23 times larger. The coefficient c 5 will have positive values for all vibration periods up to 0.5 s. This indicates that deep soil sites produce amplification of vertical seismic waves for all vibration periods smaller or equal to 0.5. This is somewhat surprising because our investigation of the horizontal spectra revealed that short-period waves at the deep soil sites are de-amplified [18]. Coefficients c 5 and c 7 from Tables 2 and 3 can be used together to examine the combined influence of deep soils and deep geological sediments. When the coefficients from Table 2 (1) and the scaling coefficients shown in Table 3. As can be observed, the actual and anticipated spectra have a high degree of agreement.
Appl. Sci. 2021, 11, x FOR PEER REVIEW 7 of 16 Figure 4. Relationship between regional vertical pseudo-acceleration spectra at deep soil sites (sL = 2) and the empirical predictions defined by Equation (1) and coefficients from Table 3.

Vertical to Horizontal PSA Ratios
The ratios of vertical to horizontal PSA estimates for various distances and site characteristics are depicted in Figures 5 and 6. These ratios were computed by dividing the vertical PSA amplitudes obtained from Equation (1) and the coefficients from Table 2 by the horizontal PSA amplitudes obtained by Equation (1) and the coefficients from Table A.1 in [18]. The plots also include the corresponding vertical to horizontal PSA estimates as suggested by Eurocode 8 [1] for Type 1 and Type 2 spectra. As can be observed, the probability of empirical estimates, source-to-site distance, vibration periods, as well as both local soil and deep geology conditions, influence the vertical to horizontal PSA ratios. The vertical to horizontal spectral ratios predicted by Type 1 or Type 2 Eurocode 8 spectra differ dramatically as well.
For the 0.05 s vibration period, the vertical to horizontal PSA ratios for deep soil sites atop deep geological sediments are in the range of 1.13-0.85 for the distances between 0 and 150 km. These ratios are 2.2-1.7 times higher than the Eurocode 8 ratio of 0.51, which is obtained for Type 2 spectra for Ground Type C and the same vibration period. For Type 1 Eurocode 8 spectra (used for magnitudes of most contributing earthquakes larger than 5.5), the ratios for T = 0.05 s are 1.5 to 2.0 times higher. The empirical ratios for the 0.30 s vibration period are in the range 0.44-0.50, which is 2.0-2.3 times higher than the corresponding ratios suggested by Eurocode 8 [1] for Type 2 spectra, and almost equal to the ratios for Eurocode 8 Type 1 spectra. For the 0.50 s vibration period, this study shows 0.30 to 0.48 ratios of vertical to horizontal PSA values, which is 1.4 to 2.2 times higher than the corresponding ratios for Eurocode 8 Type 2 spectra but 1.1 to 1.7 times lower for the Type 1 spectra. Figure 4. Relationship between regional vertical pseudo-acceleration spectra at deep soil sites (S L = 2) and the empirical predictions defined by Equation (1) and coefficients from Table 3.

Vertical to Horizontal PSA Ratios
The ratios of vertical to horizontal PSA estimates for various distances and site characteristics are depicted in Figures 5 and 6. These ratios were computed by dividing the vertical PSA amplitudes obtained from Equation (1) and the coefficients from Table 2 by the horizontal PSA amplitudes obtained by Equation (1) and the coefficients from Table A.1 in [18]. The plots also include the corresponding vertical to horizontal PSA estimates as suggested by Eurocode 8 [1] for Type 1 and Type 2 spectra. As can be observed, the probability of empirical estimates, source-to-site distance, vibration periods, as well as both local soil and deep geology conditions, influence the vertical to horizontal PSA ratios. The vertical to horizontal spectral ratios predicted by Type 1 or Type 2 Eurocode 8 spectra differ dramatically as well.
For the 0.05 s vibration period, the vertical to horizontal PSA ratios for deep soil sites atop deep geological sediments are in the range of 1.13-0.85 for the distances between 0 and 150 km. These ratios are 2.2-1.7 times higher than the Eurocode 8 ratio of 0.51, which is obtained for Type 2 spectra for Ground Type C and the same vibration period. For Type 1 Eurocode 8 spectra (used for magnitudes of most contributing earthquakes larger than 5.5), the ratios for T = 0.05 s are 1.5 to 2.0 times higher. The empirical ratios for the 0.30 s vibration period are in the range 0.44-0.50, which is 2.0-2.3 times higher than the corresponding ratios suggested by Eurocode 8 [1] for Type 2 spectra, and almost equal to the ratios for Eurocode 8 Type 1 spectra. For the 0.50 s vibration period, this study shows 0.30 to 0.48 ratios of vertical to horizontal PSA values, which is 1.4 to 2.2 times higher than the corresponding ratios for Eurocode 8 Type 2 spectra but 1.1 to 1.7 times lower for the Type 1 spectra. Appl. Sci. 2021, 11, x FOR PEER REVIEW 8 of 16

Uniform Hazard Spectra Calculations for Vertical Strong Ground Motion
We now conduct a probabilistic seismic hazard analysis (hereinafter, PSHA) for the seismic microzonation of the area within the blue rectangle shown in Figure 7. Several recent regional microzonation studies [26][27][28][29][30][31] have shown that empirical predictions based on scaling equations that consider both deep geology and local soil are in excellent agreement with recorded ground motions in the region. In this study, we define the GMPEs for the probabilistic hazard analyses for Osijek using Equation (1) and the coefficients from Table 2, as was completed for the horizontal spectra [18].

Uniform Hazard Spectra Calculations for Vertical Strong Ground Motion
We now conduct a probabilistic seismic hazard analysis (hereinafter, PSHA) for the seismic microzonation of the area within the blue rectangle shown in Figure 7. Several recent regional microzonation studies [26][27][28][29][30][31] have shown that empirical predictions based on scaling equations that consider both deep geology and local soil are in excellent agreement with recorded ground motions in the region. In this study, we define the GMPEs for the probabilistic hazard analyses for Osijek using Equation (1) and the coefficients from Table 2, as was completed for the horizontal spectra [18].

Uniform Hazard Spectra Calculations for Vertical Strong Ground Motion
We now conduct a probabilistic seismic hazard analysis (hereinafter, PSHA) for the seismic microzonation of the area within the blue rectangle shown in Figure 7. Several recent regional microzonation studies [26][27][28][29][30][31] have shown that empirical predictions based on scaling equations that consider both deep geology and local soil are in excellent agreement with recorded ground motions in the region. In this study, we define the GMPEs for the probabilistic hazard analyses for Osijek using Equation (1) and the coefficients from Table 2, as was completed for the horizontal spectra [18].
European seismic source zone model, which was one of the outcomes of the "Seismic Hazard Harmonization in Europe" (SHARE) Project [42][43][44]. Figure 7 shows the boundaries of the seismic source zones used in this study's hazard calculations. The circles in Figure 7 depict the source-to-site distances of 111 km, 185 km, and 229 km, estimated to contribute to 99 percent of the overall PSHA estimations for the 475-year return period for spectral amplitudes at 0.05 s, 0.3 s, and 0.5 s, respectively. As expected, distant severe earthquakes will have a larger effect on long-period vertical spectral amplitudes, while the short-period amplitudes will be dominated by local seismic activity. As also shown in Figure 7, the distances that will contribute to the 50% accuracy of the PSHA calculations for the 0.5 s spectral amplitude are designated by a circle with a radius of 58 km.   Figure 9 shows cumulative disaggregation examples for distance and magnitude. It also shows (in the bottom-right plot) the magnitude recurrence curve for the most contributing seismic source zone around Osijek. As can be seen, the real return periods of the most contributing earthquakes differ for different vibration periods and are much smaller than the so-called return period Tr, which is equal to the reciprocal value of the mean yearly rate of occurrence, N(PSA), of seismic events that will cause a PSA amplitude that will exceed the expected value-see the related discussions in [17,18]. This should be noted because Eurocode 8 [1] does not define Types 1 and 2 spectra in terms of conservative upper limits for magnitudes, but rather in terms of the magnitudes that "contribute most to the seismic hazard defined for the site for the purpose of probabilistic hazard assessment" (see Eurocode 8 [1]  We assumed deep soil sites and deep geological deposits for the entire case study area. The REASSES V2.0 program was used to create the PSHA estimates [38]. Although the seismic source zone can be defined based on available seismological data [39][40][41], to compare the results to our other studies for the same region [17,18] we utilized the pan-European seismic source zone model, which was one of the outcomes of the "Seismic Hazard Harmonization in Europe" (SHARE) Project [42][43][44]. Figure 7 shows the boundaries of the seismic source zones used in this study's hazard calculations. The circles in Figure 7 depict the source-to-site distances of 111 km, 185 km, and 229 km, estimated to contribute to 99 percent of the overall PSHA estimations for the 475-year return period for spectral amplitudes at 0.05 s, 0.3 s, and 0.5 s, respectively. As expected, distant severe earthquakes will have a larger effect on long-period vertical spectral amplitudes, while the short-period amplitudes will be dominated by local seismic activity. As also shown in Figure 7, the distances that will contribute to the 50% accuracy of the PSHA calculations for the 0.5 s spectral amplitude are designated by a circle with a radius of 58 km. Figure 8 shows examples of hazard disaggregation [45] for the coordinates 45 • 32 N, 18 • 23 E. Figure 9 shows cumulative disaggregation examples for distance and magnitude. It also shows (in the bottom-right plot) the magnitude recurrence curve for the most contributing seismic source zone around Osijek. As can be seen, the real return periods of the most contributing earthquakes differ for different vibration periods and are much smaller than the so-called return period Tr, which is equal to the reciprocal value of the mean yearly rate of occurrence, N(PSA), of seismic events that will cause a PSA amplitude that will exceed the expected value-see the related discussions in [17,18]. This should be noted because Eurocode 8 [1] does not define Types 1 and 2 spectra in terms of conservative upper limits for magnitudes, but rather in terms of the magnitudes that "contribute most to the seismic hazard defined for the site for the purpose of probabilistic hazard assessment" (see Eurocode    Osijek (see Figure 8) and cumulative disaggregation for distances and magnitudes for the coordinates 45°32′ N, 18°23′ E. Figure 10 shows the seismic microzonation maps for the study area for the return periods Tr of 95, 475, 975, and 2475 years. Finally, uniform hazard spectra (hereinafter, UHS) for four distinct probability levels were created for the coordinates 45°32′ N, 18°23′ E, and compared to the Eurocode 8 [1] Type 2 spectra for Ground Type C, as shown in Figure 11. The UHS can be easily constructed for any point within the examined area by reading and linking the values for the provided coordinates in the hazard maps of Figure 11 with the same probability of exceedance.
The resulting UHS spectra for various deep geology and local soil types were compared to the Eurocode 8 spectra for Ground Type C, scaled by the PGA values at  . The magnitude recurrence curve for the most contributing seismic source zone surrounding Osijek (see Figure 8) and cumulative disaggregation for distances and magnitudes for the coordinates 45°32′ N, 18°23′ E. Figure 10 shows the seismic microzonation maps for the study area for the return periods Tr of 95, 475, 975, and 2475 years. Finally, uniform hazard spectra (hereinafter, UHS) for four distinct probability levels were created for the coordinates 45°32′ N, 18°23′ E, and compared to the Eurocode 8 [1] Type 2 spectra for Ground Type C, as shown in Figure 11. The UHS can be easily constructed for any point within the examined area by reading and linking the values for the provided coordinates in the hazard maps of Figure 11 with the same probability of exceedance.
The resulting UHS spectra for various deep geology and local soil types were compared to the Eurocode 8 spectra for Ground Type C, scaled by the PGA values at  Figure 10 shows the seismic microzonation maps for the study area for the return periods Tr of 95, 475, 975, and 2475 years. Finally, uniform hazard spectra (hereinafter, UHS) for four distinct probability levels were created for the coordinates 45 • 32 N, 18 • 23 E, and compared to the Eurocode 8 [1] Type 2 spectra for Ground Type C, as shown in Figure 11. The UHS can be easily constructed for any point within the examined area by reading and linking the values for the provided coordinates in the hazard maps of Figure 11 with the same probability of exceedance. deep soil atop deep geological sediments (s = 0, sL = 2), and exists only for the 95-year and 475-year return periods. For the so-called stiff soils (those with a 15-75 m thick soil layer above the layer with VS,30 > 800 m/s [12,13]) atop deep geological sediments (s = 0, sL = 1), the SPGA is up to 23% smaller the 3.0 factor. On the other hand, for the deep soils atop geological rocks (s = 2, sL = 1), SPGA is up to 33% larger.  The resulting UHS spectra for various deep geology and local soil types were compared to the Eurocode 8 spectra for Ground Type C, scaled by the PGA values at deep soil locations (S L = 2) and for three different deep geology types ( Figure 11). As can be observed, the ratios of maximum spectral amplitudes to PGA values, designated in Figure 11 as S PGA , diverge from the 3.0 factor recommended by Eurocode 8 [1] for vertical spectra. The difference between the S PGA and the 3.0 factor is the smallest (up to 7%) for deep soil atop deep geological sediments (S = 0, S L = 2), and exists only for the 95-year and 475-year return periods. For the so-called stiff soils (those with a 15-75 m thick soil layer above the layer with V S,30 > 800 m/s [12,13]) atop deep geological sediments (S = 0, S L = 1), the S PGA is up to 23% smaller the 3.0 factor. On the other hand, for the deep soils atop geological rocks (S = 2, S L = 1), S PGA is up to 33% larger. Based on the calculated UHS amplitudes, vertical to horizontal spectral ratios were calculated for four distinct probability levels, as shown in Figure 12. These ratios were also compared to the corresponding ratios of Eurocode 8 [1] for Type 2 spectra and Ground Type C and two sets of empirical GMPEs that consider the simultaneous impacts of deep geology and local soil [15,16]. The first set, labeled USC-CA in Figure 12, was created for California [15], whereas the second set, labeled USC-ExYU, was created for the former Yugoslavia [16]. While both empirical attenuation equations ignore the effect of local soil on the vertical ground motion, the Californian equation, which was based on 1482 acceleration components, indicated that deep geology affects seismic waves in the vertical direction [15]. The ratio shown in Figure 12 for the Californian GMPE corresponds to the vertical ground motion predictions at deep geological sediments. In USC-CA [15] and USC-ExYU [16] GMPEs, the product of the categorical variable and the frequency-dependent scaling coefficients is added to be able to predict either vertical or horizontal spectra. As a result, the exact vertical to horizontal PSA ratios can be determined in the same way as was completed previously for the amplification factors for various soil and deep geology types. For example, according to USC-CA [15], the vertical to horizontal PSA ratio is equal to 10 −0.168 = 0.68 for a 0.04 s vibration period and 10 −0.212 = 0.61 for 0.50 s. According to USC-ExYU [16], the vertical to horizontal spectral ratio would be equal to 10 −0.141 = 0.72 for T = 0.04 s and 10 −0.266 = 0.54 for T = 0.60 s. As can be seen in Figure 12, these ratios are in fair agreement with our empirical predictions.
Employing the Eurocode 8 ratios for Type 2 spectra and Ground Type C results however in a severe underestimation of the vertical spectra. This underestimation exists for all the analyzed vibration periods and slightly decreased with the return period. For the vibration period of 0.5 s, the Eurocode 8 ratios are 1.6 to 1.9 times lower than the empirically determined ratios for deep soils atop deep geological sediments (s = 0, sL = 2). The maximum difference can be observed for T = 2.0 s, for which the empirical ratios are 2.6 to 3.0 times higher than the ratios predicted by Eurocode 8. We did not display Eurocode 8 ratios for Type 1 spectra in Figure 12 because they would be substantially higher than any of the empirical ratios. For a 0.05 s vibration period, for example, the empirical ratios would be 1.9 to 2.0 times lower than the ones suggested by Eurocode 8 Type 1 spectra. Based on the calculated UHS amplitudes, vertical to horizontal spectral ratios were calculated for four distinct probability levels, as shown in Figure 12. These ratios were also compared to the corresponding ratios of Eurocode 8 [1] for Type 2 spectra and Ground Type C and two sets of empirical GMPEs that consider the simultaneous impacts of deep geology and local soil [15,16]. The first set, labeled USC-CA in Figure 12, was created for California [15], whereas the second set, labeled USC-ExYU, was created for the former Yugoslavia [16]. While both empirical attenuation equations ignore the effect of local soil on the vertical ground motion, the Californian equation, which was based on 1482 acceleration components, indicated that deep geology affects seismic waves in the vertical direction [15]. The ratio shown in Figure 12 for the Californian GMPE corresponds to the vertical ground motion predictions at deep geological sediments. In USC-CA [15] and USC-ExYU [16] GMPEs, the product of the categorical variable and the frequency-dependent scaling coefficients is added to be able to predict either vertical or horizontal spectra. As a result, the exact vertical to horizontal PSA ratios can be determined in the same way as was completed previously for the amplification factors for various soil and deep geology types. For example, according to USC-CA [15], the vertical to horizontal PSA ratio is equal to 10 −0.168 = 0.68 for a 0.04 s vibration period and 10 −0.212 = 0.61 for 0.50 s. According to USC-ExYU [16], the vertical to horizontal spectral ratio would be equal to 10 −0.141 = 0.72 for T = 0.04 s and 10 −0.266 = 0.54 for T = 0.60 s. As can be seen in Figure 12, these ratios are in fair agreement with our empirical predictions.
Employing the Eurocode 8 ratios for Type 2 spectra and Ground Type C results however in a severe underestimation of the vertical spectra. This underestimation exists for all the analyzed vibration periods and slightly decreased with the return period. For the vibration period of 0.5 s, the Eurocode 8 ratios are 1.6 to 1.9 times lower than the empirically determined ratios for deep soils atop deep geological sediments (S = 0, S L = 2). The maximum difference can be observed for T = 2.0 s, for which the empirical ratios are 2.6 to 3.0 times higher than the ratios predicted by Eurocode 8. We did not display Eurocode 8 ratios for Type 1 spectra in Figure 12 because they would be substantially higher than any of the empirical ratios. For a 0.05 s vibration period, for example, the empirical ratios would be 1.9 to 2.0 times lower than the ones suggested by Eurocode 8 Type 1 spectra. Figure 12. Vertical to horizontal spectral ratios for four distinct probability levels, calculated based on UHS and Eurocode 8 amplitudes from Figure 11. Empirical estimates of vertical to horizontal spectral ratios in California [15] and former Yugoslavia [16] are labeled USC-CA and USC-ExYU, respectively.

Discussion and Conclusions
The motive for this research was the fact that most design codes, including Eurocode 8 [1], neglect the influence of the local soil on the vertical spectra. Moreover, while Eurocode 8 recognizes the potential importance of deeper geological site surroundings, stating in Clause 3.1.2(1) that the National Annex may establish a classification scheme that accounts for deep geology [1], Croatia has not included deep geology in its ground classification scheme to date.
In this study, we analyzed vertical spectra for the deep soil sites atop deep geological sediments, in low to medium seismicity regions. The city of Osijek, located on the right bank of the Drava River and in the southern part of the Pannonian Basin, was chosen as the case study region, as it was in our accompanying paper on the horizontal spectra [18]. Underneath Osijek, geological sediments reach depths of up to 2.7 km, with a total thickness of 150-180 m above layers with an s-wave velocity greater than 800 m/s [18]. Recent research studies have suggested that the city of Osijek's building stock is particularly vulnerable [46,47], necessitating the need for more reliable seismic action estimates for Osijek, especially after the devastating 29 December 2020, Mw 6.4 earthquake in central Croatia [48]. This is also the case with many other cities with similar building typologies and similar seismic activity [49][50][51][52][53][54]. New empirical regional PSA attenuation equations were constructed using only regional strong motion data. A total of 218 vertical components of strong strong motion accelerograms recorded in the north-western Balkans from 112 earthquakes with magnitudes ranging from 3 to 6.8 were included in the strong motion database used to build the predictive equations. In these equations, both the effects of local soil and deep geological conditions were taken into account. For a 0.3 s vibration period, the maximum amplification of 1.48 times was found at deep soil atop deep geological sediments compared to the rock sites. Vertical PSA amplitudes for vibration periods of 0.10 s to 0.15 s, as well as periods longer than 1.50 s, were shown to be de-amplified by up to 26%.
As was also shown, using Eurocode 8 ratios of vertical to horizontal spectra results in a significant underestimation of vertical spectra. Vertical to horizontal PSA ratios that are suggested by Eurocode 8 [1] for Type 2 spectra and Ground Type C were 2.0 to 2.3 times lower than the empirically determined ratios for PSAvert at deep soil atop deep Figure 12. Vertical to horizontal spectral ratios for four distinct probability levels, calculated based on UHS and Eurocode 8 amplitudes from Figure 11. Empirical estimates of vertical to horizontal spectral ratios in California [15] and former Yugoslavia [16] are labeled USC-CA and USC-ExYU, respectively.

Discussion and Conclusions
The motive for this research was the fact that most design codes, including Eurocode 8 [1], neglect the influence of the local soil on the vertical spectra. Moreover, while Eurocode 8 recognizes the potential importance of deeper geological site surroundings, stating in Clause 3.1.2(1) that the National Annex may establish a classification scheme that accounts for deep geology [1], Croatia has not included deep geology in its ground classification scheme to date.
In this study, we analyzed vertical spectra for the deep soil sites atop deep geological sediments, in low to medium seismicity regions. The city of Osijek, located on the right bank of the Drava River and in the southern part of the Pannonian Basin, was chosen as the case study region, as it was in our accompanying paper on the horizontal spectra [18]. Underneath Osijek, geological sediments reach depths of up to 2.7 km, with a total thickness of 150-180 m above layers with an s-wave velocity greater than 800 m/s [18]. Recent research studies have suggested that the city of Osijek's building stock is particularly vulnerable [46,47], necessitating the need for more reliable seismic action estimates for Osijek, especially after the devastating 29 December 2020, Mw 6.4 earthquake in central Croatia [48]. This is also the case with many other cities with similar building typologies and similar seismic activity [49][50][51][52][53][54]. New empirical regional PSA attenuation equations were constructed using only regional strong motion data. A total of 218 vertical components of strong strong motion accelerograms recorded in the north-western Balkans from 112 earthquakes with magnitudes ranging from 3 to 6.8 were included in the strong motion database used to build the predictive equations. In these equations, both the effects of local soil and deep geological conditions were taken into account. For a 0.3 s vibration period, the maximum amplification of 1.48 times was found at deep soil atop deep geological sediments compared to the rock sites. Vertical PSA amplitudes for vibration periods of 0.10 s to 0.15 s, as well as periods longer than 1.50 s, were shown to be de-amplified by up to 26%.
As was also shown, using Eurocode 8 ratios of vertical to horizontal spectra results in a significant underestimation of vertical spectra. Vertical to horizontal PSA ratios that are suggested by Eurocode 8 [1] for Type 2 spectra and Ground Type C were 2.0 to 2.3 times lower than the empirically determined ratios for PSA vert at deep soil atop deep geological sediments and with a vibration period of 0.3 s. The vertical to horizontal PSA ratios for a single soil type were shown to be influenced by the source-to-site distance, the probability of empirical estimates, as well as deep geology conditions. Hence, they cannot be captured by a constant that depends only on earthquake magnitude, as suggested by Eurocode 8 [1].
Seismic microzonation maps were constructed for the examined area using the proposed scaling equations. Vertical UHS amplitudes were also calculated and compared to the Eurocode 8 [1] spectra for Ground Type C, scaled by PGA values at deep soil sites, and for three different types of deep geology. The results show that the maximum spectral amplitudes to PGA ratios, referred to as S PGA in Figure 11, do not deviate too much from the 3.0 factor indicated by Eurocode 8 [1] for vertical spectra at deep soils above deep geological deposits. However, for deep soils atop deep geological rocks, the S PGA was equal to 3.9 for 95-year and 4.0 for 475-year, 975-year, and 2475-year return periods. In comparison to the 3.0 factor indicated by Eurocode 8, this is a 30-33% increase.
Vertical-to-horizontal UHS ratios were also calculated and compared to Eurocode 8 recommendations and two sets of empirical GMPEs that were developed by other authors for California [15] and former Yugoslavia [16]. These two sets of GMPEs were the only ones we could find that considered the effects of both local soil and deep geology on vertical spectral amplitudes. While our results are in good accordance with those of two other GMPEs, it was shown that the vertical to horizontal PSA ratios for all the vibration periods are underestimated by Eurocode 8 Type 2 spectra. The biggest divergence may be seen at T = 2.0 s, where the empirical ratios were 2.6 to 3.0 times higher than what Eurocode 8 predicted. The usage of Eurocode 8 Type 1 spectra, on the other hand, would result in an overestimation of vertical spectra. This overestimation would be twice the empirical ratios for T = 0.05 s.
All the results of this study indicate that the vertical design spectra should be defined directly as uniform hazard spectra calculated by using the empirical GMPEs for vertical ground motion. Even though the number of ground motion accelerograms recorded at deep soil sites in the north-western Balkans is not large, the results that are presented in this paper can be considered a first step toward defining more accurate PSA vert estimates for the case study area of Osijek. As more strong motion data for the vertical direction becomes available, it will be simple to update the GMPEs and recalculate the seismic hazard.