Sea and Lake Breezes: A Review of Numerical Studies

Numerical studies of sea and lake breezes are reviewed. The modelled dependence of sea-breeze and lake-breeze characteristics on the land surface sensible heat flux, ambient geostrophic wind, atmospheric stability and moisture, water body dimensions, terrain height and slope, Coriolis parameter, surface roughness length, and shoreline curvature is discussed. Consensus results on the influence of these geophysical variables on sea and lake breezes are synthesized as well as current gaps in our understanding. A brief history of numerical modelling, an overview of recent high-resolution simulations, and suggestions for future research related to sea and lake breezes are also presented. The results of this survey are intended to be a resource for numerical modelling, coastal air quality, and wind power studies.     

Regional climate model simulation of the AMMA Special Observing Period #3 and the pre-Helene easterly wave

The study examines results of dynamic downscaling of two global analyses: the National Center for Environmental Prediction/National Center for Atmospheric Research (NCEP/NCAR) reanalysis II and the Global Forecast System final analysis (FNL). Downscaling to a 0.5° grid over West Africa and the adjacent Atlantic Ocean is accomplished by each of two regional models, the Regional Model, version 3 (RM3) of the Center for Climate Systems Research and the Weather, Research and Forecasting model (WRF). Simulations are for September 2006, the African Monsoon Multidisciplinary Analysis (AMMA) Special Observing Period #3 (SOP-3). The aim of this study is to exploit the increased spatial detail in the simulations and representations of climate fields by the regional models to analyze meteorological systems within the SOP-3 area of interest and time frame. In particular, the paper focuses on the regional models’ representations of the structure and movement of a prominent easterly wave during September 10–13th, the precursor of Tropical Storm/Hurricane Helene. It describes the RM3 simulated structure of the developing storm in terms of circulation, precipitation, vertical motion, cumulus heating rates, and cross-sections of wind and geopotential height anomalies. Simulated cumulus heating rates within the wave’s main precipitation area imply a lowering of the bases of active cumulus in the transition from the African continent to the Atlantic, indicating that the ocean environment promotes greater upward latent heat flux that in turn intensifies overlying storms. RM3 circulation, precipitation patterns, and storm trajectory are reasonably consistent with observational evidence. Experiments show that precipitation rates near 6°N over the eastern North Atlantic are sensitive to vertical thermal stability, such that they are enhanced by warmer in situ sea-surface temperatures (SSTs) and diminished by colder SSTs. However, prescribing colder SST causes increases in precipitation north of 9°N within areas of large scale upward vertical motion where rainfall rates are less sensitive to in situ SSTs. The evaluation of WRF indicates that its storm propagation is too fast over West Africa, where associated WRF precipitation rates are exaggerated, but its performance is improved over the Atlantic.     

Low-speed impacts between rubble piles modeled as collections of polyhedra, 2

We present the results of additional calculations involving the collisions of km-scale rubble piles. In new work, we used the Open Dynamics Engine (ODE), an open-source library for the simulation of rigid-body dynamics that incorporates a sophisticated collision-detection and resolution routine. We found that using ODE resulted in a speed-up of approximately a factor of 30 compared with previous code. In this paper we report on the results of almost 1200 separate runs, the bulk of which were carried out with 1000-2000 elements. We carried out calculations with three different combinations of the coefficients of friction η and (normal) restitution ?: low (η=0,?=0.8), medium (η=0,?=0.5), and high (η=0.5,?=0.5) dissipation.For target objects of ∼1 km in radius, we found reduced critical disruption energy values in head-on collisions from 2 to 100 J kg⁻¹ depending on dissipation and impactor/target mass ratio. Monodisperse objects disrupted somewhat more easily than power-law objects in general. For oblique collisions of equal-mass objects, mildly off-center collisions (b/b₀=0.5) seemed to be as efficient or possibly more efficient at collisional disruption as head-on collisions. More oblique collisions were less efficient and the most oblique collisions we tried (b/b₀=0.866) required up to ∼200 J kg⁻¹ for high-dissipation power-law objects. For calculations with smaller numbers of elements (total impactor or 200 elements) we found that collisions were more efficient for smaller numbers of more massive elements, with values as low as for low-dissipation cases. We also analyzed our results in terms of the relations proposed by Stewart and Leinhardt [Stewart, S.T., Leinhardt, Z.M., 2009. Astrophys. J. 691, L133-L137] where where Q_R is the impact kinetic energy per unit total mass m_i+m_T. Although there is a significant amount of scatter, our results generally bear out the suggested relation.     

Idealized Large-Eddy Simulations of Sea and Lake Breezes: Sensitivity to Lake Diameter, Heat Flux and Stability

Idealized large-eddy simulations of lake and sea breezes are conducted to determine the sensitivity of these thermally-driven circulations to variations in the land-surface sensible heat flux and initial atmospheric stability. The lake-breeze and sea-breeze metrics of horizontal wind speed, horizontal extent, and depth are assessed. Modelled asymmetries about the coastline in the horizontal extent of the low-level onshore flow are found to vary as a function of the heat flux and stability. Small lake breezes develop similarly to sea breezes in the morning, but have a significantly weaker horizontal wind-speed component and a smaller horizontal extent than sea breezes in the afternoon.     

Quality of geological CO2 storage to avoid jeopardizing climate targets

We explore allowable leakage for carbon capture and geological storage to be consistent with maximum global warming targets of 2.5 and 3 °C by 2100. Given plausible fossil fuel use and carbon capture and storage scenarios, and based on modeling of time-dependent leakage of CO₂, we employ a climate model to calculate the long-term temperature response of CO₂ emissions. We assume that half of the stored CO₂ is permanently trapped by fast mechanisms. If 40?% of global CO₂ emissions are stored in the second half of this century, the temperature effect of escaped CO₂ is too small to compromise a 2.5 °C target. If 80?% of CO₂ is captured, escaped CO₂ must peak 300?years or later for consistency with this climate target. Due to much more CO₂ stored for the 3 than the 2.5 °C target, quality of storage becomes more important. Thus for the 3 °C target escaped CO₂ must peak 400?years or later in the 40?% scenario, and 3000?years or later in the 80?% scenario. Consequently CO₂ escaped from geological storage can compromise the less stringent 3 °C target in the long-run if most of global CO₂ emissions have been stored. If less CO₂ is stored only a very high escape scenario can compromise the more stringent 2.5 °C target. For the two remaining combinations of storage scenarios and climate targets, leakage must be high to compromise these climate targets.     

Geochemical distribution of arsenic,cadmium, lead and zinc in river sediments affected by tailings in Zimapán,a historical polymetalic mining zone of México

High arsenic (As) groundwater is widely distributed in northwestern Hetao Plain, an arid region with sluggish groundwater flow. Observed As concentration in groundwater from wells ranges from 76 to 1,093 μg/l. Most water samples have high total dissolved solids, with Cl and HCO₃ as the dominant anions and Na as the dominant cation. The major hydrochemical types of most saline groundwaters are Na–Mg–Cl–HCO₃ and Na–Mg–Cl. By contrast, fresh groundwaters generally belong to the Na–Mg–HCO₃ type. High concentrations of arsenic in shallow aquifers are associated with strongly reducing conditions, as evidenced by high concentrations of dissolved organic carbon, ammonium, as well as dissolved sulfide and Fe, dominance of arsenite, relatively low concentrations of nitrate and sulfate, and occasionally high content of dissolved methane (CH₄). High As groundwaters from different places at Hetao Plain experienced different redox processes. Fluoride is also present in high As groundwater, ranging between 0.40 and 3.36 mg/l. Although fluorosis poses an additional health problem in the region, it does not correlate well with As in spatial distribution. Geochemical analysis indicates that evapotranspiration is an important process controlling the enrichment of Na and Cl, as well as trace elements such as As, B, and Br in groundwater. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Similar-sized collisions and the diversity of planets

It is assumed in models of terrestrial planet formation that colliding bodies simply merge. From this the dynamical and chemical properties (and habitability) of finished planets have been computed, and our own and other planetary systems compared to the results of these calculations. But efficient mergers may be exceptions to the rule, for the similar-sized collisions (SSCs) that dominate terrestrial planet formation, simply because moderately off-axis SSCs are grazing; their centers of mass overshoot. In a “hit and run” collision the smaller body narrowly avoids accretion and is profoundly deformed and altered by gravitational and mechanical torques, shears, tides, and impact shocks. Consequences to the larger body are minor in inverse proportion to its relative mass. Over the possible impact angles, hit-and-run is the most common outcome for impact velocities v_imp between $～1.2$ and 2.7 times the mutual escape velocity v_esc between similar-sized planets. Slower collisions are usually accretionary, and faster SSCs are erosive or disruptive, and thus the prevalence of hit-and-run is sensitive to the velocity regime during epochs of accretion. Consequences of hit-and-run are diverse. If barely grazing, the target strips much of the exterior from the impactor—any atmosphere and ocean, much of the crust—and unloads its deep interior from hydrostatic pressure for about an hour. If closer to head-on $(～30\text{–}45°)$ a hit-and-run can cause the impactor core to plow through the target mantle, graze the target core, and emerge as a chain of diverse new planetoids on escaping trajectories. A hypothesis is developed for the diversity of next-largest bodies (NLBs) in an accreting planetary system—the bodies from which asteroids and meteorites derive. Because nearly all the NLBs eventually get accreted by the largest (Venus and Earth in our terrestrial system) or by the Sun, or otherwise lost, those we see today have survived the attrition of merger, evolving with each close call towards denser and volatile-poor bulk composition. This hypothesis would explain the observed density diversity of differentiated asteroids, and of dwarf planets beyond Neptune, in terms of episodic global-scale losses of rock or ice mantles, respectively. In an event similar to the Moon-forming giant impact, Mercury might have lost its original crust and upper mantle when it emerged from a modest velocity hit and run collision with a larger embryo or planet. In systems with super-Earths, profound diversity and diminished habitability is predicted among the unaccreted Earth-mass planets, as many of these will have be stripped of their atmospheres, oceans and crusts.     

Subducted seamounts in an eclogite-facies ophiolite sequence: the Andean Raspas Complex, SW Ecuador

The metamorphic Raspas Complex of southwest Ecuador consists of high-pressure mafic, ultramafic, and sedimentary rocks. The Lu–Hf ages of a blueschist, a metapelite, and an eclogite overlap at around 130 Ma and date high-pressure garnet growth. Peak metamorphic conditions in the eclogites reached 1.8 GPa at 600°C, corresponding to a maximum burial depth of ~60 km. The geochemical signatures of the eclogites suggest that their protoliths were typical mid-ocean ridge basalts (MORB), whereas the blueschists exhibit seamount-like characteristics, and the eclogite-facies peridotites seem to represent depleted, MORB-source mantle. That these rocks were subjected to similar peak PT conditions contemporaneously suggests that they were subducted together as an essentially complete section within the slab. We suggest that this section became dismembered from the slab during burial at great depth—perhaps as a consequence of scraping off the seamounts. The spatially close association of MORB-type eclogite, seamount-type blueschist, serpentinized peridotite, and metasediments points to an exhumed high-pressure ophiolite sequence.     

Timing of eclogite facies metamorphism in the southernmost Scandinavian Caledonides by Lu–Hf and Sm–Nd geochronology

The Scandinavian Caledonides contain several non-cogenetic eclogite-bearing terranes that were metamorphosed before the main orogenic stage in Scandian time (430–395 Ma). Although petrological and geochronological data from these terranes have provided essential information on the geodynamic history of the Caledonian orogenic cycle, the general picture is still patchy. To refine existing geodynamic models, we have dated the eclogite occurrence in the Jæren nappe, SW Norway, by Lu–Hf and Sm–Nd geochronology. Five out of the six studied samples provide a weighted mean Lu–Hf age of 469.9 ± 1.2 Ma (±2σ). One sample provided a significantly younger age of 457.9 ± 2.4 Ma. Garnet from the younger sample grew exclusively at eclogite facies conditions. In contrast, garnet from the other samples comprises prograde cores and peak metamorphic rims. Age estimates that take Lu-contributions of each zone into account provide an age of 471.0 ± 0.9 Ma for the cores and suggest a ca. 455 Ma age for the rims, which is identical to the bulk-garnet age of the younger sample. The same pattern is indicated by Sm–Nd ages, although these are relatively imprecise and reflect isotopic disturbance during thermal overprinting upon exhumation. The data define a new high-pressure age population for the Scandinavian Caledonides, which allows more detailed insight into the subduction history that affected the Baltoscandian margin before Scandian continental collision. Furthermore, this study highlights the potential complexities involved in garnet geochronology and shows the strength of Lu–Hf dating for unraveling the geochronological record of HP rocks.     

Sulfur isotope fractionation during the May 2003 eruption of Anatahan volcano, Mariana Islands: Implications for sulfur sources and plume processes

Sulfur isotope compositions of pumice and adsorbed volatiles on ash from the first historical eruption of Anatahan volcano (Mariana arc) are presented in order to constrain the sources of sulfur erupted during the period 10-21 May, 2003. The isotopic composition of S extracted from erupted pumice has a narrow range, from δ³⁴S_V-CDT +2.6‰ to +3.2‰, while the composition of sulfur adsorbed onto ash has a larger range (+2.8‰ to +5.3‰). Fractionation modeling for closed and open system scenarios suggests that degassing of SO₂ raised the δ³⁴S_V-CDT value of S dissolved in the melt from an initial composition of between +1.6‰ and +2.6‰ for closed-system degassing, or between −0.5‰ and +1.5‰ for open-system degassing, however closed-system degassing is the preferred model. The calculated values for the initial composition of the magma represent a MORB-like (δ³⁴S_V-CDT ∼ 0‰) mantle source with limited contamination by subducted seawater sulfate (δ³⁴S_V-CDT +21‰). Modeling also suggests that the δ³⁴S_V-CDT value of SO₂ gas in closed-system equilibrium with the degassed magma was between +0.9‰ and +2.5‰. The δ³⁴S_V-CDT value of sulfate adsorbed onto ash in the eruption plume (+2.8‰ to +5.1‰) is consistent with sulfate formation by oxidation of magmatic SO₂ in the eruption column. The sulfur isotope composition of sulfate adsorbed to ash changes from lower δ³⁴S values for ash erupted early in the eruption to higher δ³⁴S values for ash erupted later in the eruption. We interpret the temporal/stratigraphic change in sulfate isotopic composition to primarily reflect a change in the isotopic composition of magmatic SO₂ released from the progressively degassing magma and is attributed to the expulsion of an accumulated gas phase at the beginning of the eruption. More efficient oxidation of magmatic SO₂ gas to sulfate in the early water-rich eruption plume probably contributed to the change in S isotope compositions observed in the ash leachates.