港湾中海啸流的模拟及其时空演变特征对输入条件不确定性的响应分析

一直以来,海啸波特征作为表征海啸潜在破坏性的参数指标得到了广泛应用,特别是针对近场极端海啸事件造成的灾害来说,这种表征具有较好的适用性.然而总结分析历史海啸事件造成的损失发现：在远场近岸及港湾系统中,海啸诱导的强流却是造成损失的主要原因.陆架或港湾振荡导致海啸波幅快速升降诱发强流,可能促使港工设施受到威胁及损害,进而对海啸预警服务及海事应急管理提出了新的挑战.因此,全面理解与评估海啸在港湾中诱发的灾害特征,探索港湾中海啸流的数值模拟方法,发展针对港湾尺度的海啸预警服务指导产品尤为迫切.受限于海啸流验证数据的缺乏及准确模拟海啸流技术方法的诸多不确定性,大部分海啸数值模拟研究工作主要是针对水位特征的研究及验证,可能导致对港湾中海啸灾害危险性认识的曲解与低估.本研究基于非线性浅水方程,针对夏威夷群岛三个典型港湾建立了精细化海啸数值模型（空间分辨率达到10 m）,并联合有限断层破裂模型计算分析了日本东北地震海啸在三个港湾及其邻近区域的海啸特征,波、流计算结果与实测结果吻合较好,精细化的海啸港湾模型模拟结果可信.模拟发现港湾中较小的波幅,同样可以产生强流.综合分析日本东北地震海啸波、流特征对输入条件不确定性的响应结果发现：港湾中海啸波-流能量的空间分布特征差异较大,这与港湾系统中海啸波的驻波特性相关;相比海啸波幅空间特征,海啸流特征具有更强的空间敏感性;海啸流时空分布特征对输入条件的不确定性响应比海啸波幅对这些不确定性的响应更强,海啸流的模拟与预报更有挑战性;不确定性对海啸流计算精度的影响会进一步传导放大港湾海啸流危险性的评估及对港工设施产生的应力作用的误差,合理的输入条件对海啸流的精确模拟至关重要.最后,希望通过本文的研究可以从海啸波-流特征角度更加全面认识近岸海啸灾害特征,拓展海啸预警服务的广度与深度,从而为灾害应急管理部门提供更加科学合理的辅助决策产品.     

Sub-inertial characteristics of the surface flow field over the shelf of the central Mid-Atlantic Bight

Observations of surface velocity data from August 2002 to February 2004 were collected by a series of four long-range high-frequency (HF) radars along the coast of New Jersey. The shelf observations of the central Mid-Atlantic Bight (MAB) were compared to historical observations of surface flow characteristics in the area. The time-averaged spatial mean velocity of 4 cm/s in the down-shelf along-shelf direction and 3 cm/s in the offshore across-shelf direction compared very well to historical surface measurements in the study region. However, as the spatial resolution of the data set revealed, this simple measure masked significant spatial variations in the overall and seasonal mean flow structures. Three regions – the south bank of the Hudson Shelf Valley, the southern New Jersey inner shelf (LEO-15) region, and the region offshore of the Delaware Bay mouth (southwest corner) – had mean flows that favor offshore transport of surface water. In terms of temporal variability, maps of the principle axes showed that the across-shelf (minor) axis contribution was not insignificant in the surface layer ranging from 0.3 to 0.9 of along-shelf (major) axis and that there were seasonal differences in orientation and ellipticity. Analysis of the spatial changes in the temporal and spatial correlation scales over the shelf showed that shelf position, in addition to site separation, contributed to the differences in these properties. Furthermore, observations over the Hudson Shelf Valley region suggested that this was a region of transition in which the orientation of along- and across-shelf components begin to change.     

日本东北M_W9.0地震海啸在港池及邻近区域诱发的涡流危险性计算与评估分析

海啸造成的灾害与损失并非都与淹没有关,特别是港口中海啸诱导的强流会对船只及海事设施产生重要的影响及损害.由于海啸流观测数据稀缺及海啸诱导涡流机制的不确定性,过去60年海啸科学主要集中于对海啸波特征及淹没过程的研究与分析,海啸流模拟及验证工作开展较少,导致对海啸流基本特征及其造成灾害现象的曲解.开展海啸诱导的涡流研究及预警服务显得尤为重要及紧迫.考虑快速海啸预警需要,综合对比海啸诱导涡流的物理框架及模型方法,探索兼顾效率与计算精度的海啸流模拟方法是本文的核心工作及出发点.通过分析浅层湍流相干结构(TCS)产生的主要物理耗散机制,确定了考虑2D水平耗散机制的非线性浅水方程可用于海啸涡流的模拟分析.基于高精、高分辨率有限体积模型Geoclaw建立了三个精细化的港口海啸流模型,模型分辨率为5m.利用基于海啸浮标反演的海啸源模型作为初始条件,模拟分析了日本东北地震海啸在远场的海啸波流特征.海啸波流特征模拟结果与观测吻合较好,结果可信.对比发现:波驱动的自由表面流,小的位相或波幅误差就会导致大的流速误差,流的模拟和预报相对波幅来说更具挑战性.研究了海啸波流能量在港池中的分布特征,得到:港池入口及防波堤两端常被强流控制,具有极高的危险性;相对于波幅的空间变化,海啸流具有更强的空间敏感性;所建立的高分辨率海啸模型模拟再现了日本海啸在近场的涡旋结构,给出了与观测基本一致的涡流特征.最后,引入海啸流危险等级标准,分析了港口海啸流危险性等级分布、船只疏散的安全深度及回港的时间周期.针对港口、海湾同时考虑海啸波流特征的海啸预警与评估对于港口应急管理者科学决策具有重要意义.     

Common-volume measurements of mesospheric winds: 2. Small scale structure implications

Two common volume experiments were conducted in August 1996 and July 1997 between the Durham meteor wind radar (43.1°N, 70.9°W) and the Millstone Hill incoherent scatter radar, (42.6°N, 71.5°W) to compare the techniques in measuring neutral winds at an altitude of 100 km. For this comparison the vertical winds are assumed to be zero and only the horizontal components of the radar line of sight velocities are used. Analysis of the data reveals overall general agreement, but some large deviations in the wind components are observed at specific times and these are examined closely. Error analysis of the radar measurements is presented here, and emphasis is placed on the careful delineation of the effect of spatial variations in the wind field observed by the two radars. Since the spatial resolution of both radars is<3 km both horizontally and vertically, some of the three dimensional properties of the horizontal wind component can be estimated. For the incoherent scatter radar with its narrow steerable beam, the spatial location of the sampling points could be chosen; however, finer sampling of the wind field results in more temporal smearing due to the fixed measurement time for each point. For the meteor radar the spatial sample points occur randomly within the system beam since they depend on the chance location of observed meteor trails. Both systems spatially undersample the wind field in most cases, but with careful consideration of the system errors for both radars, it is shown that small scale (∼10 km) wind variations must exist at these altitudes with rms velocity differences of ∼25 m/s.     

Tsunami detection by high-frequency radar in British Columbia: performance assessment of the time-correlation algorithm for synthetic and real events

The authors recently proposed a new method for detecting tsunamis using high-frequency (HF) radar observations, referred to as “time-correlation algorithm” (TCA; Grilli et al. Pure Appl Geophys 173(12):3895–3934, 2016a, 174(1): 3003–3028, 2017). Unlike standard algorithms that detect surface current patterns, the TCA is based on analyzing space-time correlations of radar signal time series in pairs of radar cells, which does not require inverting radial surface currents. This was done by calculating a contrast function, which quantifies the change in pattern of the mean correlation between pairs of neighboring cells upon tsunami arrival, with respect to a reference correlation computed in the recent past. In earlier work, the TCA was successfully validated based on realistic numerical simulations of both the radar signal and tsunami wave trains. Here, this algorithm is adapted to apply to actual data from a HF radar installed in Tofino, BC, for three test cases: (1) a simulated far-field tsunami generated in the Semidi Subduction Zone in the Aleutian Arc; (2) a simulated near-field tsunami from a submarine mass failure on the continental slope off of Tofino; and (3) an event believed to be a meteotsunami, which occurred on October 14th, 2016, off of the Pacific West Coast and was measured by the radar. In the first two cases, the synthetic tsunami signal is superimposed onto the radar signal by way of a current memory term; in the third case, the tsunami signature is present within the radar data. In light of these test cases, we develop a detection methodology based on the TCA, using a correlation contrast function, and show that in all three cases the algorithm is able to trigger a timely early warning.     

Assessing the fidelity of surface currents from a coastal ocean model and HF radar using drifting buoys in the Middle Atlantic Bight

The rapid expansion of urbanization along the world’s coastal areas requires a more comprehensive and accurate understanding of the coastal ocean. Over the past several decades, numerical ocean circulation models have tried to provide such insight, based on our developing understanding of physical ocean processes. The systematic establishment of coastal ocean observation systems adopting cutting-edge technology, such as high frequency (HF) radar, satellite sensing, and gliders, has put such ocean model predictions to the test, by providing comprehensive observational datasets for the validation of numerical model forecasts. The New York Harbor Observing and Prediction System (NYHOPS) is a comprehensive system for understanding coastal ocean processes on the continental shelf waters of New York and New Jersey. To increase confidence in the system’s ocean circulation predictions in that area, a detailed validation exercise was carried out using HF radar and Lagrangian drifter-derived surface currents from three drifters obtained between March and October 2010. During that period, the root mean square (RMS) differences of both the east–west and north–south currents between NYHOPS and HF radar were approximately 15 cm s^?1. Harmonic analysis of NYHOPS and HF radar surface currents shows similar tidal ellipse parameters for the dominant M₂ tide, with a mean difference of 2.4 cm s^?1 in the semi-major axis and 1.4 cm s^?1 in the semi-minor axis and 3° in orientation and 10° in phase. Surface currents derived independently from drifters along their trajectories showed that NYHOPS and HF radar yielded similarly accurate results. RMS errors when compared to currents derived along the trajectory of the three drifters were approximately 10 cm s^?1. Overall, the analysis suggests that NYHOPS and HF radar had similar skill in estimating the currents over the continental shelf waters of the Middle Atlantic Bight during this time period. An ensemble-based set of particle tracking simulations using one drifter which was tracked for 11 days showed that the ensemble mean separation generally increases with time in a linear fashion. The separation distance is not dominated by high frequency or short spatial scale wavelengths suggesting that both the NYHOPS and HF radar currents are representing tidal and inertial time scales correctly and resolving some of the smaller scale eddies. The growing ensemble mean separation distance is dominated by errors in the mean flow causing the drifters to slowly diverge from their observed positions. The separation distance for both HF radar and NYHOPS stays below 30 km after 5 days, and the two technologies have similar tracking skill at the 95 % level. For comparison, the ensemble mean distance of a drifter from its initial release location (persistence assumption) is estimated to be greater than 70 km in 5 days.     

Surface circulation downstream of the Point Arena upwelling center

Measurements from recently installed 5 MHz high-frequency radar (CODAR) stations south of Point Arena, California, are used to describe surface current patterns during the upwelling season (June-August 2007). The systems provide hourly current maps on a 5-km grid, covering a region from approximately 10 to 150 km offshore (the continental shelf into the deep ocean). These HF-radar observations provide an unprecedented view of circulation in this “coastal transition zone”, between the wind-driven circulation over the shelf and the California Current circulation offshore. Circulation patterns include: (1) bifurcation of the coastal upwelling jet downstream of Point Arena into an along-shelf (down-coast) branch and an offshore branch, and (2) a large-scale anticyclonic meander that often develops into an eddy-like recirculation south of the bifurcation. The “recirculation” feature extends well offshore, with surface currents 50-100 km from the coast consistently opposing the wind stress. The spatial and temporal evolution of the surface current features during upwelling events affects surface transport from Point Arena to areas in the south, increasing the travel time of a substantial fraction of newly upwelled water from a few days to roughly two weeks. Thus, surface currents even far offshore influence coastal transport of nutrients, phytoplankton and larvae on ecologically relevant timescales, with resultant connectivity patterns very different than implied by a simple examination of the mean flow.     

Wind-speed inversion from HF radar first-order backscatter signal

Land-based high-frequency (HF) radars have the unique capability of continuously monitoring ocean surface environments at ranges up to 200 km off the coast. They provide reliable data on ocean surface currents and under slightly stricter conditions can also give information on ocean waves. Although extraction of wind direction is possible, estimation of wind speed poses a challenge. Existing methods estimate wind speed indirectly from the radar derived ocean wave spectrum, which is estimated from the second-order sidebands of the radar Doppler spectrum. The latter is extracted at shorter ranges compared with the first-order signal, thus limiting the method to short distances. Given this limitation, we explore the possibility of deriving wind speed from radar first-order backscatter signal. Two new methods are developed and presented that explore the relationship between wind speed and wave generation at the Bragg frequency matching that of the radar. One of the methods utilizes the absolute energy level of the radar first-order peaks while the second method uses the directional spreading of the wind generated waves at the Bragg frequency. For both methods, artificial neural network analysis is performed to derive the interdependence of the relevant parameters with wind speed. The first method is suitable for application only at single locations where in situ data are available and the network has been trained for while the second method can also be used outside of the training location on any point within the radar coverage area. Both methods require two or more radar sites and information on the radio beam direction. The methods are verified with data collected in Fedje, Norway, and the Ligurian Sea, Italy using beam forming HF WEllen RAdar (WERA) systems operated at 27.68 and 12.5 MHz, respectively. The results show that application of either method requires wind speeds above a minimum value (lower limit). This limit is radar frequency dependent and is 2.5 and 4.0 m/s for 27.68 and 12.5 MHz, respectively. In addition, an upper limit is identified which is caused by wave energy saturation at the Bragg wave frequency. Estimation of this limit took place through an evaluation of a year long database of ocean spectra generated by a numerical model (third generation WAM). It was found to be at 9.0 and 11.0 m/s for 27.68 and 12.5 MHz, respectively. Above this saturation limit, conventional second-order methods have to be applied, which at this range of wind speed no longer suffer from low signal-to-noise ratios. For use in operational systems, a hybrid of first- and second-order methods is recommended.     

Validation of the CUTLASS HF radar gravity wave observing capability using EISCAT CP-1 data

Quasi-periodic fluctuations in the returned ground-scatter power from the SuperDARN HF radars have been linked to the passage of medium-scale gravity waves. We have applied a technique that extracts the first radar range returns from the F-region to study the spatial extent and characteristics of these waves in the CUTLASS field-of-view. Some ray tracing was carried out to test the applicability of this method. The EISCAT radar facility at Tromsø is well within the CUTLASS field-of-view for these waves and provides a unique opportunity to assess independently the ability of the HF radars to derive gravity wave information. Results from 1st March, 1995, where the EISCAT UHF radar was operating in its CP-1 mode, demonstrate that the radars were in good agreement, especially if one selects the electron density variations measured by EISCAT at around 235 km. CUTLASS and EISCAT gravity wave observations complement each other; the former extends the spatial field of view considerably, whilst the latter provides detailed vertical information about a range of ionospheric parameters.     

Characterization of ocean surface current properties from single site HF/VHF radar

Surface current mapping from HF/VHF coastal radars traditionally requires at least two distant sites. Vector velocities are estimated by combining the radial velocity components measured by the radars. In many circumstances (e.g., failures, interferences, logistics constraints), such a combination is not possible by lack of data from one station. Two methods are evaluated to get information on surface circulation from a single site radar: the Vectorial Reconstruction Method (VRM) for current vector mapping and the Vortex Identification Method (VIM) for detecting eddy-like structures. The VRM assumes a non-divergent horizontal surface current, and the VIM analyzes radial velocities and their radial and orthoradial gradients. These two methods are tested on modeled and measured data sets in the Northwestern Mediterranean Sea where both high-resolution ocean circulation model and radar campaigns are available. The VRM performance is strongly limited by the divergence-free hypothesis which was not satisfied in our real data. The VIM succeeded in detection of vortex in the Gulf of Lions and from an operating single site radar located on the Provence coasts in summer.     

Erratum to Identification of tsunami-induced deposits using numerical modeling and rock magnetism techniques: A study case of the 1755 Lisbon tsunami in Algarve,Portugal: [Physics of the Earth and Planetary Interiors Volume 182, Issues 3–4, Pages 187–198 (October 2010)]

Storm- and tsunami-deposits are generated by similar depositional mechanisms making their discrimination hard to establish using classic sedimentologic methods. Here we propose an original approach to identify tsunami-induced deposits by combining numerical simulation and rock magnetism. To test our method, we investigate the tsunami deposit of the Boca do Rio estuary generated by the 1755 earthquake in Lisbon which is well described in the literature. We first test the 1755 tsunami scenario using a numerical inundation model to provide physical parameters for the tsunami wave. Then we use concentration (MS, SIRM) and grain size (χ_ARM, ARM, B1/2, ARM/SIRM) sensitive magnetic proxies coupled with SEM microscopy to unravel the magnetic mineralogy of the tsunami-induced deposit and its associated depositional mechanisms. In order to study the connection between the tsunami deposit and the different sedimentologic units present in the estuary, magnetic data were processed by multivariate statistical analyses. Our numerical simulation show a large inundation of the estuary with flow depths varying from 0.5 to 6 m and run up of ～7 m. Magnetic data show a dominance of paramagnetic minerals (quartz) mixed with lesser amount of ferromagnetic minerals, namely titanomagnetite and titanohematite both of a detrital origin and reworked from the underlying units. Multivariate statistical analyses indicate a better connection between the tsunami-induced deposit and a mixture of Units C and D. All these results point to a scenario where the energy released by the tsunami wave was strong enough to overtop and erode important amount of sand from the littoral dune and mixed it with reworked materials from underlying layers at least 1 m in depth. The method tested here represents an original and promising tool to identify tsunami-induced deposits in similar embayed beach environments.     

High spatial and temporal resolution observations of an impulse-driven field line resonance in radar backscatter artificially generated with the Tromsø heater

The CUTLASS Finland HF radar has been operated in conjunction with the EISCAT Tromsø RF ionospheric heater facility to examine a ULF wave characteristic of the development of a field line resonance (FLR) driven by a cavity mode caused by a magnetospheric impulse. When the heater is on, striating the ionosphere with field-aligned ionospheric electron density irregularities, a large enough radar target is generated to allow post-integration over only 1 second. When combined with 15 km range gates, this gives radar measurements of a naturally occurring ULF wave at a far better temporal and spatial resolution than has been achieved previously. The time-dependent signature of the ULF wave has been examined as it evolves from a large-scale cavity resonance, through a transient where the wave period was latitude-dependent and the oscillation had the characteristics of freely ringing field lines, and finally to a very narrow, small-scale local field line resonance. The resonance width of the FLR is only 60 km and this is compared with previous observations and theory. The FLR wave signature is strongly attenuated in the ground magnetometer data. The characterisation of the impulse driven FLR was only achieved very crudely with the ground magnetometer data and, in fact, an accurate determination of the properties of the cavity and field line resonant systems challenges the currently available limitations of ionospheric radar techniques. The combination of the latest ionospheric radars and facilities such as the Tromsø ionospheric heater can result in a powerful new tool for geophysical research.     

High spatial and temporal resolution observations of the ionospheric cusp

The Halley PACE HF radar has been operated in a new mode to provide very high time (10 s) and space (15 km) resolution measurements of the iono-spheric signatures of the cusp and the low-latitude boundary layer. The first data show that the iono-spheric signature of flux transfer events occur up to 300 km equatorward of regions showing the HF characteristics of the ionospheric cusp. Whilst larger flux transfer events are seen, on average, every 7 min, many much smaller and short-duration events have been identified. On one occasion DMSP data have been used to show that at least four flux transfer events are occurring simultaneously at the edge of the cusp over 2 h of MLT. There is strong, but not conclusive evidence, that reconnection at the magnetopause is both intermittent and patchy. These data also suggest that flux transfer events can be a significant contributor to the cross-polar cap potential.     

CUTLASS observations of a high-m ULF wave and its consequences for the DOPE HF Doppler sounder

The CUTLASS (Co-operative UK Twin Located Auroral Sounding System) Finland HF radar, whilst operating in a high spatial and temporal resolution mode, has measured the ionospheric signature of a naturally occurring ULF wave in scatter artificially generated by the Tromsø Heater. The wave had a period of 100 s and exhibited curved phase fronts across the heated volume (about 180 km along a single radar beam). Spatial information provided by CUTLASS has enabled an m-number for the wave of about 38 to be determined. This high-m wave was not detected by the IMAGE (International Monitor for Auroral Geomagnetic Effects) network of ground magnetometers, as expected for a wave of a small spatial scale size. These observations offer the first independent confirmation of the existence of the ground uncorrelated ULF wave signatures previously reported in measurements recorded from an HF Doppler sounder located in the vicinity of Tromsø. These results both demonstrate a new capability for geophysical exploration from the combined CUTLASS-EISCAT ionospheric Heater experiment, and provide a verification of the HF Doppler technique for the investigation of small scale ULF waves.     

Field Measurements and Numerical Simulations of the 2004 Tsunami Impact on the East Coast of Sri Lanka

This paper employs a numerical model of tsunami propagation together with documented observations and field measurements of the evidence left behind by the tsunami in December 2004, to identify and interpret the factors that have contributed to the significant spatial variability of the level of tsunami impact along the coastal belt of the eastern province of Sri Lanka. The model results considered in the present analysis include the distribution of the amplitude of the tsunami and the pattern of wave propagation over the continental shelf off the east coast, while the field data examined comprise the maximum water levels measured at or near the shoreline, the horizontal inundation distances and the number of housing and other buildings damaged. The computed maximum amplitude of the tsunami at water points nearest the shoreline along the east coast shows considerable variation ranging from 2.2 m to 11.4 m with a mean value of 5.7 m; moreover, the computed amplitudes agree well with the available field measurements. We also show that the shelf bathymetry off the east coast, particularly the submarine canyons at several locations, significantly influences the near-shore transformation of tsunami waves, and consequently, the spatial variation of the maximum water levels along the coastline. The measured values of inundation also show significant variation along the east coast and range from 70 m to 4560 m with a median value of 700 m. Our analyses of field data also show the dominant influence of the coastal topography and geomorphology on the extent of tsunami inundation. Furthermore, the measured inundation distances indicate no apparent correlation with the computed tsunami heights at the respective locations. We also show that both the computed tsunami heights and the measured inundation distances for the east coast closely follow the log-normal statistical distribution.     

Investigation of the relationship between optical auroral forms and HF radar E region backscatter

The SuperDARN HF radars have been employed in the past to investigate the spectral characteristics of coherent backscatter from L-shell aligned features in the auroral E region. The present study employs all-sky camera observations of the aurora from Husafell, Iceland, and the two SuperDARN radars located on Iceland, Þykkvibær and Stokkseyri, to determine the optical signature of such backscatter features. It is shown that, especially during quiet geomagnetic conditions, the backscatter region is closely associated with east-west aligned diffuse auroral features, and that the two move in tandem with each other. This association between optical and radar aurora has repercussions for the instability mechanisms responsible for generating the E region irregularities from which radars scatter. This is discussed and compared with previous studies investigating the relationship between optical and VHF radar aurora. In addition, although it is known that E region backscatter is commonly observed by SuperDARN radars, the present study demonstrates for the first time that multiple radars can observe the same feature to extend over at least 3 h of magnetic local time, allowing precipitation features to be mapped over large portions of the auroral zone.     

Surface currents and winds at the Delaware Bay mouth

Knowledge of the circulation of estuaries and adjacent shelf waters has relied on hydrographic measurements, moorings, and local wind observations usually removed from the region of interest. Although these observations are certainly sufficient to identify major characteristics, they lack both spatial resolution and temporal coverage. High-resolution synoptic observations are required to identify important coastal processes at smaller scales. Long observation periods are needed to properly sample low-frequency processes that may also be important. The introduction of high-frequency (HF) radar measurements and regional wind models for coastal studies is changing this situation. Here we analyze synoptic, high-resolution surface winds and currents in the Delaware Bay mouth over an 8-month period (October 2007 through May 2008). The surface currents were measured by two HF radars while the surface winds were extracted from a data-assimilating regional wind model. To illustrate the utility of these monitoring tools we focus on two 45-day periods which previously were shown to present contrasting pictures of the circulation. One, the low-outflow period is from 1 October through 14 November 2007; the other is the high-outflow period from 3 March through 16 April 2008. The large-scale characteristics noted by previous workers are clearly corroborated. Specifically the M2 tide dominates the surface currents, and the Delaware Bay outflow plume is clearly evident in the low frequency currents. Several new aspects of the surface circulation were also identified. These include a map of the spatial variability of the M2 tide (validating an earlier model study), persistent low-frequency cross-mouth flow, and a rapid response of the surface currents to a changing wind field. However, strong wind episodes did not persist long enough to set up a sustained Ekman response.     

The Liverpool Bay Coastal Observatory

A pilot Coastal Observatory has been established in Liverpool Bay which integrates (near) real-time measurements with coupled models and whose results are displayed on the web. The aim is to understand the functioning of coastal seas, their response to natural forcing and the consequences of human activity. The eastern Irish Sea is an apt test site, since it encompasses a comprehensive range of processes found in tidally dominated coastal seas, including near-shore physical and biogeochemical processes influenced by estuarine inflows, where both vertical and horizontal gradients are important. Applications include hypernutrification, since the region receives significantly elevated levels of nutrient inputs, shoreline management (coastal flooding and beach erosion/accretion), and understanding present conditions to predict the impact of climate change (for instance if the number and severity of storms, or of high or low river flows, change). The integrated measurement suite which started in August 2002 covers a range of space and time scales. It includes in situ time series, four to six weekly regional water column surveys, an instrumented ferry, a shore-based HF radar system measuring surface currents and waves, coastal tide gauges and visible and infra-red satellite data. The time series enable definition of the seasonal cycle, its inter-annual variability and provide a baseline from which the relative importance of events can be quantified. A suite of nested 3D hydrodynamic, wave and ecosystem models is run daily, focusing on the observatory area by covering the ocean/shelf of northwest Europe (at 12-km resolution) and the Irish Sea (at 1.8 km), and Liverpool Bay at the highest resolution of 200 m. The measurements test the models against events as they happen in a truly 3D context. All measurements and model outputs are displayed freely on the Coastal Observatory website () for an audience of researchers, education, coastal managers and the public.     

CUTLASS/IMAGE observations of high-latitude convection features during substorms

The CUTLASS Finland HF radar has been operational since February 1995. The radar frequently observes backscatter during the midnight sector from a latitude range 70–75° geographic, latitudes often associated with the polar cap. These intervals of backscatter occur during intervals of substorm activity, predominantly in periods of relatively quiet magnetospheric activity, with Kp during the interval under study being 2-and K_P for the day being only 8-. During August 1995 the radar ran in a high time resolution mode, allowing measurements of line-of-sight convection velocities along a single beam with a temporal resolution of 14s, and measurement of a full spatial scan of line-of-sight convection velocities every four minutes. Data from such scans reveal the radar to be measuring return flow convection during the interval of substorm activity. For three intervals during the period under study, a reduction in the spatial extent of radar backscatter occurred. This is a consequence of D region HF absorption and its limited extent in the present study is probably a consequence of the high latitude of the substorm activity, with the electrojet centre lying between 67° and 71° geomagnetic latitude. The high time resolution beam of the radar additionally demonstrates that the convection is highly time dependent. Pulses of equatorward flow exceeding 600 m s^–1 are observed with a duration of 5 min and a repetition period of 8 min. Their spatial extent in the CUTLASS field of view was 400–500 km in longitude, and 300–400 km in latitude. Each pulse of enhanced equatorward flow was preceded by an interval of suppressed flow and enhanced ionospheric Hall conductance. The transient features are interpreted as being due to ionospheric current vortices associated with field aligned current pairs. The relationship between these observations and substorm phenomena in the magnetotail is discussed.     

Numerical study of surface water circulation around Sekisei Lagoon, southwest Japan

This paper discusses the variability of surface currents around Sekisei Lagoon using a nested grid ocean circulation model. We developed a triple-nested grid system that consists of a coarse-resolution (1/60° or ∼1.85 km) model off Taiwan, an intermediate-resolution (1/300° or ∼370 m) model around the Yaeyama Islands, and a fine-resolution (1/900° or ∼123 m) model of Sekisei Lagoon. The nested grid system was forced by wind and heat flux calculated from six-hourly atmospheric reanalysis data and integrated over the period from May to July 2003. The coarse-resolution model was driven by lateral boundary conditions calculated from daily ocean reanalysis data to include realistic variation of the Kuroshio and mesoscale eddies with spatial scales of ∼500–700 km in the open ocean. The tidal forcing was included in the intermediate-resolution model by interpolating sea level data obtained from a data-assimilative tidal model. The results were then used to drive the fine-resolution model to simulate the surface water circulation around Sekisei lagoon. Model results show that (1) currents inside the lagoon are mainly driven by tide and wind; (2) there exists a strong southwestward current along the bottom slope in the southeast portion of the lagoon; the current is mainly driven by remote mesoscale eddies and at times intensified by the local wind; (3) the flow relaxation scheme is effective in reducing biases along the open boundaries. The simulated currents were used to examine the retention and dispersion of passive particles in the surface layer. Results show that the surface dispersion in the strong open ocean current region is significantly higher than that inside the lagoon.