Effects of the geometry of two‐dimensional fractures on their hydraulic aperture and on the validity of the local cubic law

Flow through rough fractures is investigated numerically in order to assess the validity of the local cubic law for different fracture geometries. Two‐dimensional channels with sinusoidal walls having different geometrical properties defined by the aperture, the amplitude, and the wavelength of the walls' corrugations, the corrugations asymmetry, and the phase shift between the two walls are considered to represent different fracture geometries. First, it is analytically shown that the hydraulic aperture clearly deviates from the mean aperture when the walls' roughness, the phase shift, and/or the asymmetry between the fracture walls are relatively high. The continuity and the Navier–Stokes equations are then solved by means of the finite element method and the numerical solutions compared to the theoretical predictions of the local cubic law. Reynolds numbers ranging from 0.066 to 66.66 are investigated so as to focus more particularly on the effect of flow inertial effects on the validity of the local cubic law. For low Reynolds number, typically less than 15, the local cubic law properly describes the fracture flow, especially when the fracture walls have small corrugation amplitudes. For Reynolds numbers higher than 15, the local cubic law is valid under the conditions that the fracture presents a low aspect ratio, small corrugation amplitudes, and a moderate phase lag between its walls.     

Bacterial productivity in the Prydz Bay and its adjacent waters,Antarctic

1　IntroductionBacteriaandtheiractivitiesplayanimportantroleintheelementalbiogeochemicalcyclesandenergytransformingintheocean (Zhenetal.1 997) .DortchandPackard(1 989) proposedthatfoodwebsintheeutrophicwatersaredominatedbythebiomassofprimaryproducerswhilefoodwebsintheoligotrophicwatersaredominatedbythebiomassofmicrobes.Heterotrophicbacteriahadbeenshowntoplayanimportantroleinthedecompositionoflarge ,rapidlysinkingorganicparticleswithinandbelowtheeuphot iczone ,andfurthertoaffecttheelementaldyn…     

一种从经纬度网格资料获得正方形网格资料的插值方案

本文提出了一种把经纬度网格点上的气象要素值转化到正方形网格点上的多项式插值方案。对插值方案作了稳定性讨论和定量及定性分析。结果表明,这种插值方案方便可行,而且有广泛适用性。     

Multiple longshore sand bars in the upper Chesapeake Bay

In the upper Chesapeake Bay (Maryland, U.S.A.) field surveys were conducted at 18 multiple longshore sand bar sites. The multiple bar systems were found in water depths less than approximately 2 m (mean sea level), and exhibited mild bottom slopes of 0·0052 or less. The number of bars composing each system ranged from four to 17 and the spacing between the crests typically increased in the offshore direction, ranging from 12 to 70 m. Bar height also typically increased with distance offshore and ranged from 0·03 to 0·61 m. A grain size analysis of crest and trough sediment did not reveal any significant differences and the sediment was categorized as ‘fine sand’. A review of the literature data indicated that the Chesapeake Bay multiple bars possessed similar characteristics to those found in Gelding Bay (Baltic Sea); similarities in fetch, wave height and tidal range between the two bays may account for this finding. The surf-scaling parameter indicated that the multiple bar systems were extremely dissipative with regard to wave energy, and wave height appeared to be an important factor in controlling bar spacing and bar height. A multiple wave break point hypothesis was discussed as a possible mechanism for the formation of Chesapeake Bay multiple longshore bars, and limited observational evidence appeared to support such a mechanism.     

有限元法与伪谱法混合求解弹性波动方程

在地震波场数值模拟中,有限差分法、有限元法和伪谱法都是常用的基本方法,但它们各有不同的适应性和优缺点,如有限差分法、有限元法都存在减弱网格频散和提高计算效率的矛盾,而伪谱法的网格频散小且计算效率高.有限差分法和伪谱法在处理地表结构复杂或地表剧烈起伏以及地下结构复杂的情况时存在较大的难度,而有限元法可较为理想地拟合起伏地表和任意弯曲界面,且可方便地处理自由边界条件和界面边界条件.尝试将有限元法和伪谱法相结合,形成地震波场数值模拟的一种混合方法,利用二者的优点,克服二者的缺点,达到既减弱网格频散又提高计算精度和效率的目的.并采用所谓的‘过度区域‘技术解决两种不同算法的衔接问题.模拟实例表明,给出的混合模拟方法不失为弹性波场数值模拟的一种有效方法.     

Tidal effects on temperature front in the Yellow Sea

Temperature front (TF) is one of the important features in the Yellow Sea, which forms in spring, thrives in summer, and fades in autumn as thermocline declines. TF intensity ⋎S _T ⋎ is defined to describe the distribution of TF. Based on the MASNUM wave-tide-circulation coupled model, temperature distribution in the Yellow Sea was simulated with and without tidal effects. Along 36°N, distribution of TF from the simulated results are compared with the observations, and a quantitative analysis is introduced to evaluate the tidal effects on the forming and maintaining processes of the TF. Tidal mixing and the circulation structure adapting to it are the main causes of the TF. Supported by the National Basic Research Program of China (No. G1999043809) and the National Science Foundation of China (No. 49736190).     

Investigating the mechanisms leading to the deglaciation of past continental northern hemisphere ice sheets with the CLIMBER–GREMLINS coupled model

A coupling procedure between a climate model of intermediate complexity (CLIMBER-2.3) and a 3-dimensional thermo-mechanical ice-sheet model (GREMLINS) has been elaborated. The resulting coupled model describes the evolution of atmosphere, ocean, biosphere, cryosphere and their mutual interactions. It is used to perform several simulations of the Last Deglaciation period to identify the physical mechanisms at the origin of the deglaciation process. Our baseline experiment, forced by insolation and atmospheric CO₂, produces almost complete deglaciation of past northern hemisphere continental ice sheets, although ice remains over the Cordilleran region at the end of the simulation and also in Alaska and Eastern Siberia. Results clearly demonstrate that, in this study, the melting of the North American ice sheet is critically dependent on the deglaciation of Fennoscandia through processes involving switches of the thermohaline circulation from a glacial mode to a modern one and associated warming of the northern hemisphere. A set of sensitivity experiments has been carried out to test the relative importance of both forcing factors and internal processes in the deglaciation mechanism. It appears that the deglaciation is primarily driven by insolation. However, the atmospheric CO₂ modulates the timing of the melting of the Fennoscandian ice sheet, and results relative to Laurentide illustrate the existence of threshold CO₂ values, that can be translated in terms of critical temperature, below which the deglaciation is impeded. Finally, we show that the beginning of the deglaciation process of the Laurentide ice sheet may be influenced by the time at which the shift of the thermohaline circulation from one mode to the other occurs.     

CATASTROPHIC ECO-ENVIRONMENTAL CHANGE IN THE SONGNEN PLAIN, NORTHEASTERN CHINA SINCE 1900S

Although the Songnen Plain in the northeastern China was developed relatively late in the temperate zone of the world, its eco-environment has changed greatly. This paper analyzes the changes of land cover and the rates and trends ofdesertification during the past 100 years in the Songnen Plain. According to the macroscopic analysis, we find that the eco-environment in the plain has reached to the threshold of catastrophic change since the 1950s. The Thorn Needle Catastrophic Model was used to determine and validate this conclusion. Human activities, including large-scale construction projects, such as huge dams and dikes, and excessive grazing were the primary factors contributing to regional eco-environmental catastrophe. And irrational reclamation of the wilderness also affected the eco-environmental change. The results reveal the complex human-land interactions.     

One zone modeling of irregular variability of stellar convective envelope

One zone modeling of the irregular variability of red super-giants is intended with regard to the nonlinear coupling of finite amplitude pulsation with convection. The nonlocal mixing length is employed for the evaluation of the convective flux, the turbulent pressure and the turbulent power of temperature fluctuations. The radial pulsation and the Boussinesq convection are assumed for simplicity. The one zone is defined as the layer having the entropy maximum and the minimum at the bottom and at the top, respectively. The quasi-adiabatic approximation is consistent with this definition in fixing the zone to the same mass range. The spatial derivatives are evaluated under the assumption of homologous changes with the equilibrium homologous parameters. Then, a set of 6 simultaneous first order nonlinear ordinary differential equations are obtained as the one zone representation of the irregular variability of the convective envelope.     

A MESO-β SCALE SIMULATION OF THE EFFECTS OF BOREAL FOREST ECOSYSTEM ON THE LOWER ATMOSPHERE

Based on the Intensive Field Campaign(IFC-1)data of Boreal Ecosystem-Atmosphere Study(BOREAS).a three-dimensional meso-β scale model is used to simulate the effect of boreal forests onthe lower atmosphere.A fine horizontal resolution of 2 km×2 km is used in order to distinguish thevegetative heterogeneity in the boreal region.A total of 20×25 grid points cover the entire sub-modeling area in BOREAS' South Study Area(SSA).The ecosystem types and their coverage ineach grid square are extracted from the North American Land Cover Characteristics Data Base(NALCCD)generated by the U.S.Geographical Survey(USGS)and the University of Nebraska-Lincoln(UNL).The topography of the study area is taken from the Digital Elevation Map(DEM)of USGS.The model outputs include the components of the energy balance budget within the canopyand at the ground.the turbulence parameters in the atmospheric boundary layer and the wind.temperature and humidity profiles extending up to a height of 1500 m.In addition to the fine timeand spatial step,the unique feature of the present model is the incorporation of both dynamic andbiological effects of the Boreal forest into the model parameterization scheme.The model resultscompare favorably with BOREAS' IFC-1 data in 1994 when the forest was in the luxuriant growingperiod.