Site Response Analysis of Anchorage, Alaska Using Generalized Inversions of Strong-Motion Data (2004–2019)

Anchorage, Alaska, is located in one of the most active tectonic settings in the world. The city and region were significantly impacted by the MW 9.2 Great Alaska Earthquake in 1964, and they were recently shaken by a MW 7.1 event in 2018. The city was developed in an area underlain by complex soil deposits of varied geological origins and stiffnesses, with the deposits’ thicknesses increasing east to west. Situated at the edge of the North American Plate, with the actively subducting Pacific Plate below, Anchorage is susceptible to both intraslab and interface earthquakes, along with crustal earthquakes. Strong-motion stations were installed across the city in an attempt to capture the variability in site response. Several previous studies have been performed to evaluate that variability but have not included larger magnitude events and have not benefited from the current density of instrumentation. The work presented here provides background information on the geology and tectonic setting of Anchorage and presents details related to the dataset and methods used to perform the site-response analysis. This study has collected strong-motion recordings from 35 surface stations across Anchorage for 95 events spanning from 2004 to 2019, including the MW 7.1 Anchorage Earthquake in 2018. The more than 1700 three-component recordings from those 95 events with moment magnitudes ranging from 4.5 to 7.1 were used to evaluate site response variability across the city. Using the Generalized Inversion Technique and a reference rock site, spectral amplifications were calculated and analyzed for frequencies between 0.25 and 10 Hz for each strong-motion station. The study results were used to develop contour maps at 1 Hz and 5 Hz, using logarithmic-band averages, to describe the variability of spectral amplifications at these two frequencies of interest. The results were also compared to geologic conditions across Anchorage, and the overlaying of different soil deposits can be seen to have an impact on the spectral amplification at the sites. The results of this study provide improvements on past microzonation studies and, using sensitivity analyses, offer support for the use of small and moderate earthquakes to evaluate spectral amplifications.


Introduction
Southcentral Alaska is one of the most seismically active regions of the world. It is located at the convergence of the Pacific and North American tectonic plates, which results in a variety of seismic sources, including deep subduction intraslab and interface earthquakes, as well as crustal earthquakes (Wesson et al., 2007). Southcentral Alaska also has approximately half of the state's population, primarily living and working in Anchorage. The 1964 Great Alaska Earthquake (M W 9.2), which is the second-largest earthquake recorded in modern history (USGS.gov, 2020), affected southcentral Alaska, including Anchorage, and resulted in significant ground failure. While Anchorage's population in the 1960s was approximately 44,000 (U.S. Bureau of the Census, 1963), the current population has grown to 291,000 (Census.gov, 2020); since the earthquake hazard has not reduced, the risk remains high.
The Anchorage area sits in a geologically complex area of southcentral Alaska. Mountain building at the edge of the subduction zone and episodic glaciation have affected the geology and created a basin of soft sediments composed of sands, gravels, and clay overlying a steeply dipping Tertiary rock (Combellick, 1999). An initial microzonation study was started in the 1990s by Biswas et al. (2003) within Anchorage, and several digital seismic strongmotion stations were installed across Anchorage to measure the variability of ground motions across the city Dutta et al., 2003;EERI, 2019;Franke et al., 2019;Martirosyan et al., 2002). Since then, several strong-motion stations have been added to the network by a variety of working groups. The network now consists of more than 30 surface strong-motion stations and a downhole array. Additionally, several buildings and bridges were instrumented to measure their structural response (e.g., Xiong & Yang, 2008;Yang et al., 2004).
This study presents the findings from the site response analysis of 35 surface strong-motion stations located in Anchorage, including the development of the database, methodology for performing the site response analysis for the selected stations, and a summary of the main results. This study represents a substantial improvement of the previous seismic microzonation studies for Anchorage, Alaska, as it uses data from more seismic stations representing various geological units and more and larger-amplitude strong-motion records than used in previous studies.
Microzonation studies like this consider the varied site response across Anchorage and allow the development of more resilient infrastructure that will suffer less damage and will recover more quickly after the next large earthquake. Anchorage serves a strategic community for a variety of services and operations such as the fifth busiest cargo airport in the world; United States military installations; and a deep-water port that receives 75% of the goods arriving in Alaska. Impacts on infrastructure due to earthquakes can be significant, and mitigation of the earthquake risk is critical to the continued operation of these facilities.

Tectonic Setting, Seismicity and Near-Surface Geology of Southcentral Alaska
Southcentral Alaska is on the North American Plate, with the Pacific Plate subducting below at a rate of 55 mm/yr (Haeussler, 2008). The depth of the interface between the two plates ranges from 25 to 40 km in the region around Anchorage (Hayes et al., 2018;Wesson et al., 2007). Both intraslab and interface earthquakes occur in this region. The most notable interface earthquake was the 1964 M W 9.2 Great Alaska Earthquake, which ruptured approximately 800 km of the Alaska-Aleutian megathrust (Freymueller et al., 2008). An example of a notable recent intraslab earthquake is the 2018 M W 7.1 Anchorage Earthquake, which caused extensive damage to infrastructure in southcentral Alaska (West et al., 2020). The 2018 Anchorage Earthquake occurred at a depth of more than 40 km, representing a common observation that the subduction-zone earthquakes are quite deep in this region.
In addition to the subduction-zone earthquakes, additional seismic sources include shallow crustal faults and folds (presented in Fig. 1 of Thornley et al., 2021a) in the Cook Inlet west and southwest of Anchorage and those in the northwest of Anchorage, including the Castle Mountain fault (Koehler et al., 2012). While the Cook Inlet crustal faults and folds are not likely to generate earthquakes of similar magnitude to the interface and intraslab earthquakes, they are much shallower than the subduction events and can potentially cause significant damage to the built environment. The Castle Mountain fault complex, which includes the Lake Clark fault to the west, can produce earthquakes up to M W 7.1-7.5 (Haessler et al.,2002;Wesson et al., 2007).
The geology of Anchorage varies significantly from east to west (Fig. 1). The Chugach Mountains, an accreted and lightly metamorphosed greywacke, border the city to the east (Wilson et al., 2012). Glacial valleys trend through the mountains in a northwest orientation and glacial outwash materials are found at the base of the mountains. The Chugach Mountains dip steeply to the northwest, and the sedimentary soil thickness reaches a depth of 500 m overlying Tertiary sandstone at the western edge of the city (Combellick, 1999;Glass, 1988;Schmoll & Barnwell, 1984). The soil overlying bedrock consists of a range of soils, from dense glacial outwash and till with shear wave velocities greater than 1,000 m/s (Thornley et al., 2019) to soft, cohesive lacustrine soil with shear wave velocities of 150 m/s (Updike et al. 1988). Erosional events related to several glaciation events have affected the thickness and lateral deposition of these different soils (Combellick, 1999;Ulery & Updike, 1983). Such marine transgressions and glacial advance and retreat cause soil heterogeneities, and hence significant variability in the amplification of earthquake ground motions across Anchorage.
The surficial geology consists of several soil units, including dense glacial till, glacial outwash, alluvium, and lacustrine deposits (Schmoll & Dobrovolny, 1972). The dense glacial till extends across the city and is near the surface in the eastern portion of the city. In the northern part of the city, the glacial till is overlain by glacial outwash deposits. Overlying the glacial till in the central and western portions of the city are alluvial deposits of varying thickness. The main soil unit found to have a significant impact on ground motions and site response is the Bootlegger Cove Formation (BCF). The BCF has several facies, including sand, silt, and clay. One of the sensitive clay facies is responsible for the significant ground failures in the 1964 Great Alaska Earthquake. The BCF is generally centered in the middle portion of the city, with the more sensitive clay facies located in the north and overlain by silt and sand in the south (Combellick, 1999;Updike & Ulery, 1986).

Anchorage Strong Motion Network and Strong Motion Data Used
At the time of the Great Alaska Earthquake, March 27, 1964, there were only two seismograph stations located within Alaska. These included one station in Sitka and one in Fairbanks. The Sitka Observatory was established in 1904, located 950 km southeast of Anchorage. In Fairbanks, the University of Alaska College Observatory was established in 1935. Fairbanks was the closer of the two stations to Anchorage, but it was over 400 km north of Anchorage. As a result, no local strong-motion records of the Great Alaska Earthquake exist (Hansen, 1965). In the years following this monumental event, numerous seismograph stations were installed and operated by the U.S. Geological Survey, U.S. Army Corps of Engineers, and others. By 1976 there were more than 40 stations across Alaska, including five stations in Anchorage (Procella, 1979).
As mentioned previously, in 1993, an effort was launched to provide a further understanding of the variability in Anchorage's seismicity through a multiagency seismic microzonation study. Twenty-two digital surface strong-motion stations were established across the Anchorage metropolitan area . Several studies have been performed using data from those stations, including those by Biswas et al. (2003) and Dutta et al. (2003).
In comparison to the current work, these previous studies were achieved with fewer stations (Alaska Science and Technology Foundation stations identified as K2XX Stations) and strong-motion records of small (M W 3.3) to moderate earthquakes (M W 6.5). The recent growth of the network, including the stations of the USGS National Strong Motion Program (NSMP Stations) and the Delaney Park Downhole Array (DPDA), dating from 2004, allows better characterization of the variability in response due to Anchorage's complex geology. In the early 2000s, the strong-motion data from the K2XX and NSMP stations started being collected at a single location, the Alaska Earthquake Center (AEC). Figure 1 presents the stations included in this analysis, and the latitude and longitude of each of the stations are included in ''Appendix A''. This study uses ground motions recorded by this network to update and refine the microzonation work performed previously by others.
The strong-motion data used in this study are of earthquakes between 2004 and 2019, including the November 30, 2018 M W 7.1 Anchorage Earthquake, which was recorded by 28 of the stations used in this study. Except for the DPDA data, the strong-motion records were primarily provided by the AEC. The AEC provided earthquake date, time, hypocentral location, and magnitude for the strong-motion records. The DPDA data were downloaded from the University of California, Santa Barbara (http://www. nees.ucsb.edu/), while the data from the M W 7.1 Anchorage Earthquake and larger aftershocks were downloaded from IRIS (https://www.iris.edu/hq/). The study includes recorded earthquakes with local magnitudes (M L ) greater than 4.5 and at epicentral distances less than 300 km. The cut-off distance of 300 km was chosen to include the M W 7.1 Iniskin Earthquake (January 24, 2016). A summary of the events used in this study is provided in ''Appendix B''.
The database included records that were characterized in local magnitude (M L ), body-wave magnitude (m b ), and moment magnitude (M W ) scales. The magnitude scales were unified to M W . Magnitudes in m b and M L scales were considered equal because of their relatively low magnitudes (typically below M 6) and based on research by Ruppert and Hansen (2010). These magnitudes were then converted to the M W scale using the correlation by Dutta et al. (2003) which was developed for southcentral Alaska (Eq. 1). Larger magnitude events in the database were characterized in the M W scale and did not require conversion.

Time History Processing
The strong-motion stations located across Anchorage primarily consist of Kinemetrics sensors. Many of these stations started with Altus K2 sensors and most have been upgraded over time to Basalt sensors. Most stations record at a sampling rate of 200 Hz, and, with the exception of the DPDA, all of the stations used in the study provide continuous data and are monitored by the AEC. The data used in this study was in the SAC format with a typical acceleration time history record length of around 300 s, including, in most cases, at least 20 s of pre-event (noise) data. The Seismic Analysis Code (SAC) (Goldstein & Snoke, 2005) was selected to process the strong-motion data.
The acceleration data was processed by removing the mean and trend using a least-squares fit (Goldstein & Snoke, 2005). The SAC function that removes glitches or irregularities in the records caused during data acquisition was also applied, where glitches refer to artificial spikes caused by samples that are orders of magnitude higher than the surrounding data (Goldstein et al., 2003;Helffrich et al., 2013). A fourth-order Butterworth bandpass filter was applied to the acceleration time histories, and various high and low-frequency limits were evaluated. Methods presented by Akkar and Bommer (2006) were utilized to verify the low-frequency (long-period) filter cut-off selection. Considerations for high-frequency (short-period) cutoff filters used the methods described by Douglas and Boore (2011).
The frequency range of interest for this engineering study is between 0.25 and 10 Hz, which matches the range of fundamental frequencies of structures within the city. Visual screening of records was performed to evaluate a reasonable bandwidth filter that could be applied to most records and maintain the frequency range of interest for this study. The quality of the recordings was evaluated using the signal-to-noise ratio (SNR), where the Fourier amplitude spectrum of the recorded signal and the preevent noise were calculated. An example time history and its signal and noise Fourier amplitude spectra are shown in Fig. 2. Each record was evaluated and those with a SNR of less than 3 within the frequency range of interest (0.25-10 Hz) after filtering between 0.1 and 30 Hz were excluded from this study. This practice reduced the number of available threecomponent records to 1727 from 95 events, which is 70% of the original dataset. The unused records were generally from small deep ([ 40 km) subduction earthquakes and were typically more than 100 km from Anchorage, which resulted in a SNR less than 3 within the frequency range of interest.
Surface waves can increase the uncertainty in the results of site amplification studies like this one because they have different characteristics to S waves (Bonilla et al., 1997;Oth et al., 2009). To reduce the impact of surface waves on the results, a window of 10 s was used starting at the S-wave arrival. The 10-s S-wave window was selected based on a visual evaluation of the records to capture the high-amplitude portion of the records. With the selection of a general window, there is a possibility of including some of the surface waves of some records, especially for distant lower magnitude events. Varying the window length between 5 and 15 s did not, however, have an observable effect on site amplification when only the lower magnitude events were included. The noise window was selected as the portion of the record before the P-wave arrival and was generally longer than 10 s, depending on the record. The P-and S-wave arrival times were chosen by visual inspection of the three-component records used during the time-history processing. An example of a horizontal acceleration time history with noise and S-wave windows is presented in Fig. 2. It should be noted that large low-frequency site responses can be associated with surface waves, particularly in sedimentary basins such as Anchorage, and hence the focus of this study on only S waves may be underestimating the low frequency response at some sites. Thornley et al. (2021a) presents the distribution of earthquakes used in this study and identifies the crustal, interface, and intraslab events, as determined by using the depth and geometry of the subduction zone defined by Hayes et al. (2018). Figure 3 presents the epicentral distance versus the magnitude of the events in the database. Note that the DPDA Station 8040, with latitude and longitude given in ''Appendix A'', was chosen as the point in Anchorage used to define a general epicentral distance from the city, although site-specific distances were used in the analyses. The hypocentral distance by azimuth is presented in Fig. 4 showing a concentration of earthquake events in two azimuthal ranges, 225°-250°and 330°-360°. A sensitivity analysis is presented below, evaluating the azimuth effects on the results of this study. The database consists primarily of M W 4.5-5.5 events (Fig. 5a). Additional information describing the distribution of data across epicentral distance, depth, and azimuth from Anchorage is presented in Fig. 5b-d. The station with the most recordings (70 in total) was K213 and the station with the fewest (only 6) was K207. The 35 stations had an average of 48 records each. The number of events recorded at each station and used in this study is presented in Fig. 6.

Site Response Evaluation
Not all events were recorded at every station either due to a low SNR or because the station was not functioning at the time. To include as many strong-motion stations as possible the Generalized Inversion Technique (GIT) was used to evaluate the spectral amplifications (where it is assumed that spectral amplifications provide an indicator of site response) for a site (Andrews, 1986;Castro et al., 2004;Priolo et al., 2019). The GIT methodology allows for gaps in the completeness of the dataset, providing a way to incorporate as many sites as practicable.

Figure 2
Example time history plot a. where the pre-event noise window is indicated in red and the S-wave window is in grey and b. the Fourier amplitude spectrum of the filtered signal (in black) and the noise (in red)

GIT Background Information
The use of GIT in site response analysis has become common since the original article by Andrews (1986), e.g., Parolai et al. (2000), Dutta et al. (2003), Oth et al. (2009), Bindi et al. (2017, and Laurenzano et al. (2019). Constraints applied to the source, path, and site terms in the spectral domain allows the indeterminate matrix to be solved. One constraint, such as the selection of a reference site, provides a simple way to constrain the matrix within microzonation studies (Klin et al., 2018). The main assumption of GIT is that in the frequency domain, the logarithmic value of the observed amplitude spectrum at any site due to a source is the logarithmic summation of the source, site, and path spectra. For this study, GITANES (Version 1.3) developed in MATLAB by Klin (2019) was used. Selected SAC format three-component time histories described in previous sections are loaded into GITANES. The program uses the data in the headers to select the Swave portion of the time history (10 s as shown in Fig. 2) and performs a fast Fourier transform of the data. The logarithmic value of the spectral amplitude of the S-wave data can be expressed as: where the amplitude spectrum of the recorded S-wave data, A, is a function of the hypocentral distance, r, and frequency, (f). So i (f) is the source for the ith source term, SI j (f) is the site term for the jth site, and P is the path effect describing the attenuation of the source hypocentral distance r ij for the ith earthquake event to the jth site. The path term (P) can be written as: Taking the logarithm of the path term: where r ij and f are the distance and frequency, respectively, and c is the geometric spreading coefficient. The parameters V S , Q 0 , and g are the average shear-wave velocity, quality factor at 1 Hz, and exponent, respectively. One of the tradeoffs when using GITANES is that the path terms defined in Eq. 4 are constrained. For this study, the average shear wave velocity is assumed to be 3.2 km/s (Table S1 of Silwal & Tape, 2016), due to the range of shallow to deep earthquakes. The quality factor, Q 0 , of 150 was used for this study and g was set at one. Several regional studies, including Boore (2013), McNamara (2000), and Stachnik et al. (2004) show the values of Q 0 often to range between approximately 100-300 for various models of surface, coda, and shear waves. A sensitivity analysis was performed to select the Q 0 and g values and it was found that varying them had an insignificant effect on spectral amplification. The results of the sensitivity analysis for Station 8040 are presented in Fig. 7, as an example of a typical result. It has been observed that the choice of Q 0 does not significantly affect the results of site response studies such as this one. Still, it does impact source spectra estimates (Parolai et al., 2000), which is supported by the results of the sensitivity analysis performed for this study.
The geometric spreading coefficient, c, was assumed equal to 1 (spherical waves) for distances less than 100 km and 0.5 (cylindrical waves) for distances greater than 100 km, which is similar to the assumption of Dutta et al. (2003) and used for earthquakes with epicentral distances of 75-500 km in a southcentral Alaska study by McNamara (2000). The geometric spreading coefficients are also recommended by Havskov and Ottemoller (2010) when using simplified geometric spreading parameters, such as those used in GITANES. It is noted that the change in the geometric spreading coefficient should theoretically also be reflected in the time-window chosen to capture the strongest portion of ground motions (fixed at 10 s to capture the direct S-waves). As shown in the various sensitivity studies, the impact of making such a change would be negligible so it was not done for simplicity. The distance, r ij , was collected from header information for each event. A sensitivity analysis, using a range for each coefficient and visual evaluation on the effects of the coefficient values on the site results, was performed before settling on these values.
Rewriting Eq. 2 as: logA And substituting Eq. 4: For a particular frequency we denote: where d ij is the weighted vector with path spectrum Vol. 179, (2022) Site Response Analysis of Anchorag 507 added to the amplitude spectrum and log So i (f) = so i and log SI j (f) = si j . Equation 6 can be expressed as: The matrix form of Eq. 8 can be expressed as: where G is the sparse matrix that relates m, the model vector of unknown source and site terms (logarithmic values), to d, defined in Eq. 7; and S is the row matrix that is appended to matrix G with reference site terms. These matrices are solved independently for each frequency of interest. As mentioned above the path term is constrained in GITANES but does not have an appreciable effect on the site response results. While the source spectra results are not the primary focus of this study, the source spectra have been evaluated to understand the effects of path terms Q 0 and g and to verify their appropriateness and physical meaning. The sensitivity analysis results of Q 0 and g on the source spectrum for the M W 7.1 Anchorage Earthquake are presented in Fig. 8. The variability of the source spectrum with respect to Q 0 , holding g at unity with similar results for Q 0 ranging from 150 to 250 (Fig. 8a). Figure 8b shows that the variation of g using values between 0.5 and 1.2, while maintaining Q 0 at 150, produces little impact for g ranging between 0.75 and 1.2. These sensitivity results support the use of Q 0 and g values described in the site amplification discussion above.
As shown in Fig. 9a, source spectra from the 95 earthquakes have been plotted. The source spectrum in black at the top of Fig. 9a is replotted in Fig. 9b using Q 0 of 150 and three values of g, including 0.5, 0.75, and 1.0. Figure 9b also presents Brune's omega square model (Brune, 1970) with stress drops of 2 and 5 MPa, and a seismic moment of 5.02 9 10 19 Nm to evaluate the source spectral shape. Estimates of stress drop range from 2.7 to 3.9 MPa and a seismic moment of 5.02 9 10 19 Nm by Liu et al. (2019) for the M W 7.1 Anchorage Earthquake. The results presented in Fig. 9b indicate the Brune's omega square model fits the source spectral shape for this event well, further supporting the path parameters selected for this study, given the range of earthquakes in the database. Further discussion related to calculated source spectra results has not been included in this paper for brevity.
The results of the GIT provide spectral amplification functions (SAFs) for the two horizontal components of each station from the two orthogonal components based on a selected reference station. In this study, the K216 surface station was used as the reference station. This station is located at the eastern edge of Anchorage and is in the Chugach Mountain Range on a rock outcrop. The site is underlain by glacially-carved and metamorphosed greywacke rock and is within the vicinity of the other stations. It is the best reference station available for the network as it is the stiffest site with respect to its shear-wave velocity. The other strong-motion stations are underlain by soil and were not deemed to be more effective than K216 as a reference station. The frequencies of interest for this project range from 0.25 to 10 Hz for reasons noted earlier. When evaluating the range of frequencies, Zhu et al. (2021) suggests evaluating the site-to-site variability (u S2S0 ) which identifies the parametric and modeling uncertainty and is independent of the reference station selection. Site-to-site variability is the standard deviation of the natural logarithm of the spectral amplifications. As shown in Fig. 10, u S2S0 rises sharply and the variability becomes more significant after 10 Hz indicating lower uncertainty between 0.1 and 10 Hz. Additionally, the spectral amplification at several stiffer soil sites shows ratios below unity at frequencies greater than 7 Hz, suggesting that the K216 station amplifies seismic waves above 7 Hz, which means it may not be an ideal reference site for the entire range of frequencies used in this study. Similar findings were reported by Martirosyan et al. (2002). The horizontal-to-vertical spectral ratio (HVSR) of the strong-motion stations were calculated in GITANES and shown in Fig. 11. Further discussion of the HVSR calculations for Anchorage strongmotion stations is presented in Thornley et al. (2021b). Lower horizontal site amplifications are observed at a site when HVSR is near unity. The K216 station provides generally the lowest bound HVSR ratio among the stations (black line in Fig. 11). As an example, the blue line in Fig. 11 presents the HVSR ratio for Station 8040. There are significant horizontal amplifications at Station 8040, showing why it would not be an effective reference station. However, at frequencies greater than 7 Hz, K216 diverges more significantly from unity and has the highest HVSR of all the stations, indicating site amplification complexities at this site above 7 Hz. The mean spectral amplification values for the two orthogonal horizontal component SAFs and their standard errors are provided as output from GITANES. An averaging technique, as proposed by Goulet et al. (2018), was applied to the resulting orthogonal horizontal site amplification ratios to calculate a single site amplification ratio at each station that is independent of instrument orientation. The equation: in which the EAF is the average site amplitude spectrum calculated at each frequency of interest (f ), and SAF E-W and SAF N-S are the site amplification functions for the east-west and north-south orthogonal horizontal components, respectively. This is similar to what is done for Fourier amplitude spectra (FAS) for orthogonal components within groundmotion models (Bayless & Abrahamson, 2019). From this point forward the spectral amplifications for each site have been calculated using Eq. 10.

Sensitivity of Model
A wide variety of earthquakes are used in this study, ranging in magnitude, epicentral distance, depth, and source type. An effort has been made to better understand the impacts of these elements in the dataset on the results of the site response analysis. As an example, more than 70% of the earthquakes in the database are smaller than M W 5.0 (Fig. 5a). The dataset was subdivided into earthquakes smaller than M W 5.0 and those equal to or greater than M W 5.0. The two subdivided sets were evaluated using GIT, and the results were compared. An analysis of the residuals, R k,l , was performed to evaluate the bias one dataset has compared to the other by using Eq. 11, where AmpA is the larger dataset, which in this case is the events less than M W 5.0, and AmpB is defined as the smaller dataset, with k and l being the number of the recording station and frequency, respectively. The results of the magnitude sensitivity analysis are presented in Fig. 12, with the average and the bounds of ± 1 standard deviation indicated.
The results of the magnitude analysis suggest that there is some influence of the lower magnitude data set on the results; however, zero is within one standard deviation for most of the frequencies of interest. Only for 1.8-3 Hz and 4.1-5.6 Hz is the mean more than one standard deviation above zero.
The residual analysis was also performed for epicentral distance where the division is between earthquakes less than 100 km (36%) from Anchorage versus those at distances between 100 and 300 km (64%). Sensitivity analysis considering the focal depth was also undertaken, especially given the range of earthquake depths is from 10 to 189 km. A division point was selected at a depth of 50 km, where 61% of the events were at depths greater than 50 km, and 39% are less than 50 km. It was found that, in general, there was no significant impact of epicentral distance or depth.
An additional analysis considered the impact of azimuth on the results. As shown in Figs. 4 and 5d, there are two predominant azimuthal zones in the dataset. The first is the earthquakes to the southwest of Anchorage with azimuths between 225°and 250°, and the second includes the earthquakes northnorthwest of Anchorage with an azimuthal range of 330°-360°. It was found that there was minimal impact related to whether the earthquake originated in either of these locations. Figure 13 presents the results of the residual analysis for the northern azimuthal range versus the events not in that range. Similar results were observed for the southern azimuthal range. The potential biases generally correspond to less than 30% because common (base 10) logarithms are used (10 0.1 = 1.26). It is noted that, although this potential bias is relatively minor, more data from a variety of azimuths, magnitudes and distances would better constrain the results. While our conclusions are not affected by this analysis, these effects should be considered for future studies if a range of azimuths cannot be achieved. Additional sensitivity analyses were performed to evaluate the impact on the results of event type, the November 2018 M W 7.1 earthquake, and other factors. These analyses suggest that the results are not sensitive to these factors. It is interesting to note that there were likely nonlinear site effects at some of the stations that recorded the November 2018 M W 7.1 earthquake (Thornley et al., 2021a). To evaluate the effects of this event on the overall results of the study, the GIT was repeated on the full database and the full database without the data from the M W 7.1 event. The two sets of results were compared for each strongmotion station where the M W 7.1 event was recorded and it was found that there was very little change (less than 4%) in the overall results, which is attributed to the large number of events in the database. Figure  Site Class D sites showed nonlinear response from the M W 7.1 event; however, none of the stations show any significant difference between the database with and without the M W 7.1 event (less than four percent difference at any frequency).

Site Response
As described above, there is significant geologic variability across Anchorage, which as shown by others including Souriau et al. (2007) can impact site effects across a region. That geologic variability results in nonuniform site response, where different areas of the city behave differently during strong shaking. There is a need to understand and model this variability so that better risk-mitigation decisions can be made by policymakers, engineers, and emergency response teams. The results from the GIT provide insights into those differences. For example, as shown in Fig. 15, there are significant differences in site response, as indicated by spectral amplification, between a site (K209) located in east Anchorage where the near-surface soils are dense glacial till overlying rock versus the DPDA site (8040), located in the center-north portion of the city, where there are more than 45 m of alluvium and BCF silts and clays overlying glacial till (Thornley et al., 2019). The spectral amplitude results for each station are included in ''Appendix C''. Twelve of the 35 stations presented in this study were evaluated by Dutta et al. (2003) using the GIT methodology, but with a different database of earthquakes, as described earlier.
The results of that study for those stations are compared with the spectral amplification results from this study in ''Appendix C''. In general, there is a good fit between the results from Dutta et al. (2003) and this study, where the average spectral amplification results of the past study stay within the standard error of the results of this study. There are three stations, K203, K220, and K221 where there is a difference at one frequency, but the results across the other frequencies at these stations also fit well within the standard error of the results from this study.
While acknowledging the geologic variability across Anchorage, it is important to identify areas of similarity among specific frequencies of interest. The focus of this study was to evaluate the spectral amplitudes for a range of frequencies between 0.25 and 10 Hz, with the understanding that frequencies above 7 Hz may be artificially low because of potential site amplifications at K216. To more clearly show similarities and differences across Anchorage the 1 Hz and 5 Hz frequencies were considered for further evaluation. These two frequencies of interest were selected because of their use in engineering studies; these frequencies are often used in engineering studies to define design response spectra (American Society of Civil Engineers, 2017). As shown in Fig. 16, there is a chance that the value at exactly 1 Hz may not be representative of the values above and below 1 Hz. A logarithmic-band average of the computed spectral amplifications was, therefore, used between 0.5 and 2.5 Hz for the 1 Hz band and logarithmic-band average between 4 and 6.5 Hz for the 5 Hz band, which is also shown in Fig. 16.

Spectral Amplification Variability Across Anchorage
The strong-motion stations are roughly uniformly distributed across Anchorage, from shallow soil over the rock in the eastern portion of the city to the deep soil deposits in the mid to western portions of the city. To highlight the changes in spectral amplification across the city for the specific frequency ranges, two contour maps are drawn. The contours were created using GIS software by applying an inverse distance weighting with a power of three. No smoothing was applied. The contours have not been modified to account for known changes in surficial geology between stations. Due to the high density of the strong-motion stations across Anchorage the variability of subsurface geology is generally accounted for by these contours. Trends in the contours can be seen, and initial conclusions are drawn below; however, it is acknowledged that location-specific responses may be different due to variability of surficial geology. Contour maps for 1 Hz and 5 Hz are presented in Figs. 17 and 18.
There is a general trend across Anchorage for the 1 Hz band-averaged spectral amplification. The spectral amplification of the stations compared to the reference station is less than 1.5 in the eastern portion of the city and increases to the west. Combellick (1999) presents a line that estimates the eastern extent where the BCF (primarily silt and clay) becomes less than 10 m thick. That estimated line is shown in Figs. 17 and 18. The thickness increases to the west and then decreases in thickness again at the western edge of the city as indicated by Ulery and Updike (1986). In general, the spectral amplification at 1 Hz appears to increase as the BCF increases with depth and then decreases again at the western edge of the city.
Boreholes near Stations K221 and K208 indicate cohesive soil depths of 25 and 20 m, respectively. The amplifications observed at 1 Hz for Stations K208 and K221 are similar in both this study and that of Dutta et al. (2003). Similarities like this show that, despite two different datasets, the estimated lowfrequency amplification tends to be similar. However, Station 8039, a newer station not included in Dutta et al. (2003), in southwest Anchorage, has the highest spectral amplification (3.9) at 1 Hz. Ulery and Updike (1983) indicate the BCF has deeper regions of cohesive soil in the vicinity of this site. Station 8039 is located in an area that has approximately 55 m of cohesive BCF soil below it. Interestingly, Station 8041 is located in another area with a deeper cohesive soil deposit of 45 m and a band-average spectral amplification of 3.5 at 1 Hz. The results at the western edge of the city tend to show lower amplifications (e.g., K220 in this study has a spectral amplification of 2.1 while Dutta et al. (2003) shows a spectral amplification of 2.7), which is also where the soft portions of the BCF become thinner and stiffer soil becomes more predominant (Updike & Ulery, 1986). Based on the findings from this study, the thickness of the BCF directly impacts the spectral amplification at 1 Hz, indicating that for sites located in central to western Anchorage there is a spectral amplification of more than a factor of two.
For the higher frequency band surrounding 5 Hz, there is less of an east-west contrast (Fig. 18). On the east side of Anchorage, the spectral amplification is similar to the reference site. There are two areas where the spectral amplification is above 2. The northern area is located in the Chester Creek basin (identified in Fig. 1) where the stations are located in an area with a mixture of alluvium and glacio-fluvial surficial soil. The second area is found in southwestern Anchorage where the surficial soils tend to consist of lacustrine or eolian silt, fine sand, and clay related to glacio-estuarine deposition (Combellick, 1999). The additional strong-motion stations that have been included in this study offer a significant improvement in the 5 Hz results compared to Dutta et al. (2003). This is particularly true in southwestern Anchorage where there are more significant amplifications (greater than 4) than were estimated by the previous study (approximately 2). The differences presented in this study are primarily related to the additional strong-motion stations installed in key geologic conditions, further underscoring the importance of establishing and maintaining dense networks in urban areas with complex geologic conditions.

Geological impacts on Spectral Amplification
To further evaluate the impact of geology, the spectral amplifications are plotted and compared with geologic cross-sections developed by Combellick (1999). Figures 19 and 20 present the spectral amplification for several stations along a north-south and east-west geologic section, respectively. The locations of the cross-sections are indicated as dashed lines in Fig. 1. The glacial till is represented in blue, BCF soil in green and pink indicates the estimated bedrock. Other colors represent the mixture of alluvial and other surficial soil deposits. The spectral amplification plots are centered at the locations of the strong-motion stations located along the cross-section. In areas where glacial till is shallow there is a distinct difference in the shape of the spectral amplification, especially when compared to areas where the glacial till is deep.
A third cross-section is presented in Fig. 21. This represents the NW/SE cross-section by Combellick (1999), with its location identified in Fig. 1. It offers an excellent example of the impact of soil thickness and change in response based on the type of soil. The spectral amplification starts increasing significantly in the northwest portion of the cross-section at stations located over the BCF (green) subsurface soil. This contrasts with the stations in the southeast, where the stations are located with dense to very dense glacial till close to the surface.

Conclusions
The results of this study build on past studies and provide further insight into the variability of site response through an evaluation of the spectral amplification of ground motions at strong-motion stations in Anchorage. The dataset used in this study is independent of previous studies and includes various events from magnitudes M W 4.5-7.1. There is a strong geological effect on the response of different areas in Anchorage. As seen in Figs. 17 and 18, there will be distinct impacts on structures with varying frequencies of response in future earthquakes. The results provide additional insight into which areas in Anchorage are susceptible to higher ground motions, due to the subsurface (and potentially other factors, such as topography) conditions present.
The findings of this study further indicate that the BCF substantially impacts site response, especially at frequencies around 1 Hz. The results fit well with the geologic observations related to the extent and variable thickness of that deposit. Furthermore, the results from the 5 Hz study also match well with the surficial geology, where in the southern portion of the city, the BCF is overlain by silt and sand and was in a different depositional environment compared to the northern and eastern portions of the city. As presented in the discussion of the model sensitivity, the inclusion of a wide variety of data appears not to bias the study results. While it is acknowledged that the dataset would be better if larger and closer earthquakes were included, no such data exists at this time. The sensitivity study indicates that, for several different parameters, impacts on the results were not observed. This gives confidence that the results can be used to evaluate the potential response related to larger and more damaging events.
In general, the current study agrees well with previous studies, including Dutta et al. (2003). This and the previous work are based on very different earthquake datasets and use different software, so in that way, this study validated the work performed by the previous study. This study does provide additional refinement and identifies areas where higher amplifications have been recently measured through the inclusion of additional strong-motion stations. Further investigation will be required to evaluate the impact of the site amplifications on the built environment-this study provides support to perform those analyses.

Funding
Not applicable.

Availability of data and material
Referenced in manuscript.

Code availability
Referenced in manuscript.

Conflicts of interest Not applicable.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.