Search NASA⌕ Search

SEARCH · Search NASA

Results for “Reservoir Modeling”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 37 records · Page 2

Consequences of Insufficient Production of Temporary NOx Reservoirs on Ozone Production in Models During KORUS-AQ

The advent of geostationary satellites, starting with GEMS over South Korea, offer the potential for additional insight into urban air quality. In preparation for these satellites, the joint NIER-NASA KORUS-AQ field campaign in 2016 obtained a large dataset for understanding local air quality and the contribution of local emissions to ozone and PM2.5 pollution. Models partitioned NOx primarily into inorganic nitrate (HNO3 + nitrate aerosol) while observations suggest that acyl peroxy nitrates (PNs) should be an important sink. Models tended to produce ozone too efficiently, which may be due to insufficient alkyl/multifunctional nitrate (ANs) production. Models are unable to simulate the transport of NOx away from the urban center and its impact on ozone production downwind.

Katherine R Travis↗

NASA X-HAB Water Delivery System

As part of the eXploration Systems and Habitation (X-HAB) Academic Innovation Challenge program of National Aeronautics and Space Administration (NASA), The Ohio State University was selected to improve the performance of NASA’s existing Vegetable Production System (VEGGIE), which is a deployable plant growth unit for International Space Station (ISS). During the academic year 2015-2016, The Ohio State University student team developed a passive water delivery system using capillary water transport principle (Jenson et al. 2016). The major design improvement made was directly connecting the water reservoir to the plant-rooting pillows using a single-interface capillary cord design. Harvestable plants were successfully grown from seeds using the single interface system. In addition, Nomex®, a fabric material composed of short nylon based fibers, was identified as the material for wicks. Finally, the water reservoir was modeled as a propellant management device (PMD) to ensure consistent and long term watering of the VEGGIE system. The PMDs are made of materials that utilize surface tension and adhesive forces to improve stability and fluid delivery. The team recommended using a sponge PMD in order to mitigate bubble obstruction, decrease system weight, and ensure reliable water delivery to the capillary interface.

Peter Ling↗

Lead Isotope Compositions of Acid Residues from Olivine-Phyric Shergottite Tissint: Implications for Heterogeneous Shergottite Source Reservoirs

Geochemical studies of shergottites suggest that their parental magmas reflect mixtures between at least two distinct geochemical source reservoirs, producing correlations between radiogenic isotope compositions and trace element abundances. These correlations have been interpreted as indicating the presence of a reduced, incompatible element- depleted reservoir and an oxidized, incompatible- element-enriched reservoir. The former is clearly a depleted mantle source, but there is ongoing debate regarding the origin of the enriched reservoir. Two contrasting models have been proposed regarding the location and mixing process of the two geochemical source reservoirs: (1) assimilation of oxidized crust by mantle derived, reduced magmas, or (2) mixing of two distinct mantle reservoirs during melting. The former requires the ancient Martian crust to be the enriched source (crustal assimilation), whereas the latter requires isolation of a long-lived enriched mantle domain that probably originated from residual melts formed during solidification of a magma ocean (heterogeneous mantle model). This study conducts Pb isotope and trace element concentration analyses of sequential acid-leaching fractions (leachates and the final residues) from the geochemically depleted olivine-phyric shergottite Tissint. The results suggest that the Tissint magma is not isotopically uniform and sampled at least two geochemical source reservoirs, implying that either crustal assimilation or magma mixing would have played a role in the Tissint petrogenesis.

Moriwaki, R.↗

Preliminary Report on U-Th-Pb Isotope Systematics of the Olivine-Phyric Shergottite Tissint

Geochemical studies of shergottites suggest that their parental magmas reflect mixtures between at least two distinct geochemical source reservoirs, producing correlations between radiogenic isotope compositions, and trace element abundances.. These correlations have been interpreted as indicating the presence of a reduced, incompatible-element- depleted reservoir and an oxidized, incompatible-element-rich reservoir. The former is clearly a depleted mantle source, but there has been a long debate regarding the origin of the enriched reservoir. Two contrasting models have been proposed regarding the location and mixing process of the two geochemical source reservoirs: (1) assimilation of oxidized crust by mantle derived, reduced magmas, or (2) mixing of two distinct mantle reservoirs during melting. The former clearly requires the ancient martian crust to be the enriched source (crustal assimilation), whereas the latter requires a long-lived enriched mantle domain that probably originated from residual melts formed during solidification of a magma ocean (heterogeneous mantle model). This study conducts Pb isotope and U-Th-Pb concentration analyses of the olivine-phyric shergottite Tissint because U-Th-Pb isotope systematics have been intensively used as a powerful radiogenic tracer to characterize old crust/sediment components in mantle- derived, terrestrial oceanic island basalts. The U-Th-Pb analyses are applied to sequential acid leaching fractions obtained from Tissint whole-rock powder in order to search for Pb isotopic source components in Tissint magma. Here we report preliminary results of the U-Th-Pb analyses of acid leachates and a residue, and propose the possibility that Tissint would have experienced minor assimilation of old martian crust.

Moriwaki, R.↗

Imaging a Free-Geometry Magma Reservoir Using FEM-Based Linear Inverse Modeling of InSAR Data: Application to Rabaul Caldera, Papua New Guinea

We test an innovative inversion scheme using Green’s functions from an array of pressure sources embedded in finite-element method (FEM) models to image, without assuming an a-priori geometry, the composite and complex shape of a volcano deformation source. We invert interferometric synthetic aperture radar (InSAR) data to estimate the pressurization and shape of the magma reservoir of Rabaul caldera, Papua New Guinea. The results image the extended shallow magmatic system responsible for a broad and long-term subsidence of the caldera between 2007 February and 2010 December. Elastic FEM solutions are integrated into the regularized linear inversion of InSAR data of volcano surface displacements in order to obtain a 3-D image of the source of deformation. The Green’s function matrix is constructed from a library of forward line-of-sight displacement solutions for a grid of cubic elementary deformation sources. Each source is sequentially generated by removing the corresponding cubic elements from a common meshed domain and simulating the injection of a fluid mass flux into the cavity, which results in a pressurization and volumetric change of the fluidfilled cavity. The use of a single mesh for the generation of all FEM models avoids the computationally expensive process of non-linear inversion and remeshing a variable geometry domain. Without assuming an a-priori source geometry other than the configuration of the 3-D grid that generates the library of Green’s functions, the geodetic data dictate the geometry of the magma reservoir as a 3-D distribution of pressure (or flux of magma) within the source array. The inversion of InSAR data of Rabaul caldera shows a distribution of interconnected sources forming an amorphous, shallow magmatic system elongated under two opposite sides of the caldera. The marginal areas at the sides of the imaged magmatic system are the possible feeding reservoirs of the ongoing Tavurvur volcano eruption of andesitic products on the east side and of the past Vulcan volcano eruptions of more evolved materials on the west side. The interconnection and spatial distributions of sources correspond to the petrography of the volcanic products described in the literature and to the dynamics of the single and twin eruptions that characterize the caldera. The ability to image the complex geometry of deformation sources in both space and time can improve our ability to monitor active volcanoes, widen our understanding of the dynamics of active volcanic systems and improve the predictions of eruptions.

Inverse theory↗

Imaging the Complex Geometry of a Magma Reservoir Using Fem-Based Linear Inverse Modeling of InSAR Data: Application to Rabaul Caldera, Papua New Guinea

We test an innovative inversion scheme using Green’s functions from an array of pressure sources embedded in finite-element method (FEM) models to image, without assuming an a-priori geometry, the composite and complex shape of a volcano deformation source. We invert interferometric synthetic aperture radar (InSAR) data to estimate the pressurization and shape of the magma reservoir of Rabaul caldera, Papua New Guinea. The results image the extended shallow magmatic system responsible for a broad and long-term subsidence of the caldera between 2007 February and 2010 December. Elastic FEM solutions are integrated into the regularized linear inversion of InSAR data of volcano surface displacements in order to obtain a 3-D image of the source of deformation. The Green’s function matrix is constructed from a library of forward line-of-sight displacement solutions for a grid of cubic elementary deformation sources. Each source is sequentially generated by removing the corresponding cubic elements from a common meshed domain and simulating the injection of a fluid mass flux into the cavity, which results in a pressurization and volumetric change of the fluid-filled cavity. The use of a single mesh for the generation of all FEM models avoids the computationally expensive process of non-linear inversion and remeshing a variable geometry domain. Without assuming an a-priori source geometry other than the configuration of the 3-D grid that generates the library of Green’s functions, the geodetic data dictate the geometry of the magma reservoir as a 3-D distribution of pressure (or flux of magma) within the source array. The inversion of InSAR data of Rabaul caldera shows a distribution of interconnected sources forming an amorphous, shallow magmatic system elongated under two opposite sides of the caldera. The marginal areas at the sides of the imaged magmatic system are the possible feeding reservoirs of the ongoing Tavurvur volcano eruption of andesitic products on the east side and of the past Vulcan volcano eruptions of more evolved materials on the west side. The interconnection and spatial distributions of sources correspond to the petrography of the volcanic products described in the literature and to the dynamics of the single and twin eruptions that characterize the caldera. The ability to image the complex geometry of deformation sources in both space and time can improve our ability to monitor active volcanoes, widen our understanding of the dynamics of active volcanic systems and improve the predictions of eruptions.

Inverse theory↗

Dynamics of Possible Late Heavy Bombardment Impactor Populations

The existence of the later lunar basins implies the existence of a massive dynamical reservoir that can store small bodies for some 6000 Myr after the Moon formed. We will discuss four recent models for such reservoirs. Additional information is contained in the original extended abstract.

Dones, Luke↗

The SMAP Level 4 Surface and Root-zone Soil Moisture (L4_SM) Product

The Soil Moisture Active and Passive (SMAP) mission is being developed by NASA for launch in 2013 as one of four first-tier missions recommended by the U.S. National Research Council Committee on Earth Science and Applications from Space in 2007. The primary science objectives of SMAP are to enhance understanding of land surface controls on the water, energy and carbon cycles, and to determine their linkages. Moreover, the high resolution soil moisture mapping provided by SMAP has practical applications in weather and seasonal climate prediction, agriculture, human health, drought and flood decision support. In this paper we describe the assimilation of SMAP observations for the generation of the planned SMAP Level 4 Surface and Root-zone Soil Moisture (L4_SM) product. The SMAP mission makes simultaneous active (radar) and passive (radiometer) measurements in the 1.26-1.43 GHz range (L-band) from a sun-synchronous low-earth orbit. Measurements will be obtained across a 1000 km wide swath using conical scanning at a constant incidence angle (40 deg). The radar resolution varies from 1-3 km over the outer 70% of the swath to about 30 km near the center of the swath. The radiometer resolution is 40 km across the entire swath. The radiometer measurements will allow high-accuracy but coarse resolution (40 km) measurements. The radar measurements will add significantly higher resolution information. The radar is however very sensitive to surface roughness and vegetation structure. The combination of the two measurements allows optimal blending of the advantages of each instrument. SMAP directly observes only surface soil moisture (in the top 5 cm of the soil column). Several of the key applications targeted by SMAP, however, require knowledge of root zone soil moisture (approximately top 1 m of the soil column), which is not directly measured by SMAP. The foremost objective of the SMAP L4_SM product is to fill this gap and provide estimates of root zone soil moisture that are informed by and consistent with SMAP observations. Such estimates are obtained by merging SMAP observations with estimates from a land surface model in a soil moisture data assimilation system. The land surface model component of the assimilation system is driven with observations-based surface meteorological forcing data, including precipitation, which is the most important driver for soil moisture. The model also encapsulates knowledge of key land surface processes, including the vertical transfer of soil moisture between the surface and root zone reservoirs. Finally, the model interpolates and extrapolates SMAP observations in time and in space. The L4_SM product thus provides a comprehensive and consistent picture of land surface hydrological conditions based on SMAP observations and complementary information from a variety of sources. The assimilation algorithm considers the respective uncertainties of each component and yields a product that is superior to satellite or model data alone. Error estimates for the L4_SM product are generated as a by-product of the data assimilation system.

Reichle, Rolf↗

Radicals and Reservoirs in the GMI Chemistry and Transport Model: Comparison to Measurements

We have used a three-dimensional chemistry and transport model (CTM), developed under the Global Modeling Initiative (GMI), to carry out two simulations of the composition of the stratosphere under changing halogen loading for 1995 through 2030. The two simulations differ only in that one uses meteorological fields from a general circulation model while the other uses meteorological fields from a data assimilation system. A single year's winds and temperatures are repeated for each 36-year simulation. We compare results from these two simulations with an extensive collection of data from satellite and ground-based measurements for 1993-2000. Comparisons of simulated fields with observations of radical and reservoir species for some of the major ozone-destroying compounds are of similar quality for both simulations. Differences in the upper stratosphere, caused by transport of total reactive nitrogen and methane, impact the balance among the ozone loss processes and the sensitivity of the two simulations to the change in composition.

Douglass, Anne R.↗

Radicals and Reservoirs in the GMI Chemistry and Transport Model: Comparison to Measurements

The most important use of atmospheric chemistry and transport models is to predict the future composition of the atmosphere. The amounts of gases like chlorofluorcarbons, methyl bromide, nitrous oxide and methane are changing and the stratospheric ozone layer will change because these gases are changing. Methyl bromide, nitrous oxide and methane all have natural sources, and also change because of human activity. Chlorofluorcarbons are man-made gases; these are known to decrease stratospheric ozone and future production is banned. They are long-lived gases, and many decades will pass before they are insignificant in the atmosphere. The models are used to predict changes in ozone and other gases; this is a straightforward application. The models must be also tested using observations for the present day atmosphere. This is a challenging task, because the model contains more than 50 species and more than 150 chemical reactions. Data from satellites, ground stations, aircraft and balloons are used to evaluate the model. Different models that are used in international assessments produce different results; in the most recent assessment some predict that ozone will return to 1980 levels by 2025 and others predict that this will not happen until 2050. Since all the parts of the models are conceptually the same, there must be differences in implementation that produce these differences, This work takes a single model, two different sets of winds and temperatures, and repeats the same prediction for the future. Here we compare the results for these two simulations with many observations. The purpose is to identify differences in the model results for the present atmosphere that will lead to different predictions. This sort of controlled comparison will reduce uncertainty in the predictions for stratospheric ozone.

Douglas, Anne R.↗

Water Content of Lunar Alkali Fedlspar

Detection of indigenous hydrogen in a diversity of lunar materials, including volcanic glass, melt inclusions, apatite, and plagioclase suggests water may have played a role in the chemical differentiation of the Moon. Spectroscopic data from the Moon indicate a positive correlation between water and Th. Modeling of lunar magma ocean crystallization predicts a similar chemical differentiation with the highest levels of water in the K- and Th-rich melt residuum of the magma ocean (i.e. urKREEP). Until now, the only sample-based estimates of water content of KREEP-rich magmas come from measurements of OH, F, and Cl in lunar apatites, which suggest a water concentration of < 1 ppm in urKREEP. Using these data, predict that the bulk water content of the magma ocean would have <10 ppm. In contrast, estimate water contents of 320 ppm for the bulk Moon and 1.4 wt % for urKREEP from plagioclase in ferroan anorthosites. Results and interpretation: NanoSIMS data from granitic clasts from Apollo sample 15405,78 show that alkali feldspar, a common mineral in K-enriched rocks, can have approx. 20 ppm of water, which implies magmatic water contents of approx. 1 wt % in the high-silica magmas. This estimate is 2 to 3 orders of magnitude higher than that estimated from apatite in similar rocks. However, the Cl and F contents of apatite in chemically similar rocks suggest that these melts also had high Cl/F ratios, which leads to spuriously low water estimates from the apatite. We can only estimate the minimum water content of urKREEP (+ bulk Moon) from our alkali feldspar data because of the unknown amount of degassing that led to the formation of the granites. Assuming a reasonable 10 to 100 times enrichment of water from urKREEP into the granites produces an estimate of 100-1000 ppm of water for the urKREEP reservoir. Using the modeling of and the 100-1000 ppm of water in urKREEP suggests a minimum bulk silicate Moon water content between 2 and 20 ppm. However, hydrogen loss was likely very significant in the evolution of the lunar mantle. Conclusions: Lunar granites crystallized between 4.3-3.8 Ga from relatively wet melts that degassed upon crystallization. The formation of these granites likely removed significant amounts of water from some mantle source regions, e.g. later mare basalts predicting derivation from a mantle with <10 ppm water. However, this would have been a heterogeneous pro-cess based on K distribution. Thus some, if not most of the mantle may not have been devolatilized by this process; as seen by water in volcanic glasses and melt inclusions.

Mills, R. D.↗

How Does Water Delivery System Design Impact the Microbial Load of Salad Crops?

In a microgravity setting, such as the environment aboard the International Space Station (ISS), an ideal plant water delivery system is one that can grow edible crops with minimal resource consumption and minimal risk to crew members. There are also concerns associated with the ability to control fluid escape and biofilm formation resulting in potential dangers to systems, crops, or crewmembers. To identify an appropriate system, candidate systems were assembled and operated under simulated ISS environmental conditions (T,CO2,and RH) with red romaine lettuce (Lactuca sativa cultivar 'Outredgeous') as a model crop. Fluid reservoirs and randomly selected planting sites were sampled every seven days until maturity at which point edible plant biomass and root samples were also taken. Heterotrophic bacteria and fungi growth patterns throughout each planting cycle were determined by plate counts on appropriate agar media. The candidate systems were compared to a classic hydroponics system as a control and harvested crops were compared to controls as well as Veggie-grown and market produce. Plants harvested from candidate systems yielded lower average heterotrophic bacteria and fungi per gram of plant mass levels when compared to market and Veggie samples as well as those from the control system. Additional studies to evaluate the system sanitation regimen as well as testing additional crops should be considered to aid in the selection of an ideal system.

Curry, Aaron B.↗

The mean age of mantle and crustal reservoirs

An analysis of two transport models for trace elements in mantle and crust evolution and related to the observed abundance patterns. In model I, continents are derived by melt extraction over the earth history from undepleted mantle, and the residue forms a depleted mantle which is the current source of mid-ocean ridge basalts. In model II, new additions to continents are derived from a mantle reservoir 2, which becomes more depleted by repeated extraction of melts. The isotopic composition and concentrations of trace elements are shown to reduce to simple mathematical expressions which permit calculations of basic evolutionary parameters. The mean age of the crust mass and isotopic data for the continental crust and the mantle are discussed, concluding with a consideration of their difference in compositions of newly derived crust.

Jacobsen, S. B.↗

Ion Microprobe Measurements of Comet Dust and Implications for Models of Oxygen Isotope Heterogeneity in the Solar System

The oxygen isotopic compositions of anhydrous minerals in carbonaceous chondrites reflect mixing between a O-16-rich and O-17, O18-rich reservoir. The UV photodissociation of CO (i.e. selfshielding) has been proposed as a mass-independent mechanism for producing these isotopically distinct reservoirs. Self-shielding models predict the composition for the CO gas reservoir to be O-16-rich, and that the accreting primordial dust was in isotopic equilibrium with the gaseous reservoir [1, 2]. Self-shielding also predicts that cometary water, presumed to represent the O-17, O-18-rich reservoir, should be enriched in O-17 and O-18, with compositions of 200 -1000per mille, and that the interaction with this O-17, O-18-rich H2O reservoir altered the compositions of the primordial dust toward planetary values. The bulk composition of the solar nebula, which may be an approximation to the 16O-rich gaseous reservoir, has been constrained by the Genesis results [3]. However, material representing the O-17, O-18-rich end-member is rare [4], and dust representing the original accreting primordial dust has been challenging to conclusively identify in current collections. Anhydrous dust from comets, which accreted in the distal cold regions of the nebula at temperatures below approximately 30K, may provide the best opportunity to measure the oxygen isotope composition of primordial dust. Chondritic porous interplanetary dust particles (CP-IDPs) have been suggested as having cometary origins [5]; however, until direct comparisons with dust from a known comet parent body were made, link between CP-IDPs and comets remained circumstantial. Oxygen isotope analyses of particles from comet 81P/Wild 2 collected by NASA's Stardust mission have revealed surprising similarities to minerals in carbonaceous chondrites which have been interpreted as evidence for large scale radial migration of dust components from the inner solar nebula to the accretion regions of Jupiter- family comets [6]. These studies have been largely focused on the coarse-grained terminal particles extracted from aerogel collectors; hypervelocity capture into aerogel resulted in fine-grained material that was melted and intimately mixed with the SiO2 capture medium. Hypervelocity impacts into Al foils surrounding the aerogel tiles produced impact craters that captured material from the impactor without significant oxygen contamination, allowing for analysis of both the coarse and fine-grained components of the Wild 2 dust. To date, no particles with definitive hydrated mineralogy have been observed in Stardust samples, though this may be a result of alteration due to hypervelocity capture. High-carbon hydrated CS-IDPs have been suggested as resulting from the aqueous alteration of CP-IDPs [7], and may retain evidence for interaction with O-17, O-18-enriched "cometary" water predicted by CO self-shielding. Here we present results of oxygen isotope measurements of twelve Stardust foil craters and four C-rich hydrated IDPs [8, 9], and discuss implications for models of oxygen isotope heterogeneity in the early solar system.

Snead, C. J.↗

The Sm-(146,147)-Nd-(142,143) formation interval for the lunar mantle and implications for lunar evolution

Small anomalies in the isotopic abundance of Nd-142 have been measured for two A17 high-Ti basalts, ilmenite basalt 12056, olivine-pigeonite basalt 12039, feldspathic basalt 12038, and two KREEP basalts. These anomalies correlate with Sm-147/Nd-144 for the basalt source regions as calculated from initial Nd-143/Nd-144 ratios in the basalts, and are interpreted to be from decay of Sm-146 (t sub 1/2 = 103 Ma) in distinct lunar mantle reservoirs. A three-stage model for evolution of Nd-143/Nd-144 and Nd-142/Nd-144 yields reservoir Sm-147/Nd-144 ratios which, with the Nd-142/Nd-144 ratios in the basalts, form a 'mantle isochron' giving a lunar mantle formation interval of 94+2230 Ma (2c(rho)). Calculated reservoir Sm/Nd ratios are in the range expected from some earlier models of basalt petrogenesis. The isochron value of Nd-142/Nd-144 at Sm-147/Nd-144 sub CHUR = 0.1967 is within error limits of the average Nd-142/Nd-144 measured for an L6 chondrite, an H5 chondrite, and the Orgueil carbonaceous chondrite. Evolution of Nd-143 and Nd-142 for high-Ti basalt 70135 was modeled precisely, starting from chondritic relative REE and Nd-isotopic abundances and using the initial (Sm-146/Sm-144) sub 0 ratio inferred from a previous study of angrite LEW86010 as the initial solar system value of this parameter. We infer that the initial Sm/Nd ratio in precursor lunar materials was very nearly chondritic (within approximately 8 percent) prior to lunar differentiation.

Nyquist, L. E.↗

Study of Nonclassical Fields in Phase-Sensitive Reservoirs

We show that the reservoir influence can be modeled by an infinite array of beam splitters. The superposition of the input fields in the beam splitter is discussed with the convolution laws for their quasiprobabilities. We derive the Fokker-Planck equation for the cavity field coupled with a phase-sensitive reservoir using the convolution law. We also analyze the amplification in the phase-sensitive reservoir with use of the modified beam splitter model. We show the similarities and differences between the dissipation and amplification models. We show that a super-Poissonian input field cannot become sub-Poissonian by the phase-sensitive amplification.

Kim, Myung Shik↗

The Martian climate: Energy balance models with CO2/H2O atmospheres

Progress in the development of a multi-reservoir, time dependent energy balance climate model for Mars driven by prescribed insolation at the top of the atmosphere is reported. The first approximately half-year of the program was devoted to assembling and testing components of the full model. Specific accomplishments were made on a longwave radiation code, coupling seasonal solar input to a ground temperature simulation, and conceptualizing an approach to modeling the seasonal pressure waves that develop in the Martian atmosphere as a result of sublimation and condensation of CO2 in polar regions.

Hoffert, M. I.↗

Ocean alkalinity and the Cretaceous/Tertiary boundary

A biogeochemical cycle model resolving ocean carbon and alkalinity content is applied to the Maestrichtian and Danian. The model computes oceanic concentrations and distributions of Ca(2+), Mg(2+), and Sigma-CO2. From these values an atmospheric pCO2 value is calculated, which is used to estimate rates of terrestrial weathering of calcite, dolomite, and calcium and magnesium silicates. Metamorphism of carbonate rocks and the subsequent outgassing of CO2 to the atmosphere are parameterized in terms of carbonate rock reservoir sizes, total land area, and a measure of overall tectonic activity, the sea-floor generation rate. The ocean carbon reservoir computed by the model is used with Deep Sea Drilling Project (DSDP) C-13 data to estimate organic detrital fluxes under a variety of ocean mixing rate assumptions. Using Redfield ratios, the biogenic detrital flux estimate is used to partition the ocean carbon and alkalinity reservoirs between the mixed layer and deep ocean. The calcite flux estimate and carbonate ion concentrations are used to determine the rate of biologically mediated CaCO3 titration. Oceanic productivity was severely limited for approximately 500 kyr following the K/T boundary resulting in significant increases in total ocean alkalinity. As productivity returned to the ocean, excess carbon and alkalinity was removed from the ocean as CaCO3. Model runs indicate that this resulted in a transient imbalance in the other direction. Ocean chemistry returned to near-equilibrium by about 64 mybp.

Caldeira, K. G.↗