Modelling the Solar Cycle Nonlinearities into the Algebraic Approach

Understanding and predicting solar-cycle variability requires accounting for nonlinear feedbacks that regulate the buildup of the Sun's polar magnetic field. We present a simplified but physically grounded algebraic approach that models the dipole contribution of active regions (ARs) while incorporating two key nonlinearities: tilt quenching (TQ) and latitude quenching (LQ). Using ensembles of synthetic cycles across the dynamo effectivity range $\lambda_R$, we quantify how these mechanisms suppress the axial dipole and impose self-limiting feedback. Our results show that (i) both TQ and LQ reduce the polar field, and together they generate a clear saturation (ceiling) of dipole growth with increasing cycle amplitude; (ii) the balance between LQ and TQ, expressed as $R(\lambda_R) = \mathrm{dev(LQ)}/\mathrm{dev(TQ)}$, transitions near $\lambda_R \approx 12^\circ$, with LQ dominating at low $\lambda_R$ and TQ at high $\lambda_R$; (iii) over $8^\circ \leq \lambda_R \leq 20^\circ$, the ratio follows a shallow offset power law with exponent $n \approx 0.36 \pm 0.04$, significantly flatter than the $n=2$ scaling assumed in many surface flux--transport (SFT) models; and (iv) symmetric, tilt-asymmetric, and morphology-asymmetric AR prescriptions yield nearly identical $R(\lambda_R)$ curves, indicating weak sensitivity to AR geometry for fixed transport. These findings demonstrate that nonlinear saturation of the solar cycle can be captured efficiently with algebraic formulations, providing a transparent complement to full SFT simulations. The method highlights that the LQ\--TQ balance is primarily controlled by transport ($\lambda_R$), not by active-region configuration, and statistically disfavors the SFT-based $1/\lambda_R^{2}$ dependence.