A nutrient (N), phytoplankton (P), zooplankton (Z), and detritus (D) ecosystem model coupled to an ice‐ocean model was applied to the Bering and Chukchi Seas for 2007–2008. The model reasonably reproduces the seasonal cycles of sea ice, phytoplankton, and zooplankton in the Bering–Chukchi Seas. The spatial variation of the phytoplankton bloom was predominantly controlled by the retreat of sea ice and the increased gradient of the water temperature from the south to the north. The model captures the basic structure of the measured nutrients and chl‐a along the Bering shelf during 4–23 July 2008, and along the Chukchi shelf during 5–12 August 2007. In summer 2008, the Green Belt bloom was not observed by either the satellite measurements or the model. The model‐data comparison and analysis reveal the complexity of the lower trophic dynamics in the Bering and Chukchi Seas. The complexity is due to the nature that the physical and biological components interact at different manners in time and space, even in response to a same climate forcing, over the physically distinct geographic settings such as in the Bering and North Aleutian Slopes, deep Bering basins, Bering shelf, and Chukchi Sea. Sensitivity studies were conducted to reveal the underlying mechanisms (i.e., the bottom‐up effects) of the Bering–Chukchi ecosystem in response to changes in light intensity, nutrient input from open boundaries, and air temperature. It was found that (1) a 10% increase in solar radiation or light intensity for the entire year has a small impact on the intensity and timing of the bloom in the physical–biological system since the light is not a limiting factor in the study region; (2) a 20% increase in nutrients from all the open boundaries results in an overall 7% increase in phytoplankton, with the Slope region being the largest, and the Bering shelf and Chukchi being the smallest; and (3) an increase in air temperature by 2°C over the entire calculation period can result in an overall increase in phytoplankton by 11%.
