Paleoceanographic Perturbations and the Marine Carbonate System during the Middle to Late Miocene Carbonate Crash—A Critical Review

: This study intends to review and assess the middle to late Miocene Carbonate Crash (CC) events in the low to mid latitudes of the Paciﬁc, Indian, Caribbean and Atlantic Oceans as part of the global paleoceanographic reorganisations between 12 and 9 Ma with an emphasis on record preservation and their relation to mass accumulation rates (MAR). In the Eastern Paciﬁc the accumulation changes in carbonate and opal probably reﬂect an El-Niño-like state of low productivity, which marks the beginning of the CC-event (11.5 Ma), followed by decreased preservation and inﬂux of corrosive bottom waters (10.3 to 10.1 Ma). At the same time in the Atlantic, carbonate preservation considerably increases, suggesting basin-to-basin fractionation. The low-latitude Indian Ocean, the Paciﬁc and the Caribbean are all characterised by a similar timing of preservation increase starting at ~9.6–9.4 Ma, while their MARs show drastic changes with different timing of events. The Atlantic preservation pattern shows an increase as early as 11.5 Ma and becomes even better after 10.1 Ma. The shallow Indian Ocean (Mascarene plateau) is characterised by low carbonate accumulation throughout and increasing preservation after 9.4 Ma. At the same time, the preservation in the Atlantic, including the Caribbean, is increasing due to enhanced North Atlantic deep-water formation, leading to the increase in carbonate accumulation at 10 Ma. Moreover, the shoaling of the Central American Isthmus might have helped to enhance Caribbean preservation after 9.4 Ma. Lower nannoplankton productivity in the Atlantic should have additionally contributed to low mass accumulation rates during the late CC-interval. Overall, it can be inferred that these carbonate minima events during the Miocene may be the result of decreased surface ocean productivity and oceanographically driven increased seaﬂoor dissolution. suggests that the carbon isotope pattern in the Atlantic is not related to


Introduction
The Miocene Carbonate Crash (CC) events include several sharp reductions in carbonate concentration and/or accumulation occurring in various marine systems around the world, such as the equatorial Pacific and Indian Ocean [1][2][3][4]; the Caribbean Sea [5]; the equatorial [6,7]; and South Atlantic [8][9][10]; or the Southern Ocean [11]) during the middle to late Miocene . Although global in character, the onset of this widespread phenomenon is diachronous, with the initial phases of carbonate reduction occurring at 13.2 Ma at the equatorial Indian Ocean [3] and~11.5 Ma at the Caribbean [5,12] and the eastern Pacific [13]. Moreover, beyond timing, the duration and magnitude in terms of the number of the different decline events are not consistent, and possible causes are discussed controversially. Postulations of potential mechanisms include changes in deep water circulation, shoaling of the Calcite Compensation Depth (CCD), shifts in shallow to deep carbonate fractionation, terrigenous dilution, productivity, preservation changes, evolutionary trends, and adaptations of the main calcareous plankton groups [1,3,4,[14][15][16]. During the middle to late Miocene major reorganisations on land and in the oceans induced prominent shifts in the global climate system, including important changes of plate tectonic settings (e.g., the closure of the Indonesian Seaway (IS) [17] and uplift of the Central American Seaway (CAS) [18][19][20]. The subsidence of the Greenland-Scotland Ridge is associated with secular variations in its bathymetry induced by variations in temperature and buoyancy of the Iceland mantle plume [21,22], opening and widening of the Frame Strait [23], as well as the uplift of high mountain ranges, in particular the Himalaya and Tibet Plateau as well as the Andes [24]. This, in turn, strongly influenced continental weathering [25] inducing the formation of new huge fluvial drainage systems like the Amazon [26] and therefore the increased major nutrient (e.g., phosphorous and silica) delivery to the oceans (e.g., [27][28][29]). During the Miocene strong variations in high-latitude climates are also observed. A gradual warming during the early to middle Miocene (24 Ma to 15 Ma) is followed by rapid cooling until the end of the Miocene [30][31][32][33][34][35][36][37]. Contemporaneously, the ice shield on Antarctica expanded considerably [38,39] and the first small-dimensioned ice sheets were initiated in the Northern Hemisphere [40][41][42][43][44][45][46], leading to the equivalent 40-60 m sea level drop [47]. Holbourn, et al. [48] and Tian, et al. [49] have also suggested that there was a linkage between the brief high production at~14 Ma to Atlantic cooling.
The asymmetric cooling resulted in more vigorous trade winds in the southern hemisphere and displacement of the Inter Tropical Convergence Zone (ITCZ) towards the north [50,51]. The global climate cooled stepwise as recorded in both benthic [30,35,49,52] and planktonic [49,[53][54][55] foraminiferal δ 18 O values. The Miocene isotope events (Mievents) that denote increases in ice volume and/or surface and deep water cooling [52,[56][57][58], which however are not always traceable among sites [59]. This was connected with major reorganisations of surface and deep water circulation, in particular, the initiation of modern low latitude upwelling systems in the Atlantic [8,10,60] and Pacific Oceans [61]. Wright, et al. [62] and Wright and Miller [21] place the initiation of North Component Water (NCW), a precursor to the modern North Atlantic Deep Water (NADW) and a primary component of deep-water convection, in the late early Miocene. More recently, Poore, et al. [22] stated that based on composite Atlantic, Pacific and Southern Ocean δ 13 C benthic foraminifera records there is no evidence of significant NCW overflow over the Greenland-Scotland Ridge before 12 Ma, which is consistent with evidence from Cd/Ca ratios [63].
The change in circulation mode further induced a shift in nutrient redistribution, which is manifested in the middle Miocene "silica switch" that marks the intensification of opaline deposition in the global oceans, associated with the evolutionary diversification of diatoms between 15 and 10 Ma [64][65][66][67][68]. According to Lyle and Baldauf [13], the nutrients were more accessible during that time because of a combination of oceanographic and biochemical factors, including higher organic matter degradation in the upper water column and deep-intermediate nutrient pathways that efficiently recycled these components to the upwelling zones. This opal burst was also associated with the worldwide enhancement of carbonate, phosphate and barium accumulation rates, the vertical extension of the oxygen minimum zones [8,69] and the body-size increase of many groups of large predatory marine vertebrates [70], such as fishes [71,72], marine mammals [73] and seabirds [74]. However, with the exception of species-specific (e.g., Orbulina universa [75]) and localized adaptations (e.g., [37,76,77]), the response of planktonic organisms in such oceanic environments is still poorly understood. Schmidt, et al. [78] found a shift in size of mid to low latitude foraminiferal assemblages through the Cenozoic with extreme development towards larger sizes in the late Miocene. Moreover, coccolithophores are suspected to react to oceanographic revolutions, showing extraordinary size trends during the late Miocene (e.g., Calcidicus leptoporus [79]; Reticulofenestra [80]).
Variability in the size of marine plankton is directly relevant to calcium carbonate (CaCO 3 ) and carbon dioxide (CO 2 ) production. Carbonate accumulation patterns are crucial to understand circulation patterns and climate change, since they are reflecting a balance of productivity and preservation on the sea floor. Global carbonate budgets should display a steady state of riverine influx and burial rate of calcium carbonate [81]. Because of the similarity of eustatic sea level changes and fluctuations of the CCD over the path of the Cenozoic [82], changes in the CCD were associated with shelf-to-basin fractionation. The underlying assumption is that during high sea level, carbonate accumulation would preferably take place on shelves rather than the open ocean [83,84]. In order to trace CCD changes classical approaches, certain carbonate percentages or Carbonate Mass Accumulation Rates (MAR Carbonate ) were used as cut-off values to define CCD boundaries. For instance, Roth, et al. [5] used the 0.1g/cm 2 /ky carbonate flux to reconstruct CCD changes in the deep Pacific throughout the middle to late Miocene CC-events, while Berger and Wefer [85] or Van Andel, et al. [86] used the 20 wt.%-CaCO 3 as the level of the CCD (see Figure 1A). A closer look at preservation proxies might reveal discrepancies between certain value CCD reconstructions and their relation to carbonate productivity levels or dissolution intensities. Therefore, in this study, we want to review the carbonate crash events from ODP cores of low to mid latitudes with emphasis on available carbonate preservation proxy data. The review is supplemented by original data contribution of material from ODP Site 1237 in the SE-Pacific in front of Peru (Supplementary Table S1).
Geosciences 2021, 11, 94 3 of 25 of mid to low latitude foraminiferal assemblages through the Cenozoic with extreme development towards larger sizes in the late Miocene. Moreover, coccolithophores are suspected to react to oceanographic revolutions, showing extraordinary size trends during the late Miocene (e.g., Calcidicus leptoporus [79]; Reticulofenestra [80]). Variability in the size of marine plankton is directly relevant to calcium carbonate (CaCO3) and carbon dioxide (CO2) production. Carbonate accumulation patterns are crucial to understand circulation patterns and climate change, since they are reflecting a balance of productivity and preservation on the sea floor. Global carbonate budgets should display a steady state of riverine influx and burial rate of calcium carbonate [81]. Because of the similarity of eustatic sea level changes and fluctuations of the CCD over the path of the Cenozoic [82], changes in the CCD were associated with shelf-to-basin fractionation. The underlying assumption is that during high sea level, carbonate accumulation would preferably take place on shelves rather than the open ocean [83,84]. In order to trace CCD changes classical approaches, certain carbonate percentages or Carbonate Mass Accumulation Rates (MARCarbonate) were used as cut-off values to define CCD boundaries. For instance, Roth, et al. [5] used the 0.1g/cm 2 /ky carbonate flux to reconstruct CCD changes in the deep Pacific throughout the middle to late Miocene CC-events, while Berger and Wefer [85] or Van Andel, et al. [86] used the 20 wt.%-CaCO3 as the level of the CCD (see Figure 1A). A closer look at preservation proxies might reveal discrepancies between certain value CCD reconstructions and their relation to carbonate productivity levels or dissolution intensities. Therefore, in this study, we want to review the carbonate crash events from ODP cores of low to mid latitudes with emphasis on available carbonate preservation proxy data. The review is supplemented by original data contribution of material from ODP Site 1237 in the SE-Pacific in front of Peru (Supplementary Table S1).  [83], Van Andel, et al. [86], and Hsü and Wright [87]. (B) Compiled CCD reconstructions for the mains oceans after Berger and Wefer [85].

Strategies to Detect Carbonate Crash Mechanisms: Discrimination between Dilution, Productivity, and Dissolution
Due to the technical advance of coring deep ocean floor sediments with hydraulic piston corers (e.g., during legs of the Ocean drilling program-ODP), there is the possibility to investigate a huge archive of well dated continuous sedimentary successions spanning the middle to late Miocene. A first look at a pelagic carbonate rich sedimentary sequence can give a hint to carbonate variations expressed as cyclic gradational variations Hsü & Wright (1985) van Andel (1977) North Atlantic Pacific High ------Sea Level------Low South Atlantic Miocene Oligo.

Strategies to Detect Carbonate Crash Mechanisms: Discrimination between Dilution, Productivity, and Dissolution
Due to the technical advance of coring deep ocean floor sediments with hydraulic piston corers (e.g., during legs of the Ocean drilling program-ODP), there is the possibility to investigate a huge archive of well dated continuous sedimentary successions spanning the middle to late Miocene. A first look at a pelagic carbonate rich sedimentary sequence can give a hint to carbonate variations expressed as cyclic gradational variations in sediment colour, where darker colour bands represent carbonate decreases. This method is however inappropriate for sediments of very high carbonate contents, where even severe dissolution or productivity decreases are not reflected with apparent decreases in the carbonate contents (hence no changes in their colour bands). A strong dilution signal (e.g., through a change in the ratio of opal to carbonate sedimentation or terrigenous input) can also produce low carbonate contents. Besides the observation of carbonate contents, which can be lowered due to dilution, dissolution, decreased productivity or changing ratio of opal/calcite in the rain rate, mass accumulation rates (MAR) give a measure of the component burial per surface and time and are calculated as follows: where in ODP records the Gamma ray attenuation porosity evaluator (GRAPE) data and discrete density measurements are often used to calculate dry bulk densities. Routinely MARs are reported every 1 My increment accounting for low resolution age models.
Additionally, values from composite sections have to be corrected by a growth factor to account for core expansions and disturbances [88].
MARs can help to identify whether changes in downcore carbonate contents are due to dilution or productivity/dissolution changes. The quality of these data depends on the quality of the age model and the ablated sedimentation rates, with the highest resolution gained by orbitally tuned age models. Such models interpret the sediment composition as sedimentary response to orbital cycles and is based on high resolution data of density, colour and magnetic susceptibility or geochemical measurements in the framework of classical magneto-and biostratigraphy. However, sedimentation rates might be handled with care if there are special burial circumstances involved such as the rapid sedimentation of diatom mats, which are laminated sediments that can form thick and spatial widespread deposits [89][90][91]. These sediments are rather interpreted as results of complex surface water processes than deep water dissolution phases [92][93][94]. Additionally, re-depositional events such as turbidites and slumps or phases of no deposition (hiatus) disturb the premise of slow and steady pelagic sedimentation, supposed to form an even veneer on the ocean floor topography. The process of sediment focusing and lateral advection is also under discussion. Indicative attempts to correct for this are the correction of MAR by 230 Th normalisation [95,96] or assuming a constant rain rate from extra-terrestrial 3 He [97,98].
However, MARs are the key data to observe CC-events. If decreasing accumulation rates of carbonate are accompanied by steady and low terrigenous accumulation, dilution can be excluded as a main factor. What remains then is the question of whether MAR Carbonate minima are due to lowered productivity or dissolution. Depending on the paleo-water depth of the observation, one will make a first hypothesis based on the fact that at greater depth, calcite will be dissolved more easily, and at shallower depth, productivity will be the most important factor. Lowered carbonate supply tends to increase the possibility of dissolution and will be balanced by a rising CCD. However, especially under high productivity conditions, dissolution above the lysocline contributes significantly to overall dissolution [99]. Herein, the rain ratio of organic carbon to CaCO 3 is an indicator for high productivity. Elevated proportions of organic carbon are often associated with high rates of opal accumulations [100]. The accumulation of opal in the form of diatom frustules might favour the preservation of carbonates. No evidence of enhanced dissolution of foraminifera has been observed within diatom mats [101], which might be attributed to the fast-settling velocities and/or the strong meshwork that suppresses benthic activity. The release of silica in pore waters might as well act against carbonate dissolution [102]. A higher proportion of benthic organisms can also point to higher productivity of surface waters, but might be also evidence of dissolution (e.g., [69,103]). In many cases, it is reasonable to take into account multiple proxies (e.g., sand content and/or the planktonic foraminiferal fragmentation index [104]) investigated on splits of the sand fraction. Foraminiferal fragmentation has shown to be a robust indicator of carbonate dissolution in deep sea sediments [105][106][107]. The examined Miocene tests of foraminifera appear opaque and recrystallized due to early diagenetic alterations; however, we expect that the diagenesis will not alter the whole test fragmentation index. The validity of these methods is restricted to dissolution above the CCD. Nannofossil based dissolution indices might help to find clues for dissolution intensities, despite the high potential of ecological control on abundance patterns. The application of these preservation proxies might be even better under strong dissolution conditions and for application to closely spaced inter-site comparison such as a depth transect [108]. The absence of foraminifera or a strongly reduced taphonomic assemblage in dissolved sediments often inhibits the investigation and/or interpretation of stable isotopes on benthic and planktonic faunas, which are standard proxies for circulation, productivity and climate, however here also the less well understood nannofossils might fill the gap in reflection of surface water processes in future research.
The investigations of MAR Carbonate at depth transects is especially meaningful in CCD reconstructions because site to site differences reveal the loss of carbonate from one site to another and other processes influencing carbonate MAR besides dissolution. In case of the Walvis Ridge depth transect winnowing, reworking and down-slope re-suspension seem to alter wt.% of carbonate even towards producing reversed gradients in carbonate contents [109], with higher carbonate contents at the deeper sites in comparison to the shallower ones. Another such example for a depressed carbonate depth gradient is at the Ontong Java plateau [110].
For the following comparison we have chosen several cores from low to mid latitudes of the world's main ocean representing key locations of the CC-events (see Table 1 and Figure 2). The MAR were recalculated following the method described above. Carbonate and coarse fraction data from ODP Site 1237 in the SE-Pacific represent the only new data on the middle to late Miocene. Component analysis on these samples was applied to representative splits of the > 63µm fraction with at least 300 particles counted. The fragmentation index here is the ratio of whole tests (shells comprising at least half of a whole planktonic foraminifera test) divided by the sum of whole tests and fragments. The framework of comparison of different datasets and calculations of MAR Carbonate is given by the age model. in deep sea sediments [105][106][107]. The examined Miocene tests of foraminifera appear opaque and recrystallized due to early diagenetic alterations; however, we expect that the diagenesis will not alter the whole test fragmentation index. The validity of these methods is restricted to dissolution above the CCD. Nannofossil based dissolution indices might help to find clues for dissolution intensities, despite the high potential of ecological control on abundance patterns. The application of these preservation proxies might be even better under strong dissolution conditions and for application to closely spaced inter-site comparison such as a depth transect [108]. The absence of foraminifera or a strongly reduced taphonomic assemblage in dissolved sediments often inhibits the investigation and/or interpretation of stable isotopes on benthic and planktonic faunas, which are standard proxies for circulation, productivity and climate, however here also the less well understood nannofossils might fill the gap in reflection of surface water processes in future research. The investigations of MARCarbonate at depth transects is especially meaningful in CCD reconstructions because site to site differences reveal the loss of carbonate from one site to another and other processes influencing carbonate MAR besides dissolution. In case of the Walvis Ridge depth transect winnowing, reworking and down-slope re-suspension seem to alter wt.% of carbonate even towards producing reversed gradients in carbonate contents [109], with higher carbonate contents at the deeper sites in comparison to the shallower ones. Another such example for a depressed carbonate depth gradient is at the Ontong Java plateau [110].
For the following comparison we have chosen several cores from low to mid latitudes of the world's main ocean representing key locations of the CC-events (see Table 1 and Figure 2). The MAR were recalculated following the method described above. Carbonate and coarse fraction data from ODP Site 1237 in the SE-Pacific represent the only new data on the middle to late Miocene. Component analysis on these samples was applied to representative splits of the > 63µm fraction with at least 300 particles counted. The fragmentation index here is the ratio of whole tests (shells comprising at least half of a whole planktonic foraminifera test) divided by the sum of whole tests and fragments. The framework of comparison of different datasets and calculations of MARCarbonate is given by the age model.

Age Control
The comparison of the cores from the eastern Pacific (Leg 138, Leg 202), Atlantic (Leg 154, Leg 1085) and the Caribbean (Leg 165) is advantageous because they are based on the same biostratigraphy using bioevents dated by Flores, et al. [111] from Leg 138. Age control of Leg 138 records is provided by biostratigraphical studies over the whole record, astronomical tuning of GRAPE density down to 10 Ma and magnetostratigraphical correlation in the interval younger than 13.25 Ma [112]. The astronomical tuning of Ceara Rise Site 926 based on correlation of magnetic susceptibility maxima to Northern Hemisphere insolation minima. The orbital solution of records was also compared to high resolution biostratigraphy from 5 to 14 Ma [113]. Age models from Indian Ocean Leg 115 were updated to the common standard of Berggren, et al. [114] timescale. An updated age model was also used for Site 806 from the Ontong Java Plateau [115]. Age models of the eastern Pacific Site 1010 along Baja California [116] and Site 1241 in the Eastern Equatorial Pacific (EEP) did not provide magenetostratigraphic control, however displayed the low latitude standard zonation of calcareous nannofossil. At Site 1237 off Peru magnetostratigraphy could be applied giving excellent age control with all chrons and subchrons according to GPTS between 8.4 Ma and 12.4 Ma (Chron 4r and Chron 5Ar, [117]). This scale is in good accordance with calcareous nannofossil datums.

The Carbonate Crash Events-Timing and Mechanisms
In this section, compiled data of MARs as well as preservation proxies during the CC-interval and hypotheses about trigger mechanisms will be summarized for each oceanic basin. Existing theories concerning the CC-events in the deep (>3000 m water depth) EEP from Leg 138 will be complemented by regional examples in order to compose NE and SE Pacific records into an overview picture of the Eastern Pacific. Results from Leg 202 deliver insight into intermediate water processes in the EEP (Site 1241) and (original) results of a deep-water record off Peru (Site 1237) reflect processes in the SE Pacific. Additionally, sedimentary processes from the upwelling region along the California margin (Leg 167-Site 1010) provide valuable information about the paleoceanography in the NE Pacific. In particular, proxy data from the Caribbean (Sites 998 and 999) as a key area of the CCevents will be critically reviewed and reinterpreted. The Atlantic is represented by western equatorial Ceara Rise (Site 926) and Benguela upwelling (Site 1085) records. Sites of the depth transect of the Mascarene plateau (707, 709, 710) stand here exemplarily for the low latitude Indian Ocean.
The stress in our reviewing observations lies on the available estimates of preservation of calcareous fossils along with MARs of the main biogenic components (opal and carbonate). The data and sources are listed in Table 1.

The Pacific Carbonate Crash Events
The modern equatorial Pacific surface productivity is characterised by an asymmetry between the east and the west [123]. The Western Pacific Warm water pool (WPWP) forms a stable water mass with a deep thermocline well recorded both in geochemical (isotope) and foraminiferal (assemblages) data [130,131]. In contrast to the moderate carbonate accumulation in the west, the eastern equatorial Pacific (EEP) accumulation of carbonates is mostly influenced by the upwelling intensity and an asymmetric lysocline, which is generally shallow in the east and deeper east of the East Pacific Ridge (EPR). Upwelling waters in the EEP are sourced from the Equatorial Undercurrent (EUC), which has its origin in the west Pacific and combines waters of both hemispheres (e.g., [132,133]). The eastern boundary currents (Humboldt and California current) form the prerequisite for coastal upwelling along the western boundaries of the two Americas. Oceanic productivity is furthermore distinctly influenced by the ENSO (El Niño -Southern Oscillation) phenomenon [134,135]. A strong El Niño event significantly reduces deep-water upwelling by deepening the thermocline, and finally leads to the breakdown of plankton productivity of upwelling regions in the eastern Pacific [14]. These modern patterns might have their roots in the Miocene influenced through gateway configurations [136]. The emergence of the Indonesian Seaway (IS) blocks the westward flow of warm waters to the Indian Ocean and allows warm waters to pile up to form the WPWP [115]. Generally, the strength of circulation is driven by the trade winds and boundary currents along the continents, and a slackening of these currents would tend to reduce the intensity of upwelling waters and result in a deeper thermocline in the eastern Pacific. Thus, the questions concerning the Miocene carbonate accumulation in the equatorial Pacific are: Can we find the oceanographic patterns compared to today's features and how relevant are these surface processes in comparison to deep water circulation changes?

The Eastern Equatorial Pacific
One of the first studies reporting unusually low carbonate contents during the Late Miocene (i.e., between 9.6 and 9.2 Ma) is by Vincent [137] on the section of DSDP Site 310 on the Hess Rise and in the northern central Pacific. ODP Leg 138 recovered excellent cores from the EEP. Changes in MAR Carbonate first followed a stepwise decline in various sites of the EEP from 11.2 to 9.5 Ma, with pronounced changes in the Guatemala basin (see Figure 3 Site 844) and Peru basin (Figure 3 Site 846) east of the EPR. This interval is also characterised by shifts in dominance of carbonate versus opal in sediments, indicating changing surface water ecology as well as deep circulation. At 9.5 Ma, MAR Carbonate dropped to zero reaching the carbonate crash nadir. Evidence from nannofossil preservation [4] and MAR Carbonate point to a rapid rise of the local CCD of about 800 m to water depth of about 3400-3200 m at 10 Ma [1]. After this brief dissolution event the carbonate accumulation did not recover but the preservation, as indicated by nannofossil and foraminifer indices, increase after 8.9 Ma [4]. West of the EPR in the Central Pacific basin, the MAR Carbonate also shows minima centred at 9.5 Ma. Here changes are less clear with milder influence of dissolution as indicated by better nannofossil preservation and no obvious change to higher opal/carbonate ratios like in the Guatemala basin.
enriched subsurface waters to the EEP. If the carbonate crash events in the EEP were related to the emergence of this circulation pattern one would find El Niño condition during a weak WPWP and better stratification in the west Pacific surface waters.
According to the above model, the dissolution of carbonates by corrosive deep waters during the carbonate crash event was delayed in the early phase (at about 11.5 Ma) through high productivity in a La Niña like phase in the EEP, which would explain the offset in CC-events between the Caribbean and Pacific. With regard to this hypothesis and a change from La Niña to El Niño conditions at about 10 to 9 Ma the question arises if these shifts are imprinted in the productivity along eastern boundary currents in other regions of the Eastern Pacific as found during modern El Niño events (e.g., [141,142]). In the following part, these questions will be addressed to cores from the California upwelling and from the SE Pacific.  Table 1. Color bars refer to the preservation state of planktonic foraminifera as found in the ODP Initial Site Reports.

The California Upwelling (Site 1010)
The California upwelling system was investigated during Leg 167. Sediments of Site 1010 (offshore) and 1021 (coastal upwelling) recorded a major drop in opal contents starting at 11.5 Ma while at the same time MARCarbonate dropped [61]. Despite the probable beginning of NADW circulation, which would have generally enhanced silica availability in the Pacific, certain intervals favoured offshore and others coastal sedimentation of opal. The coastal to offshore opal fractionation was enhanced through low temperature gradients among latitudes damping the strength of the California current accompanied by a deep nutricline [64]. According to this mechanism, the time interval from 12-11.5 Ma can  Table 1. Color bars refer to the preservation state of planktonic foraminifera as found in the ODP Initial Site Reports.
Another characteristic feature of the sediments is the occurrence of laminated diatom mat deposits in the near equator sites (sites 849-851) of the Central Pacific basin at the beginning and the end of the CC event ( Figure 3 indicated by stars). These sediments must have probably formed as rapid fallouts of equatorial frontal zones developed during La Niña-like events [89,115]. The compact structure of diatom mats furthermore prohibited bioturbation and thus laminae appear well preserved [94,138]. However, the diatom abundance was reduced during the crash nadir and radiolaria dominated the opal instead.
Lyle, et al. [1] concluded that the restriction of carbonate saturated deep-water flows from the Atlantic to the Pacific through the emerging CAS would trigger the dissolution in the EEP. As an alternative scenario, deep water formation in the North Atlantic would result in basin-to-basin fractionation, and displacement of corrosive Southern Component Water (SCW) towards the Pacific in turn triggering enhanced dissolution in the deep EEP. However, the comparing of timing of these events could not lead to a consistent picture concerning NADW production [4].

The West Equatorial Pacific (Site 806)
The restriction of the flow of waters through the IS and the CAS is believed to have had a major impact on the circulation system including the distribution of temperature and nutrients in surface and near surface waters in the Pacific. The development of the WPWP, the Equatorial Counter-current system and alternating hydrographic features similar to modern El Niño and La Niña patterns are placed within the middle to late Miocene [14,89,115,139,140]. These models stress the productivity component in the highly dynamic equatorial system and represent an accomplishment to the models solely based on circulation-dissolution explanations. Li, et al. [140] place the establishment of the WPWP in the South China Sea at 10 Ma as evidenced by increased abundance of mixed layer foraminifera and the decline of deeper living foraminifera (i.e., the extinction of Globoquadrina dehiscens). Similar results in the equatorial region place the initiation of the WPWP at 8 Ma [139].
Nathan and Leckie [115] found evidence for the development of a Proto-WPWP at 11.6 to 9.6 Ma based on foraminifer faunal analyses and stable isotope gradients of mixed layer and thermocline species at Site 806 (2520 m water depth) on the Ontong Java Plateau (see Figure 4). Thus, the establishment of a modern like circulation pattern influenced by the closure of the IS and the arising possibility of piling up warm waters in the west Pacific could have occurred that early. The authors mention that the equatorial undercurrent would have been strengthened by the evolution of a WPWP like today bringing nutrient enriched subsurface waters to the EEP. If the carbonate crash events in the EEP were related to the emergence of this circulation pattern one would find El Niño condition during a weak WPWP and better stratification in the west Pacific surface waters. be interpreted in terms of increased strength in current intensity favouring offshore opal sedimentation, while a slackening of the currents occurred from 11.5 to 10 Ma inducing coastal fractionation. This kind of record is associated with a general decreasing trend in opal concentration reaching its minimum at 10 Ma at both sites (shown for Site 1010, see Figure 3). Interestingly, here laminated diatom rich sediments appear as well in similar intervals in comparison to the EEP (see Figure 3). The preservation of planktonic foraminifera is poor except of an interval below 148 mbsf (approximately 11.5 Ma), where they appear well preserved within diatom rich sediments and species composition indicates an upwelling sequence [122]. The enhanced preservation of foraminifera is probably a result of rapid burial through particle loading. These results strengthen the hypothesis of a general decreasing productivity trend, as it can be associated with an El Niño phase, influencing both opal and MARCarbonate at Site 1010 at the same time when carbonate accumulation dropped in the EEP.  According to the above model, the dissolution of carbonates by corrosive deep waters during the carbonate crash event was delayed in the early phase (at about 11.5 Ma) through high productivity in a La Niña like phase in the EEP, which would explain the offset in CC-events between the Caribbean and Pacific. With regard to this hypothesis and a change from La Niña to El Niño conditions at about 10 to 9 Ma the question arises if these shifts are imprinted in the productivity along eastern boundary currents in other regions of the Eastern Pacific as found during modern El Niño events (e.g., [141,142]). In the following part, these questions will be addressed to cores from the California upwelling and from the SE Pacific.

The California Upwelling (Site 1010)
The California upwelling system was investigated during Leg 167. Sediments of Site 1010 (offshore) and 1021 (coastal upwelling) recorded a major drop in opal contents starting at 11.5 Ma while at the same time MAR Carbonate dropped [61]. Despite the probable beginning of NADW circulation, which would have generally enhanced silica availability in the Pacific, certain intervals favoured offshore and others coastal sedimentation of opal. The coastal to offshore opal fractionation was enhanced through low temperature gradients among latitudes damping the strength of the California current accompanied by a deep nutricline [64]. According to this mechanism, the time interval from 12-11.5 Ma can be interpreted in terms of increased strength in current intensity favouring offshore opal sedimentation, while a slackening of the currents occurred from 11.5 to 10 Ma inducing coastal fractionation. This kind of record is associated with a general decreasing trend in opal concentration reaching its minimum at 10 Ma at both sites (shown for Site 1010, see Figure 3). Interestingly, here laminated diatom rich sediments appear as well in similar intervals in comparison to the EEP (see Figure 3). The preservation of planktonic foraminifera is poor except of an interval below 148 mbsf (approximately 11.5 Ma), where they appear well preserved within diatom rich sediments and species composition indicates an upwelling sequence [122]. The enhanced preservation of foraminifera is probably a result of rapid burial through particle loading. These results strengthen the hypothesis of a general decreasing productivity trend, as it can be associated with an El Niño phase, influencing both opal and MAR Carbonate at Site 1010 at the same time when carbonate accumulation dropped in the EEP.

The Intermediate EEP (Site 1241) and the South East Pacific (Site 1237)
Leg 202 offers two sites to test hypotheses concerning the influence of deep and surface water processes during the CC-events in the EP. Site 1241 is positioned at 2040 m water depth above the lysocline on the Coccos Ridge in the EEP and displays the shallowest sedimentary record in comparison to the deep sites of Leg 138 (all deeper than 3000 m water depth). The calcareous microfossil tests show evidence of strong dissolution during the oldest part of recovered sediments (11-10 Ma), when opal and TOC concentrations were relatively high as well as MAR Carbonate (Figure 3). After 10 Ma, tests of foraminifera and nannofossil show moderate preservation associated with increasing carbonate content but a drop in MAR Carbonate as the site moves away from the high productivity equatorial belt. The strong dissolution event ends here at the time when the rapid drop in MAR Carbonate occurs in the Guatemala and Peru basin. Here the preservation of foraminifera is not supported by the sedimentation of diatom rich layers like in Site 1010. In this shallower site, there might have been a higher probability of organic matter arriving at the sea floor and building a benthic fluff layer in which dissolution of foraminifera fossils might have be accelerated in comparison to deeper sites [143]. In the deeper sites of the EP, therefore, the presence of diatoms might simply enhance the preservation probability by rapid settling of aggregates and particle clusters, which also in their extreme form of diatom mats hinder bioturbation and support rapid burial as positive factors in preservation of calcareous shells. Thus, dissolution at Site 1241 was likely the result of organic matter degradation; however, after 10 Ma the MAR Carbonate record shows a minimum at about 9.5 Ma contemporaneous with the nadir of EEP CC-events, suggesting a general low in carbonate productivity followed by an increase afterwards.
Site 1237 (Naszca Ridge, 3212 m water depth) is bathed by Antarctic-Circumpolar-Deep water and was mostly moving latitudinally during the last 10 Ma, coming closer to the modern high productivity belt of the Peru/Chile upwelling [121]. The MAR Carbonate record experienced a first drop at about 11 Ma and stayed low until it recovers between 10.2 and 9.8 Ma, followed by a second slight decrease. The preservation of planktonic foraminifera was preferably good until 10.3 Ma, after which the decreasing trend in the number of whole foraminifera tests and the enrichment of the coarse fraction with noncalcareous components (e.g., sponge spicules, radiolarians, ash) point to intensification of dissolution. According to findings at Site 1237, the decreasing MAR Carbonate might be assigned to generally lowered productivity until 10.3 Ma because during this period decreasing MAR Carbonate are not clearly accompanied by decreased carbonate preservation. These findings support the hypothesis of a productivity decrease in the early phase of the Miocene CC-events and point to the intensification of deep-water corrosiveness increasing at about 10.3 to 10 Ma.

The Caribbean Sea
The Caribbean basin is characterised by influx of incoming carbonate corrosive waters as a mixture of mainly AAIW and upper NADW (e.g., [19,144]). The loop current is an important part in the global thermohaline circulation as surface waters gain salinity through evaporation [145], which is believed to have enhanced deep water formation in the North Atlantic (e.g., [146]).
Roth, et al. [5] modified the alternative model of Lyle, et al. [1] calling the influx of corrosive SCW responsible for dissolution in the equatorial Pacific. The tectonic situation in the Caribbean (emergence of the CAS and opening of intrabasinal seaways like the Pedro Channel) could have allowed for a strengthening of the Caribbean Current and establishment of the loop current. Based upon the observation of different coccolith assemblages present on either side of the CAS during an interval from 10.7 to 9.4 Ma, the Pacific-Atlantic connection might have been weakened [127]. At the same time the establishment of surface water connection among sites 999 and 998 and the possible initiation of the loop current was postulated based on similar assemblages.
According to Roth, et al. [5], the changing circulation at the middle to late Miocene transition is well recorded in the contrasting carbonate preservation pattern observed in the Caribbean basin, the EEP, and the western equatorial Atlantic. In analogy to the theory of glacial interglacial preservation cycles [147] Roth, et al. [5] suggested that North Atlantic Intermediate Water (NAIW) flowed over shallow to intermediate depth sills on the Atlantic side of the Caribbean Basin during times of enhanced carbonate preservation (i.e., comparably to the glacial circulation of the Quaternary), while corrosive southern sourced intermediate water (ancient equivalent to AAIW) overflowed the sills during intervals of carbonate dissolution (i.e., comparably to interglacial circulation in the Quaternary). This would furthermore induce that during times of enhanced NCW formation the preservation becomes worse in the Caribbean (corrosive intermediate SCW), the EEP (displacement of SCW towards this region) and better in the deep Atlantic (NCW replacing SCW). This configuration was based on the correlation between carbonate MAR minima ( Figure 5), periods of more intense Northern Component Water (NCW) production [21], and associated with a closed deep water exchange between the Atlantic and the Pacific Oceans.
The study of Newkirk and Martin [15] that compared carbonate mass accumulation rate patterns with evidence from Nd isotopes from fossil fish teeth and debris at Sites 998 and 999 in the Caribbean and Sites 846 and 1241 in the eastern equatorial Pacific support the assumption that waters sourced from the Pacific dominated the flow into the Caribbean during the Miocene Caribbean CC-events. A gradual decrease in carbonate MAR Carbonate and an associated increase in εNd values at Site 999 prior to the Caribbean crash ( Figure 5) provide evidence for the introduction of a more corrosive Pacific intermediate water mass into the Caribbean as the CAS shoaled to critical depths for west to east flow. During the Caribbean carbonate crash (12-9 Ma), highly variable εNd values and MAR Carbonate record pulses of corrosive Pacific waters that filled the deep Caribbean basin. These pulses of Pacific through-flow correlate well with NCW production, suggesting that NCW production could occur with an open CAS and that flow system in the Caribbean region is linked to the global circulation patterns. After the Caribbean carbonate crash, εNd values gradually shifted to less radiogenic values, indicating a reduction in the amount of Pacific water flowing into the Caribbean coincident with the shoaling of the Isthmus of Panama.
Geosciences 2021, 11, 94 14 of 25 Ma) of the CAS, the carbonate compensation depth was deepening in the EEP rather than shoaling [153] according to the sedimentary data.

The Atlantic Ocean
The middle to late Miocene Atlantic is characterised by the fade of opal [64] and the onset of significant deep water formation in the North as well as initiated upwelling off Southwest Africa [154]. Krammer, et al. [10] and Kastanja, et al. [60], studying the Miocene records of ODP Sites 1085 and 1087 on the continental margin off Namibia, found severe reductions in carbonate contents and carbonate accumulation rates. These were attributed rather to the reduction of coccolith productivity than to dissolution. Kastanja, et al. [60] showed that a first major drop in CaCO3 concentration between 10.4 and 10.1 Ma was related mainly to changes in calcareous nannoplankton production, while another drop between 9.6 and 9 Ma is thought to have been triggered by a combination of plankton export production changes and dilution; the latter, presumably occurring in response to high shelf supply in this region during global lowering of sea level [8,60] (Figure 6). The carbonate dissolution is further supported by the increasing benthic/plankton (B/P) foraminiferal ratios, since planktic foraminiferal tests are generally thin-walled and more susceptible to dissolution than robust thick-walled benthic foraminifers [155,156]. That increase observed at Site 1085 from about 9.4 Ma is attributed to enhanced dissolution due to supply from the shelf [8] and possibly aridification with first major dust supply at 9.6 Ma [157].
Miocene South Atlantic sediments on the Walvis Ridge display short-term dissolution events that were closely related to variations in NCW circulation in the deep circulation loop of the South Atlantic, as shown in Kastanja and Henrich [9]. They registered overall good to moderate preservation in the Miocene sections evidencing persistent NCW supply to this southern location. However, some decreases of preservation at 11.6 and 10.4 Ma were found to coincide with Miocene glacial events (Mi-events), suggesting an increase of SCW influence during these intervals, which occurred as a response to the intensification of Antarctic ice sheet development. At 10.4 Ma, a change to overall better preservation points to a weakening of SCW that occurred as a response to the strengthening of NCW. The comparison of preservation records (sand fraction) with carbonate accumulation of sites 999 and 998 suggests that the carbonate preservation is persistently worse, while the MAR Carbonate already recovered at 10 Ma. The sand fraction contents of Site 999 and 998 remain relatively low until 9.4 Ma. This trend is in accordance with lower values of benthic delta 13 C at both sites and Nd isotope trends from Site 998 indicating phases of dominance of corrosive waters from probably Pacific origin and a trend to better preservation afterwards. The absence of a particular abrupt excursion in the carbon isotopic curve towards lighter values does not implicate any global rise in greenhouse gas concentrations, such as methane release [148], which was proposed to explain reductions in seafloor carbonate by dissolution in the geologic past [149]. Another interesting pattern is the high abundance of discoaster nannoliths during the carbonate crash at sites 998 and 999 [127]. Although their abundance has originally been interpreted as an ecological signal, their apparent enrichment within the sediment may well be due to dissolution of other material [108], and may thus serve as another preservation proxy indicating that the CC-ended at 9.4 Ma. This would have the consequence that the establishment of the loop current cannot be inferred from assemblages. The meaning of such high abundances as an ecologic signal might point alternatively to an interval of low plankton productivity [111,150].
The complete closure of the Isthmus might not necessarily be a precondition to enhance carbonate preservation in the deep Caribbean basin. The comparison of the carbonate mass accumulation and preservation patterns of the Caribbean sites [125] and those of the western equatorial Atlantic (i.e., Ceara Rise; Figure 4) reveals that the anti-correlation postulated by Roth, et al. [5] for these two contrasting regions is not evident from the records (see Figure 5). Instead, the records of Site 999 and 926 seem to be in phase, however much worse carbonate preservation was prevailing at Site 999 compared to the deeper Site 926, indicating the presence of more corrosive bottom waters. The similar phase trends in the preservation records might suggest that the NADW formation did not trigger a return flow and influx of corrosive waters but seem to weaken/replace/dilute the flux of corrosive waters into the Caribbean. Modelling results showed that NCW formation is possible with an open CAS [151]. However, NCW% estimates and preservation records from the Caribbean and Ceara Rise differ in this important interval as discussed in the next section. Simulation for a restricted exchange through the CAS predicts a drastic shoaling of the lysocline in the EEP [152]. However, at about the time of final closure (3 Ma) of the CAS, the carbonate compensation depth was deepening in the EEP rather than shoaling [153] according to the sedimentary data.

The Atlantic Ocean
The middle to late Miocene Atlantic is characterised by the fade of opal [64] and the onset of significant deep water formation in the North as well as initiated upwelling off Southwest Africa [154]. Krammer, et al. [10] and Kastanja, et al. [60], studying the Miocene records of ODP Sites 1085 and 1087 on the continental margin off Namibia, found severe reductions in carbonate contents and carbonate accumulation rates. These were attributed rather to the reduction of coccolith productivity than to dissolution. Kastanja, et al. [60] showed that a first major drop in CaCO3 concentration between 10.4 and 10.1 Ma was related mainly to changes in calcareous nannoplankton production, while another drop between 9.6 and 9 Ma is thought to have been triggered by a combination of plankton export production changes and dilution; the latter, presumably occurring in response to high shelf supply in this region during global lowering of sea level [8,60] (Figure 6). The carbonate dissolution is further supported by the increasing benthic/plankton (B/P) foraminiferal ratios, since planktic foraminiferal tests are generally thin-walled and more susceptible to dissolution than robust thick-walled benthic foraminifers [155,156]. That increase observed at Site 1085 from about 9.4 Ma is attributed to enhanced dissolution due to supply from the shelf [8] and possibly aridification with first major dust supply at 9.6 Ma [157].
Miocene South Atlantic sediments on the Walvis Ridge display short-term dissolution events that were closely related to variations in NCW circulation in the deep circulation loop of the South Atlantic, as shown in Kastanja and Henrich [9]. They registered overall good to moderate preservation in the Miocene sections evidencing persistent NCW supply to this southern location. However, some decreases of preservation at 11.6 and 10.4 Ma were found to coincide with Miocene glacial events (Mi-events), suggesting an increase of SCW influence during these intervals, which occurred as a response to the intensification of Antarctic ice sheet development. At 10.4 Ma, a change to overall better preservation points to a weakening of SCW that occurred as a response to the strengthening of NCW.
The western equatorial Atlantic (Ceara Rise, ODP Leg 154) sedimentation patterns are as well influenced by increasing preservation and distinct lows in high frequent carbonate accumulation occurring at 11.8 to 11 Ma and 10.4 to 10.1 Ma at Site 926 [6,125] accompanied by low sand fraction contents [126]. The comparison with Caribbean Site 999 shows similarity in MAR Carbonate and sand fraction, suggesting dissolution and/or productivity pattern in phase [125]. The general trend in preservation based on interpretation of sand fraction contents shows increasing preservation from 11.5 Ma and onwards. The preservation minima are associated with Mi events 5 and 6 suggesting a causal relation to influx of corrosive SCW acting as the dominant deep water source during pulsed cooling. The relation to NCW% estimates by Wright and Miller [21] and Poore, et al. [22] shows that the preservation record is roughly in phase until 10.1 Ma, when the increasing preservation is not reflected in increasing NCW%. The interpretation of carbon isotopes are the basis of NCW% estimates through mixing calculations [158]. The estimates might be confused by the overall low productivity during this time interval (EP, [2,EP,121]), which might prevent formation of reliable benthic δ 13 C gradients among ocean basins. Furthermore, a recent modelling study suggests that the carbon isotope pattern in the Atlantic is not related to the variations in sill depth of the Greenland Scotland ridge, questioning a major tectonic influence on circulation patterns [159].
The preservation minima are associated with Mi events 5 and 6 suggesting a causal rela-tion to influx of corrosive SCW acting as the dominant deep water source during pulsed cooling. The relation to NCW% estimates by Wright and Miller [21] and Poore, et al. [22] shows that the preservation record is roughly in phase until 10.1 Ma, when the increasing preservation is not reflected in increasing NCW%. The interpretation of carbon isotopes are the basis of NCW% estimates through mixing calculations [158]. The estimates might be confused by the overall low productivity during this time interval (EP, [2,EP,121]), which might prevent formation of reliable benthic δ 13 C gradients among ocean basins. Furthermore, a recent modelling study suggests that the carbon isotope pattern in the Atlantic is not related to the variations in sill depth of the Greenland Scotland ridge, questioning a major tectonic influence on circulation patterns [159].

The Indian Ocean
The modern Indian Ocean can be divided into a northern high productivity region (Arabian Sea and the pelagic ocean north of 10°N), which is influenced by the monsoonal

The Indian Ocean
The modern Indian Ocean can be divided into a northern high productivity region (Arabian Sea and the pelagic ocean north of 10 • N), which is influenced by the monsoonal gyre system and the lower productive region of the subtropical gyre to the south. The Indian Ocean was drilled through early ODP Legs 115, 117, 121 and 122, and recovered cores mostly shallower 2500 m water depth with the exception of Leg 115, which offers a depth transect. Peterson, et al. [2] noted that despite the completely different regional settings the carbonate sedimentation and gradients among sites were low during most of the Miocene followed by a marked step in the late Miocene accompanied by reappearing opal components in the sediment (Figure 7).
The initiation or intensification of the monsoon in connection to the uplift of the Himalaya and Tibetian Plateau as a major cause of pronounced oceanic productivity has been rejected (e.g., [160]) as well as the hypothesis of Tethyan outflow water dominating at intermediate depth during the middle to late Miocene [161,162]. The overall similarity suggests a global pattern assigned to sea level fluctuations and shelf-basin fractionation rather than a regional influence like the monsoon [2]. The increase in both carbonate and opal accumulation was named biogenic bloom and reported as well in the Pacific and Atlantic within a broad time interval between~10 and 3.8 Ma [4,8,13,163]. In the Indian Ocean it was accompanied by reappearing opal components and benthic foraminifera indicative of high productivity [160,162,164] in the sediment as indicator for high productivity. How-ever, the shift in depth gradients in the Late Miocene implies changes in circulation as well [2]. Unfortunately, winnowing gives a strong imprint in shallow sediments (Site 707) causing foraminiferal sand contributions of about 50% (Figure 7). In comparison to Site 707, the deeper Site 710 (3812 m) shows about the same MAR Carbonate in accordance with aforementioned overall reduced gradients along the depth transect from 12 to 9.4 Ma right above the LO datum of D. hamatus.
Geosciences 2021, 11, 94 16 of 25 gyre system and the lower productive region of the subtropical gyre to the south. The Indian Ocean was drilled through early ODP Legs 115, 117, 121 and 122, and recovered cores mostly shallower 2500 m water depth with the exception of Leg 115, which offers a depth transect. Peterson, et al. [2] noted that despite the completely different regional settings the carbonate sedimentation and gradients among sites were low during most of the Miocene followed by a marked step in the late Miocene accompanied by reappearing opal components in the sediment (Figure 7). The initiation or intensification of the monsoon in connection to the uplift of the Himalaya and Tibetian Plateau as a major cause of pronounced oceanic productivity has been rejected (e.g., [160]) as well as the hypothesis of Tethyan outflow water dominating at intermediate depth during the middle to late Miocene [161,162]. The overall similarity suggests a global pattern assigned to sea level fluctuations and shelf-basin fractionation rather than a regional influence like the monsoon [2]. The increase in both carbonate and opal accumulation was named biogenic bloom and reported as well in the Pacific and Atlantic within a broad time interval between ~10 and 3.8 Ma [4,8,13,163]. In the Indian Ocean it was accompanied by reappearing opal components and benthic foraminifera in- The age model of Site 710 is based on magnetostratigraphical correlations and resulted in a much earlier placement of the FO datum of D. hamatus than usual. However, the reduced depth gradient might be because of redepositional events in the deeper records but is not purely an artefact of the age model differences and evident among the other sites as described by Peterson, et al. [2]. In the deep sites, the comparison of sand fraction contents leads to interesting results: At Site 710 and 709, the record of coarse fraction suggests a change to better preservation at 9.5 Ma as well as the abundance and preservation data on planktonic foraminifera [165]. The winnowed coarse fraction record of 707 is not shown but despite high coarse fraction contents the preservation of foraminifera here is changing from moderate to good. The reappearance of opal in form of radiolaria date to 9.7 at Site 707 and to 10.3 Ma at site 709.
The proportion of the sand fraction might be controlled by other processes besides dissolution. A size change in foraminifera became evident from 10 Ma on [78]. This trend in low latitudes could have had an effect in the records shown. Winnowing could enhance the coarse fraction or dilution by finer grain sizes (nannofossils, terrigenous fraction) could decrease the coarse fraction. However, the investigation of multiple proxies or the observation of increasing sand fraction along with better preservation and higher abundance of planktonic foraminifera might confirm the interpretation as dissolution index.

Conclusions
The middle to late Miocene carbonate crash events mark a period of major perturbations in the marine carbonate system, which obviously were associated with several steps in reorganization of global deep and intermediate water circulation, affecting various parts of the global ocean basins differently in time and space. Since Miocene atmospheric pCO 2 concentrations are known to be close to preindustrial levels [166] and with no particular excursions during the time of the CC events [167], the observed reductions in seafloor carbonates may mostly be the result of decreased productivity and reorganization of the corrosive deep-water masses rather than dissolution by other exogenic causes.
A review and comparison of the eastern equatorial records seem to strengthen the hypothesis of alternating El Niño/La Niña-like states, that influenced both opal and carbonate accumulation in the EEP on the California margin and off Peru. The emergence of IS probably allowed for the development of a temporary Western Pacific Warm water pool during a La Niña-like state associated with the sedimentation of diatom mats. The alternating productivity in opaline and calcareous plankton were the main cause rather than dissolution at least in the early phase of the carbonate crash events. The deepwater corrosiveness in the Pacific was pronounced at 10.5 to 9.4 Ma and losing influence afterwards. The shallow water record of Site 1241 in the EEP showed severe supralysoclinal dissolution through phases of high productivity.
Phases of reduced calcareous nannoplankton productivity were evident both in the Atlantic (Benguela upwelling site 1085) and Pacific and were probably related to reorganizations of the upper water column and surface circulation changes. In the deep Atlantic the preservation change can be recognized at 10.5 to 10.1 Ma already. Thus 10.5 Ma might be the starting point of prolonged basin to basin fractionation through the onset of the deep convection in the North Atlantic.
Sand fraction contents and preservation of planktonic foraminifera from the deeper Caribbean and the Indian Ocean and Pacific show an increase at 9.3 to 9.5 Ma, which we interpret as the first widespread signal of better preservation indicating a lysocline turning point possibly related to CCD deepening ( Figure 8). The Caribbean CC-events ends at the same time based on reinterpreted coarse fraction contents, δ 13 C signal and Nd isotope evidence as well as nannofossil indicators independent from the step to higher carbonate accumulation at 10 Ma, which was as well recognized at Ceara Rise. The similarity of coarse fraction and mass accumulation records of Ceara Rise and Caribbean, with however worse preservation in the Caribbean, leads to the conclusion that the influence of corrosive Pacific sourced waters in the Caribbean was dampened by NCW formation and not enhanced as previously assumed. The shoaling of the CAS might have contributed at 9.4 Ma to the Atlantic type preservation pattern in the Caribbean. evidence as well as nannofossil indicators independent from the step to higher carbonate accumulation at 10 Ma, which was as well recognized at Ceara Rise. The similarity of coarse fraction and mass accumulation records of Ceara Rise and Caribbean, with however worse preservation in the Caribbean, leads to the conclusion that the influence of corrosive Pacific sourced waters in the Caribbean was dampened by NCW formation and not enhanced as previously assumed. The shoaling of the CAS might have contributed at 9.4 Ma to the Atlantic type preservation pattern in the Caribbean. While the mass of carbonate in all studies is constituted by coccoliths and other nannoliths, the use of the better-known foraminifera dominates the scientific output. However, changes in the fine fraction stable isotopes in the Atlantic, Caribbean and Pacific show drastic decreases from 11.5 to 10.5 Ma [119,168,169], but are still hardly understood and interpretations are contrasting. Understanding these patterns and the further development of high-resolution age models would help to gain better insight into the CCevents regarding considerations on the timing and carbonate budgets.
Supplementary Materials: The following are available online at www.mdpi.com/2076-3263/11/2/94/s1, Table S1: Generated carbonate mass accumulation and foraminifera fragmentation index data for core ODP 1237.  Data Availability Statement: Generated carbonate mass accumulation and foraminifera fragmentation index data for core ODP 1237 are given as supplementary material. While the mass of carbonate in all studies is constituted by coccoliths and other nannoliths, the use of the better-known foraminifera dominates the scientific output. However, changes in the fine fraction stable isotopes in the Atlantic, Caribbean and Pacific show drastic decreases from 11.5 to 10.5 Ma [119,168,169], but are still hardly understood and interpretations are contrasting. Understanding these patterns and the further development of high-resolution age models would help to gain better insight into the CC-events regarding considerations on the timing and carbonate budgets.
Supplementary Materials: The following are available online at https://www.mdpi.com/2076-326 3/11/2/94/s1, Table S1: Generated carbonate mass accumulation and foraminifera fragmentation index data for core ODP 1237.  Data Availability Statement: Generated carbonate mass accumulation and foraminifera fragmentation index data for core ODP 1237 are given as supplementary material.