Coupled climate models integrate atmospheric, oceanic, and land submodels, while the uncertainty of model parameters from different parameterization schemes or empirically derived parameters inevitably introduces systematic biases. Coupled parameter optimization (CPO) can reduce these biases to improve weather forecast and climate prediction, but must address strong nonlinearities inherent in coupled models. The analytical four-dimensional ensemble variational (A-4DEnVar) data assimilation method retains the nonlinear processing capability of the four-dimensional variational (4D-Var) data assimilation method but gets rid of the dependence on the adjoint model. In this study, a novel dynamic independent point (DIP) scheme is introduced to the improved A-4DEnVar, which reduces computational dimensionality and further explores a broader parameter space of dimensionality reduction through the outer loop. Based on the improved A-4DEnVar, a series of geographic-dependent CPO experiments with an idealized 2D coupled model are carried out. Results show that A-4DEnVar accurately captures the geographical characteristics of parameters and effectively optimizes cross-component parameters despite strong nonlinearity. Additionally, the DIP scheme presents significant advantages compared to the static independent point scheme, especially with fewer independent points. This work is offering a new perspective for parameter optimization in coupled general circulation models used for climate estimation and prediction.