Sea spray aerosols (SSA) are important contributors to the global aerosol load, and they have a significant impact on the global climate system and chemical processes. To better assess the environmental impact of SSA emissions, it is important to quantitatively characterize the parametric dependence of SSA production fluxes on other environmental variables, and assess their role in global climate models (GCMs) and chemical transport models (CTMs). There are currently significant differences in the simulation of global SSA production fluxes in various models, and these differences are mainly attributed to the differences in the sea spray source functions (SSSFs). Research in the literature over the last decade has greatly broadened the knowledge of SSSFs and improved our understanding of SSA production fluxes and how they are applied in models. This paper clarifies the derivation process of different SSSFs, and provides a critical review of the application progress of different SSSFs in different models, which plays an important role in improving the accuracy of model simulations to evaluate the impact of SSA. More specifically, we review the studies related to the parameterization of sea spray aerosol production fluxes in the last decade, mainly focusing on the influencing factors of SSA production, the research progress of SSSFs, and the application of SSSFs in the models. The main environmental factors that alter the production flux in SSA include wind speed, sea surface temperature, wave state, salinity, surface active organic matter, etc. To better parameterize the SSA production flux, different expressions of SSSFs have been proposed. This paper mainly summarizes from different variables such as particle size distribution dependence, whitecap coverage, sea surface temperature dependence and wave state dependence. Moreover, many studies have incorporated different SSSFs in different GCMs and CTMs to evaluate the reasonableness of SSSFs and the climate impact of SSA emissions. Finally, we summarize the current results and inadequacy of SSA production flux parameterization, and put forward the prospect of future research direction.