An iterative method for the construction of equilibrium N-body models for stellar disks

One widely used technique for the construction of equilibrium models of stellar disks is based on the Jeans equations and the moments of velocity distribution functions computed using these equations. Stellar disks constructed using this technique are shown to be “not entirely” in equilibrium. Our attempt to abandon the epicyclic approximation and the approximation of infinite isothermal layers, which are commonly adopted in this technique, failed to improve the situation substantially. We conclude that the main drawback of techniques based on the Jeans equations is that the system of equations employed is not closed and, therefore, requires adopting an essentially ad hoc additional closure condition. A new iterative approach to constructing equilibrium N-body models with a given density distribution is proposed. The main idea behind this approach is that a model is first constructed using some approximation method, and is then allowed to adjust to an equilibrium state with the specified density distribution—if necessary, with the required parameters of the velocity distribution remaining fixed in the process. This iterative approach was used to construct isotropic, spherically symmetric models and models of stellar disks embedded in an external potential. The numerical models constructed prove to be close to equilibrium. It is shown that the commonly adopted assumption that the profile of the radial-velocity dispersion is exponential may be wrong. The technique proposed can be applied to a wide range of problems involving the construction of models of stellar systems with various geometries.     

Temporal Variability of Total Ozone in the Asian Region Inferred from Ground-Based and Satellite Measurement Data

This paper reports investigation data on the temporal variability of total ozone content (TOC) in the Central Asian and Tibet Plateau mountain regions obtained by conventional methods, as well as by spectral, cross-wavelet, and composite analyses. The data of ground-based observation stations located at Huang He, Kunming, and Lake Issyk-Kul, along with the satellite data obtained at SBUV/SBUV2 (SBUV merged total and profile ozone data, Version 8.6) for 1980–2013 and OMI (Ozone Monitoring Instrument) and TOU (Total Ozone Unit) for 2009–2013 have been used. The average relative deviation from the SBUV/SBUV2 data is less than 1% in Kunming and Issyk-Kul for the period of 1980–2013, while the Huang He Station is characterized by an excess of the satellite data over the ground-based information at an average deviation of 2%. According to the Fourier analysis results, the distribution of amplitudes and the periods of TOC oscillations within a range of over 14 months is similar for all series analyzed. Meanwhile, according to the cross-wavelet and composite analyses results, the phase relationships between the series may considerably differ, especially in the periods of 5–7 years. The phase of quasi-decennial oscillations in the Kunming Station is close to the 11-year oscillations of the solar cycle, while in the Huang He and Issyk-Kul stations the TOC variations go ahead of the solar cycle.     

Water Resources Implications of Global Warming: A U.S. Regional Perspective

The implications of global warming for the performance of six U.S. water resource systems are evaluated. The six case study sites represent a range of geographic and hydrologic, as well as institutional and social settings. Large, multi-reservoir systems (Columbia River, Missouri River, Apalachicola-Chatahoochee-Flint (ACF) Rivers), small, one or two reservoir systems (Tacoma and Boston) and medium size systems (Savannah River) are represented. The river basins range from mountainous to low relief and semi-humid to semi-arid, and the system operational purposes range from predominantly municipal to broadly multi-purpose. The studies inferred, using a chain of climate downscaling, hydrologic and water resources systems models, the sensitivity of six water resources systems to changes in precipitation, temperature and solar radiation. The climate change scenarios used in this study are based on results from transient climate change experiments performed with coupled ocean-atmosphere General Circulation Models (GCMs) for the 1995 Intergovernmental Panel on Climate Change (IPCC) assessment. An earlier doubled-CO₂ scenario from one of the GCMs was also used in the evaluation. The GCM scenarios were transferred to the local level using a simple downscaling approach that scales local weather variables by fixed monthly ratios (for precipitation) and fixed monthly shifts (for temperature). For those river basins where snow plays an important role in the current climate hydrology (Tacoma, Columbia, Missouri and, to a lesser extent, Boston) changes in temperature result in important changes in seasonal streamflow hydrographs. In these systems, spring snowmelt peaks are reduced and winter flows increase, on average. Changes in precipitation are generally reflected in the annual total runoff volumes more than in the seasonal shape of the hydrographs. In the Savannah and ACF systems, where snow plays a minor hydrological role, changes in hydrological response are linked more directly to temperature and precipitation changes. Effects on system performance varied from system to system, from GCM to GCM, and for each system operating objective (such as hydropower production, municipal and industrial supply, flood control, recreation, navigation and instream flow protection). Effects were generally smaller for the transient scenarios than for the doubled CO₂ scenario. In terms of streamflow, one of the transient scenarios tended to have increases at most sites, while another tended to have decreases at most sites. The third showed no general consistency over the six sites. Generally, the water resource system performance effects were determined by the hydrologic changes and the amount of buffering provided by the system's storage capacity. The effects of demand growth and other plausible future operational considerations were evaluated as well. For most sites, the effects of these non-climatic effects on future system performance would about equal or exceed the effects of climate change over system planning horizons.     

Intergranular diffusion rates from the analysis of garnet surfaces: Implications for metamorphic equilibration

Novel approaches to garnet analysis have been used to assess rates of intergranular diffusion between different matrix phases and garnet porphyroblasts in a regionally metamorphosed staurolite‐mica‐schist from the Barrovian‐type area in Scotland. X‐ray maps and chemical traverses of planar porphyroblast surfaces reveal chemical heterogeneity of the garnet grain boundary linked to the nature of the adjacent matrix phase. The garnet preserves evidence of low temperature retrograde exchange with matrix minerals and diffusion profiles documenting cation movement along the garnet boundaries. Garnet–quartz and garnet–plagioclase boundaries preserve evidence of sluggish Mg, Mn and Fe diffusion at comparable rates to volume diffusion in garnet, whereas diffusion along garnet–biotite interfaces is much more effective. Evidence of particularly slow Al transport, probably coupled to Fe³⁺ exchange, is locally preserved on garnet surfaces adjacent to Fe‐oxide phases. The Ca distribution on the garnet surface shows the most complex behaviour, with long‐wavelength heterogeneities apparently unrelated to the matrix grain boundaries. This implies that the Ca content of garnet is controlled by local availability and is thought likely to reflect disequilibrium established during garnet growth. Geochemical anomalies on the garnet surfaces are also linked to the location of triple junctions between the porphyroblasts and the matrix phases, and imply enhanced transport along these channels. The slow rates of intergranular diffusion and the characteristics of different boundary types may explain many features associated with the prograde growth of garnet porphyroblasts. Thus, minerals such as quartz, Fe‐oxides and plagioclase whose boundaries with garnet are characterized by slow intergranular diffusion rates appear to be preferentially trapped as inclusions within porphyroblasts. As such grain boundary diffusion rates may be a significant kinetic impediment to metamorphic equilibrium and garnet may struggle to maintain chemical and textural equilibrium during growth in pelites.     

Quaternary sediments in the axial zone of the Brazil Basin

Lithological and micropaleontological studies of sediments were carried out along the meridional profile across the Brazil Basin. Based on nannoplankton and diatom assemblages, the sediments represented by oxidized miopelagic clays, differently reduced hemipelagic clays, and clayey-siliceous (Ethmodiscus) and calcareous (coccolith-foraminiferal) oozes are attributed to the Pleistocene-Holocene. Specific features of sedimentary material indicate its redeposition by Antarctic bottom waters, mudflows, slumps, and rockslides. It is shown that sedimentation in the Brazil Basin corresponds to the incomplete pelagic (miopelagic) type of oceanic lithogenesis transitional to the near-continental one.     

Revisiting the morphological method in two-dimensions to quantify bed-material transport in braided rivers

Research in the 1990s showed that bed-material transport rates could be estimated at the reach scale in both one-dimension and, over small spatial scales (10s of m), in two-dimensions. The limit on the latter was the spatial scale over which it was possible to obtain distributed data on morphological change. Here, we revisit the morphological method given progress in both topographical data acquisition and hydraulic modelling. The bed-material transport needed to conserve mass is calculated in both one and two dimensions for a 1600 m × 300 m Alpine braided river “laboratory”. High-resolution topographical data were acquired by laser scanning to quantify Digital Elevation Models (DEMs), and morphological changes caused by the flushing of the water intake were derived from repeated surveys. Based on DEMs of differences, 1D bed-material transport rates were calculated using the morphological method. Then, a 2D hydraulic model was combined with a topographic correction to route sediment through the network of braided channels and to obtain a spatially variable estimate of transport in both downstream and cross-stream directions. Monte Carlo simulation was applied to the routing model parameters, allowing identification of the most probable parameter values needed to minimize negative transport. The results show that within-section spatial compensation of erosion and deposition using the 1D treatment leads to substantial local errors in transport rate estimates, to a degree related to braiding intensity. Even though the 2D application showed that a large proportion of the total transport was actually concentrated into one main channel during the studied low flow event, the proportion of transport in secondary anabranches is substantial when the river starts braiding. Investigations of the effects of DEM resolution, competent flow duration and survey frequency related to ‘travelling bedload’ and sequential erosion-deposition emphasized the critical importance of careful data collection in the application of the morphological method. © 2019 John Wiley & Sons, Ltd.     

Earthquake swarms preceding the 2000 eruption of Miyakejima volcano, Japan

The first sign of magma accumulating beneath Miyakejima, an island volcano in the northern Izu islands, Japan, came at around 18:00 on 26 June 2000, when a swarm of earthquakes was detected by a volcano seismic network on the island. Earthquakes occurred initially beneath the southwest flank near the summit and gradually migrated west of the island, where a submarine eruption occurred the next morning. Earthquakes then migrated further to the northwest between Miyakejima and Kozushima, another volcanic island and developed to the most intense earthquake swarm ever observed in and around Japanese archipelago. To better image how the initial magma intrusion occurred, we relocated hypocenters by using a station-correction method and a double-difference method. The relocated epicenters are generally concentrated near the upper bound of dyke intrusions inferred from geodetic studies throughout the initial stages of the 2000 eruption at Miyakejima from 26 to 27 June 2000. As for seismic activity westward off Miyakejima in the morning on 27 June, hypocenters from both a nationwide seismic network that were relocated by the double-difference method, and those from the volcano seismic network relocated by the station-correction method, formed a very shallow cluster that ascended slowly with time as it propagated northwestward from Miyakejima. This suggests that the dykes have both a radial and upward component of movement.Editorial responsibility: S. Nakada, T. Druitt     

Review: The Yucat��n Peninsula karst aquifer, Mexico

The Yucatán Peninsula karst aquifer is one of the most extensive and spectacular karst aquifer systems on the planet. This transboundary aquifer system extends over an area of approximately 165,000? km² in México, Guatemala and Belize. The Triassic to Holocene Yucatán limestone platform is located in the vicinity of the North American/Caribbean plate boundary and has been reshaped by a series of tectonic events over its long geologic history. At the end of the Cretaceous period, the Yucatán Peninsula was hit by a large asteroid, which formed the Chicxulub impact crater. The Yucatán Peninsula karst aquifer hosts large amounts of groundwater resources which maintain highly diverse groundwater-dependent ecosystems. Large parts of the aquifer are affected by seawater intrusion. Anthropogenic pollution of the aquifer has been increasing over the past few decades, owing to relentless economic development and population growth on the Peninsula. This review summarizes the state of knowledge on the Yucatán Peninsula karst aquifer and outlines the main challenges for hydrologic research and practical groundwater-resources management on the Peninsula.