A GEOCLIM simulation of climatic and biogeochemical consequences of Pangea breakup

[1] Large fluctuations in continental configuration occur throughout the Mesozoic. While it has long been recognized that paleogeography may potentially influence atmospheric CO2 via the continental silicate weathering feedback, no numerical simulations have been done, because of the lack of a spatially resolved climate-carbon model. GEOCLIM, a coupled numerical model of the climate and global biogeochemical cycles, is used to investigate the consequences of the Pangea breakup. The climate module of the GEOCLIM model is the FOAM atmospheric general circulation model, allowing the calculation of the consumption of atmospheric CO2 through continental silicate weathering with a spatial resolution of 7.5 degrees long x 4.5 degrees lat. Seven time slices have been simulated. We show that the breakup of the Pangea supercontinent triggers an increase in continental runoff, resulting in enhanced atmospheric CO2 consumption through silicate weathering. As a result, atmospheric CO2 falls from values above 3000 ppmv during the Triassic down to rather low levels during the Cretaceous ( around 400 ppmv), resulting in a decrease in global mean annual continental temperatures from about 20 degrees C to 10 degrees C. Silicate weathering feedback and paleogeography both act to force the Earth system toward a dry and hot world reaching its optimum over the last 260 Myr during the Middle-Late Triassic. In the super continent case, given the persistent aridity, the model generates high CO2 values to produce very warm continental temperatures. Conversely, in the fragmented case, the runoff becomes the most important contributor to the silicate weathering rate, hence producing a CO2 drawdown and a fall in continental temperatures. Finally, another unexpected outcome is the pronounced fluctuation in carbonate accumulation simulated by the model in response to the Pangea breakup. These fluctuations are driven by changes in continental carbonate weathering flux. Accounting for the fluctuations in area available for carbonate platforms, the simulated ratio of carbonate deposition between neritic and deep sea environments is in better agreement with available data.