Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Statistical forecasting”

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 289 records · Page 16

Physics-Informed Gaussian Process Regression for States Estimation and Forecasting in Power Grids

Real-time state estimation and forecasting are critical for the efficient operation of power grids. In this paper, a physics-informed Gaussian process regression (PhI-GPR) method is presented and used for forecasting and estimating the phase angle, angular speed, and wind mechanical power of a three-generator power grid system using sparse measurements. In standard data-driven Gaussian process regression (GPR), parameterized models for the prior statistics are fit by maximizing the marginal likelihood of observed data. In the PhI-GPR method, we propose to compute the prior statistics offline by solving stochastic differential equations (SDEs) governing the power grid dynamics. The short-term forecast of a power grid system dominated by wind generation is complicated by the stochastic nature of the wind and the resulting uncertainty in wind mechanical power. Here, we assume that the power grid dynamics are governed by swing equations, with the wind mechanical power fluctuating randomly in time. We solve these equations for the mean and covariances of the power grid states using the Monte Carlo simulation method. We demonstrate that the proposed PhI-GPR method can accurately forecast and estimate observed and unobserved states. For the considered problem, PhI-GPR has computational advantages over the ensemble Kalman filter (EnKF) method: In PhI-GPR, ensembles are computed offline and independently of the data acquisition process, whereas for EnFK, ensembles are computed online with data acquisition, rendering real-time forecast more challenging. We also demonstrate that the PhI-GPR forecast is more accurate than the EnKF forecast when the random mechanical wind power is non-Markovian. In contrast, the two methods produce similar forecasts for the Markovian mechanical wind power. For observed states, we show that PhI-GPR provides a forecast comparable to the standard data-driven GPR; both forecasts are significantly more accurate than the autoregressive integrated moving average (ARIMA) forecast. We also show that the ARIMA forecast is more sensitive to observation frequency and measurement errors than the PhI-GPR forecast.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Investigating the potential of a global precipitation forecast to inform landslide prediction

Extreme rainfall events within landslide-prone areas can be catastrophic, resulting in loss of property, infrastructure, and life. A global Landslide Hazard Assessment for Situational Awareness (LHASA) model provides routine near-real time estimates of landslide hazard using Integrated Multi-Satellite Precipitation Retrievals for the Global Precipitation Mission (IMERG). However, it does not provide information on potential landslide hazard in the future. Forecasting potential landslide events at a global scale presents an area of open research. This study compares a global precipitation forecast provided by NASA's Goddard Earth Observing System (GEOS) with near-real time satellite precipitation estimates. The Multi-Radar Multi-Sensor gauge corrected (MRMS-GC) reference is used to assess the performance of both satellite and model-based precipitation products over the contiguous United States (CONUS). The forecast lead time of 24hrs is considered, with a focus on extreme precipitation events. The performance of IMERG and GEOS-Forecast products is assessed in terms of the probability of detection, success ratio, critical success index and hit bias as well as continuous statistics. The results show that seasonality influences the performance of both satellite and model-based precipitation products. Comparison of IMERG and GEOS-Forecast globally as well as in several event case studies (Colombia, southeast Asia, and Tajikistan) reveals that GEOS-Forecast detects extreme rainfall more frequently relative to IMERG for these specific analyses. For recent landslide points across the globe, the 24hr accumulated precipitation forecast >100 mm corresponds well with near-real time daily accumulated IMERG precipitation estimates. GEOS-Forecast and IMERG precipitation match more closely for tropical cyclones than for other types of storms. The main intention of this study is to assess the viability of using a global forecast for landslide predictions and understand the extent of the variability between these products to inform where we would expect the landslide modeling results to most prominently diverge. Results of this study will be used to inform how forecasted precipitation estimates can be incorporated into the LHASA model to provide the first global predictive view of landslide hazards.

S. Khan↗

Investigating the Potential of A Global Precipitation Forecast to Inform Landslide Prediction

Extreme rainfall events within landslide-prone areas can be catastrophic, resulting in loss of property, infra-structure, and life. A global Landslide Hazard Assessment for Situational Awareness (LHASA) model provides routine near-real time estimates of landslide hazard using Integrated Multi-Satellite Precipitation Retrievals for the Global Precipitation Mission (IMERG). However, it does not provide information on potential landslide hazard in the future. Forecasting potential landslide events at a global scale presents an area of open research. This study compares a global precipitation forecast provided by NASA’s Goddard Earth Observing System (GEOS) with near-real time satellite precipitation estimates. The Multi-Radar Multi-Sensor gauge corrected (MRMS-GC) reference is used to assess the performance of both satellite and model-based precipitation products over the contiguous United States (CONUS). The forecast lead time of 24hrs is considered, with a focus on extreme precipitation events. The performance of IMERG and GEOS-Forecast products is assessed in terms of the probability of detection, success ratio, critical success index and hit bias as well as continuous statistics. The results show that seasonality influences the performance of both satellite and model-based precipitation products. Comparison of IMERG and GEOS-Forecast globally as well as in several event case studies (Colombia, southeast Asia, and Tajikistan) reveals that GEOS-Forecast detects extreme rainfall more frequently relative to IMERG for these specific analyses. For recent landslide points across the globe, the 24hr accumulated precipitation forecast >100 mm corresponds well with near-real time daily accumulated IMERG precipitation estimates. GEOS-Forecast and IMERG precipitation match more closely for tropical cyclones than for other types of storms. The main intention of this study is to assess the viability of using a global forecast for landslide predictions and understand the extent of the variability between these products to inform where we would expect the landslide modeling results to most prominently diverge. Results of this study will be used to inform how forecasted precipitation estimates can be incorporated into the LHASA model to provide the first global predictive view of landslide hazards.

S. Khan↗

Intercomparison of Dynamically and Statistically Downscaled Climate Change Projections over the Midwest and Great Lakes Region

Downscaling of global climate model (GCMs) simulations is a key element of regional-to-local-scale climate change projections that can inform impact assessments, long-term planning, and resource management in different sectors. Here, we conduct an intercomparison between statistically and dynamically downscaled GCMs simulations using the hybrid delta (HD) and the Weather Research and Forecast (WRF) Model, respectively, over the Midwest and Great Lakes region to 1) validate their performance in reproducing extreme daily precipitation (P) and daily maximum temperature (T max ) for summer and winter and 2) evaluate projections of extremes in the future. Our results show the HD statistical downscaling approach, which includes large-scale bias correction of GCM inputs, can reproduce observed extreme P and T max reasonably well for both summer and winter. However, raw historical WRF simulations show significant bias in both extreme P and T max for both seasons. Interestingly, the convection-permitting WRF simulation at 4-km grid spacing does not produce better results for seasonal extremes than the WRF simulation at 12 km using a parameterized convection scheme. Despite a broad similarity for winter extreme P projections, the projected changes in the future summer storms are quite different between downscaling methods; WRF simulations show substantial increases in summer extreme precipitation, while the changes projected by the HD approach exhibit moderate decreases overall. The WRF simulations at 4 km also show a pronounced decoupling effect between seasonal totals and extreme daily P for summer, which suggests that there could be more intense summer extremes at two different time scales, with more severe individual convective storms combined with longer summer droughts at the end of the twenty-first century.

54 ENVIRONMENTAL SCIENCES↗

Photometric redshift uncertainties in weak gravitational lensing shear analysis: models and marginalization

ABSTRACT Recovering credible cosmological parameter constraints in a weak lensing shear analysis requires an accurate model that can be used to marginalize over nuisance parameters describing potential sources of systematic uncertainty, such as the uncertainties on the sample redshift distribution n(z). Due to the challenge of running Markov chain Monte Carlo (MCMC) in the high-dimensional parameter spaces in which the n(z) uncertainties may be parametrized, it is common practice to simplify the n(z) parametrization or combine MCMC chains that each have a fixed n(z) resampled from the n(z) uncertainties. In this work, we propose a statistically principled Bayesian resampling approach for marginalizing over the n(z) uncertainty using multiple MCMC chains. We self-consistently compare the new method to existing ones from the literature in the context of a forecasted cosmic shear analysis for the HSC three-year shape catalogue, and find that these methods recover statistically consistent error bars for the cosmological parameter constraints for predicted HSC three-year analysis, implying that using the most computationally efficient of the approaches is appropriate. However, we find that for data sets with the constraining power of the full HSC survey data set (and, by implication, those upcoming surveys with even tighter constraints), the choice of method for marginalizing over n(z) uncertainty among the several methods from the literature may modify the 1σ uncertainties on Ωm–S8 constraints by ∼4 per cent, and a careful model selection is needed to ensure credible parameter intervals.

Zhang, Tianqing (ORCID:000000025596198X)↗

A novel conditional generative model for efficient ensemble forecasts of state variables in large-scale geological carbon storage

Integrating monitoring data to efficiently update reservoir pressure and CO 2 plume distribution forecasts presents a significant challenge in geological carbon storage (GCS) applications. Inverse modeling techniques are commonly used to fuse observational data and refine reservoir model parameters, thereby improving state variable forecasts. However, these techniques often rely on linear or Gaussian assumptions, which can limit their effectiveness in accurately predicting state variables. Moreover, simulating large-scale three-dimensional (3D) GCS problems is computationally expensive, making iterative runs in inverse problems prohibitive. To address these challenges, we propose a conditional generative model utilizing the score-based diffusion method for real-time 3D pressure and saturation field distribution predictions. Our approach involves solving the score function with a mini-batch-based Monte Carlo estimator to generate labeled data. This data is subsequently employed to train a fully connected neural network, enabling it to learn the conditional sample generator within a supervised learning framework. This method enables the rapid generation of a large ensemble of predictions, facilitating comprehensive uncertainty quantification of state variables. Here we applied our method to forecast the dynamic 3D distributions of pressure and saturation fields over a 30-year injection period. The statistical assessment with low root mean square error (RMSE) values demonstrates that our method can accurately predict the spatiotemporal distributions of both pressure and saturation fields. Moreover, the developed conditional generative model shows high computational efficiency by generating 100 ensemble forecasts of 3D state variables in less than 10 min. The consistency between ensemble averages and ground truth values further illustrates the model’s capability to capture state variable dynamics during the CO 2 plume injection process. Notably, the ground truth values fall within the ensemble forecasts, indicating that our uncertainty quantification effectively captures variability and potential noise in the observations. Thus, the developed conditional generative model proves to be a more efficient, accurate, and practical tool for GCS applications, facilitating timely risk analysis and informed decision-making.

58 GEOSCIENCES↗

Coordinated Ramping Product and Regulation Reserve Procurements in CAISO and MISO using Multi-Scale Probabilistic Solar Power Forecasts (Pro2R)

How can probabilistic solar forecasts lower costs and improve reliability for independent system operator (ISO) markets? We tackle this question in three steps. First, we enhance an existing solar forecasting system to provide well-calibrated hours-ahead probabilistic forecasts. We then relate the degree of uncertainty in those forecasts to error distributions for net load ramps for the California ISO (CAISO) using statistical and machine learning methods. Projected net load errors conditioned on solar uncertainty are translated into flexible ramp requirements that therefore reflect real-time meteorological and solar conditions, improving on typical ISO procedures. Finally, a multi-period look-ahead production cost model quantifies how conditional ramp requirements can a) decrease operating costs by lowering requirements compared to often conservative unconditional methods, and b) reduce generation scarcity events and consequently improve reliability by increasing flexibility requirements at times when unconditional forecast-based requirements understate actual ramp uncertainty. In addition to the products just described (quantification of solar uncertainty, its translation into requirements for ramp capability product, and quantification of the benefits of more accurate ramp requirements), this project also developed a visualization system that alerts system operators of ramp and uncertainty conditions within the network based on solar forecasts. The system is called Resource Forecast and Ramp Visualization for Situational Awareness (RaVIS). These four products represent significant advances in the state-of-the-art of probabilistic solar forecasting, development of weather-informed reserve requirements, production costing methods for estimating the benefits of more accurate reserve requirements, and visualization of system status, respectively. Yet the products are also practical and can be immediately implemented, potentially enabling system operators to save millions of dollars in ramp product procurement costs per year.

14 SOLAR ENERGY↗

Cross-Market Price Difference Forecast Using Deep Learning for Electricity Markets

Price forecasting is in the center of decision making in electricity markets. Many researches have been done in forecasting energy prices while little research has been reported on forecasting price difference between day-ahead and realtime markets due to its high volatility, which however plays a critical role in virtual trading. To this end, this paper takes the first attempt to employ novel deep learning architecture with Bidirectional Long-Short Term Memory (LSTM) units to forecast the price difference between day-ahead and real-time markets for the same node. The raw data is collected from PJM market, processed and fed into the proposed network. The Root Mean Squared Error (RMSE) and customized performance metric are used to evaluate the performance of the proposed method. Case studies show that it outperforms the traditional statistical models like ARIMA, and machine learning models like XGBoost and SVR methods in both RMSE and the capability of forecasting the sign of price difference. Additionally, to cross-market price difference forecast, the proposed approach has the potential to be applied to solve other forecasting problems such as price spread forecast in DA market for Financial Transmission Right (FTR) trading purpose.

DA/RT price difference↗

Developing Passenger Demand Models for International Aviation from/to Egypt: A Case Study of Cairo Airport and Egyptair

This research is concerned with developing passenger demand models for international aviation from/to Egypt. In this context, aviation sector in Egypt is represented by the biggest and main airport namely Cairo airport as well as by the main Egyptian international air carrier namely Egyptair. The developed models utilize two variables to represent aviation demand, namely total number of international flights originating from and attracted to Cairo airport as well as total number of passengers using Egyptair international flights originating from and attracted to Cairo airport. Such demand variables were related, using different functional forms, to several explanatory variables including population, GDP and number of foreign tourists. Finally, two models were selected based on their logical acceptability, best fit and statistical significance. To demonstrate usefulness of developed models, these were used to forecast future demand patterns.

Abbas, Khaled A.↗

NASA Tech Briefs, November 2013

Topics include: Cryogenic Liquid Sample Acquisition System for Remote Space Applications; 5 Spatial Statistical Data Fusion (SSDF); GPS Estimates of Integrated Precipitable Water Aid Weather Forecasters; Integrating a Microwave Radiometer into Radar Hardware for Simultaneous Data Collection Between the Instruments; Rapid Detection of Herpes Viruses for Clinical Applications; High-Speed Data Recorder for Space, Geodesy, and Other High-Speed Recording Applications; Datacasting V3.0; An All-Solid-State, Room-Temperature, Heterodyne Receiver for Atmospheric Spectroscopy at 1.2 THz; Stacked Transformer for Driver Gain and Receive Signal Splitting; Wireless Integrated Microelectronic Vacuum Sensor System; Fabrication Method for LOBSTER-Eye Optics in <110> Silicon; Compact Focal Plane Assembly for Planetary Science; Fabrication Methods for Adaptive Deformable Mirrors; Visiting Vehicle Ground Trajectory Tool; Workflow-Based Software Development Environment; Mobile Thread Task Manager; A Kinematic Calibration Process for Flight Robotic Arms; Magnetostrictive Alternator; Bulk Metallic Glasses and Composites for Optical and Compliant Mechanisms; Detection of Only Viable Bacterial Spores Using a Live/Dead Indicator in Mixed Populations; and Intravenous Fluid Generation System.

Source record↗

Recent Updates in the SMAP Level-4 Soil Moisture Algorithm

The NASA Soil Moisture Active Passive (SMAP) mission generates, among other data sets, the Level-4 Soil Moisture (L4_SM) product. The L4_SM data are published with a mean latency of ~2.5 days from the time of observation and provide global, three-hourly, 9-km resolution estimates of surface and root-zone soil moisture and related land surface states and fluxes. The L4_SM algorithm is based on the assimilation of SMAP radiometer brightness temperature observations into the NASA Catchment land surface model using a spatially distributed ensemble Kalman filter (EnKF). In 2018, the L4_SM algorithm was upgraded from Version 3 to Version 4. Underlying the new version is a revised modeling system that includes improved input parameter datasets for land cover, topography, and vegetation height that are based on recent, high-quality, space-borne remote sensing observations. Additionally, SMAP Level-2 soil moisture retrievals and in situ soil moisture measurements were used to calibrate a particular Catchment model parameter that governs the recharge of surface soil moisture from below under non-equilibrium conditions, which brings the model's surface soil moisture more in line with the SMAP Level-2 and in situ soil moisture. Moreover, the calibration of the assimilated SMAP brightness temperatures changed substantially from Version 3 to Version 4, and the "catchment deficit" model variable was removed from the EnKF state vector to avoid degrading the model's groundwater estimates.Considerable effort went into the version upgrade, creating an expectation that the new version is improved over the old version. Indeed, some aspects of the new version are clearly better. However, other aspects are not. In this presentation we summarize the skill of the new and old versions vs. independent in situ measurements and in terms of data assimilation diagnostics, including, for example, the statistics of the (soil moisture) analysis increments and the observation-minus-forecast (brightness temperatures) residuals. We share our experience with trying to improve to the L4_SM product and the lessons learned from the effort.

Reichle, Rolf H.↗

Impact of AIRS Thermodynamic Profiles on Precipitation Forecasts for Atmospheric River Cases Affecting the Western United States

This project is a collaborative activity between the NASA Short-term Prediction Research and Transition (SPoRT) Center and the NOAA Hydrometeorology Testbed (HMT) to evaluate a SPoRT Advanced Infrared Sounding Radiometer (AIRS: Aumann et al. 2003) enhanced moisture analysis product. We test the impact of assimilating AIRS temperature and humidity profiles above clouds and in partly cloudy regions, using the three-dimensional variational Gridpoint Statistical Interpolation (GSI) data assimilation (DA) system (Developmental Testbed Center 2012) to produce a new analysis. Forecasts of the Weather Research and Forecasting (WRF) model initialized from the new analysis are compared to control forecasts without the additional AIRS data. We focus on some cases where atmospheric rivers caused heavy precipitation on the US West Coast. We verify the forecasts by comparison with dropsondes and the Cooperative Institute for Research in the Atmosphere (CIRA) Blended Total Precipitable Water product.

Zavodsky, Bradley T.↗

Ensemble PM2.5 Forecasting during the 2018 Camp Fire Event Using the HYSPLIT Transport and Dispersion Model

Biomass burning releases a vast amount of aerosols into the atmosphere, often leading to severe air quality and health problems. Prediction of the air quality effects from biomass burning emissions is challenging due to uncertainties in fire emission, plume rise calculation, and other model inputs/processes. Ensemble forecasting is increasingly used to represent model uncertainties. In this paper, an ensemble forecast was conducted to predict surface PM2.5 during the 2018 California Camp Fire event using the NOAA HYSPLIT dispersion model at 0.1- degree horizontal resolution. Different combinations of four satellite-based fire emission datasets (FEER, FLAMBE, GBBEPx and GFAS), two plume rise schemes (Briggs and Sofiev), various meteorology inputs and model setup options were used to create the forecast ensemble, for a total of 112 experiments. The performance of each ensemble member and the ensemble mean were evaluated using ground-based observations, with four statistical metrics and an overall rank. The ensemble spread of the 112 members reached 1000 μg/m3, highlighting the large uncertainty in wildfire forecast. The ensemble mean displayed the best performance. Each fire emission product contributed to one or more members among the top ten performers, revealing the forecasting dependence on both the quality of fire emissions data and model representation of emission, transport and removal processes. In addition, an ensemble size reduction technique was introduced. With the help of this technique, the ensemble size was reduced from 112 to 28 members and still produced an ensemble mean that yielded comparable or even better performance to that of the full ensemble.

Biomass Burning↗

Effects of Lightning and Other Meteorological Factors on Fire Activity in the North American Boreal Forest: Implications for Fire Weather Forecasting

The effects of lightning and other meteorological factors on wildfire activity in the North American boreal forest are statistically analyzed during the fire seasons of 2000-2006 through an integration of the following data sets: the MODerate Resolution Imaging Spectroradiometer (MODIS) level 2 fire products, the 3-hourly 32-kin gridded meteorological data from North American Regional Reanalysis (NARR), and the lightning data collected by the Canadian Lightning Detection Network (CLDN) and the Alaska Lightning Detection Network (ALDN). Positive anomalies of the 500 hPa geopotential height field, convective available potential energy (CAPE), number of cloud-to-ground lightning strikes, and the number of consecutive dry days are found to be statistically important to the seasonal variation of MODIS fire counts in a large portion of Canada and the entirety of Alaska. Analysis of fire occurrence patterns in the eastern and western boreal forest regions shows that dry (in the absence of precipitation) lightning strikes account for only 20% of the total lightning strikes, but are associated with (and likely cause) 40% of the MODIS observed fire counts in these regions. The chance for ignition increases when a threshold of at least 10 dry strikes per NARR grid box and at least 10 consecutive dry days is reached. Due to the orientation of the large-scale pattern, complex differences in fire and lightning occurrence and variability were also found between the eastern and western sub-regions. Locations with a high percentage of dry strikes commonly experience an increased number of fire counts, but the mean number of fire counts per dry strike is more than 50% higher in western boreal forest sub-region, suggesting a geographic and possible topographic influence. While wet lightning events are found to occur with a large range of CAPE values, a high probability for dry lightning occurs only when 500 hPa geopotential heights are above 5700m and CAPE values are near the maximum observed level, underscoring the importance of low-level instability to boreal fire weather forecasts-

Peterson, D.↗

Non-urban mobile radio market demand forecast

A national nonmetropolitan land mobile traffic model for 1990-2000 addresses user classes, density classes, traffic mix statistics, distance distribution, geographic distribution, price elasticity, and service quality elasticity. Traffic demands for business, special industrial, and police were determined on the basis of surveys in 73 randomly selected nonurban counties. The selected services represent 69% of total demand. The results were extrapolated to all services in the non-SMSA areas of the contiguous United States. Radiotelephone services were considered separately. Total non-SMSA mobile radio demand (one way) estimates are given. General functional requirements include: hand portability, privacy, reduction of blind spots, two way data transmission, position location, slow scan imagery.

Castruccio, P. A.↗

Status of the NASA GMAO Observing System Simulation Experiment

An Observing System Simulation Experiment (OSSE) is a pure modeling study used when actual observations are too expensive or difficult to obtain. OSSEs are valuable tools for determining the potential impact of new observing systems on numerical weather forecasts and for evaluation of data assimilation systems (DAS). An OSSE has been developed at the NASA Global Modeling and Assimilation Office (GMAO, Errico et al 2013). The GMAO OSSE uses a 13-month integration of the European Centre for Medium- Range Weather Forecasts 2005 operational model at T511/L91 resolution for the Nature Run (NR). Synthetic observations have been updated so that they are based on real observations during the summer of 2013. The emulated observation types include AMSU-A, MHS, IASI, AIRS, and HIRS4 radiance data, GPS-RO, and conventional types including aircraft, rawinsonde, profiler, surface, and satellite winds. The synthetic satellite wind observations are colocated with the NR cloud fields, and the rawinsondes are advected during ascent using the NR wind fields. Data counts for the synthetic observations are matched as closely as possible to real data counts, as shown in Figure 2. Errors are added to the synthetic observations to emulate representativeness and instrument errors. The synthetic errors are calibrated so that the statistics of observation innovation and analysis increments in the OSSE are similar to the same statistics for assimilation of real observations, in an iterative method described by Errico et al (2013). The standard deviations of observation minus forecast (xo-H(xb)) are compared for the OSSE and real data in Figure 3. The synthetic errors include both random, uncorrelated errors, and an additional correlated error component for some observational types. Vertically correlated errors are included for conventional sounding data and GPS-RO, and channel correlated errors are introduced to AIRS and IASI (Figure 4). HIRS, AMSU-A, and MHS have a component of horizontally correlated error. The forecast model used by the GMAO OSSE is the Goddard Earth Observing System Model, Version 5 (GEOS-5) with Gridpoint Statistical Interpolation (GSI) DAS. The model version has been updated to v. 5.13.3, corresponding to the current operational model. Forecasts are run on a cube-sphere grid with 180 points along each edge of the cube (approximately 0.5 degree horizontal resolution) with 72 vertical levels. The DAS is cycled at 6-hour intervals, with 240 hour forecasts launched daily at 0000 UTC. Evaluation of the forecasting skill for July and August is currently underway. Prior versions of the GMAO OSSE have been found to have greater forecasting skill than real world forecasts. It is anticipated that similar forecast skill will be found in the updated OSSE.

OSSE↗

Mental Models of Software Forecasting

The majority of software engineers resist the use of the currently available cost models. One problem is that the mathematical and statistical models that are currently available do not correspond with the mental models of the software engineers. In an earlier JPL funded study (Hihn and Habib-agahi, 1991) it was found that software engineers prefer to use analogical or analogy-like techniques to derive size and cost estimates, whereas curren CER's hide any analogy in the regression equations. In addition, the currently available models depend upon information which is not available during early planning when the most important forecasts must be made.

Cost Models↗

Single-column Model Intercomparison for a Stably Stratified Atmospheric Boundary Layer

The parameterization of the stably stratified atmospheric boundary layer is a difficult issue, which has a large impact on the medium-range weather forecasts and on climate integrations. A non-strongly stratified arctic case is simulated by nineteen single-column turbulence schemes. The statistics from the Large-eddy simulation (LES) intercomparison made for the same case by eight different models are used as a guiding reference. The single-column parameterizations include research schemes and operational schemes from major forecast and climate research centres. First order schemes, a large number of turbulence kinetic energy closures, and other proposals have submitted results. There is a large spread in the results; in general, the operational schemes mix more efficiently than the research ones, and the TKE and other higher order closures give results closer to the LES statistics. The sensitivities of the schemes to the parameters of their turbulence closures are partially explored.

Cuxart, J.↗