Seawater pH is a direct proxy of ocean acidification, and monitoring the global pH distribution and long-term series changes is critical to understanding the changes and responses of the marine ecology and environment under climate change. Owing to the lack of sufficient global-scale pH data and the complex relationship between seawater pH and related environmental variables, generating time-series products of satellite-derived global sea surface pH poses a great challenge. In this study, we solved the problem of the lack of sufficient data for pH algorithm development by using the massive underway sea surface carbon dioxide partial pressure (pCO(2)) dataset to structure a large data volume of near in situ pH based on carbonate calculation between underway pCO(2) and calculated total alkalinity from sea surface salinity and relevant parameters. The remote sensing inversion model of pH was then constructed through this massive pH training dataset and machine learning methods. After several tests of machine learning methods and groups of input parameters, we chose the random forest model with longitude, latitude, sea surface temperature (SST), chlorophyll a (Chla), and Mixed layer depth (MLD) as model inputs with the best performance of correlation coefficient (R-2 = 0.96) and root mean squared error (RMSE = 0.008) in the training set and R-2 = 0.83 (RMSE = 0.017) in the testing set. The sensitivity analysis of the error variation induced by the uncertainty of SST and Chla (SST <= +/- 0.5 degrees C and Chla <= +/- 20%; RMSESST <= 0.011 and RMSEChla <= 0.009) indicated that our sea surface pH model had good robustness. Monthly average global sea surface pH products from 2004 to 2019 with a spatial resolution of 0.25 degrees x 0.25 degrees were produced based on the satellite-derived SST and Chla products and modeled MLD dataset. The pH model and products were validated using another independent station-measured pH dataset from the Global Ocean Data Analysis Project (GLODAP), showing good performance. With the time-series pH products, refined interannual variability and seasonal variability were presented, and trends of pH decline were found globally. Our study provides a new method of directly using remote sensing to invert pH instead of indirect calculation based on the construction of massive underway calculated pH data, which would be made useful by comparing it with satellite-derived pCO(2) products to understand the carbonate system change and the ocean ecological environments responding to the global change.