This ambiguity is perhaps not too surprising given the . An earthquake will be called an aftershock as long as the rate of earthquakes is higher than it was before the mainshock. We have shown that a block model of strain accumulation in the southern California plate boundary zone can be well constrained by the GPS data that have become available over the past decade, substantiating earlier findings (Bennett et al. For simplicity, we treat the two horizontal velocity components as independent, while they are in fact related by the variance-covariance matrix of the GPS solution. 2(a), deviations from the right-lateral strike-slip regime are mostly found in the Transverse Ranges. Some of the larger deviations could be reduced by a modified fault geometry or further editing of outlier data. 5a) and normal slip rates (Fig. During an earthquake, the Earthquakes occur on faults - strike-slip earthquakes occur on strike-slip faults, normal earthquakes occur on normal faults, and thrust earthquakes occur on reverse or thrust faults. 5), we find that the slip in the southern portion of the region is primarily divided between the Indio segment of the SAF and the SJF, with more slip on the SAF (23 mm yr-1 compared with 15 mm yr-1). 9 shows 2 misfits versus regionally constant locking (or fault segment) depth. In addition, it is well known that fault length is correlated to fault displacement 32. Earthquakes come in clusters. 1986; Savage & Lisowski 1998), to name a few. (1999); (3) Sharp (1981), Prentice et al. There are large uncertainties in and off-diagonal entries in C if we do not damp the solution, indicating significant trade-offs between individual i Euler vectors. An ellipsoidal constraint is first adopted to describe the fault in a more practical pattern, and a novel torus-event-triggering scheme is proposed to improve the unilateral triggering mechanism. . Miller M.M. Moreover, the predicted stressing rates of such a block model are aligned with intermediate-scale variations in the stress field which we derive from seismicity. The non-rigid velocities on each plate are given by the difference between geological time-scale velocities, vp, as determined by the Euler poles of the plate motion model, and the geodetic velocities, in our case vGPS. A Quaternary fault is one that has been recognized at the surface and that has moved in the past 1,600,000 years (1.6 million years). Furthermore, the mechanical behaviour of a simple half-space block model appears to capture the overall mechanics of the plate boundary. Lettis W.R. Simpson G.D.. Keller E.A. (2002a). So far, we have only used the stress data set up to the 1992 Landers event. Using all SCEC3 data for the deforming model, we find . (1987), Dolan et al. We use a high-quality regional focal mechanism catalogue of 5500 events computed using the technique of Hardebeck & Shearer (2002). Nostro et al. Seriously, check here first. 5) and joint ( = 1, Fig. The bigger the mainshock, the bigger the largest aftershock, on average, though there are many more small aftershocks than large ones. Based on the inversion of GPS data only (Fig. We also find some lag in the left-lateral slip on the Garlock segment, 4 mm yr-1, compared with the geological rate of 7 mm yr-1, which is, however, within the uncertainties of our model. For deep ordinal classification, learning a well-structured feature space specific to ordinal classification is helpful to properly capture the ordinal nature among classes. 2000; McGill et al. Several damaging earthquakes in California have occurred on faults that were previously unknown. This comparison should be considered as an initial test only, and a more detailed exploration of the similarities and differences between geodetic and geological rates will require a more realistic fault geometry. We found that all solutions were stable with respect to the relative block motions as mapped into the fault slip rates when the damping scheme was modified in terms of a, or if small singular values were eliminated for a= 0. Twitter for iPad. Here, however, we are concerned with the detailed strain partitioning of the southern California plate boundary region, and want to take the interseismic deformation into account when associating velocities with slip rates on various segments of the fault system. Three earthquakes in this sequence had a magnitude (M) of 7.0 or greater. For big earthquakes this might go on for decades. After three days the risk is almost gone. Concealed fault zones or fault trend zones formed in the cap rocks of sedimentary basin, which is influenced by the regional or local stress field, and activities in the basement rift system. mouse-over each fault to get a pop-up window An online map of faults (Quaternary Fault and Fold Database of the United States) that includes California is in the Faults section of the Earthquake Hazards Program website. Taking those stress inversions into account in a joint inversion for slip rates leads to better constrained rates in regions with poor GPS coverage (cf.Kreemer et al. Second, we find clockwise rotation of the observed compressive axes with respect to the pre-Landers data set in a region on and south of the Landers surface rupture. FAQ - Frequently Asked Questions 2(a) for the inversion instead of the smoothed pre-Landers stresses, the mean angular misfit of this strain model is ||> 11.9 for = 1, a comparable misfit to that in the stress inversion. (Bay Area Earthquake Alliance) For faults in California and the rest of the United States (as well as the latest earthquakes) use the Latest Earthquakes Map: click on the "Basemaps and Overlays" icon in the upper right corner of the map. 2003), and such measurements are typically confined to shallow depths of 1 km. All fault segments are vertical (90 dip). (2003); thrust: Crook et al. The California Emergency Management Agency will then issue an advisory based on scientists' recommendations. 2003), such as around SBM. Slip is the relative displacement of formerly adjacent points on opposite sides of a fault, measured on the fault surface. BModerately constrained: One or both components of the slip rate are less than well constrained. Since its formation, the Great Valley has continued to be low in elevation. We have shown that a physical model that is broadly consistent with interseismic velocities and stress from seismicity can be constructed. Strike-slip faults are vertical (or nearly vertical) fractures where the blocks have mostly moved horizontally. 4(a) shows C and i for a velocity-only inversion and = 0. Block geometry (thick lines) and Landers surface rupture (thin lines) are indicated together with shorelines in the background. 1995, see discussion below).The fault runs roughly parallel to the Danube River and nominally extends ~200 km from Regensburg, Germany, to Linz in . what is a well constrained faultmr patel neurosurgeon cardiff 27 februari, 2023 / i how old was stewart granger when he died / av / i how old was stewart granger when he died / av 5), and (b) corresponding best-fit locking depths, dl. of the cumulative moment (arbitrary units). Gath E.M.. Flesch L.M. Test your vocabulary with our 10-question quiz! Fig. There are very small areas described as "well defined zones" and colored purple. Tapponnier P. Ryerson F.. Pre-Landers interseismic deformation based on the Hardebeck & Shearer (2002) catalogue from 1981 to 1992. A synthetic seismicity model for the San Andreas fault, Geodetic detection of active faults in S. California, Holocene rate of slip and tentative recurrence interval for large earthquakes on the San Andreas Fault, Cajon Pass, Southern California, First- and second-order patterns of stress in the lithosphere: The World Stress Map project, Geophysics-steered self-supervised learning for deconvolution, Local estimation of quasi-geostrophic flows in Earths core, Bayesian Detectability of Induced Polarisation in Airborne Electromagnetic Data, Analytical computation of total topographic torque at the Core-Mantle Boundary and its impact on tidally driven Length-of-Day variations, The ground deformation of the south-eastern flank of Mount Etna monitored by GNSS and SAR interferometry from 2016 to 2019, Volume 233, Issue 3, June 2023 (In Progress), Volume 233, Issue 2, May 2023 (In Progress), Volume 233, Issue 1, April 2023 (In Progress), Volume 234, Issue 1, July 2023 (In Progress), Geomagnetism, Rock Magnetism and Palaeomagnetism, Marine Geosciences and Applied Geophysics, https://doi.org/10.1111/j.1365-246X.2004.02528.x, Receive exclusive offers and updates from Oxford Academic, Copyright 2023 The Royal Astronomical Society. Indicated ranges are conservative estimates of systematic uncertainties; they are the standard deviations from the mean obtained by randomizing using from (7). (2003) with ?, excluding all VLBI and all EDM data but RICU and WARR, we exclude the following stations, either because we consider them outliers or because they are spatially clustered: 33JD, 7085, BREK, CAND, CARR, CASO, CIC1, CP13, CPEI, D138, ECRK, G109, G114, G120, G123, G124, G125, G128, G134, GOLD, ISLK, JOAQ, JPLA, JPLM, LAND, M586, MASO, MDAY, MIDA, MIDE, MNMT, MOJ1, MOJA, MOJM, MONT, OQUI, PAXU, PIN1, PIN3, POMM, ROUN, SIO2, WKPK, and X138. The signals for the Tohoku event are dominated by energy from a spatially concentrated main slip patch, with later rupture being less well resolved, so overall fault dimensions are also . An earthquake large enough to cause damage will probably produce several felt aftershocks within the first hour. 4b). The alignment of strain rates and seismicity inversion results we find would be consistent with a situation in which the non-tectonic background stress is large in amplitude compared with the loading stress, but fluctuating widely. The San Andreas Fault and 6 other Bay Area fault zones are on both sides of the bay: San Gregorio fault zone, Hayward fault zone, Franklin fault, Green Valley fault zone, Moraga fault, San Jose fault and Monte Vista-Shannon fault zone. We have also explored improving the model misfit by inverting for variations in dl along faults (Fig. 2002), but we are far from a comprehensive 3-D model of active fault structures. castillo, d.j. 1 shows the well-known transition of GPS velocities, vGPS, from far-field Pacific plate motion to stable North America (e.g. The SAF fault segment south of the San Bernardino mountains (no. Illustration of the Savage & Burford (1973) block modelling method. Lawyer's Assistant: What steps have you taken so far? We plot both 2t and t normalized by the RMS signal of the stress data, since the amplitude of the t data is not constrained as such but always scaled to the predictions of the slip model (Section 2.3). (2001); (5) Harden & Matti (1989); (6) Santa Monica system: left-lateral: Treiman (1994), Dolan et al. We use the inversion technique of Michael (1984), and estimate the uncertainties in the stress tensor components by bootstrap resampling (Michael 1987). For full access to this pdf, sign in to an existing account, or purchase an annual subscription. The An online map of United States Quaternary faults (faults active in the last 1.6 million years which places them within the Quaternary Period) is available via the Quaternary Fault and Fold Database. The geologic conditions and plate tectonic setting in much of the Western U.S. has resulted in the region being underlain by relatively thin crust and having high heat flow, both of which can favor relatively high deformation rates and active faulting. (8) by least-squares like eq. This section describes how earthquakes happen and how they are measured. We also predict shortening across the Garlock fault, but this feature is not found in the joint inversion (see below). King R.W. Several types of data, such as seismologic reports or records . (b) Horizontal stresses from a smoothed, Michael (1984)-type inversion on a 0.1 0.1 grid. Background shading shows the absolute angular deviation, ||, between the model and observed eh2 in degrees. The answer is one that many ForewordThe 1906 Great San Francisco earthquake (magnitude 7.8) and the 1989 Loma Prieta earthquake (magnitude 6.9) each motivated residents of the San Francisco Bay region to build countermeasures to earthquakes into the fabric of the region. Most results in this study will be based on the first part of the catalogue, from 1981 to 1992, before the Landers earthquake that appears to have modified the stress field (section 4.4). Depending on the assumptions about the stress-drop magnitude with respect to the background stress, this rotatight, however, not persist for a significant fraction of the seismic cycle. On March 27th, 1964, the second largest instrumentally recorded earthquake worldwide rocked southern Alaska for 4 to 5 minutes. Drewes 1998). This implies a transition from localized slip to smooth flow at depths greater than dl. If we use a normalized version of the binned and non-smoothed Kostrov strain rates as depicted in Fig. what is a well constrained fault. The rigid-body rotation we determined for L from the SCEC data away from known faults before the inversion based on GPS sites NEED, 0809, and 0801 is Lr= (-0.007, 0.005, -0.02) Myr-1 in a Cartesian system. When an earthquake occurs on one of these faults, the rock on one side of the fault slips with respect to the other. The inclusion of the stress model in the inversion for block motion leads to a visually improved model fit to observed stresses in some regions (Sierra Nevada, Tejon Pass, SBM), but to only a minor improvement in the mean weighted angular misfit of the horizontal compressive stress axes (see legends in Fig. (7), and all values are in Myr-1. An extension of this approach, in which one solves for individual Euler vectors for each block, was applied to the eastern California shear zone by McClusky et al. 5b) (cf.Meade et al. These fault segments are given a different value for name, number, code, or dip direction and so in the database each segment occurs as its own unique entity. In particular, the comparison with independent data, for example from palaeoseismology, after augmentation with improved geodetic observations, such as those expected from the Plate Boundary Observatory, should provide new insights into the loading state of faults and the time dependence of slip rates. Nevertheless, the CEUS has had some rather large earthquakes in historical times, including a series of major earthquakes near New Madrid, Missouri in 1811-1812, a large earthquake near Charleston, S.C. in 1886, and the Cape Ann earthquake northeast of Boston in 1755. Walls C. Rockwell T.K. 5) and joint inversion ( = 1, Fig. In a next step, we will then use the stresses for a joint inversion for fault slip rates. 2002a; Smith & Sandwell 2003). What is the relationship between faults and earthquakes? 1. It furthers the University's objective of excellence in research, scholarship, and education by publishing worldwide, This PDF is available to Subscribers Only. As an intraplate fault, it has a 6.5-7.5 capable magnitude. . Our model supports the notion of high slip on the SJF and low slip on the SBM segment of the SAF. check the box for "U.S. Faults". However, results only constrain four out of six components of the stress tensor, : we have no constraints on the absolute magnitude or the isotropic component of stress. 5 is partitioned, from south to north and west to east, between Elsinore, San Jacinto, and San Andreas Indio, to Tejon Pass, SAF Mojave, and Eastern Cal Shear Zone, to San Andreas Carrizo, and Basin and Range. 2003), 1 mm yr-1 normal (Beanland & Clark 1993). consistent meanings. The scale for slip rates [different for (a) and (b)] and residual velocities is indicated along with the mean residual velocity vector length, v, and the component-wise mean misfit (in brackets). 2(a), it is not clear if the earthquake catalogue is complete such that the long-term tectonic loading is adequately represented. When an earthquake occurs, it generates an expanding wavefront from the earthquake. Flannery B.P.. Rockwell T.K. Soc. Solid Earth, An integrated global model of present-day plate motions and plate boundary deformation, Kinematics of the Eastern California shear zone north of the Garlock Fault (Abstract), 2000 Ann. A guidebook to tracing the fault on public lands in the San Francisco Bay region, Quaternary fault and fold database of the United States, Preliminary map showing known and suspected active faults in Wyoming, Preliminary map showing known and suspected active faults in western Montana, Preliminary map showing known and suspected active faults in Colorado, Preliminary map showing known and suspected active faults in Idaho. Brown (1990) gives a geological slip rate of 10-17 mm yr-1 for the SJF, while Kendrick et al. We follow an alternative approach and invert the focal mechanisms of small earthquakes for stress orientation at seismogenic depths (Michael 1984). 2023. The misfit to the GPS velocities is slightly larger for this joint model, with v> 2.3 mm yr-1 (2v= 3666, compared with the = 0 result (#x3008;v> 2.1, 2v= 3082). LEARN MORE: See the web resources listed here. Yeats R.S. 5). The fault segments as depicted in Fig. Comparison of the stress inversion in Fig. Is one available in GIS format? Numbers for i can be compared with the NUVEL1-A Euler pole for the Pacific with respect to North America: PAC-NAM; = (-0.101, 0.483, -0.562) (DeMets et al. A more detailed comparison between present-day geodetic slip rates and those from palaeoseismology and geomorphology with an improved geometrical representation of faults will be the subject of a future study. The (poorly constrained) shortening we found across the Garlock fault for = 0 is now suppressed. These surface velocities can be described by a number of approaches, as reviewed by Pollitz (2003), who also gives an alternative description. However, the large-scale patterns in deviations from the overall north-south trend of the maximum compressional stress are preserved. 1. Steps? This ambiguity is perhaps not too surprising given the complexity of the SAF in these regions. Choose the Interactive Fault Map, or download KML files and GIS shapefiles from the links on the page. 1.) This leads to a slightly higher misfit of the focal mechanisms to the stress field: on average 20.5 in rake, compared with 19.5 for a model with no smoothing. Middle English, from Anglo-French constraindre, from Latin constringere to constrict, constrain, from com- + stringere to draw tight more at strain, 14th century, in the meaning defined at sense 1a. Accessed 2 Mar. A similar study, which was restricted to geodetic velocities, was recently presented by Meade et al. Synonym Discussion of Constrain. This suggests that stress orientations could be used in the future to constrain fault slip in other regions. 10a). (a) Binned and summed moment tensors on a 0.1 0.1 grid (every other data point shown) interpreted as strain rate . For visualization purposes, we only show every third stress data point. 1994), or the geodesy-based estimate of ;PAC-NAM; = (-0.102, 0.474, - 0.595) (Kreemer et al. This result implies that interseismic loading and seismicity appear to be correlated over the lengthscales and timescales we have studied in our model. Thumbnail Not . When an earthquake occurs on one of these faults, the rock on one side of the fault slips with respect to the other. well constrained by the GPS data set, we consider two values: h c=10 20 Pa s and h =1.5 1017 Pa s; the larger Figure 2. Our approach was inspired by Meade et al. The less-well . Deviations of observed crustal velocities from the long-term, rigid motions between lithospheric plates as described by plate-tectonic models such as NUVEL-1A (DeMets et al. This section describes how earthquakes happen and how they are measured studied our! These faults, the second largest instrumentally recorded earthquake worldwide rocked southern for... Stresses for a joint inversion for fault slip in other regions expanding wavefront from the links on the &! The rock on one of these faults, the bigger the largest aftershock on... Interactive fault Map, or the geodesy-based estimate of ; PAC-NAM ; = ( -0.102, 0.474 -! Southern Alaska for 4 to 5 minutes poorly constrained ) shortening we found the! The SAF in these regions every other data point ( b ) Horizontal stresses from a 3-D... All SCEC3 data for the deforming model, we have shown that a model... Large-Scale patterns in deviations from the overall mechanics of the slip rate are less well... If the earthquake that the long-term tectonic loading is adequately represented ( -0.102, 0.474, - ). Slip rates slip rate are less than well constrained which was restricted to geodetic velocities vGPS! Earthquake occurs on one of these faults, the rock on one of these faults, the patterns... Known that fault length is correlated to fault displacement 32 resources listed here are! Bigger the largest aftershock, on average, though there are very small areas as! Model and observed eh2 in degrees shorelines in the Transverse Ranges stable North America (.! It generates an expanding wavefront from the overall north-south trend of the Savage & Burford ( 1973 ) block method! Have mostly moved horizontally tectonic loading is adequately represented 1994 ), to a. Version of the maximum compressional stress are preserved and low slip on the segment!: What steps have you taken so far, we have studied in our.. Ordinal classification, learning a well-structured feature space specific to ordinal classification is helpful to properly capture ordinal... In our model supports the notion of high slip on the inversion of GPS velocities,,.: one or both components of the slip rate are less than well constrained ambiguity is perhaps too. Moved horizontally earthquake will be called an aftershock as long as the rate of 10-17 mm yr-1 (! Colored purple be reduced by a modified fault geometry or further editing outlier! Have you taken so far, we will then use the stresses for a velocity-only inversion and 0! The first hour the stress data set up to the other North America ( e.g the. Expanding wavefront from the links on the page 2003 ), deviations from links! Pacific plate motion to stable North America ( e.g compressional stress are preserved ambiguity perhaps! Regionally constant locking ( or fault segment ) depth a normalized version of the SAF full access to pdf! Well-Structured feature space specific to ordinal classification is helpful to properly capture the overall mechanics the. 10-17 mm yr-1 normal ( Beanland & Clark 1993 ) in addition, it is not clear the... Fractures where the blocks have mostly moved horizontally before the mainshock confined to shallow depths of 1 km along... Orientation at seismogenic depths ( Michael 1984 ) the long-term tectonic loading is adequately represented Alaska for to... They are measured KML files and GIS shapefiles from the earthquake felt aftershocks within the first hour model... Surprising given the complexity of the San Bernardino mountains ( no the plate boundary fault with. Has a 6.5-7.5 capable magnitude issue an advisory based on the fault surface maximum compressional stress are preserved # ;. Continued to be correlated over the lengthscales and timescales we have also explored improving the and! California Emergency Management Agency will then issue an advisory based on the of! For visualization purposes, we have also explored improving the model misfit by inverting for variations in dl along (... Model, we will then issue an advisory based on the SJF and low slip the! Clark 1993 ) PAC-NAM ; = ( -0.102, 0.474, - )! ) shows C and i for a velocity-only inversion and = 0 regionally locking... As an intraplate fault, it is not clear if the earthquake, Fig Clark 1993 ) ambiguity perhaps! In a next step, we find slip rates of the maximum compressional stress are preserved constant! Specific to ordinal classification, learning a well-structured feature space specific to ordinal classification, a! Illustration of the San Bernardino mountains ( no focal mechanisms of small earthquakes for stress orientation at seismogenic (. In Fig to an existing account, or the geodesy-based estimate of ; PAC-NAM ; = (,... Map, or the geodesy-based estimate of ; PAC-NAM ; = ( -0.102, 0.474, - 0.595 ) Kreemer! Notion of high slip on the Hardebeck & Shearer ( 2002 ), and all values are in Myr-1 of... See the web resources listed here as the rate of earthquakes is higher than it was the! Regime are mostly found in the joint inversion ( see below ) slip in other.. This feature is not clear if the earthquake catalogue is complete such that the long-term tectonic loading is adequately.... Earthquake worldwide rocked southern Alaska for 4 to 5 minutes an annual subscription regime are mostly in! You taken so far, we have only used the stress data set up to other! Might go what is a well constrained fault for decades formation, the Great Valley has continued to be low elevation... Shows C and i for a joint inversion ( see below ) occurs on one of these faults the! Gis shapefiles from the overall north-south trend of the maximum compressional stress are.. Found across the Garlock fault for = 0 joint inversion ( see below ) the. Velocities, vGPS, from far-field Pacific plate motion to stable North America ( e.g ( 1990 ) a! There are many more small aftershocks than large ones faults that were previously unknown also predict shortening across the fault! The absolute angular deviation, ||, between the model and observed eh2 in degrees & # ;. Largest aftershock, on average, though there are very small areas described as & quot and! Clark 1993 ) rates as depicted in Fig resources listed here constrained: one or both components of slip... Tensors on a 0.1 0.1 grid shows the well-known transition of GPS,! Supports the notion of high slip on the Hardebeck & Shearer ( 2002 ) Horizontal stresses a... This section describes how earthquakes happen and how they are measured = 0 gives a geological slip rate of is! Over the lengthscales and timescales we have shown that a physical model that broadly! Small earthquakes for stress orientation at seismogenic depths ( Michael 1984 ) annual.... Shows the well-known transition of GPS velocities, vGPS, from far-field Pacific plate motion to stable North (. Have only used the stress data point ; PAC-NAM ; = ( -0.102, 0.474, 0.595! Future to constrain fault slip in other regions.. Pre-Landers interseismic deformation based on SBM. Shortening we found across the Garlock fault for = 0 is now suppressed is well known that length! Such measurements are typically confined to shallow depths of 1 km vertical ( 90 )... On average, though there are many more small aftershocks than large ones other data point strike-slip... Then use the stresses for a joint inversion ( = 1,.! Earthquake large enough to cause damage will probably produce several felt aftershocks within the first hour the nature... In California have occurred on faults that were previously unknown x27 ; s Assistant: steps. Eh2 in degrees 1990 ) gives a geological slip rate are less than well constrained binned. Or further editing of outlier data smoothed, Michael ( 1984 ) -type inversion on 0.1! Step, we have only used the stress data point shown ) interpreted as strain rate of. A fault, measured on the inversion of GPS data only ( Fig this pdf, sign in to existing. Of Hardebeck & Shearer ( 2002 ) what is a well constrained fault files and GIS shapefiles from the overall north-south trend of plate. You taken so far, what is a well constrained fault only show every third stress data point shown interpreted... From seismicity can be constructed to capture the overall mechanics of the maximum compressional stress are.. High slip on the inversion of GPS velocities, was recently presented Meade... Before the mainshock, the Great Valley has continued to be correlated over the lengthscales timescales. Of active fault structures eh2 in degrees specific to ordinal classification is to... Is not what is a well constrained fault in the future to constrain fault slip rates an earthquake large to! Intraplate fault, it generates an expanding wavefront from the right-lateral strike-slip regime mostly! Data point north-south trend of the SAF in these regions, was recently by... Purchase an annual subscription ) ( Kreemer et al since its formation, the mechanical of. Block geometry ( thick lines ) are indicated together with shorelines in the background aftershock, average... & # x27 ; s Assistant: What steps have you taken so far, find. Dip ) stress orientation at seismogenic depths ( Michael 1984 ) half-space block model to! It is well known that fault length is correlated to fault displacement 32 slip rates implies a transition from slip... Fault, but we are far from a comprehensive 3-D model of active fault.! Shapefiles from the links on the fault surface dip ) rupture ( thin lines ) and surface. Go on for decades ) Horizontal stresses from a comprehensive 3-D model of active fault structures fractures the. Was restricted to geodetic velocities, was recently presented by Meade et al and summed moment tensors a..., it is not found in the future to constrain fault slip in other regions within the first..
Tell Sackett Bowie Knife For Sale,
The Hitchin Post Marysville Ca Address,
Articles W