To study the dynamic changes of soil organic carbon (SOC) content in coastal wetlands over the past 23 years, we use the third core area of the Dafeng Elk Nature Reserve in Jiangsu Province, China as the study area. A normalized difference soil spectral index (NDSOC), which is directly related to SOC, was constructed by using the spectral reflectance of Landsat images, and six auxiliary environmental variables related to SOC were extracted: surface temperature, surface moisture, soil salinity, the normalized difference vegetation index (NDVI), surface brightness, and soil texture information. The coastal wetland soil surface SOC inversion model was constructed by multiple linear regression, support vector machine, and particle swarm optimization-based random forest regression (PSO-RFR). By using the best inversion model, we obtain the dynamic distribution of soil surface SOC in the third core area of the Dafeng Elk Reserve for the past 23 years. The results indicate (1) that the NDSOC index constructed based on the short-wave infrared (SWIR) bands 1 and 2 of the Landsat images has a strong indication of SOC. Among the extracted auxiliary environmental variables, the NDVI, soil clay content, surface temperature, and soil salinity contribute more strongly to the SOC content of coastal wetlands, which averaged 27.9 %, 24.4 %, 22.1 %, and 18 %, respectively. (2) The PSO-RFR algorithm provided the best simulation of the surface SOC content of coastal wetlands: the model determination coefficient R2 is 0.731, the average fitting root mean square error RMSE is 0.917 g/kg, and the fitting uncertainty standard deviation was 0.362 g/kg. (3) From 2000 to 2022, the organic carbon storage and organic carbon density in the third core area of the Dafeng Elk Reserve gradually increased, the highest value was reached in 2016, they were 4.017 (104 t) and 1.718 kg/cm2, respectively.