This study aimed to analyze the incidence of work-related musculoskeletal disorders (WMSD) in Brazil from 2007 to 2019, examining the spatial, temporal, and spatiotemporal patterns of their occurrence. An ecological time series study was conducted using spatial analysis techniques. WMSD morbidity data from 2007 to 2019 were collected from the Brazilian Information System for Notificable Diseases of the Brazilian Health Informatics Department. Incidence rates were standardized and smoothed using the local empirical Bayes’ theorem. Time trends were analyzed by segmented linear regression. Spatial analysis was performed using Moran’s univariate global (I) and local (LISA) indexes. The spatiotemporal scan statistic was used to identify high-risk spatiotemporal clusters for WMSD. A total of 93,387 cases of WMSD were recorded in Brazil. Temporal trends showed an increase in all regions except the Northeast, which remained stable. The incidence of WMSD showed a spatial dependence, with spatial and space-time clusters identified, especially in the Southeast region, overlapping the largest economic-industrial center of the country. The spatiotemporal clustering observed in one region suggests the highest level of industrial and economic development. Our findings highlight the need to implement intersectoral surveillance policies, inspect working conditions, and invest in the prevention and promotion of workers’ health.