Ensemble prediction for nowcasting with a convection-permitting model - II: forecast error statistics

A 24-member ensemble of 1-h high-resolution forecasts over the Southern United Kingdom is used to study short-range forecast error statistics. The initial conditions are found from perturbations from an ensemble transform Kalman filter. Forecasts from this system are assumed to lie within the bounds of forecast error of an operational forecast system. Although noisy, this system is capable of producing physically reasonable statistics which are analysed and compared to statistics implied from a variational assimilation system. The variances for temperature errors for instance show structures that reflect convective activity. Some variables, notably potential temperature and specific humidity perturbations, have autocorrelation functions that deviate from 3-D isotropy at the convective-scale (horizontal scales less than 10 km). Other variables, notably the velocity potential for horizontal divergence perturbations, maintain 3-D isotropy at all scales. Geostrophic and hydrostatic balances are studied by examining correlations between terms in the divergence and vertical momentum equations respectively. Both balances are found to decay as the horizontal scale decreases. It is estimated that geostrophic balance becomes less important at scales smaller than 75 km, and hydrostatic balance becomes less important at scales smaller than 35 km, although more work is required to validate these findings. The implications of these results for high-resolution data assimilation are discussed.

Application of data assimilation to ocean and climate prediction

Ocean prediction systems are now able to analyse and predict temperature, salinity and velocity structures within the ocean by assimilating measurements of the ocean’s temperature and salinity into physically based ocean models. Data assimilation combines current estimates of state variables, such as temperature and salinity, from a computational model with measurements of the ocean and atmosphere in order to improve forecasts and reduce uncertainty in the forecast accuracy. Data assimilation generally works well with ocean models away from the equator but has been found to induce vigorous and unrealistic overturning circulations near the equator. A pressure correction method was developed at the University of Reading and the Met Office to control these circulations using ideas from control theory and an understanding of equatorial dynamics. The method has been used for the last 10 years in seasonal forecasting and ocean prediction systems at the Met Office and European Center for Medium-range Weather Forecasting (ECMWF). It has been an important element in recent re-analyses of the ocean heat uptake that mitigates climate change.

Numerical prediction of three-dimensional time-dependent viscoelastic extrudate swell using differential and algebraic models

This study investigates the numerical simulation of three-dimensional time-dependent viscoelastic free surface flows using the Upper-Convected Maxwell (UCM) constitutive equation and an algebraic explicit model. This investigation was carried out to develop a simplified approach that can be applied to the extrudate swell problem. The relevant physics of this flow phenomenon is discussed in the paper and an algebraic model to predict the extrudate swell problem is presented. It is based on an explicit algebraic representation of the non-Newtonian extra-stress through a kinematic tensor formed with the scaled dyadic product of the velocity field. The elasticity of the fluid is governed by a single transport equation for a scalar quantity which has dimension of strain rate. Mass and momentum conservations, and the constitutive equation (UCM and algebraic model) were solved by a three-dimensional time-dependent finite difference method. The free surface of the fluid was modeled using a marker-and-cell approach. The algebraic model was validated by comparing the numerical predictions with analytic solutions for pipe flow. In comparison with the classical UCM model, one advantage of this approach is that computational workload is substantially reduced: the UCM model employs six differential equations while the algebraic model uses only one. The results showed stable flows with very large extrudate growths beyond those usually obtained with standard differential viscoelastic models. (C) 2010 Elsevier Ltd. All rights reserved.

Analytical attitude prediction of spin stabilized spacecrafts perturbed by magnetic residual torque

An analytical approach for spin-stabilized spacecraft attitude prediction is presented for the influence of the residual magnetic torques. Assuming an inclined dipole model for the Earth's magnetic field, an analytical averaging method is applied to obtain the mean residual torque every orbital period. The orbit mean anomaly is utilized to compute the average components of residual torque in the spacecraft body frame reference system. The theory is developed for time variations in the orbital elements, and non-circular orbits, giving rise to many curvature integrals. It is observed that the residual magnetic torque does not have component along the spin axis. The inclusion of this torque on the rotational motion differential equations of a spin stabilized spacecraft yields conditions to derive an analytical solution. The solution shows that residual torque does not affect the spin velocity magnitude, contributing only for the precession and the drift of the spin axis of the spacecraft. (c) 2005 COSPAR. Published by Elsevier Ltd. All rights reserved.

Relationships and significance of lactate minimum, critical velocity, heart rate deflection and 3 000 m track-tests for running

The running velocities associated to lactate minimum (V-lm), heart rate deflection (V-HRd), critical velocity (CV), 3000 M (V-3000) and 10000 m performance (V-10km) were compared. Additionally the ability of V-lm and VHRd on identifying sustainable velocities was investigated.Methods. Twenty runners (28.5 +/- 5.9 y) performed 1) 3000 m running test for V3000; 2) an all-out 500 in sprint followed by 6x800 m incremental bouts with blood lactate ([lac]) measurements for V-lm; 3) a continuous velocity-incremented test with heart rate measurements at each 200 m for V-HRd; 4) participants attempted to 30 min of endurance test both at V-lm(ETVlm) and V-HRd(ETVHRd). Additionally, the distance-time and velocity-1/time relationships produced CV by 2 (500 m and 3000 m) or 3 predictive trials (500 m, 3000 m and distance reached before exhaustion during ETVHRd), and a 10 km race was recorded for V-10km.Results. The CV identified by different methods did not differ to each other. The results (m(.)min(-1)) revealed that V-.(lm) (281 +/- 14.8)< CV (292.1 +/- 17.5)=V-10km (291.7 +/- 19.3)< V-HRd (300.8 +/- 18.7)=V-3000 (304 +/- 17.5) with high correlation among parameters (P < 0.001). During ETVlm participants completed 30 min of running while on the ETVHRd they lasted only 12.5 +/- 8.2 min with increasing [lac].Conclusion. We evidenced that CV and Vim track-protocols are valid for running evaluation and performance prediction and the parameters studied have different significance. The V-lm reflects the moderate-high intensity domain (below CV), can be sustained without [lac] accumulation and may be used for long-term exercise while the V-HRd overestimates a running intensity that can be sustained for long-time. Additionally, V-3000 and V-HRd reflect the severe intensity domain (above CV).

Spin stabilized satellite's attitude analytical prediction

An analytical approach for spin stabilized attitude propagation is presented, considering the coupled effect of the aerodynamic torque and the gravity gradient torque. A spherical coordination system fixed in the satellite is used to locate the satellite spin axis in relation to the terrestrial equatorial system. The spin axis direction is specified by its right ascension and the declination angles and the equation of motion are described by these two angles and the magnitude of the spin velocity. An analytical averaging method is applied to obtain the mean torques over an orbital period. To compute the average components of both aerodynamic torque and the gravity gradient torque in the satellite body frame reference system, an average time in the fast varying orbit element, the mean anomaly, is utilized. Afterwards, the inclusion of such torques on the rotational motion differential equations of spin stabilized satellites yields conditions to derive an analytical solution. The pointing deviation evolution, that is, the deviation between the actual spin axis and the computed spin axis, is also availed. In order to validate the analytical approach, the theory developed has been applied for spin stabilized Brazilian satellite SCD1, which are quite appropriated for verification and comparison of the data generated and processed by the Satellite Control Center of the Brazil National Research Institute (INPE). Numerical simulations performed with data of Brazilian Satellite SCD1 show the period that the analytical solution can be used to the attitude propagation, within the dispersion range of the attitude determination system performance of Satellite Control Center of the Brazilian Research Institute.

On the prediction of pickup and saltation velocities in pneumatic conveying

This paper presents a comparative study of the critical pickup and saltation velocities of particles in horizontal pipelines for pneumatic conveying design. A comparative study is performed using different existing correlations in the literature for the determination of the minimum velocity of transport as a function of the particle and pipe diameter, particle density, solid mass flow rate and particle sphericity. Their limitations and difficulties in predicting those critical velocities are analyzed. For the pickup velocity, an experimental study was also carried out in order to support the analysis. Recommendations are presented on the use of such correlations.

Prediction of the combustion process in fluidized bed based on physical-chemical properties of biomass particles and their hydrodynamic behaviors

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Prediction of the Temperature Distribution of Partially Submersed Umbilical Cables

The objective of this work is to predict the temperature distribution of partially submersed umbilical cables under different operating and environmental conditions. The commercial code Fluent (R) was used to simulate the heat transfer and the air fluid flow of part of a vertical umbilical cable near the air-water interface. A free-convective three-dimensional turbulent flow in open-ended vertical annuli was solved. The influence of parameters such as the heat dissipating rate, wind velocity, air temperature and solar radiation was analyzed. The influence of the presence of a radiation shield consisting of a partially submersed cylindrical steel tube was also considered. The air flow and the buoyancy-driven convective heat transfer in the annular region between the steel tube and the umbilical cable were calculated using the standard k-epsilon turbulence model. The radiative heat transfer between the umbilical external surface and the radiation shield was calculated using the Discrete Ordinates model. The results indicate that the influence of a hot environment and intense solar radiation may affect the umbilical cable performance in its dry portion.

Prediction of the rate of decline in fetal hemoglobin levels between first and second transfusions in red cell alloimmune disease

Objective To determine variables that predict the rate of decline in fetal hemoglobin levels in alloimmune disease. Method Retrospective review of singleton pregnancies that underwent first and second intrauterine transfusions for treatment of fetal anemia because of maternal Rh alloimmunization in a tertiary referral center. Results Forty-one first intrauterine transfusions were performed at 26.1?weeks (standard deviation, SD, 4.6), mean volume of blood transfused was 44.4?mL (SD 23.5) and estimated feto-placental volume expansion was 51.3% (SD 14.5%). Between first and second transfusion, hemoglobin levels reduced on average 0.40?g/dl/day (SD 0.25). Stepwise multiple regression analysis demonstrated that this rate significantly correlated with hemoglobin levels after the first transfusion, the interval between both procedures, and middle cerebral artery systolic velocity before the second transfusion. Conclusion The rate of decline in fetal hemoglobin levels between first and second transfusions in alloimmune disease can be predicted by a combination of hemoglobin levels after the first transfusion, interval between both procedures, and middle cerebral artery systolic velocity before the second transfusion. (C) 2012 John Wiley & Sons, Ltd.

New Equation for Prediction of Reverse Remodeling after Cardiac Resynchronization Therapy

Objectives: To integrate data from two-dimensional echocardiography (2D ECHO), three-dimensional echocardiography (3D ECHO), and tissue Doppler imaging (TDI) for prediction of left ventricular (LV) reverse remodeling (LVRR) after cardiac resynchronization therapy (CRT). It was also compared the evaluation of cardiac dyssynchrony by TDI and 3D ECHO. Methods: Twenty-four consecutive patients with heart failure, sinus rhythm, QRS = 120 msec, functional class III or IV and LV ejection fraction (LVEF) = 0.35 underwent CRT. 2D ECHO, 3D ECHO with systolic dyssynchrony index (SDI) analysis, and TDI were performed before, 3 and 6 months after CRT. Cardiac dyssynchrony analyses by TDI and SDI were compared with the Pearson's correlation test. Before CRT, a univariate analysis of baseline characteristics was performed for the construction of a logistic regression model to identify the best predictors of LVRR. Results: After 3 months of CRT, there was a moderate correlation between TDI and SDI (r = 0.52). At other time points, there was no strong correlation. Nine of twenty-four (38%) patients presented with LVRR 6 months after CRT. After logistic regression analysis, SDI (SDI > 11%) was the only independent factor in the prediction of LVRR 6 months of CRT (sensitivity = 0.89 and specificity = 0.73). After construction of receiver operator characteristic (ROC) curves, an equation was established to predict LVRR: LVRR =-0.4LVDD (mm) + 0.5LVEF (%) + 1.1SDI (%), with responders presenting values >0 (sensitivity = 0.67 and specificity = 0.87). Conclusions: In this study, there was no strong correlation between TDI and SDI. An equation is proposed for the prediction of LVRR after CRT. Although larger trials are needed to validate these findings, this equation may be useful to candidates for CRT. (Echocardiography 2012;29:678-687)

Prediction of the temperature distribution of partially submersed umbilical cables

The objective of this work is to predict the temperature distribution of partially submersed umbilical cables under different operating and environmental conditions. The commercial code Fluent® was used to simulate the heat transfer and the air fluid flow of part of a vertical umbilical cable near the air-water interface. A free-convective three-dimensional turbulent flow in open-ended vertical annuli was solved. The influence of parameters such as the heat dissipating rate, wind velocity, air temperature and solar radiation was analyzed. The influence of the presence of a radiation shield consisting of a partially submersed cylindrical steel tube was also considered. The air flow and the buoyancydriven convective heat transfer in the annular region between the steel tube and the umbilical cable were calculated using the standard k-ε turbulence model. The radiative heat transfer between the umbilical external surface and the radiation shield was calculated using the Discrete Ordinates model. The results indicate that the influence of a hot environment and intense solar radiation may affect the umbilical cable performance in its dry portion.

Prediction of sediments toxicity using consensus-based sediment quality guidelines and benthic foraminifera in the Ria de Aveiro(Portugal)

Ria de Aveiro is a large and shallow lagoon on the west coast of Portugal (40º38’N, 8º45´W), characterized by a complex geometry. It includes large areas of intertidal flats and a network of narrow channels which are connected to the Atlantic by an artificial inlet. Tides are the main forcing of the hydrology and physical processes of the lagoon. The deeper areas near the inlet are characterized by strong marine influence through tidal inflow, with high values of current velocity (>1m/s) and tidal range (2–3 m at spring tides), while in remote shallow areas, the circulation and the sea water inflow are reduced. These remote areas are more influenced by fresh waters received from several rivers and several small streams. The Aveiro lagoon is a very important ecosystem but as been used as recipient for various kinds of anthropogenic wastes resulting from the high population density, urban activities and industrial development. One of the most important Portuguese industrial centre is located in the lagoon margins. Ria de Aveiro is a coastal lagoon under huge direct antropization. This system also suffers strong diffuse antropization. This work is related with diffuse antropization linked with chemical pollution which may lead to biological stress and collapse.

Prediction of near field overpressure from quarry blasting

This paper investigates the propagation of airblast or pressure waves in air produced by bench blasting (i.e. detonation of the explosive in a row of blastholes, breaking the burden of rock towards the free vertical face of the block). Peak overpressure is calculated as a function of blasting parameters (explosive mass per delay and velocity at which the detonation sequence proceeds along the bench) and the polar coordinates of the position of interest (distance to the source and azimuth with respect to the free face). The model has been fitted to empirical data using linear least squares. The data set is composed of 122 airblast records monitored at distances less than 400 m in 41 production blasts carried out in two quarries. The model is statistically significant and has a determination coefficient of 0.87. The formula is validated from 12 airblast measurements gathered in five additional blasts.

Scale dependence in earthquake phenomena and its relevance to earthquake prediction.

The recent discovery of a low-velocity, low-Q zone with a width of 50-200 m reaching to the top of the ductile part of the crust, by observations on seismic guided waves trapped in the fault zone of the Landers earthquake of 1992, and its identification with the shear zone inferred from the distribution of tension cracks observed on the surface support the existence of a characteristic scale length of the order of 100 m affecting various earthquake phenomena in southern California, as evidenced earlier by the kink in the magnitude-frequency relation at about M3, the constant corner frequency for earthquakes with M below about 3, and the sourcecontrolled fmax of 5-10 Hz for major earthquakes. The temporal correlation between coda Q-1 and the fractional rate of occurrence of earthquakes in the magnitude range 3-3.5, the geographical similarity of coda Q-1 and seismic velocity at a depth of 20 km, and the simultaneous change of coda Q-1 and conductivity at the lower crust support the hypotheses that coda Q-1 may represent the activity of creep fracture in the ductile part of the lithosphere occurring over cracks with a characteristic size of the order of 100 m. The existence of such a characteristic scale length cannot be consistent with the overall self-similarity of earthquakes unless we postulate a discrete hierarchy of such characteristic scale lengths. The discrete hierarchy of characteristic scale lengths is consistent with recently observed logarithmic periodicity in precursory seismicity.