On the growth of ice crystals in supercooled water and aqueous solution drops

Summary The results of a recent investigation on the free growth rate and the growth forms of ice in supercooled water and aqueous solutions are presented. The results are used to discuss the structure of frozen drops, the structure of hailstones, the mechanisms which are responsible for the glaciation of atmospheric clouds and the mechanisms which cause cloud electrification. It was found that the presence of dissolved salts and dust particles in cloud drops favor the formation of spongy and polycristalline ice, that it is unlikely for frozen cloud drops to develop into hexagonal shaped ice single-crystals, and that it is also unlikely that in atmospheric clouds freezing drops shatter and splinter. The latter result casts serious doubts on the splintering mechanism to contribute to thunderstorm electrification and to promote glaciation of clouds.This work was supported by the National Science Foundation under Grant No. GP 2922.     

Growth of volcanic ash aggregates in the presence of liquid water and ice: an experimental approach

Key processes influencing the aggregation of volcanic ash and hydrometeors are examined with an experimental method employing vibratory pan aggregation. Mechanisms of aggregation in the presence of hail and ice pellets, liquid water (≤30?wt%), and mixed water phases are investigated at temperatures of 18 and ?20?°C. The experimentally generated aggregates, examined in hand sample, impregnated thin sections, SEM imagery, and X-ray microtomography, closely match natural examples from phreatomagmatic phases of the 27?ka Oruanui and 2010 Eyjafjallaj?kull eruptions. Laser diffraction particle size analysis of parent ash and aggregates is also used to calculate the first experimentally derived aggregation coefficients that account for changing liquid water contents and subzero temperatures. These indicate that dry conditions (<5–10?wt% liquid) promote strongly size selective collection of sub-63?μm particles into aggregates (given by aggregation coefficients >1). In contrast, liquid-saturated conditions (>15–20?wt% liquid) promote less size selective processes. Crystalline ice was also capable of preferentially selecting volcanic ash <31?μm under liquid-free conditions in a two-stage process of electrostatic attraction followed by ice sintering. However, this did not accumulate more than a monolayer of ash at the ice surface. These quantitative relationships may be used to predict the timescales and characteristics of aggregation, such as aggregate size spectra, densities, and constituent particle size characteristics, when the initial size distribution and water content of a volcanic cloud are known. The presence of an irregularly shaped, millimeter-scale vacuole at the center of natural aggregates was also replicated during interaction of ash and melting ice pellets, followed by sublimation. Fine-grained rims were formed by adding moist aggregates to a dry mixture of sub-31?μm ash, which adhered by electrostatic forces and sparse liquid bridges. From this, we infer that the fine-grained outer layers of natural aggregates reflect recycled exposure of moist aggregates to regions of volcanic clouds that are relatively dry and dominated by <31?μm ash.     

Spherical tensor approach to the solution of the mantle stress problem

Summary The technique of spherical tensors is applied to the problem of stress in the mantle. An efficient method for stress computation is developed for the Newtonian mantle with spherical boundaries assuming that the density and gravitational potential in the form of harmonic expansions are known. The ways of including more complicated rheological models as well as the non-spherical shape of boundaries are outlined.

---

¶rt; m¶rt; ama anu amuu, m a a nuuu annaama uuauu m. ¶rt;naam, m au nmmu uaumau nmuaa n uu u um u m amu n¶rt;mam u¶rt;m uuuauau. ¶rt; mu u m¶rt;a a a uu ¶rt; u u au.

     

The adjustment of the meandering pattern of the lower Jordan River to change in water discharge

The annual discharge of the Lower Jordan River has decreased from about 1250 × 10⁶ m³ y^?1 to about 300 × 10⁶ m³ y^?1 due to water exploitation. The decrease in water discharge was not followed by a similar decrease in the suspended sediment. Evidence from air photographs and maps from 1920s to early 1970s show that the river is adjusting itself by shortening its course. The shortening is not distributed uniformly along the valley length. Local effects obscure the effect of decreasing discharge. Since the early 1970s the channel is becoming longer again due to a change in the dominant flow.     

冰川均衡调整对南极冰质量平衡监测的影响及其不确定性

本文研究了新的全球冰川均衡调整(GIA)模型对南极冰盖质量平衡监测的影响,考虑现有冰川负荷模型和地幔黏滞度模型的差异,完整评估了结果的不确定性,最后结合GRACE和卫星测高的结果进行了对比分析.结果表明,GIA对GRACE监测的等效水柱变化有重大影响,较大的GIA影响出现在西南极,沿罗斯冰架-卡姆布冰流-罗尼冰架-南极半岛分布,最大值在卡姆布冰流,达到29.8 mm/a;GIA对南极整体冰质量平衡的影响达到134±28 Gt/a.在不确定性的方差中,西南极和东南极分别以冰负荷模型差异和地幔黏滞度差异影响为主,对整个南极,冰模型差异影响占88.4%;在一些典型地区,GRACE监测的等效水柱在扣除GIA前后分别是,卡姆布冰流~32.8 mm/a和~6.3 mm/a,阿蒙森海湾~-95.3 mm/a和~-102.5 mm/a,Enderby Land~13.6 mm/a和~8.1 mm/a.整个南极冰盖总质量变化在扣除GIA贡献后为-82±29 Gt/a,该估计与卫星测高结果较吻合.此外,GIA对卫星测高监测的冰面高程变化的影响一般不超过8%.本研究为空间大地测量监测南极冰质量平衡提供了新的改正模型.     

Alternative analysis of transient infiltration experiment to estimate soil water repellency

The repellency index (RI) defined as the adjusted ratio between soil‐ethanol, S_e, and soil‐water, S_w, sorptivities estimated from minidisk infiltrometer experiments has been used instead of the widely used water drop penetration time and molarity of ethanol drop tests to assess soil water repellency. However, sorptivity calculated by the usual early‐time infiltration equation may be overestimated as the effects of gravity and lateral capillary are neglected. With the aim to establish the best applicative procedure to assess RI, different approaches to estimate S_e and S_w were compared that make use of both the early‐time infiltration equation (namely, the 1 min, S1, and the short‐time linearization approaches), and the two‐term axisymmetric infiltration equation, valid for early to intermediate times (namely, the cumulative linearization and differentiated linearization approaches). The dataset included 85 minidisk infiltrometer tests conducted in three sites in Italy and Spain under different vegetation habitats (forest of Pinus pinaster and Pinus halepensis, burned pine forest, and annual grasses), soil horizons (organic and mineral), postfire treatments, and initial soil water contents. The S1 approach was inapplicable in 42% of experiments as water infiltration did not start in the first minute. The short‐time linearization approach yielded a systematic overestimation of S_e and S_w that resulted in an overestimation of RI by a factor of 1.57 and 1.23 as compared with the cumulative linearization and differentiated linearization approaches. A new repellency index, RI_s, was proposed as the ratio between the slopes of the linearized data for the wettable and hydrophobic stages obtained by a single water infiltration test. For the experimental conditions considered, RI_s was significantly correlated with RI and WDPT. Compared with RI, RI_s includes information on both soil sorptivity and hydraulic conductivity and, therefore, it can be considered more physically linked to the hydrological processes affected by soil water repellency.     

Glacial isostatic adjustment at the Laurentide ice sheet margin: Models and observations in the Great Lakes region

Glacial Isostatic Adjustment (GIA) modelling in North America relies on relative sea level information which is primarily obtained from areas far away from the uplift region. The lack of accurate geodetic observations in the Great Lakes region, which is located in the transition zone between uplift and subsidence due to the deglaciation of the Laurentide ice sheet, has prevented more detailed studies of this former margin of the ice sheet. Recently, observations of vertical crustal motion from improved GPS network solutions and combined tide gauge and satellite altimetry solutions have become available. This study compares these vertical motion observations with predictions obtained from 70 different GIA models. The ice sheet margin is distinct from the centre and far field of the uplift because the sensitivity of the GIA process towards Earth parameters such as mantle viscosity is very different. Specifically, the margin area is most sensitive to the uppermost mantle viscosity and allows for better constraints of this parameter. The 70 GIA models compared herein have different ice loading histories (ICE-3/4/5G) and Earth parameters including lateral heterogeneities. The root-mean-square differences between the 6 best models and the two sets of observations (tide gauge/altimetry and GPS) are 0.66 and 1.57 mm/yr, respectively. Both sets of independent observations are highly correlated and show a very similar fit to the models, which indicates their consistent quality. Therefore, both data sets can be considered as a means for constraining and assessing the quality of GIA models in the Great Lakes region and the former margin of the Laurentide ice sheet.     

Glacial isostatic adjustment at the Laurentide ice sheet margin: Models and observations in the Great Lakes region

Glacial Isostatic Adjustment (GIA) modelling in North America relies on relative sea level information which is primarily obtained from areas far away from the uplift region. The lack of accurate geodetic observations in the Great Lakes region, which is located in the transition zone between uplift and subsidence due to the deglaciation of the Laurentide ice sheet, has prevented more detailed studies of this former margin of the ice sheet. Recently, observations of vertical crustal motion from improved GPS network solutions and combined tide gauge and satellite altimetry solutions have become available. This study compares these vertical motion observations with predictions obtained from 70 different GIA models. The ice sheet margin is distinct from the centre and far field of the uplift because the sensitivity of the GIA process towards Earth parameters such as mantle viscosity is very different. Specifically, the margin area is most sensitive to the uppermost mantle viscosity and allows for better constraints of this parameter. The 70 GIA models compared herein have different ice loading histories (ICE-3/4/5G) and Earth parameters including lateral heterogeneities. The root-mean-square differences between the 6 best models and the two sets of observations (tide gauge/altimetry and GPS) are 0.66 and 1.57 mm/yr, respectively. Both sets of independent observations are highly correlated and show a very similar fit to the models, which indicates their consistent quality. Therefore, both data sets can be considered as a means for constraining and assessing the quality of GIA models in the Great Lakes region and the former margin of the Laurentide ice sheet.     

Soliton generation in clear water and in the presence of surfactants

It is shown that the effect of wind in a circular wind-water tunnel can cause the formation of a single soliton, whose generation is preceded by a gradual increase in the length and amplitude of wind waves. In the presence of surfactants, another (long-wave) formation mechanism comes into play, resulting in the generation of a soliton on the surface free of wind waves in the low-frequency range, which will persist from the origination of waves to a completely formed soliton.     

应力调整改变介质渗透性对水位的影响

在岩石抽水-注水试验的基础上,通过应力和渗透系数的关系,用数值模拟的方法讨论了应力与渗流的过程。研究结果表明,随着应力场的调整介质的渗透性将发生改变,从而对渗流过程产生影响,这种影响是非线性的,渗流作用的存在使得水位的变化与介质的渗透性的变化不是对应的关系。     

水管倾斜仪电灵敏度的标定调试经验

根据以往维护水管倾斜仪的一些经验，介绍了在DSQ水管倾斜仪的标定过程中，调整换能器电压灵敏度的一些技巧。     

Modification of water masses in the Barents Sea and its coupling to ice dynamics: a model study

The physical and biological environment of the Barents Sea is characterised by large variability on a wide range of scales. Results from a numerical ocean model, SINMOD, are presented showing that the physical variability is partly forced by changes in annual net ice import. The mean contribution from ice import in the simulation period (1979–2007) is about 40% of the total amount of ice melted each year. The annual ice import into the Barents Sea varies between 143 and 1,236 km³, and this causes a substantial variability in the amount of annual ice melt in the Barents Sea. This in turn impacts the freshwater content. The simulated freshwater contribution from ice is 0.02 Sv on average and 0.04 Sv at maximum. When mixed into a mean net Atlantic Water (AW) inflow of 1.1 Sv with a salinity of 35.1, this freshwater addition decreases the salinity of the modified AW to 34.4 and 33.9 for the mean and maximum freshwater fluxes, respectively. Ice import may thus be important for the Barents Sea production of Arctic Ocean halocline water which has salinity of about 34.5. The changes in the ice melt the following summer due to ice import also affect the formation of dense water in the Barents Sea by changing stratification, altering the vertical mixing rates and affecting heat loss from the warm AW. The model results thus indicate that ice import from the Arctic has a great impact on water mass modification in the Barents Sea which in turn impacts the ventilation of the Arctic Ocean.     

砂泥岩储层孔隙度和含水饱和度同步反演

基于岩石物理和地震反演理论,提出了一种同步反演储层孔隙度和含水饱和度的方法.以岩石物理为基础,建立了砂泥岩储层物性和弹性参数之间定量的关系-Simon模型,以贝叶斯理论为手段,结合不同类型的砂泥岩储层,建立了多信息联合约束的物性参数反演目标函数,并通过蒙特卡罗和遗传算法相结合的思路求解该目标函数,最终得到孔隙度和含水饱和度的同步反演结果.将该方法应用于河道砂和砂砾岩两种不同的砂泥岩储层中,孔隙度和含水饱和度数据的联合应用,进一步减少了储层预测的多解性,为石油地质综合研究提供了更加丰富准确的基础数据.

     

Export fluxes of biogenic matter in the presence and absence of seasonal sea ice cover in the Chukchi Sea

Drifting sediment traps were deployed at 9 stations in May-June (ice-covered conditions) and July-August (ice-free conditions) 2004 in the Chukchi Sea to investigate the variability in export fluxes of biogenic matter in the presence and absence of sea ice cover. Measurements of chlorophyll-a (Chl-a), particulate organic carbon (POC), particulate nitrogen (PN), phytoplankton, zooplankton fecal pellets, and the stable carbon isotope composition (δ¹³C) of the sinking material were performed along Barrow Canyon (BC) and a parallel shelf-to-basin transect from East Hanna Shoal (EHS) to the Canada Basin. POC export fluxes were similarly high in the presence (378±106 mg C m⁻² d⁻¹) and in the absence of ice cover (442±203 mg C m⁻² d⁻¹) at the BC stations, while fluxes were significantly higher in the absence (129±98 mg C m⁻² d⁻¹) than in the presence of ice cover (44±29 mg C m⁻² d⁻¹) at the EHS stations. The C/N ratios and δ¹³C values of sinking organic particles indicated that POC export fluxes on the Chukchi continental shelf were mostly composed of freshly produced labile material, except at the EHS stations under ice cover where the exported matter was mostly composed of refractory material probably advected into the EHS region. Chl-a fluxes were higher under ice cover than in ice-free water, however, relatively low daily loss rates of Chl-a and similar phytoplankton carbon fluxes in ice-covered and ice-free water suggest the retention of phytoplankton in the upper water column. An increase in fecal pellet carbon fluxes in ice-free water reflected higher grazing pressure in the absence of ice cover. Elevated daily loss rates of POC at the BC stations confirmed other indications that Barrow Canyon is an important area of carbon export to the basin and/or benthos. These results support the conclusion that there are large spatial and temporal variations in export fluxes of biogenic matter on the Chukchi continental shelf, although export fluxes may be similar in the presence and in the absence of ice cover in highly productive regions.     

The problem of transfer of radionuclide pollution by sea ice

The problem of the radioactive pollution of the ice cover is discussed in this paper. An estimation of the possible transport of the radionuclides by drifting ice from the Kara Sea into the Arctic Ocean is presented. The conclusion about the importance of the contribution of drifting ice to the contamination of the Greenland coastal water is based on the comparison of estimates of the multiyear export of radioactive drifting ice from the Kara Sea with the observed concentrations of the radionuclides in the Greenland coastal water.     

Prediction of capillary air-liquid interfacial area vs. saturation function from relationship between capillary pressure and water saturation

This short communication investigates if the capillary air-liquid interfacial area vs. saturation relationship A_lv(S) can be predicted from the capillary pressure vs. saturation relationship S(h), using the theoretical sample scale model of Diamantopoulos and Durner (2013, 2015). We selected three published experimental datasets, where S(h) and A_lv(S) relationships had been measured for the same porous media. The sample scale model was fitted to the retention curve S(h) of each porous medium and then used to predict the air-liquid interfacial area A_lv(S). We also included in the analysis the thermodynamic models of Leverett (1941) and Grant and Gerhard (2007). For two sandy materials and especially for high saturation values, the model predicted the capillary A_lv(S) successfully, which was in one case given by a pore-network model simulation (Kibbey and Chen, 2012) and in the other case experimentally determined (Brusseau et al., 2006). For glass bead experiments, the contact angle needed to be fitted to properly describe the experimental A_lv(S) curve.     

A Brovar-type solution of the fixed geodetic boundary-value problem

For more than 150 years gravity anomalies have been used for the determination of geoidal heights, height anomalies and the external gravity field. Due to the fact that precise ellipsoidal heights could not be observed directly, traditionally a free geodetic boundary-value problem (GBVP) had to be formulated which after linearisation is related to gravity anomalies. Since nowadays the three-dimensional positions of gravity points can be determined by global navigation satellite systems very precisely, the modern formulation of the GBVP can be based on gravity disturbances which are related to a fixed GBVP using the known topographical surface of the Earth as boundary surface. The paper discusses various approaches into the solution of the fixed GBVP which after linearization corresponds to an oblique-derivative boundary-value problem for the Laplace equation. Among the analytical solution approaches a Brovar-type solution is worked out in detail, showing many similarities with respect to the classical solution of the scalar free GBVP.     

Numerical solution of the two-dimensional inverse geomagnetic induction problem

Summary An effective numerical approach to the solution of the two-dimensional inverse geomagnetic induction problem using the linearization method is presented. The numerical realization of the inversion is based on Marquardt's algorithm, for which the solution of the direct problem and the partial derivatives of this solution with respect to the electrical parameters of the medium are computed by the finite difference method. Theoretical models are studied and numerical results are presented.