Eddy generation and evolution in the North Pacific Subtropical Countercurrent (NPSC) zone

The seasonal generation and evolution of eddies in the region of the North Pacific Subtropical Countercurrent remain poorly understood due to the scarcity of available data. We used TOPEX/POSEIDON altimetry data from 1992 to 2007 to study the eddy field in this zone. We found that velocity shear between this region and the neighboring North Equatorial Current contributes greatly to the eddy generation. Furthermore, the eddy kinetic energy level (EKE) shows an annual cycle, maximum in April/May and minimum in December/January. Analyses of the temporal and spatial distributions of the eddy field revealed clearly that the velocity shear closely related to baroclinic instability processes. The eddy field seems to be more zonal than meridional, and the energy containing length scale shows a surprising lag of 2-3 months in comparison with the 1-D and 2-D EKE level. A similar phenomenon is observed in individual eddies in this zone. The results show that in this eddy field band, the velocity shear may drive the EKE level change so that the eddy field takes another 2-3 months to grow and interact to reach a relatively stable state. This explains the seasonal evolution of identifiable eddies.

卫星高度计资料在海洋潮汐研究中的应用

本文依据收集到的392个地面验潮站8个主要分潮（M2、S2、K1、O1、N2、K2、P1及Q1）的调和常数，对现有7个全球大洋潮汐模式的准确度进行了检验，结果显示各模式在深海区域均达到了比较高的准确度，相互之间差别也不大。经验模式GOT00和CSR4.0、同化模式NAO99、反演同化模式TPXO7.0、数值同化模式FES2002和FES2004的M2分潮均方根偏差在3 cm左右，其它分潮（S2、K1、O1、N2、K2、P1及Q1）大约在1～2 cm。本文还依据中国近海18个岛屿的调和常数对其中的5个大洋潮汐模式的准确度进行了检验，结果表明，M2分潮均方根偏差在6～14 cm，明显高于大洋部分的偏差，其中日本国家天文台的潮汐模式NAO99在中国近海的结果相对较准确。 我们利用1992年8月至2008年8月的TOPEX/POSEIDON和JASON-1(T/P-J)卫星高度计资料，对沿卫星轨道的302816个站点进行了14个分潮的潮汐调和分析，得到了全球大洋潮汐的8个主要分潮以及2个气象分潮Sa、Ssa的经验同潮图。主要结果有：（1）各分潮在卫星上升轨道与下降轨道的交叉点(约7000个)相关性分析表明：M2分潮的振幅和迟角的相关系数很高（分别为0.9965和0.9961）；S2，K1，O1和Sa分潮也有较好的相关性（相关系数为0.94～0.99）；（2）该结果与392地面个验潮站吻合较好，其中M2分潮的振幅、迟角和向量的均方根偏差分别为：1.73 cm，2.340和2.93 cm；S2，K1和O1分潮的振幅、迟角和向量的均方根偏差为1 cm左右，5.250～7.270和1.5～2.1 cm，该精度与最近几年国际上的主要大洋潮汐模式的准确度相近；（3）首次通过卫星资料获得了Sa、Ssa分潮的同潮图。周期为1年的Sa分潮与大洋105个地面站相比，振幅、迟角和向量的均方根偏差分别为1.50 cm、18.360和2.16 cm。在此基础上，进一步分析了构成Sa、Ssa气象分潮的两个主要因素（海水密度以及海面气压）在全球的分布。 在T/P-J等卫星资料无法覆盖到南大洋和北冰洋，本文利用Princeton Ocean Model（POM）进行了数值模拟，模拟结果与162个地面实测站（其中南大洋30个，北冰洋132个）的观测比较一致。基于卫星资料分析的结果和数值模拟结果合并得到了全球大洋的8个主要分潮同潮图。在此基础上通过全球潮汐能量耗散的计算得到潮能通量的分布，并得到全球M2、S2、K1和O1分潮的潮汐能量耗散率为2.431TW、0.401TW、0.336TW和0.176TW。 本文还利用卫星资料对南海潮汐进行了研究，在中国南海，获得了主要的半日潮、全日潮、四分日分潮和长周期分潮（M2，S2，N2，K2，K1，O1，P1，Q1，M4, MS4，Sa, Ssa）的经验同潮图。与南海沿岸94个地面验潮站的数据符合得比较好，M2，S2，K1及O1等4个主要分潮的平均振幅差为2～4 cm，均方根偏差分别是9～11 cm.其它4个主要分潮N2，K2，P1，Q1的平均振幅差为1～2 cm，均方根偏差为2～4 cm。此外，本文还利用卫星高度计资料潮汐分析结果沿卫星轨道进行高通滤波，分离得出中国近海的M2，S2，K1及O1分潮的内潮信息。

Position and structure of the subtropical/Azores front region from combined Lagrangian and remote sensing (IR/altimeter/SeaWIFS) measurements

The position and structure of the North Atlantic Subtropical Front is studied using Lagrangian flow tracks and remote sensing (AVHRR imagery: TOPEX/POSEIDON altimetry: SeaWiFS) in a broad region ( similar to 31 degree to similar to 36 degree N) of marked gradient of dynamic height (Azores Current) that extends from the Mid-Atlantic Ridge (MAR), near similar to 40 degree W, to the Eastern Boundary ( similar to 10 degree W). Drogued Argos buoy and ALACE tracks are superposed on infrared satellite images in the Subtropical Front region. Cold (cyclonic) structures, called storms, and warm (anticyclonic) structures of 100-300 km in size can be found on the south side of the Subtropical Front outcrop, which has a temperature contrast of about 1 degree C that can be followed for similar to 2500 km near 35 degree N. Warmer water adjacent to the outcrop is flowing eastward (Azores Current) but some warm water is returned westward about 300 km to the south (southern Counterflow). Estimates of horizontal diffusion in a Storm (D=2.2t10 super(2) m super(2) s super(-1)) and in the Subtropical Front region near 200 m depth (D sub(x)=1.3t10 super(4) m super(2) s super(-1), D sub(y)=2.6t10 super(3) m super(2) s super(-1)) are made from the Lagrangian tracks. Altimeter and in situ measurements show that Storms track westwards. Storms are separated by about 510 km and move westward at 2.7 km d super(-1). Remote sensing reveals that some initial structures start evolving as far east as 23 degree W but are more organized near 29 degree W and therefore Storms are about 1 year old when they reach the MAR (having travelled a distance of 1000 km). Structure and seasonality in SeaWiFS data in the region is examined.

Thermal front variability along the North Atlantic Current observed using satellite microwave and infrared data

Thermal fronts detected using multiple satellite sensors have been integrated to provide new information on the spatial and seasonal distribution of oceanic fronts in the North Atlantic. The branching of the North Atlantic Current (NAC) as it encounters the Mid-Atlantic Ridge (MAR) is reflected in surface thermal fronts, which preferentially occur at the Charlie Gibbs Fracture Zone (CGFZ) and several smaller fracture zones. North of the CGFZ there are few thermal fronts, contrasting with the region to the south, where there are frequent surface thermal fronts that are persistent seasonally and interannually. The alignment of the fronts confirms that the shallower Reykjanes Ridge north of the CGFZ is more of a barrier to water movements than the ridge to the south. Comparison of front distributions with satellite altimetry data indicates that the MAR influence on deep ocean currents is also frequently exhibited in surface temperature. The improved spatial and temporal resolution of the front analysis has revealed consistent seasonality in the branching patterns. These results contribute to our understanding of the variability of the NAC, and the techniques for visualising oceanic fronts can be applied in other regions to reveal details of surface currents that cannot be resolved using satellite altimetry or in situ measurements.

Seasonality and interannual variability of the European Slope Current from 20 years of altimeter data compared with in situ measurements

The European Slope Current (SC) is a major section of the warm poleward flow from the Atlantic to the Arctic, which also moderates the exchange of heat, salt, nutrients and carbon between the deep ocean and the European shelf seas. The mean structure of the geostrophic flow, seasonality, interannual variability and long-term trend of SC are appraised with an unprecedented continuous 20-year satellite altimeter dataset. Comparisons with long term in situ data showed a maximum correlation of r2=0.51 between altimeter and Acoustic Doppler Current Profilers (ADCP), with similar results for drogued buoy data. Mean geostrophic currents were appraised more comprehensively than previous attempts, and the paths of 4 branches of the North Atlantic Current (NAC) and positions of 5 eddies in the region were derived quantitatively. A consistent seasonal cycle in the flow of the SC was found at all 8 sections along the European shelf slope, with maximum poleward flow in the winter and minimum in the summer. The seasonal difference in the altimetry current speed amounted to ~8-10 cm s-1 at the northern sections, but only ~5 cm s-1 on the Bay of Biscay slopes. This extended altimeter dataset indicates significant regional and seasonal variations, and has revealed new insights into the interannual variability of the SC. It is shown that there is a peak poleward flow at most positions along a ~2000 km stretch of the continental slope from Portugal to Scotland during 1995-1997, but this did not clearly relate to the extreme negative North Atlantic Oscillation (NAO) in the winter of 1995-1996. The speed of the SC exhibited a long term decreasing trend of ~1% per year. By contrast the NAC showed no significant trend over the 20-year period. Major changes in the NAC occurred three times, and these changes followed decreases in the NAO index.

Southern Annular Mode and westerly-wind-driven changes in Indian-Atlantic exchange mechanisms

The dynamical link between the Indian Ocean and Atlantic Meridional Overturning Circulation (AMOC) remains poorly understood. This partly arises from the complex Agulhas leakage, which occurs via rings, cyclones, and non-eddy flux. Hindcast simulations suggest that leakage has recently increased but have not decomposed this signal into its constituent mechanisms. Here these are isolated in a realistic ocean model. Increases in simulated leakage are attributed to stronger eddy and non-eddy-driven transports, and a strong warming and salinification, especially within Agulhas rings. Variability in both regimes is associated with strengthening Indian Ocean westerly winds, reflecting an increasingly positive Southern Annular Mode. While eddy and non-eddy flux signals are tied through turbulent eddy dissipation, the ratio between the two varies decadally. Consequently, while altimetry suggests a recent increase in retroflection turbulence and implied leakage, non-eddy flux may also play a significant role in modulating the leakage AMOC connection.

Exploring multi-sensor satellite synergies to provide direction to high-resolution along-track altimtery currents

A new approach to compute along-track velocity components by combining altimetry-based across-track components and front directions from remote sensing maps of surface chlorophyll concentration is proposed. The analysis focuses on the South Madagascar region characterized by the strong East Madagascar Current and sharp gradients of surface tracers. The results are compared against in-situ observations from three moorings along the Jason-1 track 196. Accurate information on the total velocity direction is the key factor for obtaining accurate estimates of along-track velocities. Although with some limitations, surface tracer fronts can be successfully used to retrieve such information.

Rapid rise in Effective Sea-Level in southwest Bangladesh: Its causes and contemporary rates

Evidence is presented from three estuarine tide gauges located in the
Sundarban area of southwest Bangladesh of relative sea level rise
substantially in excess of the generally accepted rates from altimetry, as
well as previous tide-gauge analyses. It is proposed that the difference
arises from the use of relative mean sea level (RMSL) to characterise the
present and future coastal flood hazard, since RMSL can be misleading in
estuaries in which tidal range is changing. Three tide gauges, one located in
the uninhabited mangrove forested area (Sundarban) of southwest
Bangladesh, the others in the densely populated polder zone north of the
present Sundarban, show rates of increase in RMSL ranging from 2.8 mm
a-1 to 8.8 mm a-1. However, these trends in RMSL disguise the fact that high
water levels in the polder zone have been increasing at an average rate of
15.9 mm a-1 and a maximum of 17.2 mm a-1. In an area experiencing tidal
range amplification, RMSL will always underestimate the rise in high water
levels; consequently, as an alternative to RMSL, the use of trends in high
water maxima or ‘Effective Sea Level Rise’ (ESLR) is adopted as a more
strategic parameter to characterise the flooding hazard potential. The rate
of increase in ESLR is shown to be due to a combination of deltaic
subsidence, including sediment compaction, and eustatic sea level rise, but
principally as a result of increased tidal range in estuary channels recently
constricted by embankments. These increases in ESLR have been partially
offset by decreases in fresh water discharge in those estuaries connected
to the Ganges. The recognition of increases of the effective sea level in the
Bangladesh Sundarban, which are substantially greater than increases in
mean sea level, is of the utmost importance to flood management in this
low-lying and densely populated area.

Meridional mixing of water masses in the Kuroshio Extension: Origins of North Pacific Intermediate Water anomalies

Senior thesis written for Oceanography 445

Caracterização de ambientes de transição do SW Alentejano a distintas escalas temporais

Tese de doutoramento, Geologia (Geologia Económica e do Ambiente), Universidade de Lisboa, Faculdade de Ciências, 2016

Modernização da linha da Beira Baixa: troço Covilhã - Guarda

Trabalho de Projeto para obtenção do grau de Mestre em Engenharia Civil na Área de Especialização em Vias de Comunicação e Transportes

Seasonal and interannual variability of sea level and associated surface meteorological parameters at Cochin

This thesis entitled seasonal and interannual variability of sea level and associated surface meteorological parameters at cochin.The interesting aspect of studying sea level variability on different time scales can be attributed to the diversity of its applications.Study of tides could perhaps be the oldest branch of physical oceanography.The thesis is presented in seven chapters. The first chapter gives, apart from a general introduction, a survey of literature on sea level variability on different time scales - tidal, seasonal and interannual (geological scales excluded), with particular emphasis on the work carried out in the Indian waters. The second chapter is devoted to the study of observed tides at Cochin on seasonal and interannual time scales using hourly water level data for the period 1988-1993. The third chapter describes the long-term climatology of some important surface oceanographic and meteorological parameters (at Cochin) which are supposed to affect the sea level. The fourth chapter addresses the problem of seasonal forecasting of the meteorological and oceanographic parameters at Cochin using autoregressive, sinusoidal and exponentially weighted moving average techniques and testing their accuracy with the observed data for the period 1991-1993. The fifth chapter describes the seasonal cycles of sea level and the driving forces at 16 stations along the Indian subcontinent. It also addresses the observed interannual variability of sea level at 15 stations using available multi-annual data sets. The sixth chapter deals with the problem of coastal trapped waves between Cochin and Beypore off the Kerala coast using sea level and atmospheric pressure data sets for the year 1977. The seventh and the last chapter contains the summary and conclusions and future outlook based on this study.

Benchmarking 2D hydraulic models for urban flooding

This paper describes benchmark testing of six two-dimensional (2D) hydraulic models (DIVAST, DIVASTTVD, TUFLOW, JFLOW, TRENT and LISFLOOD-FP) in terms of their ability to simulate surface flows in a densely urbanised area. The models are applied to a 1·0 km × 0·4 km urban catchment within the city of Glasgow, Scotland, UK, and are used to simulate a flood event that occurred at this site on 30 July 2002. An identical numerical grid describing the underlying topography is constructed for each model, using a combination of airborne laser altimetry (LiDAR) fused with digital map data, and used to run a benchmark simulation. Two numerical experiments were then conducted to test the response of each model to topographic error and uncertainty over friction parameterisation. While all the models tested produce plausible results, subtle differences between particular groups of codes give considerable insight into both the practice and science of urban hydraulic modelling. In particular, the results show that the terrain data available from modern LiDAR systems are sufficiently accurate and resolved for simulating urban flows, but such data need to be fused with digital map data of building topology and land use to gain maximum benefit from the information contained therein. When such terrain data are available, uncertainty in friction parameters becomes a more dominant factor than topographic error for typical problems. The simulations also show that flows in urban environments are characterised by numerous transitions to supercritical flow and numerical shocks. However, the effects of these are localised and they do not appear to affect overall wave propagation. In contrast, inertia terms are shown to be important in this particular case, but the specific characteristics of the test site may mean that this does not hold more generally.

Calibration of uncertain flood inundation models using remotely sensed water levels

A traditional method of validating the performance of a flood model when remotely sensed data of the flood extent are available is to compare the predicted flood extent to that observed. The performance measure employed often uses areal pattern-matching to assess the degree to which the two extents overlap. Recently, remote sensing of flood extents using synthetic aperture radar (SAR) and airborne scanning laser altimetry (LIDAR) has made more straightforward the synoptic measurement of water surface elevations along flood waterlines, and this has emphasised the possibility of using alternative performance measures based on height. This paper considers the advantages that can accrue from using a performance measure based on waterline elevations rather than one based on areal patterns of wet and dry pixels. The two measures were compared for their ability to estimate flood inundation uncertainty maps from a set of model runs carried out to span the acceptable model parameter range in a GLUE-based analysis. A 1 in 5-year flood on the Thames in 1992 was used as a test event. As is typical for UK floods, only a single SAR image of observed flood extent was available for model calibration and validation. A simple implementation of a two-dimensional flood model (LISFLOOD-FP) was used to generate model flood extents for comparison with that observed. The performance measure based on height differences of corresponding points along the observed and modelled waterlines was found to be significantly more sensitive to the channel friction parameter than the measure based on areal patterns of flood extent. The former was able to restrict the parameter range of acceptable model runs and hence reduce the number of runs necessary to generate an inundation uncertainty map. A result of this was that there was less uncertainty in the final flood risk map. The uncertainty analysis included the effects of uncertainties in the observed flood extent as well as in model parameters. The height-based measure was found to be more sensitive when increased heighting accuracy was achieved by requiring that observed waterline heights varied slowly along the reach. The technique allows for the decomposition of the reach into sections, with different effective channel friction parameters used in different sections, which in this case resulted in lower r.m.s. height differences between observed and modelled waterlines than those achieved by runs using a single friction parameter for the whole reach. However, a validation of the modelled inundation uncertainty using the calibration event showed a significant difference between the uncertainty map and the observed flood extent. While this was true for both measures, the difference was especially significant for the height-based one. This is likely to be due to the conceptually simple flood inundation model and the coarse application resolution employed in this case. The increased sensitivity of the height-based measure may lead to an increased onus being placed on the model developer in the production of a valid model

Use of airborne remote sensing to detect riverside Brassica rapa to aid in risk assessment of transgenic crops

High resolution descriptions of plant distribution have utility for many ecological applications but are especially useful for predictive modelling of gene flow from transgenic crops. Difficulty lies in the extrapolation errors that occur when limited ground survey data are scaled up to the landscape or national level. This problem is epitomized by the wide confidence limits generated in a previous attempt to describe the national abundance of riverside Brassica rapa (a wild relative of cultivated rapeseed) across the United Kingdom. Here, we assess the value of airborne remote sensing to locate B. rapa over large areas and so reduce the need for extrapolation. We describe results from flights over the river Nene in England acquired using Airborne Thematic Mapper (ATM) and Compact Airborne Spectrographic Imager (CASI) imagery, together with ground truth data. It proved possible to detect 97% of flowering B. rapa on the basis of spectral profiles. This included all stands of plants that occupied >2m square (>5 plants), which were detected using single-pixel classification. It also included very small populations (<5 flowering plants, 1-2m square) that generated mixed pixels, which were detected using spectral unmixing. The high detection accuracy for flowering B. rapa was coupled with a rather large false positive rate (43%). The latter could be reduced by using the image detections to target fieldwork to confirm species identity, or by acquiring additional remote sensing data such as laser altimetry or multitemporal imagery.