Numerical procedure for estimating temperature in solarized soils

The objective of this work was to develop a simplified numerical procedure for the estimation of accumulated monthly hours of solarized soil temperatures. The proposed model requires monthly means of daily solar radiation and maximum air temperature as input data, and a daily pattern of temperature variation assumed to be sine-shaped. The procedure was verified using observations made during the years 1992 and 1993 in Jaguariúna, SP. The proposed procedure can predict monthly temperature hours at 10 cm depth in the solarized soil, with acceptable accuracy, in the region for which it was developed.

Direct numerical simulations of interface dynamics to link capillary pressure and total surface energy

The flow of two immiscible fluids through a porous medium depends on the complex interplay between gravity, capillarity, and viscous forces. The interaction between these forces and the geometry of the medium gives rise to a variety of complex flow regimes that are difficult to describe using continuum models. Although a number of pore-scale models have been employed, a careful investigation of the macroscopic effects of pore-scale processes requires methods based on conservation principles in order to reduce the number of modeling assumptions. In this work we perform direct numerical simulations of drainage by solving Navier-Stokes equations in the pore space and employing the Volume Of Fluid (VOF) method to track the evolution of the fluid-fluid interface. After demonstrating that the method is able to deal with large viscosity contrasts and model the transition from stable flow to viscous fingering, we focus on the macroscopic capillary pressure and we compare different definitions of this quantity under quasi-static and dynamic conditions. We show that the difference between the intrinsic phase-average pressures, which is commonly used as definition of Darcy-scale capillary pressure, is subject to several limitations and it is not accurate in presence of viscous effects or trapping. In contrast, a definition based on the variation of the total surface energy provides an accurate estimate of the macroscopic capillary pressure. This definition, which links the capillary pressure to its physical origin, allows a better separation of viscous effects and does not depend on the presence of trapped fluid clusters.

Quantitative comparison between simulations of seismic wave propagation in heterogeneous poro-elastic media and equivalent visco-elastic solids for marine-type environments

There is increasing evidence to suggest that the presence of mesoscopic heterogeneities constitutes an important seismic attenuation mechanism in porous rocks. As a consequence, centimetre-scale perturbations of the rock physical properties should be taken into account for seismic modelling whenever detailed and accurate responses of specific target structures are desired, which is, however, computationally prohibitive. A convenient way to circumvent this problem is to use an upscaling procedure to replace each of the heterogeneous porous media composing the geological model by corresponding equivalent visco-elastic solids and to solve the visco-elastic equations of motion for the inferred equivalent model. While the overall qualitative validity of this procedure is well established, there are as of yet no quantitative analyses regarding the equivalence of the seismograms resulting from the original poro-elastic and the corresponding upscaled visco-elastic models. To address this issue, we compare poro-elastic and visco-elastic solutions for a range of marine-type models of increasing complexity. We found that despite the identical dispersion and attenuation behaviour of the heterogeneous poro-elastic and the equivalent visco-elastic media, the seismograms may differ substantially due to diverging boundary conditions, where there exist additional options for the poro-elastic case. In particular, we observe that at the fluid/porous-solid interface, the poro- and visco-elastic seismograms agree for closed-pore boundary conditions, but differ significantly for open-pore boundary conditions. This is an important result which has potentially far-reaching implications for wave-equation-based algorithms in exploration geophysics involving fluid/porous-solid interfaces, such as, for example, wavefield decomposition.

Insights into the Design and Characteristics of the Sedimentation Process at Multi-Barrel Culverts in Iowa, TR-596, 2010

The present study is an integral part of a broader study focused on the design and implementation of self-cleaning culverts, i.e., configurations that prevent the formation of sediment deposits after culvert construction or cleaning. Sediment deposition at culverts is influenced by many factors, including the size and characteristics of material of which the channel is composed, the hydraulic characteristics generated under different hydrology events, the culvert geometry design, channel transition design, and the vegetation around the channel. The multitude of combinations produced by this set of variables makes the investigation of practical situations a complex undertaking. In addition to the considerations above, the field and analytical observations have revealed flow complexities affecting the flow and sediment transport through culverts that further increase the dimensions of the investigation. The flow complexities investigated in this study entail: flow non-uniformity in the areas of transition to and from the culvert, flow unsteadiness due to the flood wave propagation through the channel, and the asynchronous correlation between the flow and sediment hydrographs resulting from storm events. To date, the literature contains no systematic studies on sediment transport through multi-box culverts or investigations on the adverse effects of sediment deposition at culverts. Moreover, there is limited knowledge about the non-uniform, unsteady sediment transport in channels of variable geometry. Furthermore, there are few readily useable (inexpensive and practical) numerical models that can reliably simulate flow and sediment transport in such complex situations. Given the current state of knowledge, the main goal of the present study is to investigate the above flow complexities in order to provide the needed insights for a series of ongoing culvert studies. The research was phased so that field observations were conducted first to understand the culvert behavior in Iowa landscape. Modeling through complementary hydraulic model and numerical experiments was subsequently carried out to gain the practical knowledge for the development of the self-cleaning culvert designs.

Clinical management and burden of bipolar disorder: a multinational longitudinal study (WAVE-bd Study)

BACKGROUND: Studies in bipolar disorder (BD) to date are limited in their ability to provide a whole-disease perspective--their scope has generally been confined to a single disease phase and/or a specific treatment. Moreover, most clinical trials have focused on the manic phase of disease, and not on depression, which is associated with the greatest disease burden. There are few longitudinal studies covering both types of patients with BD (I and II) and the whole course of the disease, regardless of patients' symptomatology. Therefore, the Wide AmbispectiVE study of the clinical management and burden of Bipolar Disorder (WAVE-bd) (NCT01062607) aims to provide reliable information on the management of patients with BD in daily clinical practice. It also seeks to determine factors influencing clinical outcomes and resource use in relation to the management of BD. METHODS: WAVE-bd is a multinational, multicentre, non-interventional, longitudinal study. Approximately 3000 patients diagnosed with BD type I or II with at least one mood event in the preceding 12 months were recruited at centres in Austria, Belgium, Brazil, France, Germany, Portugal, Romania, Turkey, Ukraine and Venezuela. Site selection methodology aimed to provide a balanced cross-section of patients cared for by different types of providers of medical aid (e.g. academic hospitals, private practices) in each country. Target recruitment percentages were derived either from scientific publications or from expert panels in each participating country. The minimum follow-up period will be 12 months, with a maximum of 27 months, taking into account the retrospective and the prospective parts of the study. Data on demographics, diagnosis, medical history, clinical management, clinical and functional outcomes (CGI-BP and FAST scales), adherence to treatment (DAI-10 scale and Medication Possession Ratio), quality of life (EQ-5D scale), healthcare resources, and caregiver burden (BAS scale) will be collected. Descriptive analysis with common statistics will be performed. DISCUSSION: This study will provide detailed descriptions of the management of BD in different countries, particularly in terms of clinical outcomes and resources used. Thus, it should provide psychiatrists with reliable and up-to-date information about those factors associated with different management patterns of BD. TRIAL REGISTRATION NO: ClinicalTrials.gov: NCT01062607.

Comparison of Setting Time Measured Using Ultrasonic Wave Propagation with Saw-Cutting Times on Pavements, TR-675, 2015

At present, there is little fundamental guidance available to assist contractors in choosing when to schedule saw cuts on joints. To conduct pavement finishing and sawing activities effectively, however, contractors need to know when a concrete mixture is going to reach initial set, or when the sawing window will open. Previous research investigated the use of the ultrasonic pulse velocity (UPV) method to predict the saw-cutting window for early entry sawing. The results indicated that the method has the potential to provide effective guidance to contractors as to when to conduct early entry sawing. The aim of this project was to conduct similar work to observe the correlation between initial setting and conventional sawing time. Sixteen construction sites were visited in Minnesota and Missouri over a two-year period. At each site, initial set was determined using a p-wave propagation technique with a commercial device. Calorimetric data were collected using a commercial semi-adiabatic device at a majority of the sites. Concrete samples were collected in front of the paver and tested using both methods with equipment that was set up next to the pavement during paving. The data collected revealed that the UPV method looks promising for early entry and conventional sawing in the field, both early entry and conventional sawing times can be predicted for the range of mixtures tested.

The effect of workplace smoking bans on heart rate variability and pulse wave velocity of non-smoking hospitality workers

OBJECTIVES: To investigate the effect of a change in second-hand smoke (SHS) exposure on heart rate variability (HRV) and pulse wave velocity (PWV), this study utilized a quasi-experimental setting when a smoking ban was introduced. METHODS: HRV, a quantitative marker of autonomic activity of the nervous system, and PWV, a marker of arterial stiffness, were measured in 55 non-smoking hospitality workers before and 3-12 months after a smoking ban and compared to a control group that did not experience an exposure change. SHS exposure was determined with a nicotine-specific badge and expressed as inhaled cigarette equivalents per day (CE/d). RESULTS: PWV and HRV parameters significantly changed in a dose-dependent manner in the intervention group as compared to the control group. A one CE/d decrease was associated with a 2.3 % (95 % CI 0.2-4.4; p = 0.031) higher root mean square of successive differences (RMSSD), a 5.7 % (95 % CI 0.9-10.2; p = 0.02) higher high-frequency component and a 0.72 % (95 % CI 0.40-1.05; p < 0.001) lower PWV. CONCLUSIONS: PWV and HRV significantly improved after introducing smoke-free workplaces indicating a decreased cardiovascular risk.

IOWATER Sustaining the Wave

Summary of the IOWATER Program and workshops offered.

Simulation of surface waves in porous media

We present a novel numerical algorithm for the simulation of seismic wave propagation in porous media, which is particularly suitable for the accurate modelling of surface wave-type phenomena. The differential equations of motion are based on Biot's theory of poro-elasticity and solved with a pseudospectral approach using Fourier and Chebyshev methods to compute the spatial derivatives along the horizontal and vertical directions, respectively. The time solver is a splitting algorithm that accounts for the stiffness of the differential equations. Due to the Chebyshev operator the grid spacing in the vertical direction is non-uniform and characterized by a denser spatial sampling in the vicinity of interfaces, which allows for a numerically stable and accurate evaluation of higher order surface wave modes. We stretch the grid in the vertical direction to increase the minimum grid spacing and reduce the computational cost. The free-surface boundary conditions are implemented with a characteristics approach, where the characteristic variables are evaluated at zero viscosity. The same procedure is used to model seismic wave propagation at the interface between a fluid and porous medium. In this case, each medium is represented by a different grid and the two grids are combined through a domain-decomposition method. This wavefield decomposition method accounts for the discontinuity of variables and is crucial for an accurate interface treatment. We simulate seismic wave propagation with open-pore and sealed-pore boundary conditions and verify the validity and accuracy of the algorithm by comparing the numerical simulations to analytical solutions based on zero viscosity obtained with the Cagniard-de Hoop method. Finally, we illustrate the suitability of our algorithm for more complex models of porous media involving viscous pore fluids and strongly heterogeneous distributions of the elastic and hydraulic material properties.

Wave-height hazard analysis in Eastern Coast of Spain : Bayesian approach using generalized Pareto distribution

Standard practice of wave-height hazard analysis often pays little attention to the uncertainty of assessed return periods and occurrence probabilities. This fact favors the opinion that, when large events happen, the hazard assessment should change accordingly. However, uncertainty of the hazard estimates is normally able to hide the effect of those large events. This is illustrated using data from the Mediterranean coast of Spain, where the last years have been extremely disastrous. Thus, it is possible to compare the hazard assessment based on data previous to those years with the analysis including them. With our approach, no significant change is detected when the statistical uncertainty is taken into account. The hazard analysis is carried out with a standard model. Time-occurrence of events is assumed Poisson distributed. The wave-height of each event is modelled as a random variable which upper tail follows a Generalized Pareto Distribution (GPD). Moreover, wave-heights are assumed independent from event to event and also independent of their occurrence in time. A threshold for excesses is assessed empirically. The other three parameters (Poisson rate, shape and scale parameters of GPD) are jointly estimated using Bayes' theorem. Prior distribution accounts for physical features of ocean waves in the Mediterranean sea and experience with these phenomena. Posterior distribution of the parameters allows to obtain posterior distributions of other derived parameters like occurrence probabilities and return periods. Predictives are also available. Computations are carried out using the program BGPE v2.0

Mesoscale numerical analysis of the historical November 1982 heavy precipitation event over Andorra (Eastern Pyrenees)

From 6 to 8 November 1982 one of the most catastrophic flash-flood events was recorded in the Eastern Pyrenees affecting Andorra and also France and Spain with rainfall accumulations exceeding 400 mm in 24 h, 44 fatalities and widespread damage. This paper aims to exhaustively document this heavy precipitation event and examines mesoscale simulations performed by the French Meso-NH non-hydrostatic atmospheric model. Large-scale simulations show the slow-evolving synoptic environment favourable for the development of a deep Atlantic cyclone which induced a strong southerly flow over the Eastern Pyrenees. From the evolution of the synoptic pattern four distinct phases have been identified during the event. The mesoscale analysis presents the second and the third phase as the most intense in terms of rainfall accumulations and highlights the interaction of the moist and conditionally unstable flows with the mountains. The presence of a SW low level jet (30 m s-1) around 1500 m also had a crucial role on focusing the precipitation over the exposed south slopes of the Eastern Pyrenees. Backward trajectories based on Eulerian on-line passive tracers indicate that the orographic uplift was the main forcing mechanism which triggered and maintained the precipitating systems more than 30 h over the Pyrenees. The moisture of the feeding flow mainly came from the Atlantic Ocean (7-9 g kg-1) and the role of the Mediterranean as a local moisture source was very limited (2-3 g kg-1) due to the high initial water vapour content of the parcels and the rapid passage over the basin along the Spanish Mediterranean coast (less than 12 h).

Rapid and sudden advection of warm and dry air in the Mediterranean Basin

Rapid advection of extremely warm and dry air is studied during two events in the Mediterranean Basin. On 27 August 2010 a rapid advection of extremely warm and dry air affected the northeast Iberian Peninsula during a few hours. At the Barcelona city center, the temperature reached 39.3 ° C, which is the maximum temperature value recorded during 230 yr of daily data series. On 23 March 2008 a rapid increase of temperature and drop of relative humidity were recorded for a few hours in Heraklion (Crete). During the morning on that day, the recorded temperature reached 34 °C for several hours on the northern coastline of this island.According to the World Meteorological Organization none of these events can be classified as a heat wave, which requires at least two days of abnormally high temperatures; neither are they a heat burst as defined by the American Meteorological Society, where abnormal temperatures take place during a few minutes. For this reason, we suggest naming this type of event flash heat. By using data from automatic weather stations in the Barcelona and Heraklion area and WRF mesoscale numerical simulations, these events are analyzed. Additionally, the primary risks and possible impacts on several fields are presented.

Seismic Wave Velocity as a Means of In-Place Density Measurement, Part 1, HR-114

Velocity-density tests conducted in the laboratory involved small 4-inch diameter by 4.58-inch-long compacted soil cylinders made up of 3 differing soil types and for varying degrees of density and moisture content, the latter being varied well beyond optimum moisture values. Seventeen specimens were tested, 9 with velocity determinations made along two elements of the cylinder, 180 degrees apart, and 8 along three elements, 120 degrees apart. Seismic energy was developed by blows of a small tack hammer on a 5/8-inch diameter steel ball placed at the center of the top of the cylinder, with the detector placed successively at four points spaced 1/2-inch apart on the side of the specimen involving wave travel paths varying from 3.36 inches to 4.66 inches in length. Time intervals were measured using a model 217 micro-seismic timer in both laboratory and field measurements. Forty blows of the hammer were required for each velocity determination, which amounted to 80 blows on 9 laboratory specimens and 120 blows on the remaining 8 cylinders. Thirty-five field tests were made over the three selected soil types, all fine-grained, using a 2-foot seismic line with hammer-impact points at 6-inch intervals. The small tack hammer and 5/8-inch steel ball was, again, used to develop seismic wave energy. Generally, the densities obtained from the velocity measurements were lower than those measured in the conventional field testing. Conclusions were reached that: (1) the method does not appear to be usable for measurement of density of essentially fine-grained soils when the moisture content greatly exceeds the optimum for compaction, and (2) due to a gradual reduction in velocity upon aging, apparently because of gradual absorption of pore water into the expandable interlayer region of the clay, the seismic test should be conducted immediately after soil compaction to obtain a meaningful velocity value.