The Next-Generation Attenuation for subduction zone regions project (NGA-Sub) has developed data resources and ground motion models for global subduction zone regions. Here we describe the NGA-Sub database. To optimize the efficiency of data storage, access, and updating, data resources for the NGA-Sub project are organized into a relational database consisting of 20 tables containing data, metadata, and computed quantities (e.g. intensity measures, distances). A database schema relates fields in tables to each other through a series of primary and foreign keys. Model developers and other users mostly interact with the data through a flatfile generated as a time-stamped output of the database. We describe the structure of the relational database, the ground motions compiled for the project, and the means by which the data can be accessed. The database contains 71,340 three-component records from 1880 earthquakes from seven global subduction zone regions: Alaska, Central America and Mexico, Cascadia, Japan, New Zealand, South America, and Taiwan. These data were processed on a component-specific basis to minimize noise effects in the data and remove baseline drifts. Provided ground motion intensity measures include peak acceleration, peak velocity, and 5%-damped pseudo-spectral accelerations for a range of oscillator periods.