Search NASA⌕ Search

SEARCH · Search NASA

Results for “Model Basin”

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 325 records · Page 18

Assimilation of Terrestrial Water Storage from GRACE in a Snow-Dominated Basin

Terrestrial water storage (TWS) information derived from Gravity Recovery and Climate Experiment (GRACE) measurements is assimilated into a land surface model over the Mackenzie River basin located in northwest Canada. Assimilation is conducted using an ensemble Kalman smoother (EnKS). Model estimates with and without assimilation are compared against independent observational data sets of snow water equivalent (SWE) and runoff. For SWE, modest improvements in mean difference (MD) and root mean squared difference (RMSD) are achieved as a result of the assimilation. No significant differences in temporal correlations of SWE resulted. Runoff statistics of MD remain relatively unchanged while RMSD statistics, in general, are improved in most of the sub-basins. Temporal correlations are degraded within the most upstream sub-basin, but are, in general, improved at the downstream locations, which are more representative of an integrated basin response. GRACE assimilation using an EnKS offers improvements in hydrologic state/flux estimation, though comparisons with observed runoff would be enhanced by the use of river routing and lake storage routines within the prognostic land surface model. Further, GRACE hydrology products would benefit from the inclusion of better constrained models of post-glacial rebound, which significantly affects GRACE estimates of interannual hydrologic variability in the Mackenzie River basin.

Forman, Barton A.↗

Basin-scale analysis of Mokelumne River Formation for multi-well CO 2 injection

Large-scale carbon sequestration will likely require multiple projects injecting CO 2 into the same subsurface formation, raising concerns about safe operation and efficient use of storage capacity. This study evaluates the long-term response of the Mokelumne River Formation in California’s Sacramento Basin to multi-megaton CO 2 injection using three geologic models of the formation and the open-source simulator GEOS. The analysis focuses on three aspects of reservoir performance: (1) average pressure increase and dissolved CO 2 mass after 30 years for varying well counts and injection rates, (2) pressure interference in a multi-well configuration, and (3) dynamic storage capacity with identification of overpressure-prone regions. The results show that average formation pressure increases linearly with injected mass, while CO 2 dissolution exhibits mixed scaling: approximately linear with the number of wells but sublinear with injection rate, indicating that distributing injection across more wells enhances dissolution more effectively than increasing per-well rates. Pressure-interference effects are significant, with lower-permeability conditions delaying their onset but amplifying their magnitude at later times. Dynamic capacity, defined by the first occurrence of pressure exceeding the local overburden-based limit anywhere in the formation, varies across geologic models and assumed overburden pressure gradients. A lower fidelity geologic model predicts nearly twice the storage capacity of the two higher fidelity models, which consistently estimate approximately 1 Gt under the upper-bound overburden pressure gradient considered for the Sacramento Basin. In all model scenarios, overpressure develops away from injection wells, particularly in higher-elevation regions, highlighting the importance of basin-scale modelling for identifying risks beyond the immediate well vicinity.

Basin-scale↗

Determination of Shallow Crust Properties from Physical Predictors [Slides]

LANL pioneers advanced multivariate inversion techniques for Earth seismic structure modeling at continental (USA) and local scale (Parkfield). Successful application for reducing surface wave magnitude Ms threshold and detecting smaller events (Tarim basin case study). New models hint to new geodynamical and tectonics interpretation (EARS case study). Potential use for surrogate measurements in areas without access to seismic signatures (application to local scale).

54 ENVIRONMENTAL SCIENCES↗

Modeling the Effects of Artificial Drainage on Agriculture-dominated Watersheds using a Fully Distributed Integrated Hydrology Model: Datasets, code, models files

This model-data archive is for a modeling study focused on the effects of artificial drainage on agriculture-dominated watersheds using a fully distributed integrated hydrology model. This study introduced a novel strategy to represent the effects of tile drains and agricultural ditches in a fully distributed basin-scale integrated hydrology model. The numerical experiments in this study reveal the effects of surface and subsurface drainage on various hydrological states and fluxes. The details about the datasets can be found in the "readme" document attached with the dataset.

Agricultural Watersheds↗

The modification stage of basin formation - Conditions of ring formation

The impact cratering process may be divided into three stages, including a short high pressure phase of initial contact between the projectile and target, a longer cratering flow phase during which material is both ejected and displaced to form the transient crater, and a still longer modification phase during which the transient crater is modified into the final observed form. The present investigation is concerned with the nature of the modification stage as constrained by the initial boundary condition which comprises the transient crater and its surroundings, and the final boundary condition which comprises the structure of observed craters and basins. A brief description is provided of an interpretation of basin ring structure based on geophysical, geological, and dimensional analysis of craters and basins on the moon, Mercury, and Mars, on structural observations of terrestrial impact and explosion craters, and on an empirical crater restoration model. A scenario of basin formation is constructed, and relations of this scenario to observed morphology are discussed.

Croft, S. K.↗

Megafans-Some New Perspectives from a Global Study

A global study of megafans (greater than 100 km long) has revealed their widespread existence on all continents, with almost 200 documented, 93 in Africa where research is most thorough. The largest measures 705 km. Megafans are a major subset of "DFS" (distributive fluvial systems, a category that includes all fan-like features greater than 30 km long). 1. Many researchers now recognize megafans as different from floodplains, small coarse-grained alluvial fans, and deltas. Although smaller architectural elements in megafans are the same as those encountered in floodplains (channel, overbank, etc.), larger architectures differ because of the unconfined setting of megafans, versus the valley-confined setting of floodplains. 2. A length continuum is now documented between steep alluvial fans 10-20 km in length, and fluvial fans 30-50 km long. This implies a continuum of process from end-member alluvial fan processes (e.g. high-energy flows that emplace gravels, debris-flow units) to the relatively fine-grained channel and overbank deposits common to purely fluvial fans. Combinations of these different processes will then occur in many mid-sized fans. 3. The global distribution suggests a prima facie relationship with tectonic environment rather than climatic zones, with local controls being the slope of the formative river and the existence of a basin subsiding below the long profile of the river. But the global population has revealed that most megafans are relict. So it is possible that further research will show relationships to prior climatic regimes. 4. Megafans can have regional importance: e.g., along the east flank of the central Andes, nested megafans total approximately 750,000 km2-and 1.2m km2 if all megafans in S. America are counted. Modern megafan landscapes thus have basinal importance, orders of magnitude greater than alluvial fan bajadas. 5. Because so many aggrading basins are dominated today by DFS, it is claimed that DFS ought to be significant in the subsurface; and that existing fluvial models therefore may not apply to the majority of fluvial sedimentary units. Arguments have been raised against this view, but as modern megafan systems become better known they are rapidly being applied as a model in many fluvial basins. A small literature has arisen with apparent examples from every part of the world.

Wilkinson, M. Justin↗

Pantropical Indo-Atlantic temperature gradient modulates multi-decadal AMOC variability in models and observations

Abstract Interconnections between ocean basins are recognized as an important driver of climate variability. Recent modeling evidence suggests that the North Atlantic climate can respond to persistent warming of the tropical Indian Ocean sea surface temperature (SST) relative to the rest of the tropics (rTIO). Here, we use observational data to demonstrate that multi-decadal changes in pantropical ocean temperature gradients lead to variations of an SST-based proxy of the Atlantic Meridional Overturning Circulation (AMOC). The largest contribution to this temperature gradient-AMOCconnection comes from gradients between the Indian and Atlantic Oceans. TherTIOindex yields the strongest connection of this tropical temperature gradient to theAMOC. Focusing on the internally generated signal in three observational products reveals that an SST-basedAMOCproxy index has closely followed low-frequency changes ofrTIOtemperature with about 26-year lag since 1870. Analyzing the pre-industrial control simulations of 44 CMIP6 climate models shows that theAMOCproxy index lags simulated mid-latitudeAMOCvariations by 4 ± 4 years. These model simulations reveal the mechanism connectingAMOCvariations to pantropical ocean temperature gradients at a 27 ± 2 years lag, matching the observed time lag in 28 out of the 44 analyzed models. rTIO temperature changes affect the North Atlantic climate through atmospheric planetary waves, impacting temperature and salinity in the subpolar North Atlantic, which modifies deep convection and ultimately the AMOC. Through this mechanism, observed internalrTIOvariations can serve as a multi-decadal precursor ofAMOCchanges with important implications forAMOCdynamics and predictability.

Meteorology & Atmospheric Sciences↗

Compression of oceanic lithosphere - An analysis of intraplate deformation in the Central Indian Basin

The development of intraplate structure in the Central Indian Basin is examined using models in which deformation is due to flexural buckling and the hydrodynamic growth of instabilities. Comparison of the models reveal that in a strong viscous lithosphere, deformation of the layer occurs by flexural folding at a wavelength which corresponds to the flexural buckling theory; in a lithosphere of intermediate strength, the layer deforms by folding characterized by thickening which localizes beneath topographic heights; and in a relatively weak lithosphere, the layer incurs an even greater amount of localized thickening and deforms in the symmetric or pinch-and-swell mode by inverse boudinage. It is noted that the models in which the layer folds either flexurally or with periodic thickening correspond with observed depth distribution of seismicity in the Central Indian Basin and with experimental rock rheological data.

Zuber, Maria T.↗

Developing Methane Emissions Inventories by Fusing Airborne, Satellite, and Modeled Assessments: Comprehensive Surveys of the Anadarko and Haynesville Basin

The objective of this program is to create basin-specific methane emissions inventories of the Haynesville and Anadarko (Woodford shale) basins, based on truly basin-scale data collection and analysis. Emissions of natural gas, predominantly composed of methane, present a significant financial loss to the oil and gas sector and reduce the competitiveness of American energy sources. While large emissions, above 10 kg/h, occur at only a small percentage of oil and gas sites, they account for a disproportionately large volume of emissions. To address this, we executed a concurrent campaign combining cutting-edge bottom-up modeling with top-down aerial and satellite surveys to capture the heavy-tailed distribution of these large emissions to accurately characterize total emissions throughout both regions.

03 NATURAL GAS↗

Comparison of Machine Learning-Based Predictive Models of the Nutrient Loads Delivered from the Mississippi/Atchafalaya River Basin to the Gulf of Mexico

Predicting nutrient loads is essential to understanding and managing one of the environmental issues faced by the northern Gulf of Mexico hypoxic zone, which poses a severe threat to the Gulf’s healthy ecosystem and economy. The development of hypoxia in the Gulf of Mexico is strongly associated with the eutrophication process initiated by excessive nutrient loads. Due to the complexities in the excessive nutrient loads to the Gulf of Mexico, it is challenging to understand and predict the underlying temporal variation of nutrient loads. The study was aimed at identifying an optimal predictive machine learning model to capture and predict nonlinear behavior of the nutrient loads delivered from the Mississippi/Atchafalaya River Basin (MARB) to the Gulf of Mexico. For this purpose, monthly nutrient loads (N and P) in tons were collected from US Geological Survey (USGS) monitoring station 07373420 from 1980 to 2020. Machine learning models—including autoregressive integrated moving average (ARIMA), gaussian process regression (GPR), single-layer multilayer perceptron (MLP), and a long short-term memory (LSTM) with the single hidden layer—were developed to predict the monthly nutrient loads, and model performances were evaluated by standard assessment metrics—Root Mean Square Error (RMSE) and Correlation Coefficient (R). The residuals of predictive models were examined by the Durbin–Watson statistic. The results showed that MLP and LSTM persistently achieved better accuracy in predicting monthly TN and TP loads compared to GPR and ARIMA. In addition, GPR models achieved slightly better test RMSE score than ARIMA models while their correlation coefficients are much lower than ARIMA models. Moreover, MLP performed slightly better than LSTM in predicting monthly TP loads while LSTM slightly outperformed for TN loads. Furthermore, it was found that the optimizer and number of inputs didn’t show effects on the LSTM performance while they exhibited impacts on MLP outcomes. This study explores the capability of machine learning models to accurately predict nonlinearly fluctuating nutrient loads delivered to the Gulf of Mexico. Further efforts focus on improving the accuracy of forecasting using hybrid models which combine several machine learning models with superior predictive performance for nutrient fluxes throughout the MARB.

54 ENVIRONMENTAL SCIENCES↗

Floods due to Atmospheric Rivers along the U.S. West Coast: The Role of Antecedent Soil Moisture in a Warming Climate

Precipitation extremes are projected to become more frequent along the U.S. West Coast due to increased atmospheric river (AR) activity, but the frequency of less intense precipitation events may decrease. Antecedent soil moisture (ASM) conditions can have a large impact on flood responses, especially if prestorm precipitation decreases. Taken together with increased antecedent evaporative demand due to warming, this would result in reduced soil moisture at the onset of extreme precipitation events. In this work, we examine the impact of ASM on AR-related floods in a warming climate in three basins that form a transect along the U.S. Pacific Coast: the Chehalis River basin in Washington, the Russian River basin in Northern California, and the Santa Margarita River basin in Southern California. Additionally, we ran the Distributed Hydrology Soil Vegetation Model (DHSVM) over the three river basins using forcings downscaled from 10 global climate models (GCMs). We examined the dynamic role of ASM by comparing the changes in the largest 50, 100, and 150 extreme events in two periods, 1951–2000 and 2050–99. In the Chehalis basin, the projected fraction of AR-related extreme discharge events slightly decreases. In the Russian basin, this fraction increases, however, and more substantially so in the Santa Margarita basin. This is due to increases in AR-related extreme precipitation events, as well as the fact that the relationship of extreme precipitation to extreme discharge is strengthened by projected increases in year-to-year volatility of annual precipitation in California, which increases the likelihood of concurrent occurrence of large storms and wet ASM conditions.

54 ENVIRONMENTAL SCIENCES↗

Assessment of circulation and inter-basin transport in the Salish Sea including Johnstone Strait and Discovery Islands pathways

The Salish Sea consisting of Puget Sound and Georgia Basin in U.S and Canadian waters has been the subject of several independent data collection and modeling studies. However, these interconnected basins and their hydrodynamic interactions have not received attention as a contiguous unit. The Strait of Juan de Fuca is the primary pathway through which Pacific Ocean water enters the Salish Sea but the role played by Johnstone Strait and the complex channels northeast of Vancouver Island, connecting the Salish Sea and the Pacific Ocean, on overall Salish Sea circulation has not been characterized. In this paper we present a modeling-based assessment of the two-layer circulation and transport through the multiple interconnected sub-basins within the Salish Sea including the effect of exchange via Johnstone Strait and Discovery Islands. The Salish Sea Model previously developed using the finite volume community ocean model (FVCOM) was expanded over the continental shelf for this assessment encircling Vancouver Island, including Discovery Islands, Johnstone Strait, Broughton Archipelago and the associated waterways. A computational technique was developed to allow summation of volume fluxes across arbitrary transects through unstructured finite volume cells. Tidally averaged volume fluxes were computed at multiple transects. The results were used to validate the classic model of Circulation in Embracing Sills for Puget Sound and to provide quantitative estimates of the lateral distribution of tidally averaged transport through the system. Sensitivity tests with and without exchanges through Johnstone Strait demonstrate that it is a pathway for Georgia Basin runoff and Fraser River water to exit the Salish Sea and for Pacific Ocean inflow. However the relative impact of this exchange on circulation and flushing in Puget Sound Basin is small.

estuarine circulation↗

Reservoir Thermal Energy Storage Benchmarking (Rev. 3)

A benchmarking analysis of RTES research funded by GTO through the Beyond Batteries projects was conducted against the ESGC to see where they fit within the identified ESGC Use Cases. The projects were found to advance knowledge in multiple ESGC use cases, either directly or in some cases, indirectly as enabling technologies. This analysis is helpful to understand where RTES and associated research fits into the larger discussion around energy storage technologies. Also, a retrospective analysis of the Beyond Batteries projects was conducted to evaluate what the projects learned and how the results can be applied to advance the value of RTES. Major results of each of the studies are summarized in Table 2. Additionally, a comparative metrics analysis for RTES was completed to understand where RTES lies within the energy storage industry. Metrics for evaluation of RTES and its comparison to other storage technologies were selected and ranges of their values compiled. The selected metrics – LCOE (levelized cost of energy), capital costs, roundtrip efficiency, energy storage capacity, and storage time – were chosen based on data availability and have a particularly strong influence on the potential deployment of a storage technology. Charts which compare the metrics are presented in section 4.3 and show ranges for each of the 10 selected technologies. However, due to a lack of domestic operational facilities, values for RTES and for portions of the remaining technologies are based on theoretical modeling and studies of best-case scenarios. LCOE estimates for RTES fall within the lower reaches of Figure 15, but nevertheless amount to 2 – 5 times the ESGC Roadmap goal for LCOE, for example in the Facilitating and Evolving Grid Use Case. Capital costs for RTES sit on the higher end (Figure 16) but are expected to decrease as new projects are developed and the technology is refined. The theoretical roundtrip efficiency reported for RTES varies from mid to high percentages (Figure 17) with efficiencies upwards of 93% in modeled scenarios in the Portland Basin (Bershaw et al.,2020). RTES is also expected to have the largest energy storage capacities and longest storage times, likely matched only by lower efficiency hydrogen storage. To better assess the role that RTES could play in energy storage we examined it’s potential in the U.S. The potential depends on many factors. Recently, many researchers have started looking at deep sedimentary basins, depleted oil and gas fields, and basalt formations as potential targets for RTES development. The United States Geological Survey (USGS) has analyzed various cities and shown substantial RTES potential in the cooling sector (Pepin et al., 2021). By modeling RTES in low-quality groundwater (e.g., brackish), it is shown to be favorable across the U.S. with particular suitability in the Illinois Basin, Coastal Plains, and Basin and Range regions. Seasonal RTES operations have also been modeled in the Portland Basin by those at the USGS and Portland State University to simulate an RTES system supplying heating loads needed for the Oregon Health and Science University. Simulations suggest that high conductive heat loss in the initial years exists but tends to decrease with increasing time and development of the resource due to self-insulating nature of the basalts (Burns et al., 2020). Other national laboratory efforts are taking a close look at many of the technical issues involved with RTES (McLing et al., 2019, McLing et al., 2022). These include difficulties in understanding geochemical, hydrogeological, mechanical, and microbiological changes at such elevated temperatures and operational scenarios. Major gaps in research are identified and suggested for future work. With this increased focus to understand how to make RTES successful in the U.S., this technology could be a potential solution to many of the nation’s energy storage problems. For the energy independence of this country, the DOE should prioritize de-risking this technology by making future investments in pilot-scale demonstrations to attract potential investors.

15 GEOTHERMAL ENERGY↗

Insights into storm direction effect on flood response

In this study, we investigate the directional influence of storm movement on catchment flood peak response using the synthetic circular basin. Due to the complexity in defining storm movements that require meteorological modeling, we adopt a novel approach of combining the basin rotation method (BRM) with a circular basin construct. A systematic basin rotation approach provides a proxy for the storm direction relative to the river network topology to study its hydrologic response. Using Stage-IV rainfall with 4-km by 1-hour resolution for 16 years period from 2003 to 2018 in Iowa, U.S., we systematically analyze consequent flood peak response due to storm directions (basin rotations) using a distributed hydrologic model called the Hillslope-Link Model (HLM). HLM has demonstrated the ability to reproduce observed streamflow for real river basins studied in Iowa. The hydrologic simulation results show that the BRM can quantify the relationship between the rainstorm direction and catchment hydrologic response. Also, the maximum flood peak response is strongly dependent on the interaction of the peak of the rainstorm with the peak of the width function. The results indicate significant differences in runoff volume and runoff peak and their interannual variability due to rainstorm direction. The study has important implications for flood frequency predictions and developing flood resilience and watershed management strategies, particularly due to changing storm tracks under a warming climate.

42 ENGINEERING↗

OC5 Project Phase II: Validation of Global Loads of the DeepCwind Floating Semisubmersible Wind Turbine

This paper summarizes the findings from Phase II of the Offshore Code Comparison, Collaboration, Continued, with Correlation project. The project is run under the International Energy Agency Wind Research Task 30, and is focused on validating the tools used for modeling offshore wind systems through the comparison of simulated responses of select system designs to physical test data. Validation activities such as these lead to improvement of offshore wind modeling tools, which will enable the development of more innovative and cost-effective offshore wind designs. For Phase II of the project, numerical models of the DeepCwind floating semisubmersible wind system were validated using measurement data from a 1/50th-scale validation campaign performed at the Maritime Research Institute Netherlands offshore wave basin. Validation of the models was performed by comparing the calculated ultimate and fatigue loads for eight different wave-only and combined wind/wave test cases against the measured data, after calibration was performed using free-decay, wind-only, and wave-only tests. The results show a decent estimation of both the ultimate and fatigue loads for the simulated results, but with a fairly consistent underestimation in the tower and upwind mooring line loads that can be attributed to an underestimation of wave-excitation forces outside the linear wave-excitation region, and the presence of broadband frequency excitation in the experimental measurements from wind. Participant results showed varied agreement with the experimental measurements based on the modeling approach used. Modeling attributes that enabled better agreement included: the use of a dynamic mooring model; wave stretching, or some other hydrodynamic modeling approach that excites frequencies outside the linear wave region; nonlinear wave kinematics models; and unsteady aerodynamics models. Also, it was observed that a Morison-only hydrodynamic modeling approach could create excessive pitch excitation and resulting tower loads in some frequency bands.

17 WIND ENERGY↗

Evidence for a Stratified Upper Mantle Preserved within the South Pole - Aitken Basin

The evolution and compositional structure of the lunar mantle has been extensively modeled but insufficiently constrained by observations. Here, we identify and characterize mantle materials exposed by the Moon's largest impact basin to better understand the composition, stratigraphy, and evolution of the upper mantle. The vast South Pole-Aitken Basin (SPA) exhibits a broad, crescent-shaped thorium and potassium distribution. These incompatible elements are predicted to be concentrated in the dregs of the lunar magma ocean during end-stage crystallization. Through consideration of basin formation models convolved with subsequent geologic evolution, we demonstrate that the distribution and implied stratigraphy of Th- and K-bearing materials across SPA are consistent with an upper mantle ejecta origin. The most pristine exposures of these materials are confined to northwest SPA and also exhibit elevated Ti and Fe (relative to the farside highlands) in association with a gabbronoritic mineralogy. This is consistent with latestage magma ocean assemblages predicted by petrologic models. In contrast, SPA impact melt derived from greater depths is associated with a low-Ca pyroxene-dominated assemblage. Together, these compositional patterns are evidence for a stratified ancient upper mantle. Importantly, the incompatible-element-enriched, ilmenite-bearing ferroan gabbronoritic cumulates evidently had not participated in gravitational overturn at the time of SPA formation. Contrary to recent hypotheses invoking nearside sequestration of incompatible elements to explain hemispherical differences in crustal building and volcanic resurfacing, it follows that incompatible elements were globally distributed in the magma ocean at the time of SPA formation.

Lunar dichotomy↗

Model Inputs, Outputs, and Scripts associated with: “Spatial microbial respiration variations in the hyporheic zones within the Columbia River Basin”

This data package is associated with the publication “Spatial microbial respiration variations in the hyporheic zones within the Columbia River Basin” published in the Journal of Geophysical Research: Biogeosciences (Son et al. 2022) available at doi: 10.1029/2021JG006654. This data package includes the key model inputs/outputs of the river corridor model for the Columbia River Basin (CRB) and the model source codes, which were used in the manuscript. The model is a carbon-nitrogen-coupled river corridor model (RCM), and the model is used to quantify hyporheic zone (HZ) aerobic and anaerobic respiration at the NHDPLUS stream reach scales. The RCM used in this study combines empirical substrate models derived from observations and three microbially driven reactions to compute respiration of the HZ for each National Hydrography Dataset (NHD) reach within the CRB. The reactions in HZs of each NHD reach include anaerobic respiration and two-step anaerobic respiration via denitrification. Our HZ respiration estimates are limited to the lotic (or flowing) stream/river systems, and do not account for the respiration process in water column. Note that the RCM only simulates the HZ’s contribution to the dissolved carbon dioxide (CO2) concentrations in the streams, and the CO2 emissions to the atmosphere are not modelled. The model computes at hourly timesteps because of the fast reaction rates. The key input data of the model are exchange flux, residence time, and stream solute (dissolved organic carbon (DOC), dissolved oxygen (DO), and nitrate concentrations). These inputs are constant over time and represent long-term averaged values.This modeling framework successfully quantified HZ respiration components over multiple scales. It revealed key mechanisms driving the spatial variation of HZ aerobic and anaerobic respiration in reaches with varying hydrologic and substrate conditions. Thus, this modeling study offers a testing hypothesis in different river system (e.g., climate and biomes) for the HZ respiration processes, and can be used as a sampling design tool for large-scale HZ experimental studies.This dataset contains five folders: (1) model_inputs, (2) model_outputs, (3) Rscripts, (4) figures, and (5) model_codes. It also contains a readme, file level metadata (FLMD), and data dictionary (dd). Please see the FLMD for a list of all the files contained in this data package and descriptions for each. The model_inputs folder contains the model inputs used to drive the model simulations. The model_outputs folder contains key model output files from the river corridor model. The Rscripts folder contains the Rscripts for pre- and post- processing model results. The figures folder contains the raw figures associated with the manuscript. The model_codes folder includes key model source codes/input files. All files are .jpg, .jpeg, .out, .e, .od, .dat, .sub, .F90, .0, .R, .sbx, .cpg, .sbn, .shx, .shp, .dbf, .prj, .tfw, .tif, .xml, .pdf, or .csv.

54 ENVIRONMENTAL SCIENCES↗

Thermal cooling of the oceanic lithosphere from geoid height data

Another type of geophysical observation has proved to be very useful in the study of thermal cooling of the oceanic lithosphere. It is the geoid height derivative with respect to plate age, a quantity computed from the short wavelength geoid step across fracture zones measured by altimeter satellites. Two categories of simples models are proposed to describe cooling and contraction of the oceanic lithosphere with age. Both plate model and half space model, give almost similar results up to ages of 50 to 70 ma, but predict quite distinct behavior of seafloor depth, heat flow and other parameters in old basins. Tests of thermal models are based on heat flow and topography data. However, heat flow is not very sensitive to the form of the thermal model. Large areas of the ocean floor are particularly shallow, and as a result topography data may not be very appropriate to discriminate between plate and half space models, and no clear concensus on a preferred model yet exists.

Cazenae, A.↗