The mid-Cretaceous is well known for its ocean anoxic events. The causal mechanisms are controversial: stagnant deepwater, high biological productivity in the surface waters, and other possibilities have been suggested. Our study simulated the mid-Cretaceous ocean, using general circulation models combined with a marine biogeochemical cycle model to explore the relationship between thermohaline circulation and biogeochemical cycles and investigate the causes of ocean anoxic events. The simulated thermohaline circulation shows an unsteady inactive state. Oxygen concentrations in the deepwater decrease under the inactive state, but a horizontal gradient develops, with higher oxygen concentrations in the Tethys and lower concentrations in eastern Panthalassa. This is not due to the different ages of the deepwater but rather to the differences in biological productivity in the surface water, meaning that the relationship between thermohaline circulation and biogeochemical cycles under the inactive state is different from that in the present ocean. In the standard simulation, assuming the present level of the total amount of phosphate in the ocean, 29% of the bottom water is anoxic. The experiments increasing the amount of phosphate show its high sensitivity for extending the anoxic region with global-scale anoxia simulated under the doubled amount of phosphate. The high amount of phosphate would be reasonable because the inactive state would induce an imbalance of phosphate between riverine input and sediment output. Therefore, both the inactive thermohaline circulation and the increase in the total amount of phosphate in the ocean induce the global-scale anoxic condition in the deepwater. (c) 2008 Elsevier Ltd. All rights reserved.