Development and testing of numerical methods for two-way nested air pollution modelling

A new 3-D model REGINA (REGIonal high resolutioN Air pollution model) is under development at the National Environmental Research Institute (NERI). The model is based on pieces from several models developed over the last decades at NERI. The aim of the work is to develop a nested model which can operate with very high resolution in both space and time. To fulfill this aim the choice and implementation of accurate numerical methods is crucial. The model will be applied for studying air pollution phenomena (both monitoring, forecasting and scenarios) over Denmark. The present paper is focussed on the modification, implementation and testing of a numerical method for treating the horizontal advection in the model as well as the implementation of two-way nesting techniques. The horizontal transport in the model is solved using an accurate space derivative algorithm. This method is traditionally implemented with periodic boundary conditions, however, this is not an option for nested modelling. A new method for calculating non-periodic boundary conditions has been developed in order to overcome this problem. Extensive testing of the numerical solution of the advection and the coupling of the solution of advection and chemistry in the model using Molenkamp–Crowley rotation tests have been carried out. The results show that the model with the current implementation of numerical methods is suitable for calculating air pollution levels with high resolution.     

A two-way nested-grid ocean-circulation model for the Meso-American Barrier Reef System

A two-way nested-grid ocean-circulation model is developed for the Meso-American Barrier Reef System (MBRS), using a newly developed two-way interactive nesting technique. The unique feature of this new nesting technique is its use of the semi-prognostic method (Sheng et al. 2001) to exchange information between different grids through the model momentum equations. The nested-grid model for the MBRS has a fine-resolution inner model embedded in a coarse-resolution outer model. The outer model is the western Caribbean Sea model developed by Sheng and Tang (2003), with a horizontal resolution of roughly 19 km. The inner model domain covers the northwest Caribbean Sea (NWCS) between 79°W and 89°W and between 15.5°N and 22°N, with a horizontal resolution of roughly 6 km. The nested-grid ocean model is initialized with the January mean temperature and salinity and forced by the monthly mean COADS (comprehensive ocean-atmosphere data set) wind stress and surface heat flux. The model sea-surface salinity is restored to the monthly mean climatology. The nested-grid model is integrated for 2 years and the second-year model results are presented in this paper. The model-calculated annual-mean near-surface currents over the NWCS agree reasonably well with the time-mean near-surface currents inferred by Fratantoni (2001) from trajectories of the satellite-tracked 15-m drogued drifters in the 1990s. The two-way nested model is also used to quantify the role of local wind stress, local density gradients and boundary forcings of the outer model in driving the annual-mean circulation in the region.Responsible Editor: Phil Dyke     

台湾及台湾海峡地震活动对大陆地区的影响

根据已有和近年来补充调查的资料,阐明了台湾及台湾海峡地震对大陆东南沿海地区影响的差异,并通过对以下问题的分析,探讨了其结果在震害防御与历史地震考证方面的意义:①台湾地区的强震频度较高,间隔时间从几年到数十年,但对大陆沿海地带的最大影响烈度仅为Ⅵ度,最大有感半径达1100km,有感地带的宽度大于500km．这些地震在大陆地区可产生数毫米位移的长周期地面运动;②台湾海峡西侧滨海断裂带发生的强震,其频度不如台湾地区高,却因距离大陆较近,最大影响烈度达Ⅷ——Ⅸ度,最大有感半径大于1000km,有感地带的宽度大于400km;③我国东南沿海地带的震害防御策略,宜在着重考虑滨海断裂带强震影响的同时,还要注意台湾地震对工程结构造成的损伤,以及恐震心理引发的哄动社会效应影响;④大陆地区记载到的1517年5月19日地震事件,宜视为台湾地区的强震影响更为妥切．      

Study on the transverse structures across Taiwan Strait

The mountain ranges are almost NE- or NNE-trending in mainland and Taiwan island across the strait. But a number of correspondent transverse structures across Taiwan Strait have been identified on the basis of morphology, surface structure, deep seismic survey and gravity field etc. Moreover, it is inferred from the square or/and rectangular shape of the isoseist maps of 8 M≥ 6 earthquakes occurring in the western foothill plain of Taiwan that there is an apparent segmentation of the Taiwan Island along the N-S orientation. This has been verified from the GPS’s results of the ground motions measured during the Chi-Chi earthquake in 1999. As a matter of fact, the transverse structures within the Taiwan Island could be correlated with more salient NW- and EW-trending structures in the mainland of Fujian than in Taiwan. The investigation of the transverse structures across Taiwan Strait is conducive to realizing the geology within the strait. Hence, the tectonic sketch map of the strait as well as mainland and island drawn in this context shows NW- and EW-trending structures for the study area.     

A particle-tracking model for simulating pollutant dispersion in the Strait of Gibraltar

A particle-tracking model to simulate the dispersion of contaminants in the Strait of Gibraltar has been developed. The model solves the hydrodynamic equations off-line and tidal analysis is carried out to determine tidal constants for the two main constituents. Tidal constants and residuals are stored in files that are read by the dispersion model. A lagrangian approach is used to solve dispersion; diffusion and decay are simulated by a Monte Carlo method. A method for assessing the areas of the Strait with higher probability of being affected by contamination occurring after an accident in the shipping routes is given. Generally speaking, the fate of a pollutant discharge strongly depends on wind conditions. Winds from the east tend to retain contamination into the Strait. As a consequence, transverse mixing occurs and both Spain and Morocco coasts are affected by contamination. Under calm conditions and west winds, contaminants are flushed out of the Strait faster and transverse mixing does not occur. Thus, only part of Morocco coast has a higher probability of being affected by contamination.     

Interannual variability of summer coastal upwelling in the Taiwan Strait

This study dealt with the interannual variability of summer coastal upwelling in the Taiwan Strait, based on empirical orthogonal function (EOF) analysis. Three datasets were used for the analysis: the National Oceanic and Atmospheric Administration (NOAA), Advanced Very High Resolution Radiometer (AVHRR) sea surface temperature dataset from 1985 to 2005; hydrographic records at two coastal stations from 1970 to 2001; and cruise measurements in 1988 and 2004. The results indicated that the first mode (85.3%) of the spatial variance showed a persistent front, which was generally aligned northeast–southwestward in the western Taiwan Strait. This front separated colder water on the west side from warmer water on the east side. The eigenvector time series showed that the variability of this front with time was closely correlated with the change in the wind stress anomaly of the alongshore wind component, derived from 17 years of the European remote sensing (ERS) satellite and QuickScat wind dataset from 1992 to 2005. Records of water temperature and salinity anomaly at Pingtan Island (Is.) located in the northwestern Taiwan Strait, and Dongshan Is. located in the southwestern Taiwan Strait, showed that a negative temperature anomaly appeared along with a positive salinity anomaly in some years. This suggested a dominant influence of cold and saline upwelling water at the surface. The years of notable cooling events derived from the station measurements were generally consistent with the time series of the EOF Mode 1. The change in upwelling derived from cruise measurements further confirmed the change shown by the EOF Mode 1 time series. These 1985–2005 results indicated that for the entire western Taiwan Strait summer coastal upwelling was strong in 1987, 1993, and 1998, and that upwelling in the northwest and the southwest Taiwan Strait showed different behavior. A delayed ENSO (El Niño Southern Oscillation) effect was suggested as a major mechanism for the interannual variability of Taiwan Strait coastal upwelling.     

A coupled tide-wave model for the NW European shelf seas

Understanding the interaction of tides and waves is essential in many studies, including marine renewable energy, sediment transport, long-term seabed morphodynamics, storm surges and the impacts of climate change. In the present research, a COAWST model of the NW European shelf seas has been developed and applied to a number of physical processes. Although many aspects of wave–current interaction can be investigated by this model, our focus is on the interaction of barotropic tides and waves at shelf scale. While the COWAST model was about five times more computationally expensive than running decoupled ROMS (ocean model) and SWAN (wave model), it provided an integrated modelling system which could incorporate many wave–tide interaction processes, and produce the tide and wave parameters in a unified file system with a convenient post-processing capacity. Some applications of the model such as the effect of tides on quantifying the wave energy resource, which exceeded 10% in parts of the region, and the effect of waves on the calculation of the bottom stress, which was dominant in parts of the North Sea and Scotland, during an energetic wave period are presented, and some challenges are discussed. It was also shown that the model performance in the prediction of the wave parameters can improve by 25% in some places where the wave-tide interaction is significant.     

台湾海峡架桥铺路的层析技术问题的初步分析

关于在台湾海峡架桥和修筑隧道的问题从1987年至2002年在福省召开过三次"台湾海峡两岸地质地震和架桥与修筑隧道的研讨会",以后陆续出版了相应的论文和论文集,提出了"台湾海峡架桥和修筑隧道的可行性研究",和"台湾海峡隧道工程设计总建议"等[1-4],论述了修建该工程的可行性等问题.提出平潭岛可作为海峡北线的端点,探讨了如何具体建桥和相关工程的问题.近年来地球层析技术的发展,可对台湾海峡隧道工程的工程地质问题,海港的开发和沿海陆地的充分利用等,做出探测台湾海峡地下地层三维成像细结构和地球物理场的空间分布,作为台湾海峡架桥铺路工程基本设计的必要依据.为了让大中小轮船在台湾海峡顺利通航,利用现代的科技研讨台湾海峡两岸连接的最佳方案.为了架成大型的高桥和修建堤坝式的地基,必须用地震层析和应力波等方法探测桥梁和堤坝地基的结构和地层的坚实程度,以便决定桥墩深入海底的深度,和设计坚固而又经济的桥墩与堤坝的基础.参考开采海底石油平台施工和荷兰填海造田的工程经验,以双边海堤与大型高桥及活动的船闸相结合,从平潭岛东至台湾新竹西海岸,在124千米距离上架设4-6段桥梁和5-7段堤坝型的路基,既让大中小船舶照常航行,海洋生物照常游动,又可在其上通车,以缩短两岸距离,适应中国大陆和台湾经济快速发展的需要,争取和平沟通的愿望.由于架桥铺路的工程费用较修筑隧道的投资低得多,在福建近海地带局部填海造地,充分利用海洋资源和建筑深水码头,可得到较大的经济补偿,在现阶段较易实施.本文初步分析了有关的地球层析和工业CT技术的作用问题.     

A regional model of low flow for southern Taiwan

The primary purpose of this study is to develop regional models of the lower part of flow duration curves (LPFDCs) to synthesize low‐flow characteristics at ungauged sites in southern Taiwan. Because of the close relationship between low streamflow regimes and hydrogeological features, the model development first involved delimiting homogeneous hydrogeological regions by using two‐step cluster analysis. Each homogeneous region was then discriminated by an equation developed on the basis of its hydrogeological features, which was then used to determine which of three sets of regional LPFDC models would be appropriate for a particular ungauged site. Each of the three sets of regional LPFDC models were developed using both conventional multivariate statistical regression and fuzzy regression. Thirty‐four stream‐gauged watersheds located in southern Taiwan provide the data set. The study results reveal that the regional LPFDC models developed in this study could be applied reasonably at ungauged sites. Copyright © 2002 John Wiley & Sons, Ltd.     

A coupled physical-biochemical lake model for forecasting water quality

A new one-dimensional numerical model that includes physical and biochemical processes has been developed. The biochemical processes, influenced by the lake dynamics, are required for forecasting water quality. The model is used to investigate the effects of different internal restoration measures, such as artificial mixing, input of oxygen and drainage of deep water.The model is applied to the Northern Basin of Lake Lugano, a Swiss-Italian border lake. The lake is highly eutrophic and chemically stratified throughout the year. The model was calibrated over one year and validated over a period of several years. The results agree well with the measured data. The coupled model reproduces the observed depth dependency of conductivity even during long simulation times. Due to the predominant mixing, decoupled physical models cannot maintain such gradients. The forecasting capabilities of the model are demonstrated for different case studies. The impact of restoration measures on water quality is rather small. Best results are achieved by reducing the external nutrient loading. Caution is recommended for internal measures as these have to be studied in greater detail.     

台湾海峡南部一次5.0级地震的sPn震相分析

运用Pn震相对齐测定sPn震相的方法分析了发生在台湾海峡南部的一次5．0级地震的 sPn震相,初步了解sPn震相的一些特征,并在不同地壳模型下用sPn-Pn的走时差计算了该地震的震源深度。结果表明:在该地震中,sPn震相特征明显;用Pn震相对齐测定sPn震相的方法可以快速、可靠地测定出sPn震相;用sPn-Pn的走时差计算震源深度时,不同模型计算出的结果相差不大。     

1994年台湾海峡7.3级地震发震断层的地震学证据

通过地震学观测资料的分析和研究,并综合其他作者的研究成果,论证了台湾海峡１９９４年７．３级地震发震断层的主破裂面走向应为北西向,地震的发生是由分布于该区域的巴士系构造活动所致。     

Generation and evolution of internal solitary waves in the southern Taiwan Strait

ABSTRACT

The generation processes and potential energy sources of internal solitary waves (ISWs) in the southern Taiwan Strait are investigated by driving a high resolution non-hydrostatic numerical model with realistic background conditions. Two main types of ISWs are clarified according to their different energy sources. One is generated by the nonlinear disintegration of remote internal tides emanating from Luzon Strait, and the other type is generated by local tide-topography interaction at the continental slope. The basic properties and evolution processes differ between these two kinds of ISWs. The waves originated from the remote internal tides at Luzon Strait have amplitudes comparable to previous field observations. In contrast, the ISWs generated locally are much weaker than observed waves, even in the presence of a steady offshore background current, which intensifies the generation of onshore ISWs. The ISWs induced by remotely generated M₂ internal tides are stronger than those induced by K₁ internal tides, and the fraction of internal wave energy transmitted onto the shelf is not significantly influenced by the intensity of remotely generated internal tides.     

The role of coastal-trapped waves on the 2008 cold disaster in the Taiwan Strait

In early 2008, cold water in the Taiwan Strait (TWS) was moved sequentially by a cross-strait flow and a southward flow to the Penghu Island, causing a cold-related fish kill disaster. Except for the local wind forcing, the coastal-trapped waves (CTWs), intermittently propagating toward the TWS from north in winter, are an additional factor that could impact the flow patterns by changing cross-strait sea-level gradient during the disaster. In the first stage (January 28–February 7), the reach of a large CTW trough induced an additional northward flow, which formed a cyclone after turning around the Zhangyun Ridge. Then, the cyclone led to an additional cross flow, which enhanced an eastward (offshore) movement of cold water. In the second stage (February 7–14), the arrival of a large CTW crest triggered an additional southward flow, which intensified a southward movement of the cold water. Due to the additional eastward and southward movements caused by the CTWs, the cold water could reach Penghu Island inducing a cold disaster.     

A 2.5D coupled FE-BE model for the prediction of railway induced vibrations

Ground vibrations induced by railway traffic at grade and in tunnels are often studied by means of two-and-half dimensional (2.5D) models that are based on a Fourier transform of the coordinate in the longitudinal direction of the track. In this paper, the need for 2.5D coupled finite element-boundary element models is demonstrated in two cases where the prediction of railway induced vibrations is considered. A recently proposed novel 2.5D methodology is used where the finite element method is combined with a boundary element method, based on a regularized boundary integral equation. In the formulation of the boundary integral equation, Green's functions of a layered elastic halfspace are used, so that no discretization of the free surface or the layer interfaces is required. In the first case, two alternative models for a ballasted track on an embankment are compared. In the first model, the ballast and the embankment are modelled as a continuum using 2.5D solid elements, whereas a simplified beam representation is adopted in the second model. The free field vibrations predicted by both models are compared to those measured during a passage of the TGVA at a site in Reugny (France). A very large difference is found for the free field response of both models that is due to the fact that the deformation of the cross section of the embankment is disregarded in the simplified representation. In the second case, the track and free field response due to a harmonic load in a tunnel embedded in a layered halfspace are considered. A simplified methodology based on the use of the full space Green's function in the tunnel–soil interaction problem is investigated. It is shown that the rigorous finite element-boundary element method is required when the distance between the tunnel and the free surface and the layer interfaces of the halfspace is small compared to the wavelength in the soil.     

台湾海峡大容量气枪震源海陆联测初探

本文利用在我国台湾海峡采用大容量气枪震源开展海陆联测获得的广角地震测线HX9, 采用二维射线追踪法反演得到了HX9剖面的地壳二维速度结构和地壳界面形态, 初步探明了福建—台湾海峡海陆过渡带的深部构造． 结果表明: HX9剖面的地壳内存在两个速度间断面, 即C界面和莫霍面, 其中: C界面为上、 下地壳的分界面, 是一个小的速度不连续面, 速度变化值达0.08—0.16 km/s; 而地壳底部的莫霍面则有较大的速度反差, 变化值达1.02—1.29 km/s, 莫霍面上、 下的速度分别为6.75—6.97 km/s和8.00—8.07 km/s． 沿剖面的地壳界面形态总体起伏不大, 陆域上、 下地壳的厚度和界面变化趋势均相似, 从陆域到海域呈微倾斜变化趋势, 表现为减薄陆壳的特征． 莫霍面陆域埋深约为31.6 km, 向福建东南沿海逐渐减薄至27.4 km左右．      

台湾海峡中部障碍层日变化及其局地反馈

利用“延平2号”科学调查船于2005年6月下旬在台湾海峡中部（119．2°E,24．3°N）实施的定点连续观测资料,研究了夏季风期间障碍层的日变化过程及其局地反馈特征．障碍层开始形成时,由于非太阳辐射局限于上混合层,而太阳辐射可以穿透混合层底直至温跃层,导致障碍层内累积了最多热能,形成逆温层．随着障碍层内部热量的不断累积,障碍层与上混合层之间的热交换过程逐渐增强,上混合层显著增暖。此外,障碍层对应浮力频率极小值区,具有较弱的层结稳定性;针对海气通量的分析表明障碍层造成了显著的局地感热、潜热通量异常．进一步分析表明,海表风搅拌产生的湍动能在向下传输过程中,次盐跃层和温跃层构成了两个界面,大部分湍动能被较浅的盐度层结阻挡,限制在混合层以上,部分穿透混合层底的湍动能再次被温跃层阻挡,障碍层对应湍动能的薄弱区．     

Diurnal evolution of the barrier layer and its local feedback in the central Taiwan Strait

Diurnal evolution of the barrier layer (BL) and its local feedback features in the central Taiwan Strait (119.2oE,24.3oN) during summertime monsoon are investigated using in situ moored observations conducted by the Yan-Ping 2 research vessel in late June 2005.During the initiation phase,for the non-solar radiation tends to be trapped in the upper mixed layer,whereas the solar radiation can penetrate deeply through the mixed layer approaching the thermocline,most heat is accumulated inside the BL inducing a...     

An intermediate coupled model for the tropical ocean-atmosphere system

An intermediate ocean-atmosphere coupled model is developed to simulate and predict the tropical interannual variability. Originating from the basic physical framework of the Zebiak-Cane(ZC) model, this tropical intermediate couple model(TICM) extends to the entire global tropics, with a surface heat flux parameterization and a surface wind bias correction added to improve model performance and inter-basin connections. The model well reproduces the variabilities in the tropical Pacific and Indian basins. The simulated El Ni?o-Southern Oscillation(ENSO) shows a period of 3–4 years and an amplitude of about 2°C, similar to those observed. The variabilities in the Indian Ocean, including the Indian Ocean basin mode(IOBM) and the Indian Ocean Dipole(IOD), are also reasonably captured with a realistic relationship to the Pacific. However, the tropical Atlantic variability in the TICM has a westward bias and is overly influenced by the tropical Pacific. A 47-year hindcast experiment using the TICM for the period of 1970–2016 indicates that ENSO is the most predictable mode in the tropics. Skillful predictions of ENSO can be made one year ahead, similar to the skill of the latest version of the ZC model, while a "spring predictability barrier" still exists as in other models. In the tropical Indian Ocean, the predictability seems much higher in the west than in the east. The correlation skill of IOD prediction reaches 0.5 at a 5-month lead, which is comparable to that of the state-of-the-art coupled general circulation models. The prediction of IOD shows a significant "winter-spring predictability barrier", implying combined influences from the tropical Pacific and the local sea-air interaction in the eastern Indian Ocean. The TICM has little predictive skill in the equatorial Atlantic for lead times longer than 3 months, which is a common problem of current climate models badly in need of further investigation.