(Table S2) The reconstructed sea surface temperature (SST) over the eastern equatorial Pacific

Abundant hydroclimatic evidence from western Amazonia and the adjacent Andes documents wet conditions during Heinrich Stadial 1 (HS1, 18-15 ka), a cold period in the high latitudes of the North Atlantic. This precipitation anomaly was attributed to a strengthening of the South American summer monsoon due to a change in the Atlantic interhemispheric sea surface temperature (SST) gradient. However, the physical viability of this mechanism has never been rigorously tested. We address this issue by combining a thorough compilation of tropical South American paleorecords and a set of atmosphere model sensitivity experiments. Our results show that the Atlantic SST variations alone, although leading to dry conditions in northern South America and wet conditions in northeastern Brazil, cannot produce increased precipitation over western Amazonia and the adjacent Andes during HS1. Instead, an eastern equatorial Pacific SST increase (i.e., 0.5-1.5 °C), in response to the slowdown of the Atlantic Meridional Overturning Circulation during HS1, is crucial to generate the wet conditions in these regions. The mechanism works via anomalous low sea level pressure over the eastern equatorial Pacific, which promotes a regional easterly low-level wind anomaly and moisture recycling from central Amazonia towards the Andes.

ROCK OUTCROPS IN THE CANADIAN SHIELD: AN INVESTIGATION OF CONTAMINANT TRANSPORT FROM SURFACE SOURCES IN FRACTURED ROCK AQUIFERS

Flow, recharge and transport dynamics in fractured rock aquifers with low lying rock outcrops is a largely unexplored area of study in hydrogeology. The purpose of this thesis is to examine these topics in an agricultural area in Eastern Ontario. The study consists of a regional scale groundwater quality study, an infiltration experiment that considers bacteria transport from the ground surface to a well, and a numerical modelling study that tests the parameters that affect surface infiltration of a tracer from a rock outcrop to a deeper horizontal fracture. In the water quality study, approximately 65% of the samples contained total coliform, 16% contained E. coli, and 1% contained nitrate-N at greater than 5 mg/L. Occurrence of E. coli increased when considering seasonality, where wells were drilled on rock outcrops, and for shallow well intervals. Nitrate-N did not occur above the Guidelines for Canadian Drinking Water Quality (Health Canada, 2012) of 10 mg/L. Rapid arrival times were observed in the infiltration study for both the microspheres (30 minutes) and a dye tracer (45 minutes) in a well approximately 6.0 m in horizontal and 2.8 m in vertical distance from the tracer source. Transport velocities were approximately 38.9 m/day for the dye tracer and 115.2 m/day for the colloidal tracer. Results of the model runs indicate that overburden can provide an effective protective layer to transport in fractures, that high groundwater velocities occur in larger fracture apertures and higher gradients dilute tracer concentrations, and that lower groundwater velocities occur with smaller fracture apertures and lower gradients result in elevated tracer concentrations. Lower rainfall rates, larger fracture apertures, early tracer time, larger gradients, and lower water levels maintained unsaturated conditions for longer time periods such that tracer transport was delayed until saturated conditions were attained. The overall heterogeneity of this aquifer environment creates a source water protection conundrum where the water quality is generally good, while transport can occur very quickly in proximity to rock outcrops and in areas with limited overburden.

Modeling the impact of large-scale solid mass transfers on the gravity field

The focus of this thesis is to explore and quantify the response of large-scale solid mass transfer events on satellite-based gravity observations. The gravity signature of large-scale solid mass transfers has not been deeply explored yet; mainly due to the lack of significant events during dedicated satellite gravity missions‘ lifespans. In light of the next generation of gravity missions, the feasibility of employing satellite gravity observations to detect submarine and surface mass transfers is of importance for geoscience (improves the understanding of geodynamic processes) and for geodesy (improves the understanding of the dynamic gravity field). The aim of this thesis is twofold and focuses on assessing the feasibility of using satellite gravity observations for detecting large-scale solid mass transfers and on modeling the impact on the gravity field caused by these events. A methodology that employs 3D forward modeling simulations and 2D wavelet multiresolution analysis is suggested to estimate the impact of solid mass transfers on satellite gravity observations. The gravity signature of various submarine and subaerial events that occurred in the past was estimated. Case studies were conducted to assess the sensitivity and resolvability required in order to observe gravity differences caused by solid mass transfers. Simulation studies were also employed in order to assess the expected contribution of the Next Generation of Gravity Missions for this application.

Modeling long-term variability and change of soil moisture and groundwater level - from catchment to global scale

The water stored in and flowing through the subsurface is fundamental for sustaining human activities and needs, feeding water and its constituents to surface water bodies and supporting the functioning of their ecosystems. Quantifying the changes that affect the subsurface water is crucial for our understanding of its dynamics and changes driven by climate change and other changes in the landscape, such as in land-use and water-use. It is inherently difficult to directly measure soil moisture and groundwater levels over large spatial scales and long times. Models are therefore needed to capture the soil moisture and groundwater level dynamics over such large spatiotemporal scales. This thesis develops a modeling framework that allows for long-term catchment-scale screening of soil moisture and groundwater level changes. The novelty in this development resides in an explicit link drawn between catchment-scale hydroclimatic and soil hydraulics conditions, using observed runoff data as an approximation of soil water flux and accounting for the effects of snow storage-melting dynamics on that flux. Both past and future relative changes can be assessed by use of this modeling framework, with future change projections based on common climate model outputs. By direct model-observation comparison, the thesis shows that the developed modeling framework can reproduce the temporal variability of large-scale changes in soil water storage, as obtained from the GRACE satellite product, for most of 25 large study catchments around the world. Also compared with locally measured soil water content and groundwater level in 10 U.S. catchments, the modeling approach can reasonably well reproduce relative seasonal fluctuations around long-term average values. The developed modeling framework is further used to project soil moisture changes due to expected future climate change for 81 catchments around the world. The future soil moisture changes depend on the considered radiative forcing scenario (RCP) but are overall large for the occurrence frequency of dry and wet events and the inter-annual variability of seasonal soil moisture. These changes tend to be higher for the dry events and the dry season, respectively, than for the corresponding wet quantities, indicating increased drought risk for some parts of the world.

Modelling the Sun's small scale global photospheric magnetic field

We present a new model for the Sun's global photospheric magnetic field during a deep minimum of activity, in which no active regions emerge. The emergence and subsequent evolution of small- scale magnetic features across the full solar surface is simulated, subject to the influence of a global supergranular flow pattern. Visually, the resulting simulated magnetograms reproduce the typical structure and scale observed in quiet Sun magnetograms. Quantitatively, the simulation quickly reaches a steady state, resulting in a mean field and flux distribution that are in good agreement with those determined from observations. A potential coronal magnetic field is extrapolated from the simulated full Sun magnetograms, to consider the implications of such a quiet photospheric magnetic field on the corona and inner heliosphere. The bulk of the coronal magnetic field closes very low down, in short connections between small-scale features in the simulated magnetic network. Just 0.1% of the photospheric magnetic flux is found to be open at 2:5 Rʘ, around 10 - 100 times less than that determined for typical HMI synoptic map observations. If such conditions were to exist on the Sun, this would lead to a significantly weaker interplanetary magnetic field than is presently observed, and hence a much higher cosmic ray flux at Earth.

Detection of mesoscale thermal fronts from 4km data using smoothing techniques: Gradient-based fronts classification and basin scale application

In order to optimize frontal detection in sea surface temperature fields at 4 km resolution, a combined statistical and expert-based approach is applied to test different spatial smoothing of the data prior to the detection process. Fronts are usually detected at 1 km resolution using the histogram-based, single image edge detection (SIED) algorithm developed by Cayula and Cornillon in 1992, with a standard preliminary smoothing using a median filter and a 3 × 3 pixel kernel. Here, detections are performed in three study regions (off Morocco, the Mozambique Channel, and north-western Australia) and across the Indian Ocean basin using the combination of multiple windows (CMW) method developed by Nieto, Demarcq and McClatchie in 2012 which improves on the original Cayula and Cornillon algorithm. Detections at 4 km and 1 km of resolution are compared. Fronts are divided in two intensity classes (“weak” and “strong”) according to their thermal gradient. A preliminary smoothing is applied prior to the detection using different convolutions: three type of filters (median, average and Gaussian) combined with four kernel sizes (3 × 3, 5 × 5, 7 × 7, and 9 × 9 pixels) and three detection window sizes (16 × 16, 24 × 24 and 32 × 32 pixels) to test the effect of these smoothing combinations on reducing the background noise of the data and therefore on improving the frontal detection. The performance of the combinations on 4 km data are evaluated using two criteria: detection efficiency and front length. We find that the optimal combination of preliminary smoothing parameters in enhancing detection efficiency and preserving front length includes a median filter, a 16 × 16 pixel window size, and a 5 × 5 pixel kernel for strong fronts and a 7 × 7 pixel kernel for weak fronts. Results show an improvement in detection performance (from largest to smallest window size) of 71% for strong fronts and 120% for weak fronts. Despite the small window used (16 × 16 pixels), the length of the fronts has been preserved relative to that found with 1 km data. This optimal preliminary smoothing and the CMW detection algorithm on 4 km sea surface temperature data are then used to describe the spatial distribution of the monthly frequencies of occurrence for both strong and weak fronts across the Indian Ocean basin. In general strong fronts are observed in coastal areas whereas weak fronts, with some seasonal exceptions, are mainly located in the open ocean. This study shows that adequate noise reduction done by a preliminary smoothing of the data considerably improves the frontal detection efficiency as well as the global quality of the results. Consequently, the use of 4 km data enables frontal detections similar to 1 km data (using a standard median 3 × 3 convolution) in terms of detectability, length and location. This method, using 4 km data is easily applicable to large regions or at the global scale with far less constraints of data manipulation and processing time relative to 1 km data.

Small-scale variability of the current in the Strait of Bonifacio

Current dynamics in the Strait of Bonifacio (south Corsica) were investigated at a small scale during the STELLAMARE1 multidisciplinary cruise in summer 2012, using in situ measurements and modeling data. The Strait of Bonifacio is a particularly sensitive marine area in which specific conservation measures have been taken to preserve the natural environment and wild species. Good knowledge of the hydrodynamics in this area is essential to optimize the Marine Protected Area's management rules. Therefore, we used a high-resolution model (400 m) based on the MARS3D code to investigate the main flux exchanges and to formulate certain hypotheses about the formation of possible eddy structures. The aim of the present paper is first to synthetize the results obtained by combining Acoustic Doppler Current Profiler data, hydrological parameters, Lagrangian drifter data, and satellite observations such as MODIS OC5 chlorophyll a data or Metop-A AVHRR Sea Surface Temperature (SST) data. These elements are then used to validate the presence of the mesoscale eddies simulated by the model and their recurrence outside the cruise period. To complete the analysis, the response of the 3D hydrodynamical model was evaluated under two opposing wind systems and certain biases were detected. Strong velocities up to 1 m s(-1) were recorded in the east part due to the Venturi effect; a complementary system of vortices governed by Coriolis effect and west wind was observed in the west part, and horizontal stratification in the central part has been identified under typical wind condition.

Observations of turbulence within the surf and swash zone of a field-scale sandy laboratory beach

Current coastal-evolution models generally lack the ability to accurately predict bed level change in shallow (<~2 m) water, which is, at least partly, due to the preclusion of the effect of surface-induced turbulence on sand suspension and transport. As a first step to remedy this situation, we investigated the vertical structure of turbulence in the surf and swash zone using measurements collected under random shoaling and plunging waves on a steep (initially 1:15) field-scale sandy laboratory beach. Seaward of the swash zone, turbulence was measured with a vertical array of three Acoustic Doppler Velocimeters (ADVs), while in the swash zone two vertically spaced acoustic doppler velocimeter profilers (Vectrino profilers) were applied. The vertical turbulence structure evolves from bottom-dominated to approximately vertically uniform with an increase in the fraction of breaking waves to ~ 50%. In the swash zone, the turbulence is predominantly bottom-induced during the backwash and shows a homogeneous turbulence profile during uprush. We further find that the instantaneous turbulence kinetic energy is phase-coupled with the short-wave orbital motion under the plunging breakers, with higher levels shortly after the reversal from offshore to onshore motion (i.e. wavefront).

On Growth Rate of Wind Waves: Impact of Short-Scale Breaking Modulations

The wave generation model based on the rapid distortion concept significantly underestimates empirical values of the wave growth rate. As suggested before, inclusion of the aerodynamic roughness modulations effect on the amplitude of the slope-correlated surface pressure could potentially reconcile this model approach with observations. This study explores the role of short-scale breaking modulations to amplify the growth rate of modulating longer waves. As developed, airflow separations from modulated breaking waves result in strong modulations of the turbulent stress in the inner region of the modulating waves. In turn, this leads to amplifying the slope-correlated surface pressure anomalies. As evaluated, such a mechanism can be very efficient for enhancing the wind-wave growth rate by a factor of 2-3.

Origin and fate of surface drift in the oceanic convergence zones of the eastern Pacific

This study investigates the structure and intensity of the surface pathways connecting to and from the central areas of the large-scale convergence regions of the eastern Pacific Ocean. Surface waters are traced with numerical Lagrangian particles transported in the velocity field of three different ocean models with horizontal resolutions that range from ¼° to 1/32°. The connections resulting from the large-scale convergent Ekman dynamics agree qualitatively but are strongly modulated by eddy variability that introduces meridional asymmetry in the amplitude of transport. Lagrangian forward-in-time integrations are used to analyze the fate of particles originating from the central regions of the convergence zones and highlight specific outflows not yet reported for the southeastern Pacific when using the currents at the highest resolutions (1/12° and 1/32°). The meridional scales of these outflows are comparable to the characteristic width of the fine-scale striation of mean currents.

Seasonal changes of heat and moisture fluxes from the Caspian Sea surface and its effect on synoptic systems with heavy rainfalls in southern shores of the sea during 2005-2010

Caspian Sea with its unique characteristics is a significant source to supply required heat and moisture for passing weather systems over the north of Iran. Investigation of heat and moisture fluxes in the region and their effects on these systems that could lead to floods and major financial and human losses is essential in weather forecasting. Nowadays by improvement of numerical weather and climate prediction models and the increasing need to more accurate forecasting of heavy rainfall, the evaluation and verification of these models has been become much more important. In this study we have used the WRF model as a research-practical one with many valuable characteristics and flexibilities. In this research, the effects of heat and moisture fluxes of Caspian Sea on the synoptic and dynamical structure of 20 selective systems associated with heavy rainfall in the southern shores of Caspian Sea are investigated. These systems are selected based on the rainfall data gathered by three local stations named: Rasht, Babolsar and Gorgan in different seasons during a five-year period (2005-2010) with maximum amount of rainfall through the 24 hours of a day. In addition to synoptic analyses of these systems, the WRF model with and without surface flues was run using the two nested grids with the horizontal resolutions of 12 and 36 km. The results show that there are good consistencies between the predicted distribution of rainfall field, time of beginning and end of rainfall by the model and the observations. But the model underestimates the amounts of rainfall and the maximum difference with the observation is about 69%. Also, no significant changes in the results are seen when the domain and the resolution of computations are changed. The other noticeable point is that the systems are severely weakened by removing heat and moisture fluxes and thereby the amounts of large scale rainfall are decreased up to 77% and the convective rainfalls tend to zero.

Statistical Analysis of Sea Surface Temperature and Chlorophyll-a Concentration Patterns in the Gulf of Tadjourah (Djibouti)

The sea surface temperature (SST) and chlorophyll-a concentration (CHL-a) were analysed in the Gulf of Tadjourah from two set of 8-day composite satellite data, respectively from 2008 to 2012 and from 2005 to 2011. A singular spectrum analysis (SSA) shows that the annual cycle of SST is strong (74.3% of variance) and consists of warming (April-October) and cooling (November-March) of about 2.5C than the long-term average. The semi-annual cycle captures only 14.6% of temperature variance and emphasises the drop of SST during July-August. Similarly, the annual cycle of CHL-a (29.7% of variance) depicts high CHL-a from June to October and low concentration from November to May. In addition, the first spatial empirical orthogonal function (EOF) of SST (93% of variance) shows that the seasonal warming/cooling is in phase across the whole study area but the southeastern part always remaining warmer or cooler. In contrast to the SST, the first EOF of CHL-a (54.1% of variance) indicates the continental shelf in phase opposition with the offshore area in winter during which the CHL-a remains sequestrated in the coastal area particularly in the south-east and in the Ghoubet Al-Kharab Bay. Inversely during summer, higher CHL-a quantities appear in the offshore waters. In order to investigate processes generating these patterns, a multichannel spectrum analysis was applied to a set of oceanic (SST, CHL-a) and atmospheric parameters (wind speed, air temperature and air specific humidity). This analysis shows that the SST is well correlated to the atmospheric parameters at an annual scale. The windowed cross correlation indicates that this correlation is significant only from October to May. During this period, the warming was related to the solar heating of the surface water when the wind is low (April-May and October) while the cooling (November-March) was linked to the strong and cold North-East winds and to convective mixing. The summer drop in SST followed by a peak of CHL-a, seems strongly correlated to the upwelling. The second EOF modes of SST and CHL-a explain respectively 1.3% and 5% of the variance and show an east-west gradient during winter that is reversed during summer. This work showed that the seasonal signals have a wide spatial influence and dominate the variability of the SST and CHL-a while the east-west gradient are specific for the Gulf of Tadjourah and seem induced by the local wind modulated by the topography.

Optimization of the large scale synthesis of the LSF-20 cathode material for SOFCs

Póster presentado en: 21st World Hydrogen Energy Conference 2016. Zaragoza, Spain. 13-16th June, 2016

Catalytic CO2 hydrogenation: study at the molecular scale using field emission techniques

An accurate characterization of the rhodium specimen was obtained via FIM experiments. Reaction behaviors between H2 and CO2 were observed in FEM mode at 700 K. At this temperature, CO desorption occurs, preventing CO+H2 reaction. Surface is mainly recovered by oxygen; reaction with hydrogen occurs. Finally, we can identify the reaction as the Reverse Water Gas Shift.

The role of the Indian Ocean sector and sea surface salinity for prediction of the coupled Indo-Pacific system

The purpose of this dissertation is to evaluate the potential downstream influence of the Indian Ocean (IO) on El Niño/Southern Oscillation (ENSO) forecasts through the oceanic pathway of the Indonesian Throughflow (ITF), atmospheric teleconnections between the IO and Pacific, and assimilation of IO observations. Also the impact of sea surface salinity (SSS) in the Indo-Pacific region is assessed to try to address known problems with operational coupled model precipitation forecasts. The ITF normally drains warm fresh water from the Pacific reducing the mixed layer depths (MLD). A shallower MLD amplifies large-scale oceanic Kelvin/Rossby waves thus giving ~10% larger response and more realistic ENSO sea surface temperature (SST) variability compared to observed when the ITF is open. In order to isolate the impact of the IO sector atmospheric teleconnections to ENSO, experiments are contrasted that selectively couple/decouple the interannual forcing in the IO. The interannual variability of IO SST forcing is responsible for 3 month lagged widespread downwelling in the Pacific, assisted by off-equatorial curl, leading to warmer NINO3 SST anomaly and improved ENSO validation (significant from 3-9 months). Isolating the impact of observations in the IO sector using regional assimilation identifies large-scale warming in the IO that acts to intensify the easterlies of the Walker circulation and increases pervasive upwelling across the Pacific, cooling the eastern Pacific, and improving ENSO validation (r ~ 0.05, RMS~0.08C). Lastly, the positive impact of more accurate fresh water forcing is demonstrated to address inadequate precipitation forecasts in operational coupled models. Aquarius SSS assimilation improves the mixed layer density and enhances mixing, setting off upwelling that eventually cools the eastern Pacific after 6 months, counteracting the pervasive warming of most coupled models and significantly improving ENSO validation from 5-11 months. In summary, the ITF oceanic pathway, the atmospheric teleconnection, the impact of observations in the IO, and improved Indo-Pacific SSS are all responsible for ENSO forecast improvements, and so each aspect of this study contributes to a better overall understanding of ENSO. Therefore, the upstream influence of the IO should be thought of as integral to the functioning of ENSO phenomenon.