“…It is easy to know that scalar Nicholson's blowflies Eq. (1.1) is a special case of Nicholson's blowflies system (1.2), where x i (t) denotes the density of the ith-population at time t, a ij (t) (i = j) is the rate of the population moving from class j to class i at time t, a ii (t) is the coefficient of instantaneous loss (which integrates both the death rate and the dispersal rates of the population in class i moving to the other classes), β ij (t)x i (tτ ij (t))e -γ ij (t)x i (t-τ ij (t)) is the birth function, β ij (t) is the birth rate for the species, 1 γ ij (t) is the ith-population reproducing at its maximum rate, and τ ij (t) is the generation time of the ith-population at time t. For the feedback function xe -x and its derivative 1-x e x , the author in [8] It is worth noting that the global exponential stability of almost periodic solutions of (1.1) has been shown in [5,6] under the restriction that the almost periodic solution exists in a small interval [κ, κ] ≈ [0.7215355, 1.342276], and the global exponential stability of (1.2) has been established in [7] where the authors adopted the restraint that the almost periodic solution exists in a small domain Obviously, the above restriction and restraint do not accord with the biological significance of the population models. On the other hand, γ ij (t) ≥ 1 for all t ∈ R, i ∈ Q, j ∈ I := {1, 2, .…”