[Project page]

Ground-motion records and site information for Puerto Rico and the U.S. Virgin Islands, 2006-2024

Contents

• Mixed Effects Data Files
• Ground-Motion Models
• Format of Mixed Effects Data Files
  • Total residuals (CSV)
  • Mixed effects (JSON)
  • Mixed effects summary (JSON)

Mixed Effects Data Files

Download All

Ground-Motion Models (GMMs)

The GMMs labels ending in PRVI indicate regionally adjusted models for Puerto Rico and the U.S. Virgin Islands developed as part of the analysis of the data provided in this dataset.

*-PRVI

Aagaard BT, Moschetti MP and Withers KB (2026) Suitability of ground-motion models for seismic hazard assessments in Puerto Rico and the U.S. Virgin Islands. Earthquake Spectra XX(XX): XX-XX. DOI:XXXX.

AbrahamsonEtAl2014

Abrahamson NA, Silva WJ and Kamai R (2014) Summary of the ASK14 ground motion relation for active crustal regions. Earthquake Spectra 30(3): 1025–1055. DOI:10.1193/070913EQS198M.

BooreEtAl2014

Boore DM, Stewart JP, Seyhan E and Atkinson GM (2014) NGA-West2 equations for predicting PGA, PGV, and 5% damped PSA for shallow crustal earthquakes. Earthquake Spectra 30(3): 1057–1085. DOI:10.1193/070113EQS184M.

CampbellBozorgnia2014

Campbell KW and Bozorgnia Y (2014) NGA-West2 ground motion model for the average horizontal components of PGA, PGV, and 5% damped linear acceleration response spectra. Earthquake Spectra 30(3): 1087–1115. DOI:10.1193/062913EQS175M.

ChiouYoungs2014

Chiou BSJ and Youngs RR (2014) Update of the Chiou and Youngs NGA model for the average horizontal component of peak ground motion and response spectra. Earthquake Spectra 30(3): 1117–1153. DOI:10.1193/072813EQS219M.

AbrahamsonGulerce2022

Abrahamson NA and Gulerce Z (2020) Regionalized ground-motion models for subduction earthquakes based on the NGA-SUB database. Technical Report 2020-25, Pacific Earthquake Engineering Research Center. peer.berkeley.edu/publications/2020-25.
Abrahamson NA and Gulerce Z (2022) Summary of the Abrahamson and Gulerce NGA-SUB ground-motion model for subduction earthquakes. Earthquake Spectra 38(4): 2638–2681. DOI:10.1177/87552930221114374.

KuehnEtAl2020

Kuehn NM, Bozorgnia Y, Campbell KW and Gregor N (2020) Partially non-ergodic ground-motion model for subduction regions using the NGA subduction database. Technical Report 2020-04, Pacific Earthquake Engineering Research Center. peer.berkeley.edu/publications/2020-04.
Kuehn NM, Bozorgnia Y, Campbell KW and Gregor N (2023) A regionalized partially nonergodic ground-motion model for subduction earthquakes using the NGA-Sub database. Earthquake Spectra 39(3): 1625–1657. DOI:10.1177/87552930231180906.

ParkerEtAl2020

Parker GA, Stewart JP, Boore DM, Atkinson GM and Hassani B (2020) NGA-subduction global ground motion models with regional adjustment factors. Technical Report 2020-03, Pacific Earthquake Engineering Research Center. peer.berkeley.edu/publications/2020-03.
Parker GA, Stewart JP, Boore DM, Atkinson GM and Hassani B (2022) NGA-subduction global ground motion models with regional adjustment factors. Earthquake Spectra 38(1): 456–493. DOI:10.1177/87552930211034889.

ClarosGomez2022

Claros Gomez DF (2022) Ground motion prediction equations for Puerto Rico. Doctor of Philosophy in Civil Engineering, University of Puerto Rico, Mayaguez.

KothaEtAl2020

Kotha SR, Weatherill G, Bindi D and Cotton F (2020) A regionally-adaptable ground-motion model for shallow crustal earthquakes in Europe. Bulletin of Earthquake Engineering 18(9): 4091–4125. DOI:10.1007/s10518-020-00869-1.

MotazedianAtkinson2005USGS

Motazedian D and Atkinson G (2005) Ground-motion relations for Puerto Rico. In: Mann P (ed.) Active Tectonics and Seismic Hazards of Puerto Rico, the Virgin Islands, and Offshore Areas, volume 385. Geological Society of America. ISBN 978-0-8137-2385-3. DOI:10.1130/0-8137-2385-X.61.

Formats of Mixed Effects Data Files

Total residuals (csv)

These files contain the natural logarithm of the total residuals with respect to the given ground-motion model (GMM), ln(R_total) = ln(Y_obs) - ln(Y_gmm), where R_total is the total residual, Y_obs is the ground-motion metric from the ground-motion record, and Y_gmm is the predicted ground-motion metric from the GMM. We compute total residuals for all spectral values even those at periods outside the passband of the filtered ground-motion records. We use negative spectral acceleration values to indicate values at spectral periods outside the passband of the filtered record. The list of columns is organized by earthquake, station, and record parameters followed by ground-motion metrics; it does not match the order of the columns in the CSV files.

EarthquakeId: Event ID from ANSS ComCat
EarthquakeTime: Earthquake origin time (UTC)
EarthquakeLatitude: Earthquake latitude (decimal degrees)
EarthquakeLongitude: Earthquake longitude (decimal degrees)
EarthquakeDepth: Earthquake depth (km)
EarthquakeMagnitude: Earthquake magnitude
EarthquakeMagnitudeType: Earthquake magnitude type
Network: Network code
DataProvider: Data source provider
StationCode: Station code
StationID: Concatenated network, station, and instrument codes
StationDescription: Station description
StationLatitude: Station latitude (decimal degrees)
StationLongitude: Station longitude (decimal degrees)
StationElevation: Station elevation (m)
SamplingRate: Record sampling rate (Hz)
BackAzimuth: Site-to-source azimuth (decimal degrees)
EpicentralDistance: Epicentral distance (km)
HypocentralDistance: Hypocentral distance (km)
SourceFile: Source file
RuptureDistance: Closest distance to the rupture plane (km)
RuptureDistanceVar: Variance of rupture distance estimate (km^2)
JoynerBooreDistance: Joyner-Boore distance (km)
JoynerBooreDistanceVar: Variance of Joyner-Boore distance estimate (km^2)
GC2_rx: Distance measured perpendicular to the fault strike from the surface projection of the up-dip edge of the fault plane (km)
GC2_ry: Distance measured parallel to the fault strike from the midpoint of the surface projection of the fault plane (km)
GC2_ry0: Horizontal distance off the end of the rupture measured parallel to strike (km)
GC2_U: Strike-normal (U) coordinate, defined by Spudich and Chiou (2015; https://doi.org/10.3133/ofr20151028) (km)
GC2_T: Strike-parallel (T) coordinate, defined by Spudich and Chiou (2015; https://doi.org/10.3133/ofr20151028) (km)
H1Lowpass: H1 channel lowpass frequency (Hz)
H1Highpass: H1 channel highpass frequency (Hz)
H2Lowpass: H2 channel lowpass frequency (Hz)
H2Highpass: H2 channel highpass frequency (Hz)
DURATIONp-q: p-q percent significant duration (s)
PGA: Total residual for the peak ground acceleration (%g)
PGV: Total residual for the peak ground velocity (cm/s); not included for the AbrahamsonGulerce2022 GMMs
SA(X): Total residual for the median horizontal spectral acceleration (%g) at X seconds
SiteId: Combination of network code and station code in the form NETWORK_CODE.STATION_CODE
FocalMechanism: Focal mechanism from earthquake classification
vs30: Vs30 (m/s) from station site information
vs30_measured: Flag indicating whether Vs30 is measured (AbrahamsonEtAl2014 and ChiouYoungs2014 GMMs only)
z1pt0: Depth (m) to the 1.0 km/s shear-wave velocity horizon (AbrahamsonEtAl2014 and ChiouYoungs2014 GMMs only)

Mixed Effects Results

These JSON files contain the results of the linear mixed effects analysis for each ground-motion model. The mixed effects include groupings by EarthquakeId and SiteId. PGV is in cm/s, PGA and spectral acceleration are in %g. NaN values indicate values that are not available from the analysis, such as cases in which a ground-motion metric is outside the passband for all records.

Values in uppercase denote variable values:

M: Number of earthquakes
N: Number of sites
R: Number of records
IMT: Ground-motion intensity metric, such as PGV, PGA, SA(X) where X is the spectral period in seconds.
EQ_ID_X: Earthquake id at index X in list
COEFF_X: Coefficient at index X in list
UNC_X: Variance of coefficient at index X in list
IMT_X: Ground-motion intensity metric at index X in list
IMT_VALUE: Value of ground-motion intensity metric

JSON schema:

{
    "fixed_effect": {
        "imts": {
            "IMT": IMT_VALUE,
        }
    },
    "random_effects": {
        "EarthquakeId": {
            "ids": [EQ_ID_0, EQ_ID_1, ..., EQ_ID_M],
            "imts": {
                "IMT": {
                    "coefficients": [COEFF_0, COEFF_1, ..., COEFF_M],
                    "uncertainty": [UNC_0, UNC_1, ..., UNC_M]
                },
            }
        },
        "SiteId": {
            "ids": [SITE_ID_0, SITE_ID_1, ..., SITE_ID_N],
            "imts": {
                "{IMT}": {
                    "coefficients": [COEFF_0, COEFF_1, ..., COEFF_N],
                    "uncertainty": [UNC_0, UNC_1, ..., UNC_N]
                },
            }
        }
    },
    "residuals": {
        "ids": {
            "EarthquakeId": [EQ_ID_0, EQ_ID_1, ..., EQ_ID_R],
            "SiteId": [SITE_ID_0, SITE_ID_1, ..., SITE_ID_R]
        },
        "imts": {
            "IMT": [IMT_0, IMT_1, ..., IMT_R],
        }
    }
}

Mixed Effects Summary (json)

These JSON files contain the fixed effect (bias, often denoted by k) and standard deviation in the event term (std_EarthquakeId, often denoted tau), site term (std_SiteId, often denoted phiS2S), and remaining residual (std_residuals, often denoted phiSS) for each ground-motion intensity metric (IMT: PGV, PGA, SA(X) where X is the spectral period in seconds). PGV is in cm/s, PGA and spectral acceleration are in %g.

Values in uppercase denote variable values:

IMT: Ground-motion intensity metric, such as PGV, PGA, SA(X) where X is the spectral period in seconds.
REAL_VALUE: Real value

JSON schema:

{
    "IMT": {
        "fixed_effect": REAL_VALUE,
        "std_residuals": REAL_VALUE,
        "std_EarthquakeId": REAL_VALUE,
        "std_SiteId": REAL_VALUE,
    }
}