The position of the lowest levels in the boundary layer of atmospheric circulation models

Abstract

The error associated with the position of the lowest wind level in atmospheric boundary‐layer modelling is studied in connection with vertical resolutions typical of parametrization schemes for atmospheric circulation models. The test case is the neutrally stratified steady‐state boundary layer. Finite‐difference and finite‐element schemes of two types are used: in one case the lowest wind level is centred with respect to the surface and the lowest internal level where the shear stress is calculated, in the other case the lowest wind level is set very close to the surface (5 m). It is found that schemes of the latter type underestimate the friction force at the lowest level and consequently overestimate the wind and the surface stress. This error is largest at low resolution since it is due to the uncentring of the lowest wind level with respect to the stress levels. The error in schemes of the former type is different, and is associated with the determination of the surface stress from a wind at a height that may exceed the extent of the surface layer. For all neutral cases, this error can be made small by adding a corrective term to the traditional logarithmic formulation. This paper shows that considerably more error is generated by uncentred differencing than by deepening the surface layer.     

Displacements and transformations of nitrate-rich and nitrate-poor water masses in the tropical Pacific during the 1997 El Niño

A Lagrangian analysis was applied to the outputs of a coupled physical-biogeochemical model to describe the redistribution of nitrate-rich and nitrate-poor surface water masses in the tropical Pacific throughout the major 1997 El Niño. The same tool was used to analyze the causes of nitrate changes along trajectories and to investigate the consequences of the slow nitrate uptake in the high nutrient low chlorophyll (HNLC) region during the growth phase of the event. Three patterns were identified during the drift of water masses. The first mechanism is well known along the equator: oligotrophic waters from the western Pacific are advected eastward and retain their oligotrophic properties along their drift. The second concerns the persistent upwelling in the eastern basin. Water parcels have complex trajectories within this retention zone and remain mesotrophic. This study draws attention to the third process which is very specific to the HNLC region and to the El Niño period. During the 1997 El Niño, horizontal and vertical inputs of nitrate decreased so dramatically that nitrate uptake by phytoplankton became the only mechanism driving nitrate changes along pathways. The study shows that because of the slow nitrate uptake characteristic of the tropical Pacific HNLC system, nitrate in the pre-El Niño photic layer can support biological production for a period of several months. As a consequence, the slow nitrate uptake delays the gradual onset of oligotrophic conditions over nearly all the area usually occupied by upwelled waters. Owing to this process, mesotrophic conditions persist in the tropical Pacific during El Niño events.     

Long-term groundwater resource observatory for Southwestern Madagascar

Madagascar has one of the highest poverty rates in the world and consequently the long-term monitoring of groundwater resources is not a priority for the authorities. However, groundwater is often the only sustainable resource that has a satisfactory quality to supply the population. This is especially true in the south-west of the country, which is a semi-arid region and a global change hot spot (intense land use and climate changes). In response to the lack of data, the Groundwater Resource Observatory for Southwestern Madagascar (GROSoM) was established to monitor piezometry and meteorology over the longer term as part of a humanitarian response. The first site was setup in 2014 in a catchment located over a carbonate plateau; in 2018, a second site was installed in an alluvial setting within a crystalline basement catchment and a third site will be installed in 2020 to monitor groundwater dynamics in a coastal setting. The three sites, located between Toliara and Taolagnaro cities, are complementary and representative of various hydrogeological systems in Southwestern Madagascar. Each site includes a weather station and between 3 and 6 piezometric probes. The monitoring data indicate a strong inter-annual variability in precipitation, which induces a strong variability in aquifers recharge. One of the driest years in 2016 seems to be consistent with strong El Niño – Southern Oscillation (ENSO) effects observed at the global scale, while years with higher recharge appear to be related to cyclones such as Fundi in 2015 and Eketsang in 2019. Preliminary results of cross-disciplinary studies demonstrated a link between groundwater and health issues (i.e., admissions to basic health centres). This observatory aims to produce long-term data and has two objectives: (i) strengthening the early warning system for humanitarian crises in Madagascar; (ii) contributing to a better understanding of the effects of climate change on groundwater resources in this semi-arid region.     

Comparison of LiDAR waveform processing methods for very shallow water bathymetry using Raman,near‐infrared and green signals

Airborne light detection and ranging (LiDAR) bathymetry appears to be a useful technology for bed topography mapping of non‐navigable areas, offering high data density and a high acquisition rate. However, few studies have focused on continental waters, in particular, on very shallow waters (<2 m) where it is difficult to extract the surface and bottom positions that are typically mixed in the green LiDAR signal. This paper proposes two new processing methods for depth extraction based on the use of different LiDAR signals [green, near‐infrared (NIR), Raman] of the SHOALS‐1000T sensor. They have been tested on a very shallow coastal area (Golfe du Morbihan, France) as an analogy to very shallow rivers. The first method is based on a combination of mathematical and heuristic methods using the green and the NIR LiDAR signals to cross validate the information delivered by each signal. The second method extracts water depths from the Raman signal using statistical methods such as principal components analysis (PCA) and classification and regression tree (CART) analysis. The obtained results are then compared to the reference depths, and the performances of the different methods, as well as their advantages/disadvantages are evaluated. The green/NIR method supplies 42% more points compared to the operator process, with an equivalent mean error (?4·2 cm verusu ?4·5 cm) and a smaller standard deviation (25·3 cm verusu 33·5 cm). The Raman processing method provides very scattered results (standard deviation of 40·3 cm) with the lowest mean error (?3·1 cm) and 40% more points. The minimum detectable depth is also improved by the two presented methods, being around 1 m for the green/NIR approach and 0·5 m for the statistical approach, compared to 1·5 m for the data processed by the operator. Despite its ability to measure other parameters like water temperature, the Raman method needed a large amount of reference data to provide reliable depth measurements, as opposed to the green/NIR method. Copyright © 2010 John Wiley & Sons, Ltd.     

Evidence of an early alteration process driven by magmatic fluid in Mururoa volcano

Evidence for a deuteric alteration process induced by a magmatic fluid has been found in the feeder zone of the Mururoa volcano (French Polynesia). Within the dikes, where basaltic glass does not show any evidence of pervasive alteration, vesicles are filled with dioctahedral smectites and calcite, while olivine phenocrysts are replaced by dioctahedral smectites, ankerite and calcite.The ¹³C signature of carbonates, the carbon and H₂O content of the whole rocks and their impoverishment in deuterium are compatible with the presence of magmatic CO₂ during the crystallization of intruding lavas and exclude contamination by seawater. Mass balance calculations on selected thin sections photographs of partly filled up vesicles and replaced olivine crystals, constrain, assuming a closed system interaction, the chemical composition of the initial fluid and the respective amounts of the initial solid phases involved in the alteration process. Thermodynamic modelings using the EQ3/6 software package correctly predict the mineralogic, chemical and isotopic exchanges accompanying alteration, thus validating the closed system assumption. The model which allows prediction of the influence of CO₂ on the alteration products, shows that, above a 0.25 CO₂ mole fraction in the initial fluid, the alteration is entirely controlled by the chemical composition of the initial solid phases. The presence of CO₂ implies the precipitation of dioctahedral smectites and carbonates instead of the magnesian smectites commonly observed in CO₂-free systems.The Mururoa feeder zone shows alteration features typical of a closed system interaction between the basaltic rock and a magmatic fluid in which seawater did not take part.     

The snowball Earth aftermath: Exploring the limits of continental weathering processes

Carbonates capping Neoproterozoic glacial deposits contain peculiar sedimentological features and geochemical anomalies ascribed to extraordinary environmental conditions in the snowball Earth aftermath. It is commonly assumed that post-snowball climate dominated by CO₂ partial pressures several hundred times greater than modern levels, would be characterized by extreme temperatures, a vigorous hydrological cycle, and associated high continental weathering rates. However, the climate in the aftermath of a global glaciation has never been rigorously modelled. Here, we use a hierarchy of numerical models, from an atmospheric general circulation model to a mechanistic model describing continental weathering processes, to explore characteristics of the Earth system during the supergreenhouse climate following a snowball glaciation. These models suggest that the hydrological cycle intensifies only moderately in response to the elevated greenhouse. Indeed, constraints imposed by the surface energy budget sharply limit global mean evaporation once the temperature has warmed sufficiently that the evaporation approaches the total absorbed solar radiation. Even at 400 times the present day pressure of atmospheric CO₂, continental runoff is only 1.2 times the modern runoff. Under these conditions and accounting for the grinding of the continental surface by the ice sheet during the snowball event, the simulated maximum discharge of dissolved elements from continental weathering into the ocean is approximately 10 times greater than the modern flux. Consequently, it takes millions of years for the silicate weathering cycle to reduce post-snowball CO₂ levels to background Neoproterozoic levels. Regarding the origin of the cap dolostones, we show that continental weathering alone does not supply enough cations during the snowball melting phase to account for their observed volume.     

Identifying parameters controlling soil delayed behaviour from laboratory and in situ pressuremeter testing

The aim of this paper is to present a methodology for identifying the soil parameters controlling the delayed behaviour from laboratory and in situ pressuremeter tests by using an elasto‐viscoplastic model (EVP‐MCC) based on Perzyna's overstress theory and on the elasto‐plastic Modified Cam Clay model. The influence of both the model parameters and the soil permeability was studied under the loading condition of pressuremeter tests by coupling the proposed model equations with Biot's consolidation theory. On the basis of the parametric study, a methodology for identifying model parameters and soil permeability by inverse analysis from three levels of constant strain rate pressuremeter tests was then proposed and applied on tests performed on natural Saint‐Herblain clay. The methodology was validated by comparing the optimized values of soil parameters and the values of the same parameters obtained from laboratory test results, and also by using the identified parameters to simulate other tests on the same samples. The analysis of the drainage condition and the strain rate effect during a pressuremeter test demonstrated the coupled influence of consolidation and viscous effects on the test results. The numerical results also showed that the inverse analysis procedure could successfully determine the parameters controlling the time‐dependent soil behaviour. Copyright © 2008 John Wiley & Sons, Ltd.     

Pliocene extensional tectonics in the Eastern Central Patagonian Cordillera: geochronological constraints and new field evidence

Recent field work and review of radiometric data obtained from Neogene lavas and plutonic rocks exposed in the Eastern Central Patagonian Cordillera (46–48ºS), which overlie subducted segments of the South Chile Ridge, suggest important Late Miocene to Pleistocene morphological changes in relation to base level variations and/or tectonic events. We present new field observations from a region south of Lago General Carrera‐Buenos Aires, between the main Cordillera and the Meseta del Lago Buenos Aires, demonstrating that normal faulting controlled valley incisions and occurred during lava emplacement at 5–4 Ma and after 3 Ma. We also show that the 12 Ma basaltic flows of the Meseta del Lago Buenos Aires (∼2000 m a.s.l.) have been subjected to deep incision, with younger lavas dated at 1.2 Ma partially filling the valleys. These incisions are thought to reflect progressive eastward tilting of the entire meseta. Our new observations, together with additional features from Central Patagonia, strongly suggest that tectonic events led to a regional widespread morphological change after 5 Ma. The coincidence in time and space between the subduction of segments of the South Chile Ridge at 6 and 3 Ma causing opening of a slab window, and strong base level variations in the studied area, suggests a cause‐and‐effect relationship. In Central Patagonia, compressional tectonics ended well before extensional events reported here. Causes of uplift and further extension are probably completely disconnected. The uplift is purely tectonic in origin and occurred prior to the subduction of the South Chile Ridge. Extension should be a consequence of this subduction.