Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Reanalysis”

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 469 records · Page 26

Development and evaluation of a new 4DEnVar-based weakly coupled ocean data assimilation system in E3SMv2

The development, implementation, and evaluation of a new weakly coupled ocean data assimilation (WCODA) system for the fully coupled Energy Exascale Earth System Model version 2 (E3SMv2) utilizing the four-dimensional ensemble variational (4DEnVar) method are presented in this study. The 4DEnVar method, based on the dimension-reduced projection four-dimensional variational (DRP-4DVar) approach, replaces the adjoint model with the ensemble technique, thereby reducing computational demands. Monthly mean ocean temperature and salinity data from the EN4.2.1 reanalysis are integrated into the ocean component of E3SMv2 from 1950 to 2021 with the goal of providing realistic initial conditions for decadal predictions and predictability studies. The performance of the WCODA system is assessed using various metrics, including the reduction rate of the cost function, root mean square error (RMSE) differences, correlation differences, and model biases. Results indicate that the WCODA system effectively assimilates the reanalysis data into the climate model, consistently achieving negative reduction rates of the cost function and notable improvements in RMSE and correlation across various ocean layers and regions. Significant enhancements are observed in the upper ocean layers across the majority of global ocean regions, particularly in the north Atlantic, north Pacific, and Indian Ocean. Model biases in sea surface temperature and salinity are also substantially reduced. For sea surface temperature, cold biases in the north Pacific and north Atlantic are diminished by about 1–2 °C, and warm biases in the Southern Ocean are corrected by approximately 1.5–2.5 °C. In terms of salinity, improvements are observed with bias reductions of about 0.5–1 psu in the north Atlantic and north Pacific and up to 1.5 psu in parts of the Southern Ocean. The ultimate goal of the WCODA system is to advance the predictive capabilities of E3SM for subseasonal to decadal climate predictions, thereby supporting research on strategic energy-sector policies and planning.

54 ENVIRONMENTAL SCIENCES↗

Intense windstorms in the northeastern United States

Abstract. Windstorms are a major natural hazard in many countries. The objective of this study is to identify and characterize intense windstorms during the last 4 decades in the US Northeast and determine both the sources of cyclones responsible for these events and the manner in which those cyclones differ from the cyclone climatology. The windstorm detection is based on the spatial extent of locally extreme wind speeds at 100 m height from the ERA5 reanalysis database. During the top 10 windstorms, wind speeds exceed their local 99.9th percentile over at least one-third of land-based ERA5 grid cells in this high-population-density region of the USA. Maximum sustained wind speeds at 100 m during these windstorms range from 26 to over 43 ms−1, with wind speed return periods exceeding 6.5 to 106 years (considering the top 5 % of grid cells during each storm). Property damage associated with these storms, with inflation adjusted to January 2020, ranges from USD 24 million to over USD 29 billion. Two of these windstorms are linked to decaying tropical cyclones, three are Alberta clippers, and the remaining storms are Colorado lows. Two of the 10 re-intensified off the east coast, leading to development of nor'easters. These windstorms followed frequently observed cyclone tracks but exhibit maximum intensities as measured using 700 hPa relative vorticity and mean sea level pressure that is 5–10 times the mean values for cyclones that followed similar tracks over this 40-year period. The time evolution of wind speeds and concurrent precipitation for those windstorms that occurred after the year 2000 exhibit good agreement with in situ ground-based and remote sensing observations, plus storm damage reports, indicating that the ERA5 reanalysis data have a high degree of fidelity for large, damaging windstorms such as these. A larger pool of the top 50 largest windstorms exhibit evidence of only weak serial clustering, which is in contrast to the relatively strong serial clustering of windstorms in Europe.

58 GEOSCIENCES↗

Evaluation of a high-resolution regional climate simulation for surface and hub-height wind climatology over North America

Assessing the availability of key wind resources requires augmenting observations to support the implementation of wind energy infrastructure. However, observations are limited, necessitating the development of high-resolution, long-term gridded datasets. This study presents a robust, dynamically downscaled climatological dataset, offering 20 years of hourly wind data at a 4 km spatial resolution across North America, and evaluates its performance against observations, including meteorological towers and automated surface-observing system (ASOS) stations, as well as coarse-resolution reanalysis data (the European Centre for Medium-Range Weather Forecasts (ECMWF) reanalysis version 5 (ERA5)). Results demonstrate that the downscaled high-resolution wind data outperform ERA5 in regions of complex terrain and coastal areas, with improved overlap coefficients for wind data distributions and reduced root mean square errors (RMSEs) for hub-height and near-surface diurnal wind patterns. The downscaled simulation also captures the synoptic drivers of seasonal wind direction patterns reasonably well, indicated by high wind rose similarity indices. This study also provides an analysis of interannual variability, utilizing the dataset's full 20-year period, and model uncertainty, generated by varying model initial conditions and physics parameterizations across 1-year ensemble members, which are key considerations for wind resource assessment in wind farm development.

17 WIND ENERGY↗

Evaluation of obstacle modelling approaches for resource assessment and small wind turbine siting: case study in the northern Netherlands

Abstract. Growth in adoption of distributed wind turbines for energy generation is significantly impacted by challenges associated with siting and accurate estimation of the wind resource. Small turbines, at hub heights of 40 m or less, are greatly impacted by terrestrial obstacles such as built structures and vegetation that can cause complex wake effects. While some progress in high-fidelity complex fluid dynamics (CFD) models has increased the potential accuracy for modelling the impacts of obstacles on turbulent wind flow, these models are too computationally expensive for practical siting and resource assessment applications. To understand the efficacy of available models in situ, this study evaluates classic and commonly used methods alongside new state-of-the-art lower-order models derived from CFD simulations and machine learning approaches. This evaluation is conducted using a subset of an extensive original dataset of measurements from more than 300 operational wind turbines in the northern Netherlands. The results show that data-driven methods (e.g. machine learning and statistical modelling) are most effective at predicting production at real sites with an average error in annual energy production of 2.5 %. When sufficient data may not be available de novo to support these data-driven approaches, models derived from high-fidelity simulations show promise and reliably outperform classic methods. On average these models have 6.3 %–11.5 % error compared with 26 % for classic methods and 27 % baseline error for reanalysis data without obstacle correction. While more performant on average, these methods are also sensitive to the quality of obstacle descriptions and reanalysis inputs.

17 WIND ENERGY↗

ERA5-Land Data for LASSO-CACTI Overview Paper

The European Centre for Medium-Range Weather Forecasts (ECMWF) generated a soil reanalysis dataset for the land component of the fifth generation of European ReAnalysis (ERA5), referred to as ERA5-Land. This is a model-generated dataset, with the original version available for the period 1950 to present. The version archived in this DOE ARM product is a subset of the data is for the period of the CACTI field campaign plus several preceding months, specifically from August 1, 2018 through March 22, 2019 with hourly intervals. The ARM copy is also a sub-region of the original global product; the ARM copy is for -60 to -5 °N by -105 to -30 °W. Only variables necessary to drive the WRF-Hydro model are included, which are the 2-m temperature and specific humidity, 10-m wind components, surface pressure, rain rate, and downward surface short and longwave radiation. These data have been obtained from the Copernicus Data Store.

10m wind u-component↗

An exact zooming method for finite element analyses

An exact zooming technique which employs static condensation and exact structural reanalysis methods was developed. Successive application of static condensation reduces the system to one that is only associated with the degrees of freedom (dof) of the original model. Application of an exact static reanalysis technique permits the displacements at the dof of the original model that are contained in the zoomed portion of the structure to be obtained first. The response external to the zoom, as well as the response of additional dof within various levels of zooming, is computed. With the triangular factor of the stiffness matrix of the original system available, this approach involves only the solution of a system of equations of small order.

Hirai, I.↗

An improved search for large-scale convection cells in the solar atmosphere

A reanalysis of Mount Wilson solar velocity observations was made to search for giant cellular patterns. The reanalysis avoids several errors made in a previous search. No cells are detected with sensitivity of 3 to 12 m/s depending upon wavenumber. The observed amplitudes do not conflict with recent model predictions.

Labonte, B. J.↗

Stress constraints in optimality criteria design

Procedures described emphasize the processing of stress constraints within optimality criteria designs for low structural weight with stress and compliance constraints. Prescreening criteria are used to partition stress constraints into either potentially active primary sets or passive secondary sets that require minimal processing. Side constraint boundaries for passive constraints are derived by projections from design histories to modify conventional stress-ratio boundaries. Other procedures described apply partial structural modification reanalysis to design variable groups to correct stress constraint violations of unfeasible designs. Sample problem results show effective design convergence and, in particular, advantages for reanalysis in obtaining lower feasible design weights.

Levy, R.↗

MR Cygni revisited

New analysis tools and additional unanalyzed observations justify a reanalysis of MR Cygni. The reanalysis applied successively more restrictive physical models, each with an optimization program. The final model assigned separate first and second order limb darkening coefficients, from model atmospheres, to individual grid points. Proper operation of the optimization procedure was tested on simulated observational data, produced by light synthesis with assigned system parameters, and modulated by simulated observational error. The iterative solution converged to a weakly-determined mass ratio of 0.75. Assuming the B3 primary component is on the main sequence, the HR diagram location of the secondary from the light ratio (ordinate) and adjusted T sub eff (abscissa) was calculated. The derived mass ratio, together with a main-sequence mass for the B3 component, implies a main-sequence secondary spectral type of B4. The photometrically-determined secondary radii agree with this spectral type, in marginal disagreement with the B7 type from the HR diagram analysis. The individual masses, derived from the radial velocity curve of the primary component, the photometrically-determined i, and alternative values of derived mass ratio are seriously discrepant with main sequence objects. The imputed physical status of the system is in disagreement with representations that have appeared in the literature.

Linnell, Albert P.↗

A hierarchical structure for automatic meshing and adaptive FEM analysis

A new algorithm for generating automatically, from solid models of mechanical parts, finite element meshes that are organized as spatially addressable quaternary trees (for 2-D work) or octal trees (for 3-D work) is discussed. Because such meshes are inherently hierarchical as well as spatially addressable, they permit efficient substructuring techniques to be used for both global analysis and incremental remeshing and reanalysis. The global and incremental techniques are summarized and some results from an experimental closed loop 2-D system in which meshing, analysis, error evaluation, and remeshing and reanalysis are done automatically and adaptively are presented. The implementation of 3-D work is briefly discussed.

Kela, Ajay↗

Boundary formulations for sensitivity analysis without matrix derivatives

A new hybrid approach to continuum structural shape sensitivity analysis employing boundary element analysis (BEA) is presented. The approach uses iterative reanalysis to obviate the need to factor perturbed matrices in the determination of surface displacement and traction sensitivities via a univariate perturbation/finite difference (UPFD) step. The UPFD approach makes it possible to immediately reuse existing subroutines for computation of BEA matrix coefficients in the design sensitivity analysis process. The reanalysis technique computes economical response of univariately perturbed models without factoring perturbed matrices. The approach provides substantial computational economy without the burden of a large-scale reprogramming effort.

Kane, J. H.↗

Preliminary solutions for the lunar gravity field from analysis of lunar orbiter tracking data

Knowledge of the gravitation field, in combination with surface topography, provides one of the principal means of inferring the internal structure of a planetary body. Previous analyses of the lunar gravitational field have been based on data from the Lunar Orbiters, the Apollo subsatellites, and the low altitude passes of the Apollo spacecraft. Recently, Konopliv et al. have reanalyzed all available Lunar Orbiter and Apollo subsatellite tracking data, producing a 60th degree and order solution. In preparation for the Clementine Mission to the Moon, we have also initiated a reanalysis of the Lunar Orbiter and Apollo subsatellite data. Our reanalysis takes advantage of advanced force and measurement modeling techniques as well as modern computational facilities. We applied the least squares collocation technique which stabilizes the behavior of the solution and high degree and order. The extension of the size of the field reduces the aliasing coming from the omitted portion of the gravitational field. This is especially important for the analysis of the tracking data from the Lunar Orbiters, as the periapse heights frequently ranged from 50 to 100 km.

Lemoine, F. G.↗

Technical report series on global modeling and data assimilation. Volume 4: Documentation of the Goddard Earth Observing System (GEOS) data assimilation system, version 1

This report describes the analysis component of the Goddard Earth Observing System, Data Assimilation System, Version 1 (GEOS-1 DAS). The general features of the data assimilation system are outlined, followed by a thorough description of the statistical interpolation algorithm, including specification of error covariances and quality control of observations. We conclude with a discussion of the current status of development of the GEOS data assimilation system. The main components of GEOS-1 DAS are an atmospheric general circulation model and an Optimal Interpolation algorithm. The system is cycled using the Incremental Analysis Update (IAU) technique in which analysis increments are introduced as time independent forcing terms in a forecast model integration. The system is capable of producing dynamically balanced states without the explicit use of initialization, as well as a time-continuous representation of non- observables such as precipitation and radiational fluxes. This version of the data assimilation system was used in the five-year reanalysis project completed in April 1994 by Goddard's Data Assimilation Office (DAO) Data from this reanalysis are available from the Goddard Distributed Active Center (DAAC), which is part of NASA's Earth Observing System Data and Information System (EOSDIS). For information on how to obtain these data sets, contact the Goddard DAAC at (301) 286-3209, EMAIL daac@gsfc.nasa.gov.

Suarez, Max J.↗

Implications of a high population 2 B/Be ratio

The observed boron/beryllium ratio in extreme Population II stars has been interpreted as evidence of Be and B synthesis by early Galactic cosmic rays. However, a recent reanalysis of the boron abundance in the Population II halo star HD 140823 suggests that B/H may be larger than previously reported, by as much as a factor of 4. This would yield a B/Be ratio lying in the range 14 less than or approximately equal to B/Be less than or approximately equal to 50. The possibilty of a high Population II B/Be ratio stresses the importance of the upper limit to the B/Be ratio arising from cosmic-ray production. It is found that the limit to cosmic-ray-produced B/Be depends upon the assumed cosmic-ray spectrum. For any Population II cosmic-ray spectrum that is a single power law in either total energy per nucleon or in momentum, the B/Be ratio constrained to lie in the range 7.6 less than or approximately equal to B/Be less than or approximately equal to 14. Thus, if the new B/Be ratio is correct, it requires either a bimodal cosmic-ray flux with a large low-energy component, or, for another B source, possibly the proposed nu-process in supernovae, either of which may be helpful in explaining the observed B-11/B-10 ratio. Finally, it is noted that the boron reanalysis highlights the uncertainty in our knowledge of the B/Be ratio, and the need for additional data on Be and B abundances.

Fields, Brian D.↗

A Satellite-Derived Upper-Tropospheric Water Vapor Transport Index for Climate Studies

A new approach is presented to quantify upper-level moisture transport from geostationary satellite data. Daily time sequences of GOES-7 water vapor imagery were used to produce estimates of winds and water vapor mixing ratio in the upper-troposphere sensed by the 6.7-microns water vapor channel. The winds and mixing ratio values were gridded and then combined to produce a parameter called the Water Vapor Transport Index (WVTI) which represents the magnitude of the two dimensional transport of water vapor in the upper troposphere. Daily grids of WVTI, meridional moisture transport, mixing ratio, pressure and other associated parameters were averaged to produce monthly fields for June, July and August of 1987 and 1988 over the Americas and surrounding oceanic regions. The WVTI was used to compare upper-tropospheric moisture transport between the summers of 1987 and 1988, contrasting the latter part of the 1986/87 El Nino event and the La Nina period of 1988. A similar product derived from the NCEP reanalysis was used to help validate the index and understand interannual variability of moisture transport from the modeling perspective. Both datasets showed large upper-level water vapor transport associated with synoptic features over the Americas and with outflow from tropical convective systems. Minimal transport occurred over tropical and subtropical high pressure regions where winds were light. Index values from NCEP were 2-3 times larger than that determined from the GOES satellite This difference resulted from an over estimate of the zonal winds and upper-tropospheric moisture in the reanalysis model. A comparison of the satellite-derived monthly values between the summers of 1987 and 1988 provided some insight into the impact of the ENSO event on upper-level moisture and its transport during the period. During July 1987, a large portion of the tropics in the eastern Pacific Ocean and Caribbean Sea was dominated by strong vapor transport in excess of 5 g/kg/ms with relatively small amounts in the other months. In contrast, JJA 1988 showed an opposite trend with July 1988 being less dominated by tropical water vapor transport. The meridional transport of upper-level water vapor indicated significant poleward transport from the tropics to the higher latitudes. This transport favored the Southern Hemisphere with large transport occurring south of the ITCZ which extended across the eastern Pacific and northern South America. Zonally-averaged monthly transport values were shown to provide a simple way to quantify the monthly and interannual changes in water vapor transport. Zonally-averaged WVTI values peaked in the Southern Hemisphere subtropics during both Austral winters. In the tropics a single more pronounced peak located over the equator and south latitudes occurred in 1988 as opposed to a dual peak in 1987. The second peak around 20 deg.N latitude is consistent with findings of others where upper-tropospheric winds were noted to be stronger in this region during warm ENSO events. Zonally-averaged meridional transport was southward for all summer months and stronger in 1988. The asymmetric nature of the zonally-averaged meridional transport (more southerly water vapor transport) was enhanced during JJA 1998 thus indicating a stronger upper-level branch of the Hadley circulation during this notable strong La Nina period.

Jedlovec, Gary J.↗

Interannual Variability of the Mosaic Land-Surface Model

Recently, NASA Goddard Earth Observing System (GEOS-1) reanalysis data has been used to provide forcing for the Koster and Suarez Mosaic Land-surface Model (LSM). The LSM was integrated off-line at all global land points for the period of 1983 - 1995 by the Off-line Land- surface GEOS Assimilation system (OLGA). Here, we compare the interannual variability of OLGA, GEOS-1 and surface observing stations temperature and moisture. Particular attention is given to the United States because of the extreme seasons of 1988 and 1993. Furthermore, the comparison of OLGA is extended to include the analysis of data on the'tiles' (different surface types) in the Mosaic LSM. Results indicate that the GEOS-1 near-surface temperature and moisture reasonably represents the interannual variability in more normal years. However, OLGA also simulates the extreme drought and floods years well. The analysis of the tile information shows that the "Bare soil" surface type is most sensitive to the climate extremes. Off-line testing has provided valuable information on the performance of the Mosaic LSM prior to its incorporation into the new version of the GEOS Data Assimilation System and the integration of a new long reanalysis.

Bosilovich, Michael G.↗

Wintertime East Asian Jet Stream and Its Association with the Asian-Pacific Climate

Interannual variability of the wintertime East Asian westerly jet stream and the linkage between this variability and the Asian-Pacific climate are investigated. The study emphasizes on the variability of the jet core and its association with the Asian winter monsoon, tropical convection, upper tropospheric wave patterns, and the teleconnection of the jet with other climate systems. The relationship between the jet and North Pacific sea surface temperature pattern (SST) is also explored. NCEP/NCAR reanalysis, NASA GISS surface temperature, NASA GEOS reanalysis, NOAA reconstructed SST, GPCP precipitation, and NOAA snow cover data sets are analyzed in this study. An index of the East Asian jet has been defined by the December-February means of the 200 mb zonal winds that are averaged within a box enclosing the jet maximum, which shifts only moderately from one year to another especially in the south-north direction. The jet links to a teleconnection pattern whose major climate anomalies appear over the Asian continent and western Pacific (west of the dateline). This pattern differs distinctly from the teleconnection pattern associated with El Nino/Southern Oscillation (ENSO), which causes the Pacific/North American pattern to the east of the dateline. A strong jet is accompanied clearly by an increase in the intensity of the atmospheric circulation over Asia and the Pacific. In particular, the winter monsoon strengthens over East Asia, leading to cold climate in the region, and convection intensifies over the tropical Asia-Australia sector. Changes in the jet are associated with broad-scale modification in the upper tropospheric wave patterns that leads to downstream climate anomalies over the eastern Pacific. Through this downstream influence, the East Asian jet causes climate signals in North America as well. A strong jet gives rise to warming and less snow cover in the western United States but reverse climate anomalies in the eastern part of the country, although these signals are relatively weaker than the jet-related anomalies in East Asia. There is a strong association between the East Asian jet and the North Pacific SST (NPSST). A strong jet is accompanied by a cooling in the extratropical Pacific and a warming in the tropical-subtropical Pacific. Evidence also indicates that the extratropical NPSST pattern plays a role in modulating the intensity of the jet stream. ENSO, the jet, and the NPSST are mutually interactive on certain time scales and such an interaction links closely to the climate anomalies in the Asian-Pacific-American regions.

Yang, Song↗

Comparison of TOPEX/Poseidon Sea Level and Linear Model Results forced by Various Wind Products for the Tropical Pacific

The goal of this paper is to compare TOPEX/Posaidon (T/P) sea level with sea level results from linear ocean model experiments forced by several different wind products for the tropical Pacific. During the period of this study (October 1992 - October 1995), available wind products include satellite winds from the ERS-1 scatterometer product of [HALP 97] and the passive microwave analysis of SSMI winds produced using the variational analysis method (VAM) of [ATLA 91]. In addition, atmospheric GCM winds from the NCEP reanalysis [KALN 96], ECMWF analysis [ECMW94], and the Goddard EOS-1 (GEOS-1) reanalysis experiment [SCHU 93] are available for comparison. The observed ship wind analysis of FSU [STRI 92] is also included in this study. The linear model of [CANE 84] is used as a transfer function to test the quality of each of these wind products for the tropical Pacific. The various wind products are judged by comparing the wind-forced model sea level results against the T/P sea level anomalies. Correlation and RMS difference maps show how well each wind product does in reproducing the T/P sea level signal. These results are summarized in a table showing area average correlations and RMS differences. The large-scale low-frequency temporal signal is reproduced by all of the wind products, However, significant differences exist in both amplitude and phase on regional scales. In general, the model results forced by satellite winds do a better job reproducing the T/P signal (i.e. have a higher average correlation and lower RMS difference) than the results forced by atmospheric model winds.

Hackert, Eric C.↗