In this study, the model Riccati equation with variable coefficients as functions, as well as a derivative of a fractional variable order (VO) of the Gerasimov-Caputo type, is used to approximate the data for some physical processes with saturation. In particular, the proposed model is applied to the description of solar activity (SA), namely the number of sunspots observed over the past 25 years. It is also used to describe data from Johns Hopkins University on coronavirus infection COVID-19, in particular data on the Russian Federation and the Republic of Uzbekistan. Finally, it is used to study issues related to seismic activity, in particular, the description of data on the volumetric activity of Radon (RVA). The Riccati equation used in the mathematical model was numerically solved by constructing an implicit finite difference scheme (IFDS) and its implementation by the modified Newton method (MNM). The calculated curves obtained in the study are compared with known experimental data. It is shown that if the model parameters are chosen appropriately, the model curves will give results that correlate well with real experimental data. Moreover, with other parameters of the model, it is possible to make some prediction about the possible course of the considered processes.