An earthquake of M s7.0 struck Jiuzhaigou Country in Sichuan Province on 8 August 2017. It occurred at the northeastern boundary of the Bayan Har block. The Jiuzhaigou earthquake was one of the largest earthquakes in continental China since the 2013 M s7.0 Lushan earthquake, and it was widely felt across Sichuan and adjacent provinces. The earthquake left 25 dead and up to 500 injured. Because the Jiuzhaigou earthquake did not produce obvious surface ruptures, the causative fault of this event was much debated. The shape of the seismogenic fault, the focal depth of the mainshock and the characteristics of the aftershock sequence were also not very clear. The Sichuan earthquake administration deployed 6 temporary seismic stations around the source region within 4 d of the mainshock. In this study, we collected seismic waveform data and phase arrival data recorded by permanent and temporary stations. We present a systematic investigation about the location of the Jiuzhaigou earthquake sequence with high precision. These results provide new constraints for refining earthquake rupture process, constructing seismogenic models, and evaluating regional earthquake risks. The mainshock of the Jiuzhaigou M s7.0 earthquake was relocated using a 3D velocity model. The procedure takes station elevation, topography and the Earth’s ellipticity into consideration. The location of the mainshock was found to be 103.806°E, 33.201°N, and the focal depth to be 20.4 km. The depth of the initial rupture point of the mainshock was estimated to be below 14 km with S-P arrival time difference recorded at the nearest strong motion station. The locations of early aftershocks were corrected with latter aftershocks recorded using dense temporary seismic stations. A double difference algorithm was used to relocate 1-month aftershock sequence of the Jiuzhaigou earthquake. The locations of 3030 aftershocks were here determined. The average location error was found to be 0.16, 0.15, and 0.18 km in the E-W, N-S, and U-D directions, respectively. The relocation results showed that the aftershocks spread approximately 42 km, trending NWW. The aftershock zone was found to connect the Tazang fault to the north and to the Huya fault to the south. The mainshock is located at the center of the aftershock zone, with similar aftershock length at the both sides. There is a ~5 km sparse aftershock segment to the NW of the mainshock which is consistent with large coseismic slip area. This phenomenon can be explained due to the large amount of stress released during the mainshock rupture. The depth of the aftershock in the northwest end was relatively shallow, and the aftershock zone was about 6 km wide. The depth of aftershock in the southeast section was deeper and the area was narrower, about 4 km. The predominant distribution of the focal depth is 4–20 km. The dip angle of the seismogenic fault is steep, with an average value of about 84°, and changes obviously along the strike direction. The seismogenic fault inclined to the southwest in the shallow part and to the northeast in the deep part. The complex fault shape should be considered when constructing seismogenic model and inverting earthquake rupture process. The depth of the initial rupture point of the mainshock is deeper than the centroid depth and average depth of the aftershock sequence, indicating that the earthquake rupture spreads from deep space to shallow. We find that the aftershocks migrated in the along-strike direction with logarithmic time since the mainshock. The length of the aftershock area expanded from 25 to 42 km in 2 d. The rate of spread was approximately 1.3 km/log(s), consistent with aftershock expansion caused by propagating afterslip. There may be afterslip in the source area of the Jiuzhaigou earthquake.