“…Then we performed a sensitivity analysis to verify the generality of our results under a wider range of ecological constraints and initial conditions ( u b = [0, 0.5, 1, 1.5, 2]; u # = [1, 1.5, 2, 3, 4, 5]; f b = [3.2, 2.5, 1.9, 1.5, 1.1]; h b = [0, 5, 7, 10, 12, 13, 14]; h # = [10, 12, 14, 17, 20, 30]; N(t 0 ) = [25, 30, 35] ). Constraint values were chosen to explore the range of possible states and controls observed in our study area on lapwing nesting fields (0≤h≤30 and 0≤u≤5; [18]). As f(u) is a monotonous function of u , values of f b were thus chosen to correspond to the different thresholds on u # .…”