This project proposes an algorithm that utilizes data from the Next Generation Liquefaction Database to create geospatial meshes of ground motion intensity, and to determine intensities at sites where earthquake-induced liquefaction has been observed. This project applies a combined modeling and data-driven approach to determine residuals or model deviance from true recorded intensity measure values at recording stations. Going forward, this project serves to provide the geotechnical engineering community with a new technique for more consistent ground motion estimations at liquefaction sites.