Storm Surges in the Strait of Georgia Simulated with a Regional Model

The Strait of Georgia is a large, semi-enclosed body of water between Vancouver Island and the mainland of British Columbia connected to the Pacific Ocean via Juan de Fuca Strait at the south and Johnstone Strait at the north. During the winter months, coastal communities along the Strait of Georgia are at risk of flooding caused by storm surges, a natural hazard that can occur when a strong storm coincides with high tide. This investigation produces storm surge hindcasts using a three-dimensional numerical ocean model for the Strait of Georgia and the surrounding bodies of water (Juan de Fuca Strait, Puget Sound, and Johnstone Strait) collectively known as the Salish Sea. The numerical model employs the Nucleus for European Modelling of the Ocean architecture in a regional configuration. The model is evaluated through comparisons of tidal elevation harmonics and storm surge with observations. Important forcing factors contributing to storm surges are assessed. It is shown that surges entering the domain from the Pacific Ocean make the most significant contribution to surge amplitude within the Strait of Georgia. Comparisons between simulations and high-resolution and low-resolution atmospheric forcing further emphasize that remote forcing is the dominant factor in surge amplitudes in this region. In addition, local wind patterns caused a slight increase in surge amplitude on the mainland side of the Strait of Georgia compared with Vancouver Island coastal areas during a major wind storm on 15 December 2006. Generally, surge amplitudes are found to be greater within the Strait of Georgia than in Juan de Fuca Strait.     

Toward an early warning system for dengue prevention: modeling climate impact on dengue transmission

Dengue fever is the most prevalent mosquito-borne viral disease of humans in tropical lands. As an efficient vaccine is not yet available, the only means to prevent epidemics is to control mosquito populations. These are influenced by human behavior and climatic conditions and thus, need constant effort and are very expansive. Examples of succeeded prevention are rare because of the continuous reintroduction of virus or vector from outside, or growing resistance of mosquito populations to insecticides. Climate variability and global warming are other factors which may favour epidemics of dengue. During a pilot study in Claris EC project, a model for the transmission of dengue was built, to serve as a tool for estimating the risk of epidemic transmission and eventually forecasting the risk under climatic change scenarios. An ultimate objective would be to use the model as an early warning system with meteorological forecasts as input, thus allowing better decision making and prevention.     

Calibration framework for a Kalman filter applied to a groundwater model

The paper presents a novel approach to the setup of a Kalman filter by using an automatic calibration framework for estimation of the covariance matrices. The calibration consists of two sequential steps: (1) Automatic calibration of a set of covariance parameters to optimize the performance of the system and (2) adjustment of the model and observation variance to provide an uncertainty analysis relying on the data instead of ad-hoc covariance values. The method is applied to a twin-test experiment with a groundwater model and a colored noise Kalman filter. The filter is implemented in an ensemble framework. It is demonstrated that lattice sampling is preferable to the usual Monte Carlo simulation because its ability to preserve the theoretical mean reduces the size of the ensemble needed. The resulting Kalman filter proves to be efficient in correcting dynamic error and bias over the whole domain studied. The uncertainty analysis provides a reliable estimate of the error in the neighborhood of assimilation points but the simplicity of the covariance models leads to underestimation of the errors far from assimilation points.     

Plate Motion of India and Interseismic Strain in the Nepal Himalaya from GPS and DORIS Measurements

We analyse geodetically estimated deformation across the Nepal Himalaya in order to determine the geodetic rate of shortening between Southern Tibet and India, previously proposed to range from 12 to 21 mm yr^?1. The dataset includes spirit-levelling data along a road going from the Indian to the Tibetan border across Central Nepal, data from the DORIS station on Everest, which has been analysed since 1993, GPS campaign measurements from surveys carried on between 1995 and 2001, as well as data from continuous GPS stations along a transect at the logitude of Kathmandu operated continuously since 1997. The GPS data were processed in International Terrestrial Reference Frame 2000 (ITRF2000), together with the data from 20 International GNSS Service (IGS) stations and then combined using quasi- observation combination analysis (QOCA). Finally, spatially complementary velocities at stations in Southern Tibet, initially determined in ITRF97, were expressed in ITRF2000. After analysing previous studies by different authors, we determined the pole of rotation of the Indian tectonic plate to be located in ITRF2000 at 51.409±1.560° N and ?10.915±5.556°E, with an angular velocity of 0.483±0.015°. Myr^?1. Internal deformation of India is found to be small, corresponding to less than about 2 mm yr^?1 of baseline change between Southern India and the Himalayan piedmont. Based on an elastic dislocation model of interseismic strain and taking into account the uncertainty on India plate motion, the mean convergence rate across Central and Eastern Nepal is estimated to 19±2.5 mm yr^?1, (at the 67% confidence level). The main himalayan thrust (MHT) fault was found to be locked from the surface to a depth of about 20 km over a width of about 115 km. In these regions, the model parameters are well constrained, thanks to the long and continuous time-series from the permanent GPS as well as DORIS data. Further west, a convergence rate of 13.4±5 mm yr^?1, as well as a fault zone, locked over 150 km, are proposed. The slight discrepancy between the geologically estimated deformation rate of 21±1.5 mm yr^?1 and the 19±2.5 mm yr^?1 geodetic rate in Central and Eastern Nepal, as well as the lower geodetic rate in Western Nepal compared to Eastern Nepal, places bounds on possible temporal variations of the pattern and rate of strain in the period between large earthquakes in this region.     

Epiclastic deposits and ‘horseshoe-shaped’ calderas in Tahiti (Society Islands) and Ua Huka (Marquesas Archipelago), French Polynesia

Occurrences of debris avalanche deposits newly identified in Tahiti (Society Islands) and Ua Huka (Marquesas Archipelago) are described and interpreted here. In both islands, the breccias are located within horseshoe-shaped residual calderas. In Tahiti, the epiclastic formations, up to 500 m thick, lie on the floor of the central depression and in the valley of the northwards running Papenoo River. In Ua Huka, the breccias crop out within a depression limited by a semicircular crest in four bays along the southern coast. Their thickness is ca. 100 m. A few clasts collected in the Tahitian breccias and some rocks forming their substratum have been dated (K–Ar datings) and analysed (major and trace elements, Sr–Nd isotopes) for this study. Using these data, we show that the debris avalanche(s) occurred in Tahiti Nui at the end of the growth of the shield volcano (between 570 000 and 390 000 years ago), maybe in consequence of the emplacement of the plutonic body which occupies the central part of the caldera. In Ua Huka, the collapse took place nearly 3 Ma ago, between the construction of the shield volcano and that of the inner one. The southwards orientation of the caldera, like that of the neighbouring island Nuku Hiva, might reflect a preferential direction of weakness in the substratum of the central Marquesas.     

High-resolution modelling of a coastal harbour in the presence of strong tides and significant river runoff

Ocean Dynamics - In the context of Canada’s Ocean Protection Plan (OPP), improved coastal and near-shore modelling is needed to enhance marine safety and emergency response capacity in the...     

The planet formation imager

The Planet Formation Imager (PFI, www.planetformationimager.org) is a next-generation infrared interferometer array with the primary goal of imaging the active phases of planet formation in nearby star forming regions. PFI will be sensitive to warm dust emission using mid-infrared capabilities made possible by precise fringe tracking in the near-infrared. An L/M band combiner will be especially sensitive to thermal emission from young exoplanets (and their disks) with a high spectral resolution mode to probe the kinematics of CO and H₂O gas. In this paper, we give an overview of the main science goals of PFI, define a baseline PFI architecture that can achieve those goals, point at remaining technical challenges, and suggest activities today that will help make the Planet Formation Imager facility a reality.     

A possible link between gabbros bearing High Temperature Iddingsite alteration and huge pegmatoid intrusions: The Society Islands, French Polynesia

This paper reports new field observations and new petrological, textural and geochemical data on two gabbroic intrusions and one pegmatoid dyke from the French Polynesian islands Maupiti and Bora Bora, respectively. Olivine crystals from the Faataufi and Barque de Hiro gabbros (Maupiti) include High Temperature Iddingsite (HTI), exhibiting three distinct morphological facies. Chemical and crystallographical data performed by electron microprobe, scanning microscope, X-ray diffractometry and Near InfraRed spectroscopy show that HTI results from modifications of olivine structure through intense hydroxylation and Fe-oxidation. The HTI-free 120 m-wide Cloche de Hiro dyke (Bora Bora) exhibits a textural zonation from heterogranular/intergranular gabbro (peripheral part) to typical pegmatoid (inner area). We propose a textural tetrahedral classificatory diagram which can be used to clearly distinguish the textural characteristics of the pegmatoids from those of the other coarse-grained samples. Through comparison with a lava flow containing vesicle-rich segregation sheets, we suggest that the Cloche de Hiro dyke may correspond to a huge segregation structure derived from a Maupiti-type HTI-bearing gabbro through a vapor-differentiation process. Both intrusions might be considered as representative of two superposed levels in a vertically arranged magmatic complex. In this view, the Maupiti gabbroic bodies would represent the lower part of the postulated system. Under high oxygen fugacity, olivine crystals would have undergone important modifications of their structure, leading to HTI. The huge Bora Bora pegmatoid dyke would be the upper part of the complex. After expulsion from the HTI-bearing gabbroic solidification zone by build-up of gas pressure, vesicle-rich residual melt would have gathered into a shallow level and crystallized there as a pegmatoid body.     

Neural network based daily precipitation generator (NNGEN-P)

Daily weather generators are used in many applications and risk analyses. The present paper explores the potential of neural network architectures to design daily weather generator models. Focusing this first paper on precipitation, we design a collection of neural networks (multi-layer perceptrons in the present case), which are trained so as to approximate the empirical cumulative distribution (CDF) function for the occurrence of wet and dry spells and for the precipitation amounts. This approach contributes to correct some of the biases of the usual two-step weather generator models. As compared to a rainfall occurrence Markov model, NNGEN-P represents fairly well the mean and standard deviation of the number of wet days per month, and it significantly improves the simulation of the longest dry and wet periods. Then, we compared NNGEN-P to three parametric distribution functions usually applied to fit rainfall cumulative distribution functions (Gamma, Weibull and double-exponential). A data set of 19 Argentine stations was used. Also, data corresponding to stations in the United States, in Europe and in the Tropics were included to confirm the results. One of the advantages of NNGEN-P is that it is non-parametric. Unlike other parametric function, which adapt to certain types of climate regimes, NNGEN-P is fully adaptive to the observed cumulative distribution functions, which, on some occasions, may present complex shapes. On-going works will soon produce an extended version of NNGEN to temperature and radiation.     

Segregation structures in vapor-differentiated basaltic flows

 Vesicle cylinders represent a spectacular kind of segregation structure involving residual liquids formed in situ during the cooling of lava flows. These vertical pipes, commonly found within basalt flows typically 2–10 m thick, are interpreted as the product of a vapor-driven differentiation process. The olivine phenocrysts and the earliest generation of groundmass olivines found in cylinder-bearing basalts appear to have been generally affected by magmatic oxidation, resulting in high-temperature iddingsite (HTI) alteration. This feature is also observed within cylinder-free basalt flows which exhibit other kinds of vesicular segregation structures, such as vesicle-rich pegmatoid segregation sheets and/or segregation vesicles. Detailed textural, petrological, and geochemical characteristics of two types of cylinders, three types of vesicle sheets, and five types of segregation vesicles are described, based on the study of 12 occurrences of HTI-bearing basalt flows from oceanic shield volcanoes or continental basalt plateaus. We propose a general classification of these segregation structures likely to derive from vapor differentiation. Flow thickness is probably the main factor influencing their morphology. Finally, we suggest that the concomitant occurrence of olivine oxidation and vapor-differentiation effects results from the late persistence of water oversaturation after eruption, perhaps due to a high rate of magma ascent. Received: 27 March 1999 / Accepted: 15 February 2000