Regional recharge assessment in the crystalline bedrock aquifer of the Kenogami Uplands,Canada

Abstract

Estimates of trends of climatic changes at basin and state scales are required for developing adaptation strategies related to planning, development and management of water resources. In the present study, seasonal and annual trends of changes in maximum temperature (T _max), minimum temperature (T _min), mean temperature (T _mean), temperature range (T _range), highest maximum temperature (H _max) and lowest minimum temperature (L _min) have been examined at the basin scale. The longest available records over the last century, for 43 stations covering nine river basins in northwest and central India, were used in the analysis. Of the nine river basins studied, seven showed a warming trend, whereas two showed a cooling trend. The Narmada and Sabarmati river basins experienced the maximum warming and cooling, respectively. The majority of basins in the study area show increasing trend in T _range, H _max and L _min. Seasonal analysis of different variables shows that the greatest changes in T _max and T _mean were observed in the post-monsoon season, while T _min experienced the greatest change in the monsoon season. This analysis provides scenarios of temperature changes which may be used for sensitivity analysis of water availability for different basins, and accordingly in planning and implementation of adaptation strategies.     

Correlation between clouds at different altitudes and solar activity: Fact or Artifact?

Studies of the relation between cosmic rays (CR) (solar activity) and atmospheric cloudiness are mostly based on the satellite ISCCP cloud data. However, doubts have been cast that these relations can be an artifact of instrumental effects, i.e., of the masking/obscuring low clouds by higher clouds in the satellite view. If this is the case, most of the earlier results based on ISCCP data would be devaluated. Here, we reanalyze the ISCCP cloud coverage data and its relation with the cosmic ray-induced ionization, and show that the correlation between low clouds and CR is affected by higher clouds in some geographical regions, but not everywhere. In turn, our results show that low clouds also may affect the relation of higher clouds with CR in some regions. Accordingly, correlation analysis can be performed only when the strong relation between clouds of different types is taken into account. In particular, studies based on global or latitudinal (zonally averaged) cloud data should be revised.     

Volcanic eruption clouds and the thermal power output of explosive eruptions

The maximum height attained by a volcanic eruption cloud is principally determined by the convective buoyancy of the mixture of volcanic gas + entrained air + fine-sized pyroclasts within the cloud. The thermal energy supplied to convection processes within an eruption cloud is derived from the cooling of pyroclastic material and volcanic gases discharged by an explosive eruption. Observational data from six recent eruptions indicates that the maximum height attained by volcanic eruption clouds is positively correlated with the rate at which pyroclastic material is produced by an explosive eruption (correlation coefficient r = + 0.97). The ascent of industrial hot gas plumes is also governed by the thermal convection process. Empirical scaling relationships between plume height and thermal flux have been developed for industrial plumes. Applying these scaling relationships to volcanic eruption clouds suggests that the rate at which thermal energy is released into the atmosphere by an explosive eruption increases in an approximately linear manner as an eruption's pyroclastic production rate increases.     

Effect of atmospheric gases,surface albedo and cloud overlap on the absorbed solar radiation

Recent studies have provided new evidence that models may systematically underestimate cloud solar absorption compared to observations. This study extends previous work on this “absorption anomaly” by using observational data together with solar radiative transfer parameterisations to calculate f_s (the ratio of surface and top of the atmosphere net cloud forcings) and its latitudinal variation for a range of cloud types. Principally, it is found that (a) the zonal mean behaviour of f_s varies substantially with cloud type, with the highest values obtained for low clouds; (b) gaseous absorption and scattering can radically alter the pattern of the variation of f_s with latitude, but gaseous effects cannot in general raise f_s to the level of around 1.5 as recently determined; (c) the importance of the gaseous contribution to the atmospheric ASR is such that whilst f_s rises with surface albedo, the net cloud contribution to the atmospheric ASR falls; (d) the assumed form of the degree of cloud overlap in the model can substantially affect the cloud contribution to the atmospheric ASR whilst leaving the parameter f_s largely unaffected; (e) even large uncertainties in the observed optical depths alone cannot account for discrepancies apparent between modelled and newly observed cloud solar absorption. It is concluded that the main source of the anomaly may derive from the considerable uncertainties regarding impure droplet microphysics rather than, or together with, uncertainties in macroscopic quantities. Further, variable surface albedos and gaseous effects may limit the use of contemporaneous satellite and ground-based measurements to infer the cloud solar absorption from the parameter f_s.     

Stochastic analysis of neutrally stratified cirrus layer

Two cirrus cloud systems observed during the winter of 2001 at the Southern Great Plains site of the Atmospheric Radiation Measurements program in Oklahoma, USA are studied because of the distinct neutrally stratified layers formed within the clouds. Observations are obtained with 35 GHz millimeter-wave radar and backscattering cross-section η(t) signals within radar-reflectivity restricted sublayers of the clouds are analyzed. The neutrally stratified layers of cirrus clouds are known to be associated with the existence of generating cells. The statistics of radiative properties within the neutrally stratified layers is obtained to be non-Gaussian and time-dependent. The purpose of this research is to derive a model of the cloud-generating cells layer in cirrus based on the statistics of observations. The Fokker–Planck equation approach provides suitable framework to treat non-Gaussian, time-dependent probability density functions (pdfs) such as those found for the η(t) signals. It is shown that the deviations from Gaussianity of radiative properties of the neutrally stratified generating cells layer in cirrus can be modeled by linear stochastically perturbed dynamics with multiplicative noise statistics. Because the multiplicative noise is often identified with state-dependent variations of stochastic feedbacks from unresolved system components it is expected that derived stochastic model will be useful for parameterization of cirrus in global circulation models (GCMs).     

Characteristics of Cloud Drop Spectra at Different Levels and with Respect to Cloud Thickness

--Measurements on drop size spectra were made in cumulus clouds over Pune (inland) region on many days during the summer monsoon seasons. In this paper, the measurements in non-raining cumulus clouds made in the years 1984, 1985 and 1986 at different levels and for different cloud thickness have been studied. In general, the drop size spectra broadened with height and the concentration of drops with diameter > 50 wm (N_L), mean volume diameter (MVD), liquid water content (LWC) and dispersion increased with height while the concentration of drops with diameter < 20 wm (N_S) and the total concentration of drops (N_T) decreased with height. The average drop size distributions were unimodal at the lower levels while they were bimodal at the higher levels. High water contents were confined to drops in the size range 5-25 wm at both higher and lower levels. The average drop size spectra were broader and N_L, LWC, MVD and dispersion greater while N_T and N_S smaller for thicker clouds (range of vertical extent 1.1-2.1 km) as compared to those for thinner clouds (range of vertical extent 0.3-1.1 km). Water contents for the drops > 28 wm were higher while those for the drops > 28 wm lower in thicker clouds than in thinner clouds. The average drop size distributions were bimodal in the former case, while they were unimodal in the other case.     

Graphical analysis of enthalpy-composition relationships in mixed magmas

Enthalpy-composition diagrams (H^*-X) calculated from existing temperature-composition (T-X) phase diagrams and thermodynamic data provide a simple and effective means for evaluating the enthalpy-temperature effects of magma mixing. If athermal mixing is assumed, adiabatic mixing lines on H^*-X diagrams are straight lines because enthalpy, unlike temperature, is an extensive property of a system. Comparison of binary T-X diagrams with their derivative H^*-X diagrams shows that incorrect predictions can be obtained when T-X diagrams are used to analyze mixing problems.An H^*-X diagram calculated using experimentally determined phase equilibria for anhydrous basalt-rhyolite mixtures predicts that adiabatic mixing of basalt and rhyolite at their respective liquidus temperatures will result in small amounts of crystallization ( < 2 wt.% ). Because the phase equilibria of hydrous basalt-rhyolite mixtures have not yet been determined, an H^*-X diagram for such mixtures cannot be constructed. However, existing hydrous phase equilibria can be used to predict whether adiabatic mixtures of anhydrous basalt and hydrous rhyolite will be super- or subliquidus. Calculations at P_total = 200 MPa show that on an H^*-X diagram a mixing chord drawn between anhydrous basalt and water-saturated rhyolite at their respective liquidus temperatures lies below the enthalpy values calculated for a Paricutin andesite and Mt. St. Helens dacite at their estimated liquidus temperatures. This indicates that the liquidus for mixtures of anhydrous basalt and water-saturated rhyolite is noticeably convex upward, suggesting that larger amounts of crystallization will occur than in the anhydrous case.     

Peak and integral seismic parameters of L��Aquila 2009 ground motions: observed versus code provision values

The M _w 6.3 L’Aquila earthquake of April 6, 2009 hit a wide area of the Abruzzo region (Central Italy). The epicentre of the main shock was very close to the urban centre of L’Aquila, the regional capital, with an epicentral distance less than 10 km. It was the strongest earthquake ever recorded in Italy which has provided ground motion recordings from accelerometric stations located in close proximity to the epicentre. Because of this, several remarkable results can be achieved by analysing the strong motion recorded signals in terms of peak (PGA, PGV and PGD) and integral (Housner Intensity, I _H) seismic parameters. Additionally, an alternative time-domain representation of recorded signals has been used to furnish a rapid comparison of traces recorded at different stations and along different directions. Some comparisons between the response spectra derived from the recordings and the elastic demand spectra provided in the new seismic Italian code have also been performed. PGA recorded values are very high and generally higher than code values for seismic actions with return period T _R = 475 years. In some cases, this also happens for seismic actions with T _R up to 2,475 years. With regard to I _H, recorded values are generally higher for T _R = 475 years, whilst they are remarkably lower for T _R = 2,475 years. Accurate analyses have been carried out in the article to better understand the above differences and their significance and implications.     

Predictors of the peak width for networks with exponential links

We investigate optimal predictors of the peak (S) and distance to peak (T) of the width function of drainage networks under the assumption that the networks are topologically random with independent and exponentially distributed link lengths. Analytical results are derived using the fact that, under these assumptions, the width function is a homogeneous Markov birth-death process. In particular, exact expressions are derived for the asymptotic conditional expectations ofS andT given network magnitudeN and given mainstream lengthH. In addition, a simulation study is performed to examine various predictors ofS andT, includingN, H, and basin morphometric properties; non-asymptotic conditional expectations and variances are estimated. The best single predictor ofS isN, ofT isH, and of the scaled peak (S divided by the area under the width function) isH. Finally, expressions tested on a set of drainage basins from the state of Wyoming perform reasonably well in predictingS andT despite probable violations of the original assumptions.     

Studying the Relationship between High-Latitude Geomagnetic Activity and Parameters of Interplanetary Magnetic Clouds with the Use of Artificial Neural Networks

The cause-and-effect relations of the dynamics of high-latitude geomagnetic activity (in terms of the AL index) and the type of the magnetic cloud of the solar wind are studied with the use of artificial neural networks. A recurrent neural network model has been created based on the search for the optimal physically coupled input and output parameters characterizing the action of a plasma flux belonging to a certain magnetic cloud type on the magnetosphere. It has been shown that, with IMF components as input parameters of neural networks with allowance for a 90-min prehistory, it is possible to retrieve the AL sequence with an accuracy to ~80%. The successful retrieval of the AL dynamics by the used data indicates the presence of a close nonlinear connection of the AL index with cloud parameters. The created neural network models can be applied with high efficiency to retrieve the AL index, both in periods of isolated magnetospheric substorms and in periods of the interaction between the Earth’s magnetosphere and magnetic clouds of different types. The developed model of AL index retrieval can be used to detect magnetic clouds.     

Relations of sand trapping efficiency and migration speed of transverse dunes to wind velocity

A self‐consistent model which describes transverse dune migration in equilibrium is introduced. It shows that an equilibrium expression for dune migration speed (c _d) must take into account sand trapping efficiency (T _E), and that T _E is strongly related to the wind speedup over the windward surface. An expression for sand trapping efficiency (T _E) is analytically derived from a microscale analysis of sand grain deposition on the slip face. Sand trapping efficiency (T _E) is mainly determined by shear velocity on a level surface (u_*(−∞)), and rapidly decreases as u_*(−∞) increases. For each dune height (H), dune migration speed (c _d) first increases, and then decreases monotonically after reaching the maximum, as the shear velocity on a level surface (u_*(−∞)) increases. Dune migration speed (c _d) is not inversely proportional to dune height (H). For low dunes, small sand trapping efficiency (T _E) suppresses c _d, whereas for high dunes, wind speedup and large T _E resist the decrease of c _d. Some field data show the same tendency. The dune‐to‐plane‐bed transition observed in subaqueous and venusian bedforms could be associated with the decrease of sand trapping efficiency (T _E). Copyright © 2000 John Wiley & Sons, Ltd.     

Reconciling pyroclastic flow and surge: the multiphase physics of pyroclastic density currents

Two end-member types of pyroclastic density current are commonly recognized: pyroclastic surges are dilute currents in which particles are carried in turbulent suspension and pyroclastic flows are highly concentrated flows. We provide scaling relations that unify these end-members and derive a segregation mechanism into basal concentrated flow and overriding dilute cloud based on the Stokes number (S_T), the stability factor (Σ_T) and the dense-dilute condition (D_D). We recognize five types of particle behaviors within a fluid eddy as a function of S_T and Σ_T: (1) particles sediment from the eddy, (2) particles are preferentially settled out during the downward motion of the eddy, but can be carried during its upward motion, (3) particles concentrate on the periphery of the eddy, (4) particles settling can be delayed or ‘fast-tracked’ as a function of the eddy spatial distribution, and (5) particles remain homogeneously distributed within the eddy. We extend these concepts to a fully turbulent flow by using a prototype of kinetic energy distribution within a full eddy spectrum and demonstrate that the presence of different particle sizes leads to the density stratification of the current. This stratification may favor particle interactions in the basal part of the flow and D_D determines whether the flow is dense or dilute. Using only intrinsic characteristics of the current, our model explains the discontinuous features between pyroclastic flows and surges while conserving the concept of a continuous spectrum of density currents.     

The impact of cloud inhomogeneities on the Earth radiation budget: the 14 October 1989 I.C.E. convective cloud case study

Through their multiple interactions with radiation, clouds have an important impact on the climate. Nonetheless, the simulation of clouds in climate models is still coarse. The present evolution of modeling tends to a more realistic representation of the liquid water content; thus the problem of its subgrid scale distribution is crucial. For a convective cloud field observed during ICE 89, Landsat TM data (resolution: 30m) have been analyzed in order to quantify the respective influences of both the horizontal distribution of liquid water content and cloud shape on the Earth radiation budget. The cloud field was found to be rather well-represented by a stochastic distribution of hemi-ellipsoidal clouds whose horizontal aspect ratio is close to 2 and whose vertical aspect ratio decreases as the cloud cell area increases. For that particular cloud field, neglecting the influence of the cloud shape leads to an over-estimate of the outgoing longwave flux; in the shortwave, it leads to an over-estimate of the reflected flux for high solar elevations but strongly depends on cloud cell orientations for low elevations. On the other hand, neglecting the influence of cloud size distribution leads to systematic over-estimate of their impact on the shortwave radiation whereas the effect is close to zero in the thermal range. The overall effect of the heterogeneities is estimated to be of the order of 10 W m⁻² for the conditions of that Landsat picture (solar zenith angle 65○, cloud cover 70%); it might reach 40 W m⁻² for an overhead sun and overcast cloud conditions.     

The dimensions and dynamics of volcanic eruption columns

Eruption columns can be divided into three regimes of physical behaviour. The basal gas thrust region is characterized by large velocities and decelerations and is dominated by momentum. This region is typically a few hundred metres in height and passes upwards into a much higher convective region where buoyancy is dominant. The top of the convective region is defined by the level of neutral density (heightH _B ) where the column has a bulk density equal to the surrounding atmosphere. Above this level the column continues to ascend to a heightH _T due to its momentum. The column spreads horizontally and radially outwards between heightH _T andH _B to form an umbrella cloud. Numerical calculations are presented on the shape of eruption columns and on the relationships between the heightH _B and the mass discharge rate of magma, magma temperature and atmospheric temperature gradients. Spreading rate of the column margins increases with height principally due to the decrease in the atmospheric pressure. The relationship between column height and mass discharge rate shows good agreement with observations. The temperature inversion above the tropopause is found to only have a small influence on column height and, eruptions with large discharge rates can inject material to substantially greater heights than the inversion level. Approximate calculations on the variation of convective velocities with height are consistent with field data and indicate that columns typically ascend at velocities from a few tens to over 200 m/s. In very large columns (greater than 30 km) the calculated convective velocities approach the speed of sound in air, suggesting that compressibility effects may become important in giant columns. Radial velocities in the umbrella region where the column is forced laterally into the atmosphere can be substantial and exceed 55 m/s in the case of the May 18th Mount St. Helens eruption. Calculations on motions in this region imply that it plays a major role in the transport of coarse pyroclastic fragments.     

Impacts of forest structure on precipitation interception and run‐off generation in a semiarid region in northern China

Water resource scarcity and uneven distribution are 2 major environmental issues in China today. Forest structure is a dominant factor that influences hydrological processes, but the specific interactions remain uncertain due to the predominant use of individual or 1‐dimensional forest structure metrics in previous studies. In this study, forest structures in 8 run‐off plots on Mount Miaofeng in north China were parameterized by metrics of different dimensionalities. The relation between canopy interception and forest structure, shrub/litter interception, and forest structure as well as run‐off and forest structure were analysed by regression method and validated by leave‐one‐out cross test. The results showed that canopy interception rates ranged from less than 0.10 all the way to 0.80, affected by forest structure and precipitation, with interception rate decreasing logarithmically as precipitation increased. Forests with a larger canopy area (CA), leaf area index (LAI), and higher average height (H) had a narrow range of canopy interception rates, and forest with larger value of diameters at breath height (DBH), H, LAI, vertical heterogeneity coefficient (T), and structure complexity index (SCI) had higher interception rates. Forests with higher value of DBH, H, and horizontal heterogeneity coefficient (R) had higher shrub/litter interception rates on the forest floor. The run‐off coefficient was only significantly associated with LAI, T, and SCI. The validation test indicated that regression analysis of canopy interception rates and shrub interception are reliable and SCI is a key factor to influence the run‐off coefficient. However, the regression results of litter interception have a relatively large error. According to the results, to reduce the risks of the landslides and floods, forest managers should complicate the canopy and preserve trees with thicker stems and larger canopies. By contrast, to obtain more water resource from run‐off in arid regions, forest managers should harvest trees with large canopies and construct complex vertical structures by intermediate cutting.     

Demographic parameters and mixis of three Brachionus angularis Gosse (Rotatoria) strains fed on different algae

Rotifers are one of the smallest metazoans, served as perfect material for evolution theories and excellent food resources to larva in aquaculture. Population dynamics of rotifers help to optimize culture conditions, and to understand selection of asexual and sexual reproduction types. We described population dynamics of three Brachionus angularis strains using the life table method, with three types of food. The following demographic parameters were analyzed: intrinsic rate of increase (r), net reproductive rate (R₀), generation time (T), average lifespan (L), and the total percentage of mictic females in the offspring. The results showed that there were significant effects of strain on the parameters except for R₀. However, the food type significantly affected r, R₀, T, and the total percentage of mictic females in the offspring. Demographic parameters differed among strains, showing the possible adaptation to local niches. In comparison to Chlorella pyrenoidosa, Scenedesmus obliquus served more nutrition to enlarge population size, and more mixis was induced in rotifer populations fed on it.     

Downdraft effects in a shallow convection adjustment

Summary The paper deals with the role of penetrative downdrafts in the convective adjustment suggested by Betts[1]. Three different types of downdraft sinking levels are used in the assessment of the reference profiles of temperature and humidity in cumulus cloud layers and these are compared with the profiles with no downdraft. It is shown that the depths of downdraft penetration may significantly influence the reference profiles mainly for clouds of vertical extents larger than one kilometre.     

PHYSICAL PROBLEMS OF WEATHER MODIFICATION / Des problèmes physiques de la modification artificielle du temps

Abstract

There are many contrasting views on the success of experiments into weather modification. his paper reviews a large number of these studies commencing with a summary of some of the pertinent problems in cloud physics and those leading to cloud modification. It deals with the techniques of modification such as by the use of dry ice and silver iodide and gives details of the results of experiments with clouds of different types. The conclusion reached is that the result of methods of modification used at present are not sufficiently predictable to warrant their use commercially and that the results of most of the interesting studies have not been evaluated by sufficiently reliable methods.     

Roundhouse (RND) Mountain Top Research Site: Measurements and Uncertainties for Winter Alpine Weather Conditions

The objective of this work is to better understand and summarize the mountain meteorological observations collected during the Science of Nowcasting Winter Weather for the Vancouver 2010 Olympics and Paralympics (SNOW-V10) project that was supported by the Fog Remote Sensing and Modeling (FRAM) project. The Roundhouse (RND) meteorological station was located 1,856 m above sea level that is subject to the winter extreme weather conditions. Below this site, there were three additional observation sites at 1,640, 1,320, and 774 m. These four stations provided some or all the following measurements at 1 min resolution: precipitation rate (PR) and amount, cloud/fog microphysics, 3D wind speed (horizontal wind speed, U _h; vertical air velocity, w _a), visibility (Vis), infrared (IR) and shortwave (SW) radiative fluxes, temperature (T) and relative humidity with respect to water (RH_w), and aerosol observations. In this work, comparisons are made to assess the uncertainties and variability for the measurements of Vis, RH_w, T, PR, and wind for various winter weather conditions. The ground-based cloud imaging probe (GCIP) measurements of snow particles using a profiling microwave radiometer (PMWR) data have also been shown to assess the icing conditions. Overall, the conclusions suggest that uncertainties in the measurements of Vis, PR, T, and RH can be as large as 50, >60, 50, and >20 %, respectively, and these numbers may increase depending on U _h, T, Vis, and PR magnitude. Variability of observations along the Whistler Mountain slope (~500 m) suggested that to verify the models, model space resolution should be better than 100 m and time scales better than 1 min. It is also concluded that differences between observed and model based parameters are strongly related to a model’s capability of accurate prediction of liquid water content (LWC), PR, and RH_w over complex topography.     

Characteristics of cloud drop spectra in tropical warm clouds

Characteristics of cloud drop spectra were studied using 400 samples obtained from 120 warm cumulus clouds formed during the summer monsoon season.The total concentration of cloud drops (N _T) varied from 384 to 884 cm^–3 and the maximum concentration was observed in the layer below the cloud-top. The width of the drop spectrum was broader in the cloud-base region and in the region below the cloud-top. The spectrum was multimodal at all levels except in the cloud-top region where it was unimodal. The concentration of drops with diameter greater than 50 m (N _L) varied from 0.0 to 0.674 cm^–3.N _L was larger in the cloud-base region.N _L decreased with height up to the middle level and thereafter showed an increase. In the cloud-top region no large drops were present. The computed values of the liquid water varied between 0.132 and 0.536 g m^–3 and the mean volume diameter (MVD) varied between 8.1 and 12.0 m. The LWC and MVD showed a decrease with height except in the middle region of the cloud where the values were higher than the adjacent levels. The dispersion of the cloud drops was lower (0.65) in the cloud-top region and higher (1.01) in the cloud-base region.The observed cloud microphysical characteristics were attributed to vertical mixing in clouds induced by the cloud-top gravity oscillations (buoyancy oscillations) generated by the intensification of turbulent eddies due to the buoyant production of energy by the microscale-fractional-condensation (MFC) in turbulent eddies.