Near-surface wind field characteristics of the desert-oasis transition zone in Dunhuang,China

The desert-oasis transition zone(DOTZ)serves as a buffer area between the desert and oasis.Understanding its wind field characteristics is of great significance for the prevention and control of aeolian disasters in the oasis.In this study,we used meteorological data during 2013–2019 from the portable meteorological stations at five sites(site A on the edge of the oasis,sites B,C,and D in the DOTZ,and site O in the desert)in Dunhuang,China to analyze the near-surface wind field characteristics and their causes,as well as to reveal the key role of the DOTZ in oasis protection.The results showed that the mean wind speed,frequency of sand-driving wind,and directional variability of wind decreased from west to east within the DOTZ,and wind speed was significantly affected by air temperature.The terrain influenced the prevailing winds in the region,mainly from northeast and southwest.Only some areas adjacent to the oasis were controlled by southeasterly wind.This indicated that the near-surface wind field characteristics of the DOTZ were caused by the combined effects of local terrain and surface hydrothermal difference.At site D,the annual drift potential(DP)was 24.95 vector units(VU),indicating a low wind energy environment,and the resultant drift direction(RDD)showed obvious seasonal differences.Additionally,the DOTZ played an important buffering role between the desert and oasis.Compared with the desert,the mean wind speed in the oasis decreased by 64.98%,and the prevailing wind direction was more concentrated.The results of this study will be useful in interpreting the aeolian activity of the DOTZ in Dunhuang.

3D near-surface P-wave velocity structure imaging with Distributed Acoustic Sensing and electric hammer source

Distributed Acoustic Sensing(DAS) is an emerging technique for ultra-dense seismic observation, which provides a new method for high-resolution sub-surface seismic imaging. Recently a large number of linear DAS arrays have been used for two-dimensional S-wave near-surface imaging in urban areas. In order to explore the feasibility of three-dimensional(3D) structure imaging using a DAS array, we carried out an active source experiment at the Beijing National Earth Observatory. We deployed a 1 km optical cable in a rectangular shape, and the optical cable was recast into 250 sensors with a channel spacing of 4 m. The DAS array clearly recorded the P, S and surface waves generated by a hammer source. The first-arrival P wave travel times were first picked with a ShortTerm Average/Long-Term Average(STA/LTA) method and further manually checked. The P-wave signals recorded by the DAS are consistent with those recorded by the horizontal components of short-period seismometers. At shorter source-receiver distances, the picked P-wave arrivals from the DAS recording are consistent with vertical component recordings of seismometers, but they clearly lag behind the latter at greater distances.This is likely due to a combination of the signal-to-noise ratio and the polarization of the incoming wave. Then,we used the Tomo DD software to invert the 3D P-wave velocity structure for the uppermost 50 m with a resolution of 10 m. The inverted P-wave velocity structures agree well with the S-wave velocity structure previously obtained through ambient noise tomography. Our study indicates the feasibility of 3D near-surface imaging with the active source and DAS array. However, the inverted absolute velocity values at large depths may be biased due to potential time shifts between the DAS recording and seismometer at large source-receiver distances.

Near-surface velocity model construction based on a Monte-Carlo scheme

In tomographic statics seismic data processing, it 1s crucial to cletermme an optimum base for a near-surface model. In this paper, we consider near-surface model base determination as a global optimum problem. Given information from uphole shooting and the first-arrival times from a surface seismic survey, we present a near-surface velocity model construction method based on a Monte-Carlo sampling scheme using a layered equivalent medium assumption. Compared with traditional least-squares first-arrival tomography, this scheme can delineate a clearer, weathering-layer base, resulting in a better implementation of damming correction. Examples using synthetic and field data are used to demonstrate the effectiveness of the proposed scheme.

Bearing splitting and near-surface source ranging in the direct zone of deep water

Sound multipath propagation is very important for target localization and identification in different acoustical zones of deep water. In order to distinguish the multipath characteristics in deep water, the Northwest Pacific Acoustic Experiment was conducted in 2015. A low-frequency horizontal line array towed at the depth of around 150 m on a receiving ship was used to receive the noise radiated by the source ship. During this experiment, a beating-splitting phenomenon in the direct zone was observed through conventional beamforming of the horizontal line array within the frequency band 160 Hz- 360 Hz. In this paper, this phenomenon is explained based on ray theory. In principle, the received signal in the direct zone of deep water arrives from two general paths including a direct one and bottom bounced one, which vary considerably in arrival angles. The split bearings correspond to the contributions of these two paths. The beating-splitting phenomenon is demonstrated by numerical simulations of the bearing-time records and experimental results, and they are well consistent with each other. Then a near-surface source ranging approach based on the arrival angles of direct path and bottom bounced path in the direct zone is presented as an application of bearing splitting and is verified by experimental results. Finally, the applicability of the proposed ranging approach for an underwater source within several hundred meters in depth in the direct zone is also analyzed and demonstrated by simulations.

VALIDATION OF NEAR-SURFACE WINDS OBTAINED BY A HYBRID WRF/CALMET MODELING SYSTEM OVER A COASTAL ISLAND WITH COMPLEX TERRAIN

The results from a hybrid approach that combines a mesoscale meteorological model with a diagnostic model to produce high-resolution wind fields in complex coastal topography are evaluated.The diagnostic wind model(California Meteorological Model,CALMET) with 100-m horizontal spacing was driven with outputs from the Weather Research and Forecasting(WRF) model to obtain near-surface winds for the 1-year period from 12 September 2003 to 11 September 2004.Results were compared with wind observations at four sites.Traditional statistical scores,including correlation coefficients,standard deviations(SDs) and mean absolute errors(MAEs),indicate that the wind estimates from the WRF/CALMET modeling system are produced reasonably well.The correlation coefficients are relatively large,ranging from 0.5 to 0.7 for the zonal wind component and from 0.75 to 0.85 for the meridional wind component.MAEs for wind speed range from 1.5 to 2.0 m s-1 at 10 meters above ground level(AGL) and from 2.0 to 2.5 m s-1 at 60 m AGL.MAEs for wind direction range from 30 to 40 degrees at both levels.A spectral decomposition of the time series of wind speed shows positive impacts of CALMET in improving the mesoscale winds.Moreover,combining the CALMET model with WRF significantly improves the spatial variability of the simulated wind fields.It can be concluded that the WRF/CALMET modeling system is capable of providing a detailed near-surface wind field,but the physics in the diagnostic CALMET model needs to be further improved.

Numerical Study of Viscous Wave-Making Resistance of Ship Navigation in Still Water

The prediction of a ship＇s resistance especially the viscous wave-making resistance is an important issue in CFD applications. In this paper, the resistances of six ships from hull 1 to hull 6 with different hull forms advancing in still water are numerically studied using the solver naoe-FOAM-SJTU, which was developed based on the open source code package OpenFOAM. Different components of the resistances are computed and compared while considering two speed conditions （12 kn and 16 kn）. The resistance of hull 3 is the smallest while that of hull 5 is the largest at the same speed. The results show hull 3 is a good reference for the design of similar ships, which can provide some valuable guidelines for hull form optimization.

Numerical Research on Wave-making Resistance of Trimaran

A 3D rankine panel method was developed for calculating the linear wave-making resistance of a tri-maran with Wigley hulls. In order to calculate the normal vector and derivative of the body surface accurately, non-uniform rational B-spline (NURBS) was adopted to represent body surface and rankine source density. The radiation condition is satisfied using the numerical technology of staggered grids. Numerical results show that the linear wave-making resistance of the trimaran can be calculated effectively using this method.

Near-surface absorption compensation technology and its application in the Daqing Oilfields

High-frequency seismic data components can be seriously attenuated during seismic wave propagation in unconsolidated （low-velocity） layers, resulting in reduced seismic resolution and signal-to-noise （S/N） ratio. In this paper, first, based on Wiener filter theory, inverse filter calculations for near-surface absorption attenuation compensation were accomplished by analysis of the direct wave spectral components from different distances near the surface. The direct waves were generated by detonators in uphole shots and were acquired by receivers on the surface. The spatially varying inverse filters were designed to compensate for the frequency attenuation of 3D pre-stack CRG （common receiver-gather） data. After applying the filter to CRG data, the high frequency components were compensated with the low frequencies maintained. The seismic resolution and S/N ratio are enhanced and match better with synthetic seismograms and better meet the needs of geological interpretation.

A Havelock Source Panel Method for Near-surface Submarines

A panel method is described for calculating potential flow around near-surface submarines. The method uses Havelock sources which automatically satisfy the linearized free-surface boundary condition. Outputs from the method include pressure field, pressure drag, wave resistance, vertical force, trim moment and wave pattern. Comparisons are made with model tests for wave resistance of Series 58 and DARPA SUBOFF hulls, as well as with wave resistance, lift force and trim moment of three length-to-diameter variants of the DSTO Joubert submarine hull. It is found that the Havelock source panel method is capable of determining with reasonable accuracy wave resistance, vertical force and trim moment for submarine hulls. Further experimental data are required in order to assess the accuracy of the method for pressure field and wave pattern prediction. The method is implemented in the computer code“HullWave”and offers potential advantages over RANS-CFD codes in terms of speed, simplicity and robustness.

Quasi-plane-hypothesis of strain coordination for RC beams seismically strengthened with externally-bonded or near-surface mounted fiber reinforced plastic

The application of fiber reinforced plastic（FRP）,including carbon FRP and glass FRP,for structural repair and strengthening has grown due to their numerous advantages over conventional materials such as externally bonded reinforcement（EBR） and near-surface mounted（NSM） strengthening techniques.This paper summarizes the results from 21 reinforced concrete beams strengthened with different methods,including externally-bonded and near-surface mounted FRP,to study the strain coordination of the FRP and steel rebar of the RC beam.Since there is relative slipping between the RC beam and the FRP,the strain of the FRP and steel rebar of the RC beam satisfy the quasi-plane-hypothesis;that is,the strain of the longitudinal fiber that parallels the neutral axis of the plated beam within the scope of the effective height（h 0） of the cross section is in direct proportion to the distance from the fiber to the neutral axis.The strain of the FRP and steel rebar satisfies the equation：ε FRP =βε steel,and the value of β is equal to 1.1-1.3 according to the test results.

Improved seismoacoustic analysis model and its application to source parameter inversion of near-surface small-yield chemical explosions

The seismoacoustic analysis method has broad potential applications to source parameter estimation for near-surface explosion events such as industrial explosions and terrorist attacks.In this study,current models were improved by modifying the acoustic model and adopting the Bayesian Markov-chain-Monte-Carlo inversion method.The source parameters of near-surface small-yield chemical explosions were analyzed via the improved seismoacoustic analysis model and by the estimation accuracy of seismoacoustic joint inversion.Estimation and analysis results showed that the improved seismoacoustic analysis model considered ground shock coupling and the impact of explosion products ejecting from the surface so that the improved acoustic impulse relation was more consistent with the measured data than the Ford impulse relation.It is suitable for deep-burial,shallow-burial,and near-surface aerial explosions.Furthermore,trade-off relationships were declined through the application of the improved model to source parameter inversion for near-surface small-yield chemical explosions,and source parameter estimation accuracy was improved.

Active source seismic imaging on near-surface granite body:case study of siting a geological disposal repository for high-level radioactive nuclear waste

In order to research whether it is suitable to set a geological disposal repository for high-level radioactive nuclear waste into one target granite body,two active source seismic profles were arranged near a small town named Tamusu,Western China.The study area is with complex surface conditions,thus the seismic exploration encountered a variettraveltimey of technical difculties such as crossing obstacles,de-noising harmful scattered waves,and building complex near-surface velocity models.In order to address those problems,techniques including cross-obstacle seismic geometry design,angle-domain harmful scattered noise removal,and an acoustic wave equation-based inversion method jointly utilizing both the and waveform of frst arrival waves were adopted.The fnal seismic images clearly exhibit the target rock’s unconformable contact boundary and its top interface beneath the sedimentary and weathered layers.On this basis,it could be confrmed that the target rock is not thin or has been transported by geological process from somewhere else,but a native and massive rock.There are a few small size fractures whose space distribution could be revealed by seismic images within the rock.The fractures should be kept away.Based on current research,it could be considered that active source seismic exploration is demanded during the sitting process of the geological disposal repository for nuclear waste.The seismic acquisition and processing techniques proposed in the present paper would ofer a good reference value for similar researches in the future.

Double linear strain distribution assumption of RC beam strengthened with external-bonded or near-surface mounted fiber reinforced plastic

Rehabilitation of existing structures with fiber reinforced plastic(FRP)has been growing in popularity because they offer superior performance in terms of resistance to corrosion and high specific stiffness.The strain coordination results of 34 reinforced concrete beams(four groups)strengthened with different methods were presented including external-bonded or near-surface mounted glass or carbon FRP or helical rib bar in order to study the strain coordination of the strengthening materials and steel rebar of RC beam.Because there is relative slipping between concrete and strengthening materials(SM),the strain of SM and steel rebar of RC beam satisfies the double linear strain distribution assumption,that is,the strain of longitudinal fiber parallel to the neutral axis of plated beam within the scope of effective height(h0)of the cross section is in direct proportion to the distance from the fiber to the neutral axis.The strain of SM and steel rebar satisfies the equation εGCH=βεsteel,where the value of β is equal to 1.1-1.3 according to the test results.

Linking moisture and near-surface wind with winter temperature to reveal the Holocene climate evolution in arid Xinjiang region of China

An increasing number of palaeo-climatic records have been reported to identify the Holocene climate history in the arid Xinjiang region of northwest China. However, few studies have fully considered the internal linkages within the regional climate system, which may limit our understanding of the forcing mechanisms of Holocene climate change in this region. Here, we systematically consider three major issues of the moisture/precipitation, temperature and near-surface wind relevant to the Holocene climate history of Xinjiang. First, despite there still has debated for the Holocene moisture evolution in this region, more climatic reconstructions from lake sediments, loess, sand-dunes and peats support a long-term regional wetting trend. Second, temperature records from ice cores, peats and stalagmites demonstrate a long-term winter warming trend during the Holocene in middle-to high-latitudes of Asia. Third, recent studies of aeolian sedimentary sequences reveal that the near-surface winds in winter gradually weakened during the Holocene, whereas the winter mid-latitude Westerlies strengthened in the Tienshan Mountains. Based on this evidence, in the arid Xinjiang region we propose an early to middle Holocene relatively cold and dry interval, with strong near-surface winds;and a warmer, wetter interval with weaker near-surface winds in the middle to late Holocene during winter. Additionally,we develop a conceptual model to explain the pattern of Holocene climate changes in this region.From the early to the late Holocene, the increasing atmospheric COcontent and winter insolation,and the shrinking of high-latitude continental ice-sheets, resulted in increasing winter temperatures in middle to high latitudes in the Northern Hemisphere. Subsequently, the increased winter temperature strengthened the winter mid-latitude Westerlies and weakened the Siberian high-pressure system,which caused an increase in winter precipitation and a decrease in near-surface wind strength. This scenario is strongly supported by evidence from geological records, climate simulation results, and modern reanalysis data. Our hypothesis highlights the important contribution of winter temperature in driving the Holocene climatic evolution of the arid Xinjiang region, and it implies that the socio-economic development and water resources security of this region will face serious challenges presented by the increasing winter temperature in the future.

Temporal-spatial cross-correlation analysis of non-stationary near-surface wind speed time series

Temporal-spatial cross-correlation analysis of non-stationary wind speed time series plays a crucial role in wind field reconstruction as well as in wind pattern recognition.Firstly,the near-surface wind speed time series recorded at different locations are studied using the detrended fluctuation analysis(DFA),and the corresponding scaling exponents are larger than 1.This indicates that all these wind speed time series have non-stationary characteristics.Secondly,concerning this special feature( i.e.,non-stationarity)of wind signals,a cross-correlation analysis method,namely detrended cross-correlation analysis(DCCA) coefficient,is employed to evaluate the temporal-spatial cross-correlations between non-stationary time series of different anemometer pairs.Finally,experiments on ten wind speed data synchronously collected by the ten anemometers with equidistant arrangement illustrate that the method of DCCA cross-correlation coefficient can accurately analyze full-scale temporal-spatial cross-correlation between non-stationary time series and also can easily identify the seasonal component,while three traditional cross-correlation techniques(i.e.,Pearson coefficient,cross-correlation function,and DCCA method) cannot give us these information directly.

Effects of Near-surface Low Temperature on Freezing Injury of ‘Tunisian Soft-Seed’ Pomegranate

[Objective]The paper was to explore the influence of near-surface low temperature on cultivation of soft-seed pomegranate,and to provide guidance for planting location of soft-seed pomegranate.[Method]Taking 10 soft-seed pomegranate planting plots under different site conditions as the research objects,the near-surface low temperature of 45-50 cm was dynamically monitored from December 1,2018 to February 20,2019,and comparative analysis was made based on the local meteorological data over the same period.[Result]The near-surface low temperature of each temperature monitoring point was lower than the local meteorological data,which were all in the range of low temperature causing freezing in-jury of pomegranate trees,but the degree of freezing injury was different.The variation of near-surface low temperature was positively correlated with the altitude of terrain,but negatively correlated with the difference of topography and landform.When the local topography and landform were similar,the accumulation time of near-surface low temperature was negatively correlated with the altitude of terrain,while the duration of low tem-perature directly affected the degree of freezing injury.[Conclusion]The development of soft-seed pomegranate cultivation in Tunisia along Huang Mangling region in Henan Province refers to the local meteorological data.Meantime,it is also necessary to pay attention to the impact of regional microclimate environment,especially early monitoring of near-surface temperature to select suitable site and natural conditions.

Bias-Corrected Short-Range Ensemble Forecasts for Near-Surface Variables during the Summer Season of 2010 in Northern China

A running mean bias （RMB） correction ap- proach was applied to the forecasts of near-surface variables in a seasonal short-range ensemble forecasting experiment with 57 consecutive cases during summer 2010 in the northern China region. To determine a proper training window length for calculating RMB, window lengths from 2 to 20 days were evaluated, and 16 days was taken as an optimal window length, since it receives most of the benefit from extending the window length. The raw and 16-day RMB corrected ensembles were then evaluated for their ensemble mean forecast skills. The results show that the raw ensemble has obvious bias in all near-surface variables. The RMB correction can remove the bias reasonably well, and generate an unbiased ensemble. The bias correction not only reduces the ensemble mean forecast error, but also results in a better spreaderror relationship. Moreover, two methods for computing calibrated probabilistic forecast （PF） were also evaluated through the 57 case dates： 1） using the relative frequency from the RMB-eorrected ensemble; 2） computing the forecasting probabilities based on a historical rank histogram. The first method outperforms the second one, as it can improve both the reliability and the resolution of the PFs, while the second method only has a small effect on the reliability, indicating the necessity and importance of removing the systematic errors from the ensemble.

Unraveling the distribution patterns of near-surface temperature lapse rates in the Northwestern Kunlun Mountains

The near-surface temperature lapse rates for the core area of the Kunlun Mountains remain critically unresolved due to data scarcity.Here,we revealed the spatial and temporal patterns of nearsurface temperature lapse rate in the Kunlun Mountain regions based on both long-term meteorological records(1961-2017)and field surveys measured data(2012-2017).The results showed that(1)The near-surface temperature lapse rates(β;)has spatiotemporal distribution patterns on the Northwestern Kunlun Mountains(NWKM),and in complex terrain areas the complexity of the temperature-elevation relationship cannot be explained by the constant environmental temperature lapse rate(0.65℃/100 m)throughout alone.(2)Theβ;for the daily mean,minimum,and maximum temperature on the north slopes in the Kunlun mountain area are 0.41,0.47,and 0.37℃/100 m and on the Tiznafu River(TR)basin are 0.51,0.47 and 0.53℃/100 m,respectively.(3)The variations ofβ;for daily maximum and minimum temperature of the two regions exhibit similar monthly characteristics,which are lower in the winter and spring months than in other months.A greatest variability ofβ;for the daily mean,minimum,and maximum temperature appears in winter and a light variability ofβ;occurs in spring.The seasonal variability ofβ;for daily maximum temperature is greater than that for daily minimum temperature,and the seasonal variability ofβ;for daily average temperature has the smallest variability.(4)There is no significant trend of change occurred in theβ;of NWKM.(5)The spatial and temporal variations ofβ;for the NWKM are linked to the geographic differences and climate factors.The results of Grey Relational Analysis showed that theβ;distribution is mainly influenced by the wind speed and relative humidity of the NWKM.

IMPROVEMENT OF REGIONAL PREDICTION OF SEA FOG ON GUANGDONG COASTLAND USING THE FACTOR OF TEMPERATURE DIFFERENCE IN THE NEAR-SURFACE LAYER

The relationship between the factor of temperature difference of the near-surface layer(T_(1000 hPa)-T_(2m))and sea fog is analyzed using the NCEP reanalysis with a horizontal resolution of l°xl°(2000 to 2011) and the station observations(2010 to 2011).The element is treated as the prediction variable factor in the GRAPES model and used to improve the regional prediction of sea fog on Guangdong coastland.(1) The relationship between this factor and the occurrence of sea fog is explicit:When the sea fog happens,the value of this factor is always large in some specific periods,and the negative value of this factor decreases significantly or turns positive,suggesting the enhancement of warm and moist advection of air flow near the surface,which favors the development of sea fog.(2) The transportation of warm and moist advection over Guangdong coastland is featured by some stages and the jumping among these states.It also gets stronger over time.Meanwhile,the northward propagation of warm and moist advection is quite consistent with the northward advancing of sea fog from south to north along the coastland of China.(3) The GRAPES model can well simulate and realize the factor of near-surface temperature difference.Besides,the accuracy of regional prediction of marine fog,the relevant threat score and Heidke skill score are all improved when the factor is involved.

Spatial-temporal variations in near-surface soil freeze-thaw cycles in the source region of the Yellow River during the period 2002–2011 based on the Advanced Microwave Scanning Radiometer for the Earth Observing System(AMSR-E) data

Detecting near-surface soil freeze-thaw cycles in high-altitude cold regions is important for understanding the Earth＇s surface system, but such studies are rare. In this study, we detected the spatial-temporal variations in near-surface soil freeze-thaw cycles in the source region of the Yellow River（SRYR） during the period 2002–2011 based on data from the Advanced Microwave Scanning Radiometer for the Earth Observing System（AMSR-E）. Moreover, the trends of onset dates and durations of the soil freeze-thaw cycles under different stages were also analyzed. Results showed that the thresholds of daytime and nighttime brightness temperatures of the freeze-thaw algorithm for the SRYR were 257.59 and 261.28 K, respectively. At the spatial scale, the daily frozen surface（DFS） area and the daily surface freeze-thaw cycle surface（DFTS） area decreased by 0.08% and 0.25%, respectively, and the daily thawed surface（DTS） area increased by 0.36%. At the temporal scale, the dates of the onset of thawing and complete thawing advanced by 3.10（±1.4） and 2.46（±1.4） days, respectively; and the dates of the onset of freezing and complete freezing were delayed by 0.9（±1.4） and 1.6（±1.1） days, respectively. The duration of thawing increased by 0.72（±0.21） day/a and the duration of freezing decreased by 0.52（±0.26） day/a. In conclusion, increases in the annual minimum temperature and winter air temperature are the main factors for the advanced thawing and delayed freezing and for the increase in the duration of thawing and the decrease in the duration of freezing in the SRYR.