Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Primitive equations”

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 235 records · Page 13

Ocean Thermal and Color Evolution During the 1997/1998 ENSO Event

A reduced gravity primitive equation modeling and assimilation system is used to study the evolution of the tropical Pacific during the 1997/1998 ENSO cycle. The modeling/assimilation scheme ingests satellite altimeter data and TAO temperature profiles and uses SSM/I satellite derived winds as surface boundary forcing. The four-dimensional structure of the upper ocean circulation structure will be compared against available in situ observations across the Pacific basin. In particular, variability near the Galapagos Islands will be highlighted during the spring of 1998 when phytoplankton concentrations were observed to increase a hundred-fold over a two week period.

Rienecker, Michele↗

An Extended Kalman Filter to Assimilate Altimetric Data into a Non-Linear Model of the Tropical Pacific

A new implementation of the extended Kaman filter is developed for the purpose of assimilating altimetric observations into a primitive equation model of the tropical Pacific. Its specificity consists in defining the errors into a reduced basis that evolves in time with the model dynamic. Validation by twin experiments is conducted and the method is shown to be efficient in quasi real conditions. Data from the first 2 years of the Topex/Poseidon mission are assimilated into the Gent & Cane [1989] model. Assimilation results are evaluated against independent in situ data, namely TAO mooring observations.

Gourdeau, L.↗

Can TOPEX/Poseidon Data be Used for Monitoring the Indonesian Throughflow?

The effects of the Indonesian throughflow(ITF) are studied in a reduced gravity, primitive equation, sigma coordinate model. The surface heat fluxes are provided by coupling the ocean GCAF to an advective atmospheric mixed layer model. There is an ENSO related signal in the ITF but the correlation with the Southern Oscillation Index (SOI) is only -0.31. When the winds over the Indian Ocean are held to climatology, this correlation jumps to -0.65 indicating that the non-EASO signal in the ITF is caused by the downstream winds. On interannual time-scales the ITF can be explained in terms of sea level differences between the western Pacific and eastern Indian Oceans when appropriate representative locations are chosen as demonstrated in both model and TOPEX data. It is shown that the main climatological effect of the ITF is to warm the Indian Ocean and cool the Pacific. Mile the main ENSO indices, NIN03 and NIN04, are almost identical with and without the ITF, the total SSTs show significant ENSO dependence.

Murtugudde, Ragu↗

Modeling of Dense Water Production and Salt Transport from Alaskan Coastal Polynyas

The main significance of this paper is that a realistic, three-dimensional, high-resolution primitive equation model has been developed to study the effects of dense water formation in Arctic coastal polynyas. The model includes realistic ambient stratification, realistic bottom topography, and is forced by time-variant surface heat flux, surface salt flux, and time-dependent coastal flow. The salt and heat fluxes, and the surface ice drift, are derived from satellite observations (SSM/I and NSCAT sensors). The model is used to study the stratification, salt transport, and circulation in the vicinity of Barrow Canyon during the 1996/97 winter season. The coastal flow (Alaska coastal current), which is an extension of the Bering Sea throughflow, is formulated in the model using the wind-transport regression. The results show that for the 1996/97 winter the northeastward coastal current exports 13% to 26% of the salt produced by coastal polynyas upstream of Barrow Canyon in 20 to 30 days. The salt export occurs more rapidly during less persistent polynyas. The inclusion of ice-water stress in the model makes the coastal current slightly weaker and much wider due to the combined effects of surface drag and offshore Ekman transport.

Signorini, Sergio R.↗

Data Assimilation and Transport Modeling in Terrestrial and Planetary Atmospheres

Data assimilation is a blanket term used to describe a number of techniques for retrieving important physical parameters from observational data, subject to constraints imposed by prior knowledge (such as, in the case of meteorology, the primitive equations that govern atmospheric motion). Since these newly developed methods make efficient use of computational resources, they are of great importance in the interpretation of the voluminous datasets that are now produced by satellite missions. As proposed, these techniques have been applied to the study of the Martian and terrestrial atmospheres based on available satellite observations. In addition, a sophisticated hydrodynamic model (non-hydrostatic, and therefore applicable to the study of the interiors of the giant planets) has also been developed and successfully applied to the study of tidally induced motions in Jupiter.

Houben, Howard C.↗

Effects of Convective Asymmetries on Hurricane Intensity: A Numerical Study

The influence of the uniform large-scale flow, beta effect, and vertical shear of the environmental flow on hurricane intensity is investigated in the context of the induced convective or potential vorticity asymmetries with a hydrostatic primitive equation hurricane model. In agreement with the previous studies, imposing of one of these environmental effects can substantially weaken the simulated tropical cyclones. In response t o the environmental influence, significant asymmetries develop with a structure similar to the spiral bands in real hurricanes, which are dominated by wavenumber-one components. The tendencies of the mean radial, azimuthal winds and temperature associated with the environment-induced convective asymmetries are evaluated respectively. The resulting asymmetries can effectively reduce hurricane intensity by directly producing the negative tendency of the mean tangential wind in the vicinity of the radius of maximum wind, and by weakening the mean radial circulation. The reduction effects are closely associated with the spiral structure of the induced asymmetries. The time lag observed between the imposition of the environmental influence and the resulting rise in the minimum central pressure is the time required for developing the spiral structure. This study also confirms the axisymmetrization process associated with the induced wavenumber-one components of potential vorticity asymmetries, but it exists only within the radius of maximum wind.

Wu, Liguang↗

Optimal Observations for Variational Data Assimilation

An important aspect of Ocean state estimation is the design of an observing system that allows the efficient study of climate aspects in the ocean. A solution of the design problem is presented here in terms of optimal observations that emerge as nondimensionalized singular vectors of the modified data resolution matrix. The actual computation is feasible only for scalar quantities in the limit of large observational errors. In the framework of a lo resolution North Atlantic primitive equation model it is demonstrated that such optimal observations when applied to determining the strength of the volume and heat transport across the Greenland-Scotland ridge, perform significantly better than traditional section data. On seasonal to inter-annual time-scales optimal observations are located primarily along the continental shelf and information about heat-transport, wind stress and stratification is being communicated via boundary waves and advective processes. On time-scales of about a month, sea surface height observations appear to be more efficient in reconstructing the cross-ridge heat transport than hydrographic observations. Optimal observations also provide a tool for understanding how the ocean state is effected by anomalies of integral quantities such as meridional heat transport.

Koehl, Armin↗

Sea Surface Salinity Variability from Simulations and Observations: Preparing for Aquarius

Oceanic fresh water transport has been shown to play an important role in the global hydrological cycle. Sea surface salinity (SSS) is representative of the surface fresh water fluxes and the upcoming Aquarius mission scheduled to be launched in December 2010 will provide excellent spatial and temporal SSS coverage to better estimate the net exchange. In most ocean general circulation models, SSS is relaxed to climatology to prevent model drift. While SST remains a well observed variable, relaxing to SST reduces the range of SSS variability in the simulations (Fig.1). The main objective of the present study is to simulate surface tracers using a primitive equation ocean model for multiple forcing data sets to identify and establish a baseline SSS variability. The simulated variability scales are compared to those from near-surface argo salinity measurements.

Jacob, S. Daniel↗

Megadroughts in Southwestern North America in ECHO-G Millennial Simulations and Their Comparison to Proxy Drought Reconstructions

Simulated hydroclimate variability in millennium-length forced transient and control simulations from the ECHAM and the global Hamburg Ocean Primitive Equation (ECHO-G) coupled atmosphere-ocean general circulation model (AOGCM) is analyzed and compared to 1000 years of reconstructed Palmer drought severity index (PDSI) variability from the North American Drought Atlas (NADA). The ability of the model to simulate megadroughts in the North American southwest is evaluated. (NASW: 25deg42.5degN, 125deg-105degW). Megadroughts in the ECHO-G AOGCM are found to be similar in duration and magnitude to those estimated from the NADA. The droughts in the forced simulation are not, however, temporally synchronous with those in the paleoclimate record, nor are there significant differences between the drought features simulated in the forced and control runs. These results indicate that model-simulated megadroughts can result from internal variability of the modeled climate system rather than as a response to changes in exogenous forcings. Although the ECHO-G AOGCM is capable of simulating megadroughts through persistent La Nina-like conditions in the tropical Pacific, other mechanisms can produce similarly extreme NASW moisture anomalies in the model. In particular, the lack of low-frequency coherence between NASW soil moisture and simulated modes of climate variability like the El Nino-Southern Oscillation, Pacific decadal oscillation, and Atlantic multidecadal oscillation during identified drought periods suggests that stochastic atmospheric variability can contribute significantly to the occurrence of simulated megadroughts in the NASW. These findings indicate that either an expanded paradigm is needed to understand multidecadal hydroclimate variability in the NASW or AOGCMs may incorrectly simulate the strength and/or dynamics of the connection between NASW hydroclimate variability and the tropical Pacific.

variability↗

The Impact of Satellite Sea Surface Salinity for Prediction of the Coupled Indo-Pacific System

We assess the impact of satellite sea surface salinity (SSS) observations on seasonal to interannual variability of tropical Indo-Pacific Ocean dynamics as well as on dynamical ENSO forecasts. Our coupled model is composed of a primitive equation ocean model for the tropical Indo-Pacific region that is coupled with the global SPEEDY atmospheric model (Molteni, 2003). The Ensemble Reduced Order Kalman Filter is used to assimilate observations to constrain dynamics and thermodynamics for initialization of the coupled model. The baseline experiment assimilates satellite sea level, SST, and in situ subsurface temperature and salinity observations. This baseline is then compared with experiments that additionally assimilate Aquarius (version 4.0) and SMAP (version 2.0) SSS. Twelve-month forecasts are initialized for each month from Sep. 2011 to Dec. 2016. We find that including satellite SSS significantly improves NINO 3.4 sea surface temperature anomaly validation after 1 out to 12 month forecast lead times. For initialization of the coupled forecast, the positive impact of SSS assimilation is brought about by surface freshening near the eastern edge of the western Pacific warm pool and density changes that lead to shallower mixed layer between 10 degrees South latitude-5 degrees North latitude. SST differences at initialization force wide-spread downwelling favorable curl over most of the tropical Pacific. Over an average forecast, SST remains warmer with SSS assimilation at the eastern edge of the warm pool. This warm SST propagates into the eastern Pacific and drags westerly wind anomalies eastward into the NINO 3.4 region. In addition, salting near the ITCZ (Intertropical Convergence Zone) leads to a deepening of the mixed layer and thermocline near 8 degrees North latitude. These patterns together lead to a funneling effect that provides the background state to amplify equatorial Kelvin waves. We show that the downwelling Kelvin waves are amplified by assimilating satellite SSS and lead to significantly improved forecasts particularly for the 2015 El Nino.

El Nino↗

The Impact of Satellite Sea Surface Salinity for Prediction of the Coupled Indo-Pacific System

We assess the impact of satellite sea surface salinity (SSS) observations on seasonal to interannual variability of tropical Indo-Pacific Ocean dynamics as well as on dynamical ENSO forecasts. Our coupled model is composed of a primitive equation ocean model for the tropical Indo-Pacific region that is coupled with the global SPEEDY atmospheric model (Molteni, 2003). The Ensemble Reduced Order Kalman Filter is used to assimilate observations to constrain dynamics and thermodynamics for initialization of the coupled model. The baseline experiment assimilates satellite sea level, SST, and in situ subsurface temperature and salinity observations. This baseline is then compared with experiments that additionally assimilate Aquarius (version 4.0) and SMAP (version 2.0) SSS. Twelve-month forecasts are initialized for each month from Sep. 2011 to Dec. 2016. We find that including satellite SSS significantly improves NINO3.4 sea surface temperature anomaly validation after 1 out to 12 month forecast lead times. For initialization of the coupled forecast, the positive impact of SSS assimilation is brought about by surface freshening near the eastern edge of the western Pacific warm pool and density changes that lead to shallower mixed layer between 10S-5N. SST differences at initialization force wide-spread downwelling favorable curl over most of the tropical Pacific. Over an average forecast, SST remains warmer with SSS assimilation at the eastern edge of the warm pool. This warm SST propagates into the eastern Pacific and drags westerly wind anomalies eastward into the NINO3.4 region. In addition, salting near the ITCZ leads to a deepening of the mixed layer and thermocline near 8N. These patterns together lead to a funneling effect that provides the background state to amplify equatorial Kelvin waves. We show that the downwelling Kelvin waves are amplified by assimilating satellite SSS and lead to significantly improved forecasts particularly for the 2015 El Nino.

SST↗

An Approximate Kalman Filter for Ocean Data Assimilation; An Example with an Idealized Gulf Stream Model

A practical method of data assimilation for use with large, nonlinear, ocean general circulation models is explored. A Kalman filter based on approximations of the state error covariance matrix is presented, employing a reduction of the effective model dimension, the error's asymptotic steady-state limit, and a time-invariant linearization of the dynamic model for the error integration. The approximations lead to dramatic computational savings in applying estimation theory to large complex systems. We examine the utility of the approximate filter in assimilating different measurement types using a twin experiment of an idealized Gulf Stream. A nonlinear primitive equation model of an unstable east-west jet is studied with a state dimension exceeding 170,000 elements. Assimilation of various pseudo measurements is examined, including velocity, density, and volume transport at localized arrays, and realistic distributions of satellite altimetry and acoustic tomography observations. Results are compared in terms of their effects on the accuracies of the estimation. The approximate filter is shown to outperform a previous study that used an empirical nudging scheme. The examples demonstrate that useful approximate estimation errors can be computed in a practical manner for general circulation models.

['Kalman↗

Australian tidal currents – assessment of a barotropic model (COMPAS v1.3.0 rev6631) with an unstructured grid

While the variations of tidal range are large and fairly well known across Australia (less than 1 m near Perth but more than 14 m in King Sound), the properties of the tidal currents are not. We describe a new regional model of Australian tides and assess it against a validation dataset comprising tidal height and velocity constituents at 615 tide gauge sites and 95 current meter sites. The model is a barotropic implementation of COMPAS, an unstructured-grid primitive-equation model that is forced at the open boundaries by TPXO9v1. The mean absolute error (MAE) of the modelled M2 height amplitude is 8.8 cm, or 12 % of the 73 cm mean observed amplitude. The MAE of phase (10°), however, is significant, so the M2 mean magnitude of vector error (MMVE, 18.2 cm) is significantly greater. The root sum square over the eight major constituents is 26 % of the observed amplitude. We conclude that while the model has skill at height in all regions, there is definitely room for improvement (especially at some specific locations). For the M2 major axis velocity amplitude, the MAE across the 95 current meter sites, where the observed amplitude ranges from 0.1 to 156 cm s −1 , is 6.9 cm s −1 , or 22 % of the 31.7 cm s −1 observed mean. This nationwide average result is encouraging, but it conceals a very large regional variation. Relative errors of the tidal current amplitudes on the narrow shelves of New South Wales (NSW) and Western Australia exceed 100 %, but tidal currents are weak and negligible there compared to non-tidal currents, so the tidal errors are of little practical significance. Looking nationwide, we show that the model has predictive value for much of the 79 % of Australia's shelf seas where tides are a major component of the total velocity variability. In descending order this includes the Bass Strait, the Kimberley to Arnhem Land, and southern Great Barrier Reef regions. There is limited observational evidence to confirm that the model is also valuable for currents in other regions across northern Australia. We plan to commence publishing “unofficial” tidal current predictions for chosen regions in the near future based on both our COMPAS model and the validation dataset we have assembled.

Tidal currents↗

Theory and modeling of atmospheric turbulence, part 2

Two dimensional geostrophic turbulence driven by a random force is investigated. Based on the Liouville equation, which simulates the primitive hydrodynamical equations, a group-kinetic theory of turbulence is developed and the kinetic equation of the scaled singlet distribution is derived. The kinetic equation is transformed into an equation of spectral balance in the equilibrium and non-equilibrium states. Comparison is made between the propagators and the Green's functions in the case of the non-asymptotic quasi-linear equation to prove the equivalence of both kinds of approximations used to describe perturbed trajectories of plasma turbulence. The microdynamical state of fluid turbulence is described by a hydrodynamical system and transformed into a master equation analogous to the Vlasov equation for plasma turbulence. The spectral balance for the velocity fluctuations of individual components shows that the scaled pressure strain correlation and the cascade transfer are two transport functions that play the most important roles.

Chen, C. M.↗

Study of three-dimensional effects on vortex breakdown

The incompressible axisymmetric steady Navier-Stokes equations in primitive variables are used to simulate vortex breakdown. The equations, discretized using a second-order, central-difference scheme, are linearized and then solved using an exact LU decomposition, Gaussian elimination, and Newton iteration. Solutions are presented for Reynolds numbers, based on vortex-core radius, as high as 1500. An attempt to study the stability of the axisymmetric solutions against three-dimensional perturbations is discussed.

Salas, M. D.↗