Estimation of the High-Frequency Attenuation Parameter Kappa for the Zagreb (Croatia) Seismic Stations

Featured Application: The results from this study of the high-frequency attenuation conﬁrm the e ﬀ ect of local structures on local attenuation and of deep structures on the attenuation at long distances. Additionally, it should be considered as a step towards the better analysis of seismic hazards for wider area of the city of Zagreb. Abstract: The city of Zagreb (Croatian capital) is situated in the contact area of three major regional tectonic units: the SE Alps, NW Dinarides, and Tisza Unit in the southwestern margin of the Pannonian Basin. The Zagreb seismic zone encompasses the Medvednica Mountains and the city of Zagreb with its surrounding areas, which was struck by the strongest instrumentally recorded earthquake (M5.5) on 22 March 2020. The objective of this contribution is the estimation of the high-frequency attenuation spectral parameter kappa ( κ ) and its local site-speciﬁc component for the Zagreb (Croatia) seismic stations to which we were particularly encouraged after the scale of the damage after the Zagreb 2020 earthquake. We tested linear dependence of κ with epicentral distance using traditional linear least square regression, linear regression for data with errors, and constrained model at close distances to estimate near-site attenuation ( κ 0 ). Regression-estimated site kappa values at zero-distance are within the range of the uncertainty ( ± 1 standard deviation) with constrained κ 0 value as well within the range of existing global κ 0 and V S 30 (shear wave velocity in the top 30 m) values. Spatial distribution of κ within the Zagreb seismic zone shows that κ is not isotropic and high-frequency attenuation anisotropy is probably a ﬀ ected by local and regional geological variability, regional active faults and a complex tectonic structure in each direction.


Introduction
Attenuation of seismic waves is one of the key factors in seismic hazard assessment for earthquake-prone regions as well as being important for the quantification of earthquakes and plays a significant role in studies of seismic source and crustal structure [1]. The Fourier amplitude spectrum (FAS) of ground motion is influenced by source parameters, propagation path, and local site conditions as modelled by the Brune's [2] theoretical omega-square source model. Hanks [3] was one of the first who observed the "crashing spectrum syndrome" in the acceleration spectrum of shear waves (S-waves) on strong-motion accelerograms recorded at close distances. The high-frequency band limitation parameter (f max ) was defined as the cut-off frequency at which the spectrum starts to decay very rapidly. Hanks [3] concluded that the f max observed in acceleration FAS is controlled by the local site conditions, particularly by the subsurface geological structures below and near the site. One of the first uses of the term "site attenuation parameter" to model the high-frequency spectral attenuation was by Cormier [4]. Empirical spectral decay parameter kappa (κ) was introduced by Anderson and Hough [5] to describe the difference between the observed acceleration spectrum of S-waves and Brune's [2] omega-square source model at the high-frequency part of the FAS. Anderson and Hough [5] observed a linear relation between the κ of the ground motion and the epicentral distance (R e ) of the recording station. The zero-distance intercept κ 0 (also called the near-site attenuation or site-kappa) represents the attenuation contribution to κ from subsurface geological structure beneath the site. The slope (κ R ) is related to the regional attenuation due to the horizontal propagation of S-waves through the crust below and near the site within a few kilometres (e.g., [5][6][7]). Over the last three decades, the near-site attenuation parameter κ 0 has been used in a variety of applications. Particularly in the creation and calibration of ground-motion models (GMMs) based on stochastic simulations (e.g., [8][9][10][11]), host-to-target adjustments of GMMs (e.g., [12,13]) and site-specific ground response analysis using Random Vibration Theory (e.g., [14]) for critical facilities where the choice of seismological parameters play an important role.
On 22 March 2020 Zagreb was struck by the strongest ever instrumentally recorded earthquake (M5.5) that caused significant damage in the city itself. The seismological data were reliable enough for preliminary 3D interpretation of the rupture along a segment of an active NW verging thrust fault in the subsurface of the Medvednica Mountains [15]. The aim of this paper is a more reliable interpretation of near-site attenuation or site-kappa based on seismological data obtained up to now.
The objective of this contribution is to present the calculation of high-frequency attenuation spectral parameter κ and its local site-specific component κ 0 for two seismological stations: Puntijarka-PTJ (situated on the top of Medvednica Mountains) and Zagreb-ZAG (situated in Zagreb city), using classical Anderson and Hough [5] approaches. The results from this study of the high-frequency attenuation will contribute to a better knowledge of the effect of local structures on local attenuation and of deep structures on the attenuation at long distances. Additionally, attenuation parameters (κ, κ 0 ) should be considered as a step towards a better analysis of seismic hazard and finally towards seismic hazard assessment for the wider area of the city of Zagreb that has a population close to 1 million and the city's strategic infrastructure.

Geological Features of the Zagreb Area
The city of Zagreb (the Croatian capital) is situated in the contact area of three major regional tectonic units [16,17]: the SE Alps, NW Dinarides, and Tisza Unit in the southwestern margin of the Pannonian Basin (Figure 1a). Interaction of the units causes earthquakes in the upper crust [18] that are distributed predominantly along the regional active faults ( [19], Figure 1b). The NW verging reverse fault system striking along the Medvednica Mountains is probably related to a major seismogenic thrust fault [15] in the wider city of Zagreb region, and is dissected by active transversal strike-slip faults (Figure 1b). A system of south-dipping longitudinal reverse faults has been mapped on the surface along the Medvednica Mountains ( [20][21][22], Figure 1c), and interpreted in one of the latest tectonic models as the youngest (Pliocene to Quaternary) inverse structure ( [19], Figure 1c). A master thrust could represent a basal detachment (décollement) of the system of steeper reverse faults striking along the mountain range. Although the frontal part of the thrust could be steeply inclined towards the surface along the northern margin of the Medvednica Mountains, or it could be buried below the Neogene Hrvatsko Zagorje basin ( [15], Figure 1c). Due to the Pliocene-Quaternary N-S directed shortening, Medvednica Mountains experienced approx. up to 2000 m of differential uplift and now exposes pre-Neogene and Neogene tectonic and stratigraphic units that are surrounded by Pliocene-Quaternary sediments [23]. According to the DEM-based morphometric analysis the ongoing tectonic uplift in Medvednica Mountains area is concentrated along the SW corner of the NW mountain front and along a narrow boundary zone that divides the Medvednica Mountains into the NE and SW regions. According to [24] this area is considered to be tectonically active, with low uplift rates of approx. 0.023 mm/y which is in good agreement with the 0.17-0.4 mm/y Quaternary uplift rate reported by [25] for the central part of Medvednica Mountains. During the recent rupture (March 2020) along a segment of an active NW verging thrust fault in the subsurface, the central part of the mountain was uplifted by up to 3 cm [15].
The mountain is considered as uplifted bedrock, composed of Paleozoic and Triassic metamorphic and volcanic rocks in the core, obducted ofiolitic mélange and overlying Cretaceous carbonate and flysch rocks. Miocene transgressive successions of the Pannonian basin unconformably cover the eroded Phanerozoic basement all around the mountain ( [20][21][22], Figure 1c). Predominantly unconsolidated superficial Pliocene-Quaternary sediments (gravel, sand and clay) are up to 2000 m thick in the subsurface of Sava basin on the south and the Hrvatsko Zagorje basin on the north (Figure 1c), and pinch-out towards Medvednica Mt. ( [20][21][22]. Thus, the seismic station Puntijarka (PTJ) is situated on the bedrock that can be considered as a firmly consolidated basement, while the seismic station Zagreb (ZAG) is situated on the Quaternary alluvial sediments of the Sava river that are mostly unconsolidated materials (Figure 1a).

Seismicity of the Zagreb Area
Several seismic zones (Zagreb, Novo Mesto-Krško, Karlovac-Metlika and Pokupsko-Petrinja) are situated on the wider Zagreb area as can be seen on Figure 2 (marked with coloured shadowing).
The most important one is the Zagreb seismic zone which encompasses the Medvednica Mountains and the city of Zagreb with its surrounding areas. Epicentre locations of the strongest earthquakes are mostly situated on the north-western slopes of the Medvednica Mountains. The strongest earthquake was on 9 November 1880, close to Kašina, in the village Planina, and is known as the "Great Zagreb earthquake". This was the first Croatian earthquake for which the value of the focal depth (16 km) was determined based on macroseismic investigations. The analysis of the earthquake effects enabled the determination of the epicentral intensity (VIII • MCS). The earthquake itself is extremely well documented [26] due to extensive material damage. Of the 3670 buildings (Zagreb had around 30,000 inhabitants at the time) all were damaged and about 13% were destroyed. After 1900 the strongest event that happened in Zagreb area occurred on 22 March, 2020, 140 years after the "Great Zagreb earthquake" [15]. Its epicentre was 7 km north of the centre of Zagreb, in the vicinity of Markuševac andČučerje, and the hypocentre was at a depth of 10 km. This M5.5 earthquake was felt with a maximum intensity of VII-VIII • MSK. Regarding other recent seismic events, prior to the Zagreb event of 2020 there was only a single strong earthquake in this area. The M5 event occurred in 1990, 3.5 km south of Kraljev Vrh, where the worst of the damage was observed.
The strongest earthquakes in the Novo Mesto-Krško seismic zone occurred mainly in the period up to 1900. The maximum intensity of these events did not exceed VIII • MCS. The strongest earthquakes (with the maximum intensity VIII • MCS) occurred in 1634 (Novo Mesto area) and 1917 (around Brežice). The strongest earthquake after 1900 was in 1974. It occurred approximately 15 kilometres north of Krško with the magnitude of 5.1.
Seismic zone Karlovac-Metlika encompasses the areas of Žumberačka gora and Samoborska gora together with Ozalj and Jastrebarsko. It is characterized by the occurrence of a large number of small earthquakes and smaller number of stronger ones. The majority of earthquakes (95%) that have occurred in this area have a magnitude of less than 3.0. Some authors have mentioned two historical earthquakes southeast of Metlika, within this seismic zone, with the intensity of IX • MCS, in 567 and 1097. However, the information that is available on these earthquakes is very unreliable. We will emphasize the earthquakes whose locations and intensities are known with a higher degree of reliability. These are two earthquakes of intensity VIII • MCS, in 1645, 6 km west of Ozalj and in 1699, with the epicentre in Metlika. In addition to these two events, four more occurred with intensities greater or equal to VII • MCS (in 1645, 1881, 1887 and 1974).

Figure 2.
Spatial distribution of earthquake locations in the investigated area (373BC-2019, according to the Croatian Earthquake Catalogue-CEC, the updated version first described in [27]. Seismic zones are marked as: Zagreb-green, Novo Mesto-Krško-red, Karlovac-Metlika-blue and Pokupsko-Petrinja-yellow. Faults are marked with black lines [28,29]. The locations of epicentres of strong earthquakes in the Pokupsko-Petrinja seismic zone stretch in the N-S direction. The strongest event here is one of the best-known earthquakes, not only in Croatia but also globally. It is the earthquake that occurred on 8 October, 1909, with the epicentre 9 km north of Pokupsko. Therefore, it is known as the Pokupsko earthquake. The magnitude of the earthquake was 5.8, with the focal depth of 14 km, and the effects that it caused were rated with a maximum intensity of VIII • MCS. The Pokupsko earthquake is one of the few historical earthquakes that are regularly referred to in seismological textbooks and other books on the history of science due to it's study leading to important discoveries that have marked a turning point in the understanding of earthquakes and their effects. By analysing the earthquake, the great Croatian scientist and geophysicist Andrija Mohorovičić proved the existence of the discontinuity of velocity that separates the crust from the mantle of the Earth (the Mohorovičić discontinuity). The main Pokupsko earthquake was followed by a series of more than 50 aftershocks before the end of 1910. The strongest of which had a magnitude of 5.3, with no other earthquakes after 1910 exceeding M4.5.

Methods and Data
We applied Anderson and Hough's [5] well-known method (AH84) for estimating κ by analysing the S-wave spectra in this study following Ktenidou et al. [7,11] guidelines. The total path attenuation of S-waves within the crust is separated into two attenuation parameters: the frequency-dependent quality factor Q(f ) and the near-surface attenuation parameter κ (also termed site diminution parameter, κ 0 ) as shown in Equation (1) where A(f,t) is acceleration spectrum that contains the effects of source, distance and local site effects. In Equation (1), t * is attenuation time [4], f is frequency, R is distance, β 0 is shear wave velocity of the Earth's medium: Anderson and Hough [5] observed a linear dependence between the calculated κ of the ground motion and the epicentral distance (R e ) of the recording station and proposed a mathematical formula (Equation (2)) that treats κ as a function of R e . In Equation (2), the zero-distance intercept κ 0 represents the attenuation contribution to κ from geological structure beneath the site (near-site attenuation or site-kappa), and the distance dependent part in term of the slope (κ R ) is related to the regional attenuation: Value of κ gradually increases with distance, and degree of the slope κ R depends on the local and regional geology. Since its introduction, general observation of Anderson and Hough [5] study that κ is a linear function of R e has proven to be good approximation as shown in the recent studies that estimated κ (e.g., [1,6,7,[30][31][32][33]). In the last decade relationship between estimated κ 0 and measured V S30 values were proposed; κ 0 has lower values for sites on harder rocks (higher V S30 values) and higher κ 0 values are obtained on softer rocks (e.g., [11]).
Two seismological stations included in this study are: Puntijarka (PTJ), located at the top of Medvednica Mountains and Zagreb (ZAG), located in the city area ( Figure 1). A compiled ground motion data set is selected within M L ≥ 3.0, R e ≤ 150 km and a focal depth of h ≤ 30 km. Figure 3 evaluates the compiled ground motion dataset distribution of the epicentres (R e ) and focal depth with magnitudes (M L ) at each station. The number of recordings at each station depends not only on local seismicity and data selection that is the same for both stations, but also on the operational period of each station (from 2000) and usable frequency sampling rate. In this study we used 50 Hz seismograms recorded from 2004 to 2019 (PTJ, 127 recordings) and from 2007 to 2019 (ZAG, 92 recordings). We calculated the parameter κ by using the classical Anderson and Hough [5] approach (Equations (1) and (2)) for all three components (north-south-NS, east-west-EW and vertical-Z) of weak motion seismograms with frequency sampling rate of 50 Hz. Figure 4 represents two examples of step-by-step κ calculation for seismograms recorded on ZAG and PTJ stations. The S-wave window (with a part of coda tail that cannot be avoided) is selected for each record with a minimum duration of 3 s and transformed via FFT to obtain FAS. The dropout of FAS at frequencies above 20 Hz was due to anti-alias filter. Therefore, FAS above 20 Hz (80% of the Nyquist frequency) was not considered in the κ calculation. We have used case-by-case evaluation in order to visually inspect the possible reasons of excluding data from the analysis (e.g., missing data, poor recording, strong resonance peaks, noise effects) that can significantly affect calculation of κ [7]. The high-frequency decay parameter κ is calculated from the slope of FAS in the linear-logarithmic space for the high-frequency range between f 1 and f 2 . In most cases, f 1 is picked as the lower bound of the high-frequency slope before the FAS starts to decrease rapidly (f max ) and is always picked to be higher than the corner frequency f c to exclude possible source contribution (e.g., [7,33]), whereas f 2 is fixed to 20 Hz in all cases. We calculate corner frequency shown in Figure 4 based on theoretical FAS (black line, [10]) using known seismological parameters for each recording based on EQ catalogue (updated version first described in [27]). We adopted Q(f) value from [34] as Q c (f ) = 78f c 0.69 for PTJ and Q c (f ) = 45f c 0.92 for ZAG stations with preselected kappa value equal to 0.0 s. Afterwards, we used an empirically estimated kappa value using AH84 method (magenta fit line) to determine the Brune [2] theoretical FAS for each recording in order to compare how the high-frequency attenuation parameter kappa would fit within the full theoretical spectrum fit (black dashed line). This spectrum can be called the semi-empirical spectrum as most known empirical earthquake parameters are used to plot this spectrum. An exact match is impossible due to differences and errors in estimated seismological parameters. However, a good match between the semi-empirical spectrum and high-frequency kappa fitting is observed in most cases.

Results
In the present work we applied original the Anderson and Hough [5] linear formulation (Equation (2)) to a pair of individual κ values with epicentral distances R e for all records at each station.
We estimated the value of site-specific (near-site) attenuation parameter κ 0 based on proposed κ model as a function of epicentral distance. Calculated individual κ values for EW and NS components we combined into an average value of κ hor for each recording. In certain cases where they differ significantly from each other (difference > 25%), the recording is excluded for further regression analysis to obtain the site kappa value (e.g., [7,11]). 81 usable pairs of κ and R e data for κ models were finally chosen (from the initial 92) for the ZAG station and 118 (from the initial 127) for the PTJ station.
In this study, both κ models (κ hor and κ ver ) as a function of epicentral distance are proposed to estimate site-specific parameters κ 0hor and κ 0ver using linear least-square regression (e.g., [7]). Both variables, κ and R e are imperfectly known and errors in R e could have impact on the final κ 0 value and the slope of regression κ R . Therefore, instead of using the traditional linear least-squares regression, we performed linear regression suitable for data with errors, following the method by York et al. [35]. This allowed us to test how the existence and correlation for the observational errors in two coordinates (R e and κ) affect values of κ 0 and κ R , and if there exists a significant difference regarding the standard linear regression.
Typically, we set the standard error for R e to vary ±2, 5 and 10 km and error in κ set to be ±1 or 2 standard deviations based on observed empirical data. A regression test showed that differences between standard and error-in-variables linear regression κ 0 and κ R values vary mostly between 1-5% (Table 1). Highest differences up to 9% are within standard error in R e ± 10 km and smallest below 1% within standard error in R e ±2 km with respect to ±1 or 2 standard deviations in κ. However, ±1 standard deviations in κ is too small to be consider representative as can be seen by data scatter in Figure 5. We choose to evaluate final κ 0 and κ R values using errors for epicentral distances to be in order of 5 km with errors in κ to be ±2 standard deviations as shown in Figure 5. Summarized results and differences between standard linear regression and errors-in-variables regression are shown in Table 1.   [35]) compared to linear least-square (LSE) regression with differences (in %). Estimation of individual kappa values is very sensitive as can be observed on Figure 4. This is particularly important to consider when kappa is estimated by different users and different computational procedures as described in [7]. For this reason, taking 5% uncertainty as a difference between standard least-square regression without consideration of errors in variables and error-in-variables, yield to site kappa value that is more conveniently estimated as shown in Table 1. With less data, large data scatter and lack of data at shorter epicentral distances, differences between two regression methods could be significant, and that can affect final κ 0 and κ R values. Less data and lack of data at shorter epicentral distances are related to the seismicity of an area and operational years of seismological stations.

HOR VER
However, large data scatter in κ and R e is not only sensitive to differences in methods and data, and to the local and regional characteristics of study area (the effect of distance and local site effects), but also to complex Q structure, different ray paths across faulting structure for different earthquakes [36,37] and anisotropic attenuation properties [1]. Therefore, we choose an area within about 20 km of the epicentre as near-field attenuation [38] to constrain κ 0 value to κ value at shortest R e ( Figure 6) using nearby records (hockey-stick model proposed by [36]). This enabled us to determine site kappa value affected primarily by upper local site condition that would not be influenced by changes in slope of κ R related to Q; in order to minimize influence of regional effects (deeper local site effect + source and path effects, [32]) on κ 0 . Additionally, by restricting the distance range for kappa we ensure that local site attenuation at zero distance is overwhelmed by the path component at larger distances [36,37].  Figure 6 shows comparison of results from different methods to estimate κ 0 for the PTJ and ZAG stations. It can be observed that κ 0 values are constrained for shorter distances to catch mainly upper local site conditions and data ranges are similar within errors. κ 0 values estimated from the extrapolation to zero distance (preferably using error-in-variables regression) are within the range of the uncertainty (±1 standard deviation) with a constrained κ 0 value ( Figure 6, Table 2). At larger distances, variations between empirical kappa values are scattered among similar distances due to heterogeneity along the paths between different sources and stations.  [39,40]. The range of the uncertainty with constrained κ 0 value is set to ±1 standard deviation. It is possible to compare estimated values of frequency-dependent Q(f) for the high-frequency range with frequency-independent Q estimated from κ [5,7,11,32]. The Q value is inversely proportional to slope κ R as a regional attenuation contribution, Q(κ R ) = 1/β·κ R (typical value of β = 3.5 km/s is used as an average shear wave velocity of the crust). To compare two independent attenuation studies, frequency-dependent attenuation, the Q(f ) value, is adopted from [34] as Q c (f) = 78f c 0.69 for the PTJ and Q c (f) = 45f c 0.92 for ZAG stations for a time lapse equal to 30 s. For the high-frequency range 10-20 Hz this yields to Q c (f = 10-20 Hz) = 382-616 for PTJ and Q c (f = 10-20 Hz) = 374-708 for ZAG. Common regional effects Q(κ R ) averaged for high frequency attenuation and azimuthal changes, estimated from the slope κ R for both components (horizontal and vertical) are: Q(κ R ) = 913-1120 for PTJ and Q(κ R ) = 1143-1348 for ZAG. Differences in the frequency independent Q(κ R ) are probably the result of higher local site attenuation effects at ZAG station compared to PTJ station, azimuthal differences and data quality between stations. Regional Q(κ R ) slope variation within ±2 standard deviations ( Figure 6) are shown based on Q c (f ) increase with frequency for a different lapse time t L [34]. Comparing high-frequency coda attenuation values with frequency independent, Q c (f ) and Q(κ R ) for both stations it can be seen that they are comparable (difference is 32-67%). However, different azimuthal regional effects are neglected within average Q values under the assumption that they are affected by the same regional characteristics. Q(κ R ) slope variation shown in Figure 6 could indicate influence of regional azimuthal data subsets attenuation effects. But also differences in the average shear wave velocity of the crustal structure at the contact of the Dinarides and the Pannonian Basin affect Q(κ R ) values. Large differences among different Q approaches can be explained by inherent errors of Q and κ measurement (often in the order of ±50%) and different attenuation methods in Q and κ calculation (different data, band-pass filters, time window length). This shows that frequency-dependent attenuation cannot be ruled out for κ [1,33].
As expected, site kappa values for PTJ and ZAG are different as these stations are situated in different local site characteristics (Figure 1, Table 2). In terms of local site conditions, estimated κ 0 values for both stations compared to V S30 values ( [39,40]; for PTJ: site type A (800 m/s) by EC8 [41] and for ZAG: site type B (360-800 m/s) by EC8) are within the range of existing κ 0 -V S30 correlations [11,32]. It confirms that lower values of κ 0 are observed for sites on harder rocks (lower attenuation) and higher κ 0 values for sites on softer rocks (higher attenuation). Evaluation of the site-specific κ 0 is very sensitive since the AH84 method requires high-magnitude events, something that is difficult to acquire when the local site is situated in a low-to-moderate seismicity area. Correlations between κ 0 -V S30 [11,32] could be helpful to estimate κ 0 for a certain target local site. However, both variables κ 0 and V S30 have some uncertainties. Especially, the errors in V S30 may even be significantly higher than the error in κ 0 (percentage-wise) due to numerous reasons (different geophysical methods, correlations, analysis, regions, etc.), as observed in [36]. The preferred approach should be estimating the site kappa value from earthquake recordings, which is occasionally difficult in low seismicity areas. The scarcity in short distances affects the estimation of site kappa and constraining to fixed value within errors could be helpful. In contrast, greater distances recordings control the slope of kappa regression where potential uncertainties in location could exists (error-in-variable regression is preferred) as well as stronger influences of regional attenuation effects [36].

Discussion
Local and regional geological and tectonic characteristics around each station are important to define primary effects on the κ distribution within their azimuthal area subsets [1,7,42] estimated κ from the vertical component of ground motion and observed that the κ ver values are slightly lower than the average κ hor , something that is also observed in our study. Generally, the vertical component of the ground motion is mainly controlled by the source effect and exhibits relatively less sediment-induced amplification than horizontal components (e.g., [43]). Under the assumption that the local site effects are more pronounced in the horizontal components than in the vertical components of ground motion, we present the spatial distribution of horizontal κ values (Figure 7) to determine if there exist different trends in the high-frequency regional attenuation between different azimuthal area subsets. The lowest κ values are spatially distributed within a few kilometres around the stations due to near-site effects and their values matching the κ 0 value for horizontal components (PTJ: 0.0283-0.0360 s, ZAG: 0.0411-0.0475 s), particularly using constrained κ 0 value based on hockey-stick model [36]. Station PTJ is situated on the top of the Medvednica Mountains as a part of the Pannonian Basin-External Dinarides-Alps transition zone. The station site area is composed of bedrock Palaeozoic-Triassic ortometamorphites and parametamorphites. In the lowland Zagreb area, where station ZAG is situated, Miocene semi-consolidated rocks to recent unconsolidated Pliocene and Quaternary sediments of the rivers Sava to the south and Krapina are present (Figure 1). Compared to [1], κ 0 values for Zagreb area are higher than those for Kvarner region stations (Brijuni: 0.0249-0.0362 s, Rijeka: 0.0212-0.0239 s and Novalja: 0.0225-0.0235 s). The Kvarner stations sites and surrounding area are composed mainly of carbonates and the lowest κ values are regionally stabilized for hard rocks [11,32,37].
In an opposite way, the addition of softer rocks, within a few kilometres around the stations increases κ values [11,30,32]. The κ regionalization model schematically illustrated by [32] proposed division of near-site attenuation effects into upper local site effects (correlation with V S30 ), deeper site effects along near-field regional path (constrained model) and regional path effects from source. The κ 0 values for Zagreb stations (PTJ and ZAG) are similar to those of Ozalj station (0.0377-0.0412 s) that is situated in a zone consisting mainly of consolidated rocks that are covered by thin alluvial sediments along the near Kupa River valley [1]. According to the model illustrated by [32], surrounding areas of Hrvatsko Zagorje and Sava Basins, with sediments depth up to 2 km (Figure 1), have influence on the increase of the high frequency attenuation and local site attenuation due to the addition of the deeper site attenuation. The older and more competent carbonate rocks in the top few kilometres of the crust are expected to have a lower kappa [37] as the situation is in the Kvarner region [1], than younger and less competent Miocene softer rocks in the uppermost crust surrounding the Medvednica Mountains.
If we take a look at the gradual increase of κ with distance from the stations, the regional path effect from the source attributed to the κ is described by the slope regression κ R . It can be observed that the spatial κ distribution is not isotropic and different trends between estimated lower and higher attenuation of ground motion subsets are observed. According to [1], the source of anisotropy is most likely attributed to orientations of the cracks and fractures under the local tectonic stress field, waveguides along major faults and attenuation within fault zones, S-wave reflections from different parts of the shallower Moho discontinuity. Variations in observed κ distribution between lower and higher attenuation subset areas is clearly under influence of local and regional geological structures; to the north influenced by the Ivanščica Mt. and Kalnik Mt., to the south by Žumberak Mt., and in its centre by the Medvednica Mountains, that are distributed predominantly within regional active faults ( [19], Figure 1) and main seismic zones (Figure 2). We limit spatial distribution of κ ( Figure 7) within distance of 80-100 km from Zagreb area to avoid additional path complications, since ground motion ray paths from source to station show a wider range of Q ( Figure 6, [37]), and to derive conclusions within seismic zones Novo Mesto-Krško, Karlovac-Metlika and Pokupsko-Petrinja (Figure 2), that influence seismic hazard of the Zagreb city area.
As pointed in [1], observed lower attenuation properties may be caused by S-wave reflection from different parts of the shallower Mohorovičić discontinuity, which is less than 25 km in the Pannonian basin and thickens towards the Dinarides [44]. However, lower values of kappa along the SW-NE striking Medvednica structure (Figure 7) are probably caused by the shallow subsurface to surface position of the metamorphic and volcanic core of the mountain (Figure 1b). The core is situated below the major detachment and, thus, probably has not undergone thin-skin deformations that additionally cracked the rocks above the detachment.
The attenuation of the earthquake ground motion is complex in nature, and high-frequency kappa physics is still unresolved [32,33,36]. The regional propagation difference in the transition zone between Dinarides and the Pannonian basin is observed from the kappa spatial distribution. It can be concluded that the regional difference in attenuation properties of the rocks in the transition zone between Dinarides and the Pannonian basin in this area are far from isotropic. In comparison to [1], the obtained results are similar and consistent. However, some differences are evident, pointing to effects such as local intrinsic attenuation anisotropy from different causes, e.g., scattering due to heterogeneity, fracturing, flow of fluids in rocks. The results from both studies taken together indicate a change in the attenuation properties from the Dinarides to the Pannonian basin.
The main question is how the high-frequency attenuation observed along the propagation path affects seismic hazard of the Zagreb area. As shown in Figure 2, Zagreb city is under the influence of the near-field Zagreb seismic zone and close seismic zones Novo Mesto-Krško, Karlovac-Metlika and Pokupsko-Petrinja. It is also situated in the contact area of three major regional tectonic units ( Figure 1). On 22 March 2020, Zagreb was hit by the strongest earthquake since 1880 [15], 5.5 on the Richter scale, with a shallow focal depth of <10 km. It occurred at 06:24 local time and hit the Croatian capital Zagreb as well as two neighbouring counties, namely Zagreb County (Sava Basin) and Krapina-Zagorje County (Hrvatsko Zagorje Basin). The maximal intensity of VII-VIII • EMS-98 scale was observed close to the earthquake epicentral area, just 7 km north of the centre of Zagreb, a city with a population of close to one million. Earthquake was felt throughout Croatia (Figure 8), as well as in neighbouring countries (Slovenia, Bosnia and Herzegovina, Hungary). Attenuation of macroseismic field (Figure 8) can be possibly linked to the kappa spatial distribution (Figure 7) if the same physical mechanisms are considered [1]. If particular paths from Zagreb city are looked comparing Figures 7 and 8 (Varaždin, Ivanščica Mt.-Ivanec, Križevci-Kalnik Mt., Hrvatsko Zagorje Basin, Sava Basin, Karlovac, Sisak), Zagreb M5.5 earthquake was felt more widely due to different high-frequency attenuation anisotropy effects (Figure 7). This shows that the local attenuation and the attenuation at long distances are influenced by local and regional geological variability, regional active faults and complex tectonic structure in each direction. When all this is taken together, they certainly plays significant role in seismic hazard of the Zagreb city area. Thus, this is a step towards a better definition of the input parameters for seismic hazard and risk assessment.

Conclusions
This study presents the estimation of high-frequency attenuation of the spectral parameter κ and its local site-specific component κ 0 for two seismological stations within the Zagreb seismic zone: Puntijarka-PTJ (situated on the peak of the Medvednica Mountains) and Zagreb-ZAG (situated in the city of Zagreb).
Near-site attenuation kappa (or site kappa, κ 0 ) we estimated using traditional linear least-square regression and the linear regression suitable for data with errors (their differences is taken as 5% uncertainty). Site kappa is estimated by restricting the distance range for kappa, to ensure that local site attenuation at zero distance is not overwhelmed by the path component for larger distances.
The main results and conclusions are: Comparison of the high-frequency coda attenuation values with frequency independent, Q c (f ) and Q(κ R ) for both stations can be considered comparable when looking at average values with the assumption that they are affected by the same regional characteristics; • The lowest κ values are spatially distributed within a few kilometres around the stations due to near-site effects; • Variations in observed κ distribution between lower and higher attenuation subset areas is clearly under the influence of local and regional geological structures that are distributed predominantly within regional active faults and complex tectonic structure.
Attenuation of seismic waves is one of the key factors in seismic hazard assessment. Additionally, it is important for quantification of earthquakes and plays a significant role in studies of seismic source and crustal structure. Thus, this study is a step towards the better defining of input parameters for seismic hazard and risk assessment for the wider area of the city of Zagreb.