Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “ocean model”

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 1,297 records · Page 72

Vertical ocean loading amplitudes from VLBI measurements

Vertical ocean loading amplitudes are determined by analysis of IRIS geodetic Very Long Baseline Interferometry (VLBI) data. The 4 nearly diurnal (K(sub 1), P(sub 1), O(sub 1), Q(sub 1)) and 4 nearly semidiurnal (K(sub 2), S(sub 2), M(sub 2), and N(sub 2)) component amplitudes can be inferred from the data with accuracies of 1-2 mm. Uncertainties of total displacements can approach 1 cm. Empirically determined total displacements are considerably larger than the values calculated from two geophysical models. The Scherneck model is found to give a better representation of VLBI delay data than the model of Pagiatakis by about 3 mm in residuals at 6 sites. Empirical estimation of the ocean loading amplitudes reduces Chi squared by 3067 for the 96 additional degrees of freedom in a fit to 273,000 IRIS VLBI observations, and reduces the RMS residuals by 3 mm relative to a fit using the fixed Scherneck model. Vertical ocean loading amplitudes can thus be inferred from VLBI data at a level which improves the overall model, but detailed assessment of individual tidal components is presently obscured by incomplete modeling at the tidal frequencies.

Sovers, O. J.↗

Evolution of the 1997-98 El Nino and La Nina: the Role of Altimetry in ENSO Studies and Prediction

El Nino and the Southern Oscillation (ENSO) are known to be a quasi-periodic coupled ocean-atmosphere phenomenon occurring roughly every 3-7 years. Its largest oceanic manifestations are found in the equatorial Pacific, but it is of global climatic significance. During the flight of the TOPEX]Poseidon altimeter, the equatorial Pacific was in the midst of an unusual period of prolonged ENSO warming in the mid-1990's. In 1997, the tropical Pacific witnessed the development of a major El Nino event, rivaling in strength the 1982-83 El Nino. The intensity of this El Nino, and its rapid intensification early in the calendar year, caught the scientific community by surprise. None of the forecast models, statistical or dynamical, captured these features of the 1997 El Nino, although those models that assimilated in situ ocean observations did better than most in predicting that the event would occur in 1997. The onset of the El Nino coincided with the occurrence of several westerly wind events in the western Pacific rather than to delayed oscillator physics. However, western Pacific sea levels were anomalously elevated a year and a half prior to the onset, which may have helped precondition the system to a particularly strong ENSO episode. The coming of La Nina was foreshadowed by large-scale drops in sea level in the western tropical Pacific and cyclonic wind stress curl anomalies far beyond the +\- 8 degrees latitude bounds of the TAO mooring array. In past El Nino events, the evolution of the coupled system could only be studied from the perspective of in situ observations. The 1997-1998 El Nino event will be the first time a major El Nino event will have been observed from start to finish via remotely-sensed measurements of sea surface topography, sea surface temperature, sea surface winds, ocean color, and precipitation. These observations have considerable potential for advancing coupled climate forecasts on time scales from seasons to decades. This presentation will address questions along the lines of: (1) What are the scientific problems in ENSO that can benefit most from altimetry data? (2) Are there any precursors in the global ocean-atmosphere system for the 1997-98 extreme event? (3) Is there a unique role of altimetry in detecting potential precursors?(4) Has Topex/Poseidon data been adequately used to improve ENSO forecasts? Is there any firm evidence? (6) What is the status of the transition of El Nino into a La Nina?

Busalacchi, Antonio J.↗

Using 3-D Numerical Weather Data in Piloted Simulations

This report describes the process of acquiring and using 3-D numerical model weather data sets in NASA Langley's Research Flight Deck (RFD). A set of software tools implement the process and can be used for other purposes as well. Given time and location information of a weather phenomenon of interest, the user can download associated numerical weather model data. These data are created by the National Oceanic and Atmospheric Administration (NOAA) High Resolution Rapid Refresh (HRRR) model, and are then processed using a set of Mathworks' Matlab(TradeMark) scripts to create the usable 3-D weather data sets. Each data set includes radar re ectivity, water vapor, component winds, temperature, supercooled liquid water, turbulence, pressure, altitude, land elevation, relative humidity, and water phases. An open-source data processing program, wgrib2, is available from NOAA online, and is used along with Matlab scripts. These scripts are described with sucient detail to make future modi cations. These software tools have been used to generate 3-D weather data for various RFD experiments.

Daniels, Taumi S.↗

Eastern Equatorial Pacific Warming Delayed By Aerosols and Thermostat Response to Co2 Increase

Understanding the tropical Pacific response to global warming remains challenging. Here, we use a range of CMIP6 greenhouse warming experiments to assess the recent and future evolution of the equatorial Pacific east-west temperature gradient and corresponding Walker circulation. In abrupt CO2-increase scenarios many models generate an initial strengthening of this gradient resembling an ocean thermostat (OT), followed by a small weakening; other models generate an immediate weakening that becomes progressively stronger establishing a pronounced eastern equatorial Pacific (EP) warming pattern. The initial response in these experiments is a strong predictor for the future EP pattern simulated in both abrupt and realistic warming scenarios, but not in historical simulations showing no multi-model trend. The likely explanation is that the recent CO2-driven changes in the tropical Pacific are masked by aerosol effects and a potential OT-related delay, while the EP warming pattern will emerge as greenhouse gases overcome aerosol forcing.

Climate and Earth system modelling↗

Direct numerical simulations of stably-stratified sheared turbulence: Implications for oceanic mixing

Direct numerical simulations of the time evolution of homogeneous stably stratified turbulent shear flows have been performed for several Richardson numbers Ri and Reynolds numbers R(sub lambda) in earlier works. The results show excellent agreement with length scale models developed from laboratory experiments to characterize oceanic turbulence. When the Richardson number Ri is less than the stationary value Ri(sub s), the turbulence intensity grows at all scales, and the growth rate appears to be a function of Ri. The size of the vertical density inversions also increases. On the other hand, when Ri is greater than or equal to Ri(sub s) the largest turbulent eddies become vertically constrained by buoyancy when the Ellison (turbulence) scale L(sub E) and the Ozmidov (buoyancy) scale L(sub O) are equal. At this point, the mixing efficiency is maximal and corresponds to a flux Richardson number R(sub f) = 0.20. The vertical mass flux becomes counter-gradient when epsilon = 19(nu)N(exp 2) and vertical density overturns are suppressed in less than half a Brunt-Vaisala period. The results of the simulations were also recast in terms of the Hydrodynamic Phase Diagram introduced in fossil turbulence models. The so-called point of fossilization occurs when epsilon = 4DCN(exp 2); Gibson proposed 13DCN(exp 2). This value is in agreement with indirect laboratory observations and field observations. Finally, the validity of the steady-state models to estimate vertical eddy diffusivities in the oceanic thermocline is discussed.

Itsweire, E. C.↗

The effects of atmospheric cloud radiative forcing on climate

In order to isolate the effects of atmospheric cloud radiative forcing (ACRF) on climate, the general circulation of an ocean-covered earth called 'Seaworld' was simulated using the Colorado State University GCM. Most current climate models, however, do not include an interactive ocean. The key simplifications in 'Seaworld' are the fixed boundary temperature with no land points, the lack of mountains and the zonal uniformity of the boundary conditions. Two 90-day 'perpetual July' simulations were performed and analyzed the last sixty days of each. The first run included all the model's physical parameterizations, while the second omitted the effects of clouds in both the solar and terrestrial radiation parameterizations. Fixed and identical boundary temperatures were set for the two runs, and resulted in differences revealing the direct and indirect effects of the ACRF on the large-scale circulation and the parameterized hydrologic processes.

Randall, David A.↗

Chemistry of Venus’ Recent Basalts as Clues to its Ancient Past

Perhaps the most important question about Venus is whether it ever had a hydrosphere: liquid water, oceans, and thus an environment suitable for life as we know it. The D/H of Venus’ atmosphere suggests extensive water loss, and climate models may be consistent with global oceans. Rocks dating from that epoch may be preserved in Venus’ tesserae (and especially Ishtar Terra) – inferred to of felsic or silicic rock, which would suggest abundant water. However, the exact elemental compositions of tesserae (silicic or not) will be difficult to retrieve and most of Venus’ surface is basalt flows and volcanic constructs, inferred to be much younger than the oceanic epoch. Remotely sensed data are (and will be) ambiguous about the specific rock types that make up tesserae, and lander spacecraft are almost certain (in the near term) to avoid the rough and precipitous topography of the tesserae and touch down instead on flat, safe basaltic plains or volcanic rises. ‘Recent’ basalts at Venus’ surface could preserve chemical tracers of an ancient aqueous past, if their source regions (material that was re-melted) had been affected by water. This scenario occurs on Earth in Island Arc Basalts (IAB), where their compositions are thought to reflect aqueous alteration of parental midocean ridge basalts (MORB) and incorporation of oceanic sediments. Ancient Venus might have supported plate tectonic and thus have had IAB equivalents; absent plate tectonics, basalt affected by aqueous alteration could have been cycled into its mantle by burial under thick sections of later basalt. Our purpose is to suggest specific chemical clues in current basalts that would permit recognition of those with a history of aqueous interactions from those that did not. We compare Earth’s Ocean Island Basalts (OIB) that involved little water with IAB that show chemical effects ascribable to water. It is not clear that Venus’ basalts can be mapped into terrestrial tectonic settings; however, it seems reasonable that aqueous geochemical processes could have (had) similar effects on both planets.

A H Treiman↗

Provenance of the K/T boundary layers

An array of chemical, physical and isotopic evidence indicates that an impact into oceanic crust terminated the Cretaceous Period. Approximately 1500 cu km of debris, dispersed by the impact fireball, fell out globally in marine and nonmarine environments producing a 2 to 4 mm thick layer (fireball layer). In North American locales, the fireball layer overlies a 15 to 25 mm thick layer of similar but distinct composition. This 15 to 25 mm layer (ejecta layer) may represent approximately 1000 cu km of lower energy ejecta from a nearby impact site. Isotopic and chemical evidence supports a mantle provenance for the bulk of the layers. The extraordinary REE pattern of the boundary clays was modelled as a mixture of oceanic crust, mantle, and approximately 10 percent continental material. The results are presented. If the siderophiles of the ejecta layer were derived solely from the mantle, a test may be available to see if the siderophile element anomaly of the fireball layer had an extraterrestrial origin. Radiogenic Os-187 is depleted in the mantle relative to an undifferentiated chondritic source. Os-187/Os-186 ratios of 1.049 and 1.108 were calculated for the ejecta and fireball layers, respectively.

Hildebrand, A. R.↗

Using a laser aureole to invert lidar return

An aureole generated by a laser beam was studied. The strength of the signal redirected towards a sensor high above the surface by a combination of one scattering event in the marine boundary layer (mbl) and one single reflection event from the ocean surface was estimated. A model of mbl aerosol size distributions was used to estimate Mie scattering for a wide range of meteorolocial conditions. The sea surface reflection was determined from a Gaussian model of the wave slopes. These laser aureoles which were estimated over the wide range of conditions and were normalized by the reflected laser light were found to be highly correlated with the optical depth of the boundary layer. By estimating optical depth from the aureole, the Bernoulli-Riccati inversion of lidar return could be constrained and the inversion accuracy improved. A Monte Carlo program was developed to study the laser aureole generated by up to 8 orders of reflection and scattering. The aureole was generated by a narrow, 10 nsec laser pulse at 1.06 microns and measured by a receiver 10 km above the ocean surface. The original theoretical computation compared well with the Monte Carlo method. When multiple scattereffects were included, the normalized aureole was still highly correlated with the mbl optical depth over the range of conditions.

Hooper, W. P.↗

Theory for synthetic aperture radar imaging of the ocean surface - With application to the Tower Ocean Wave and Radar Dependence experiment on focus, resolution, and wave height spectra

A one-dimensional model for simulating azimuthal SAR imaging of the ocean surface is developed which can admit both the 'distributed surface' and 'velocity bunching' approaches. Computer simulations demonstrate that the time-dependent modulation patterns due to the radar cross section variation and the velocity bunching effects provide optimum focusing around half the phase velocity of the long wave. The results indicate that in the Tower Ocean Wave and Radar Dependence experiment, SAR imaging at L band is approximately linear.

Kasilingam, Dayalan P.↗

Constraining Remote Oxidation Capacity with ATom Observations

The global oxidation capacity, defined as the tropospheric mean concentration of the hydroxyl radical (OH), controls the lifetime of reactive trace gases in the atmosphere such as methane and carbon monoxide (CO). Models tend to underestimate the methane lifetime and CO concentrations throughout the troposphere, which is consistent with excessive OH. Approximately half the oxidation of methane and non-methane volatile organic compounds (VOCs) is thought to occur over the oceans where oxidant chemistry has received little validation due to a lack of observational constraints. We use observations from the first two deployments of the NASA ATom aircraft campaign during July–August 2016 and January–February 2017 to evaluate the oxidation capacity over the remote oceans and its representation in the GEOS-Chem chemical transport model. The model successfully simulates the magnitude and vertical profile of remote OH within the measurement uncertainties. Comparisons against the drivers of OH production (water vapor, ozone, and NOy concentrations, ozone photolysis frequencies) also show minimal bias with the exception of wintertime NOy, for which a model overestimate may indicate insufficient wet scavenging and/or missing loss on seasalt aerosol but large uncertainties remain that require further studies of NOy partitioning and removal in the troposphere. During the ATom-1 deployment, OH reactivity (OHR) below 3 km is significantly enhanced, and this is not captured by the sum of its measured components (cOHRobs) or by the model (cOHRmod). This enhancement could suggest missing reactive VOCs but cannot be explained by new estimates of ocean VOC sources and additional modeled reactivity in this region would be difficult to reconcile with the full suite of ATom measurement constraints. The model generally reproduces the magnitude and seasonality of cOHRobs but underestimates the contribution of oxygenated VOC, mainly acetaldehyde, which is severely underestimated throughout the troposphere despite its calculated lifetime of less than a day. Missing model acetaldehyde in previous studies was attributed to measurement uncertainties that have been largely resolved. Observations of peroxyacetic acid (PAA) provide new support for remote levels of acetaldehyde. The underestimate in modeled acetaldehyde and PAA is present throughout the year in both hemispheres and peaks during Northern Hemisphere summer. The addition of ocean VOC sources in the model increases annual surface cOHRmod by 10 % and improves model-measurement agreement for acetaldehyde particularly in winter but cannot resolve the model summertime bias. Doing so would require a 100 Tg yr−1 source of a long-lived unknown precursor throughout the year with significant additional emissions in the Northern Hemisphere summer. Improving the model bias for remote acetaldehyde and PAA is unlikely to fully resolve previously reported model global biases in OH and methane lifetime, suggesting that future work should examine the sources and sinks of OH over land.

Travis, Katherine R.↗

Coupled Data Assimilation for Integrated Earth System Analysis and Prediction: Goals, Challenges, and Recommendations

The purpose of this report is to identify fundamental issues for coupled data assimilation (CDA), such as gaps in science and limitations in forecasting systems, in order to provide guidance to the World Meteorological Organization (WMO) on how to facilitate more rapid progress internationally. Coupled Earth system modeling provides the opportunity to extend skillful atmospheric forecasts beyond the traditional two-week barrier by extracting skill from low-frequency state components such as the land, ocean, and sea ice. More generally, coupled models are needed to support seamless prediction systems that span timescales from weather, subseasonal to seasonal (S2S), multiyear, and decadal. Therefore, initialization methods are needed for coupled Earth system models, either applied to each individual component (called Weakly Coupled Data Assimilation - WCDA) or applied the coupled Earth system model as a whole (called Strongly Coupled Data Assimilation - SCDA). Using CDA, in which model forecasts and potentially the state estimation are performed jointly, each model domain benefits from observations in other domains either directly using error covariance information known at the time of the analysis (SCDA), or indirectly through flux interactions at the model boundaries (WCDA). Because the non-atmospheric domains are generally under-observed compared to the atmosphere, CDA provides a significant advantage over single-domain analyses. Next, we provide a synopsis of goals, challenges, and recommendations to advance CDA: Goals: (a) Extend predictive skill beyond the current capability of NWP (e.g. as demonstrated by improving forecast skill scores), (b) produce physically consistent initial conditions for coupled numerical prediction systems and reanalyses (including consistent fluxes at the domain interfaces), (c) make best use of existing observations by allowing observations from each domain to influence and improve the full earth system analysis, (d) develop a robust observation-based identification and understanding of mechanisms that determine the variability of weather and climate, (e) identify critical weaknesses in coupled models and the earth observing system, (f) generate full-field estimates of unobserved or sparsely observed variables, (g) improve the estimation of the external forcings causing changes to climate, (h) transition successes from idealized CDA experiments to real-world applications. Challenges: (a) Modeling at the interfaces between interacting components of coupled Earth system models may be inadequate for estimating uncertainty or error covariances between domains, (b) current data assimilation methods may be insufficient to simultaneously analyze domains containing multiple spatiotemporal scales of interest, (c) there is no standardization of observation data or their delivery systems across domains, (d) the size and complexity of many large-scale coupled Earth system models makes it is difficult to accurately represent uncertainty due to model parameters and coupling parameters, (e) model errors lead to local biases that can transfer between the different Earth system components and lead to coupled model biases and long-term model drift, (e) information propagation across model components with different spatiotemporal scales is extremely complicated, and must be improved in current coupled modeling frameworks, (h) there is insufficient knowledge on how to represent evolving errors in non-atmospheric model components (e.g. as sea ice, land and ocean) on the timescales of NWP.

CDA↗

Mapping Aerosol Lidar Ratios Over Ocean for CALIPSO Using Modis AOD-Constrained Retrievals and A Global Aerosol Model

The CALIPSO aerosol algorithms currently assign one lidar ratio value for each of the seven tropospheric aerosol types. In a future data products release, the CALIPSO project aims to improve these algorithms by developing regional and seasonal lidar ratio climatologies. In this study, aerosol lidar ratios are inferred through Cloud-Aerosol Lidar with Orthogonal Polarization (CALIOP) backscatter profiles constrained by collocated aerosol optical depth from Aqua Moderate Resolution Imaging Spectroradiometer (MODIS) datasets. The analysis is subsampled for only those profiles that are cloud-free and contain one aerosol type (based on CALIOP feature classification). In addition, the CALIOP profiles are collocated with aerosol volume fractions simulated by the Goddard Chemistry Aerosol Radiation and Transport (GOCART) model. In this talk, the twelve-year (2006-2017) mean spatial distributions of inferred aerosol lidar ratios for CALIOP-classified marine aerosols and the corresponding modeled sea salt volume fractions (SSVF) will be shown. Model-assisted climatological lidar ratio maps on a seasonal scale will also be presented, developed from the empirical relationship found between the modeled SSVF and lidar ratios. A comparison with past studies will be provided as well as results of a sensitivity study regarding the variability of retrieved lidar ratios as a function of horizontal averaging resolution. While the majority of this talk will focus on lidar ratios for CALIOP-classified marine aerosols, preliminary results will be shown for other aerosol types over ocean, such as dust and elevated smoke.

Travis Toth↗

Synthesis and Assimilation Systems - Essential Adjuncts to the Global Ocean Observing System

Ocean assimilation systems synthesize diverse in situ and satellite data streams into four-dimensional state estimates by combining the various observations with the model. Assimilation is particularly important for the ocean where subsurface observations, even today, are sparse and intermittent compared with the scales needed to represent ocean variability and where satellites only sense the surface. Developments in assimilation and in the observing system have advanced our understanding and prediction of ocean variations at mesoscale and climate scales. Use of these systems for assessing the observing system helps identify the strengths of each observation type. Results indicate that the ocean remains under-sampled and that further improvements in the observing system are needed. Prospects for future advances lie in improved models and better estimates of error statistics for both models and observations. Future developments will be increasingly towards consistent analyses across components of the Earth system. However, even today ocean synthesis and assimilation systems are providing products that are useful for many applications and should be considered an integral part of the global ocean observing and information system.

Rienecker, Michele M.↗

Dust Transport, Deposition and Radiative Effects Observed from MODIS

Carlson (1977) used satellite (AVHRR) observation of dust episodes 3 estimate that 90 tg of dust are emitted from Africa (0-30 N) to the Atlantic Ocean between June and August. MODIS systematic measurements of aerosol optical thickness (AOT) and the fraction of the AOT (f) due to the fine mode (see Remer et al abstract), are used to derive the column concentration, flux and deposition of African dust over the Atlantic Ocean. The main data set is for 2001 but the results are consistent with MODIS measurements from 2002. The analysis first determines the properties of maritime baseline aerosol (AOT=0.06, f=0.5); followed by linear scaling of the dust AOT and the anthropogenic AOT, based on MODIS measured values of the fraction "f" being 0.9 for anthropogenic aerosol and 0.5 for dust. NCEP winds are used in the analysis and are evaluated against observed dust movements between the Terra and Aqua passes (see Koren et al. abstract). Monthly values of dust transport and deposition are calculated. Preliminary results show that 280 tg of dust are emitted annually from Africa to the Atlantic Ocean between 20s and 30N, with 40 tg returning to Africa and Europe between 30N and 50N. 85 tg reach the Americas, with 130-150 tg are deposited in the Atlantic Ocean. The results are compared with dust transport models that indicate 110-230 tg of dust being deposited in the Ocean. It is interesting to note that the early estimates of Carlson (1977) and Carlson & Prosper0 (1972) are very close to our estimate from MODIS of 100 tg for the same latitude range and monthly period.

Kaufman, Y. J.↗

Antarctic Ice Shelf Thickness Change from Multimission Lidar Mapping

We calculate rates of ice thickness change and bottom melt for ice shelves in West Antarctica and the Antarctic Peninsula from a combination of elevation measurements from NASA–CECS Antarctic ice mapping campaigns and NASA Operation IceBridge corrected for oceanic processes from measurements and models, surface velocity measurements from synthetic aperture radar, and high-resolution outputs from regional climate models. The ice thickness change rates are calculated in a Lagrangian reference frame to reduce the effects from advection of sharp vertical features, such as cracks and crevasses, that can saturate Eulerian-derived estimates. We use our method over different ice shelves in Antarctica, which vary in terms of size, repeat coverage from airborne altimetry, and dominant processes governing their recent changes. We find that the Larsen-C Ice Shelf is close to steady state over our observation period with spatial variations in ice thickness largely due to the flux divergence of the shelf. Firn and surface processes are responsible for some short-term variability in ice thickness of the Larsen-C Ice Shelf over the time period. The Wilkins Ice Shelf is sensitive to short-timescale coastal and upper-ocean processes, and basal melt is the dominant contributor to the ice thickness change over the period. At the Pine Island Ice Shelf in the critical region near the grounding zone, we find that ice shelf thickness change rates exceed 40m/yr, with the change dominated by strong submarine melting. Regions near the grounding zones of the Dotson and Crosson ice shelves are decreasing in thickness at rates greater than 40m/yr, also due to intense basal melt. NASA–CECS Antarctic ice mapping and NASA Operation IceBridge campaigns provide validation datasets for floating ice shelves at moderately high resolution when coregistered using Lagrangian methods.

Sutterley, Tyler C.↗