In taxi management, taxi-driver shift behaviors play a key role in arranging the operation of taxis, which affect the balance between the demand and supply of taxis and the parking spaces. At the same time, these behaviors influence the daily travel of citizens. An analysis of the distribution of taxi-driver shifts, therefore, contributes to transportation management. Compared to the previous research using the real shift records, this study focuses on the spatiotemporal analysis of taxi-driver shifts using big trace data. A two-step strategy is proposed to automatically identify taxi-driver shifts from big trace data without the information of drivers’ identities. The first step is to pick out the frequent spatiotemporal sequential patterns from all parking events based on the spatiotemporal sequence analysis. The second step is to construct a Gaussian mixture model based on prior knowledge for further identifying taxi-driver shifts from all frequent spatiotemporal sequential patterns. The spatiotemporal distribution of taxi-driver shifts is analyzed based on two indicators, namely regional taxi coverage intensity and taxi density. Taking the city of Wuhan as an example, the experimental results show that the identification precision and recall rate of taxi-driver shift events based on the proposed method can achieve about 95% and 90%, respectively, by using big taxi trace data. The occurrence time of taxi-driver shifts in Wuhan mainly has two high peak periods: 1:00 a.m. to 4:00 a.m. and 4:00 p.m. to 5:00 p.m. Although taxi-driver shift behaviors are prohibited during the evening peak hour based on the regulation issued by Wuhan traffic administration, experimental results show that there are still some drivers in violation of this regulation. By analyzing the spatial distribution of taxi-driver shifts, we find that most taxi-driver shifts distribute in central urban areas such as Wuchang and Jianghan district.