Data assimilation using the particle filter for identifying the elasto‐plastic material properties of geomaterials

A computational method, incorporating the finite element model (FEM) into data assimilation using the particle filter, is presented for identifying elasto‐plastic material properties based on sequential measurements under the known changing traction boundary conditions to overcome some difficulties in identifying the parameters for elasto‐plastic problems from which the existing inverse analysis strategies have suffered. A soil–water coupled problem, which uses the elasto‐plastic constitutive model, is dealt with as the geotechnical application. Measured data on the settlement and the pore pressure are obtained from a synthetic FEM computation as the forward problem under the known parameters to be identified for both the element tests and the ground behavior during the embankment construction sequence. Parameter identification for elasto‐plastic problems, such as soil behavior, should be made by considering the measurements of deformation and/or pore pressure step by step from the initial stage of construction and throughout the deformation history under the changing traction boundary conditions because of the embankment or the excavation because the ground behavior is highly dependent on the loading history. Thus, it appears that sequential data assimilation techniques, such as the particle filter, are the preferable tools that can provide estimates of the state variables, that is, deformation, pore pressure, and unknown parameters, for the constitutive model in geotechnical practice. The present paper discusses the priority of the particle filter in its application to initial/boundary value problems for elasto‐plastic materials and demonstrates a couple of numerical examples. Copyright © 2012 John Wiley & Sons, Ltd.     

On the Contraction of Protostellar Clouds with Different Metallicities

We examine the thermal and chemical evolution of gravitationally collapsing protostellar clouds with metallicity 0≤Z/Z _⊙≤1.During the first collapse stage, the temperatures are higher for lower metallicity clouds. However, in the course of the adiabatic contraction of transient cores, the evolutionary trajectories of the clouds converge to a curve that is determined only by fundamental physical constants. The trajectories coincide each other thereafter. The size of the stellar core at formation is the same regardless of metallicity. This revised version was published online in July 2006 with corrections to the Cover Date.     

Origin of near-surface high-salinity water observed in the Kuroshio Extension region

Hydrographic data in the Kuroshio Extension (KE) region from 2008 to 2010 show large year-to-year variability in near-surface salinity, including a very large anomalous event in February 2010. During this event, the deep winter mixed layer in the southern KE region had higher salinity than had existed during the previous summer in September 2009. Our analysis shows that advection from the Philippine Sea along the western branch of the North Pacific subtropical gyre, taking approximately 9 months, resulted in this large salinity anomaly in February 2010 and contributes to the interannual salinity variability in the southern KE region.     

Regional climatic effects according to different estimations of biogenic volatile organic compounds during the asian summer monsoon

A series of 60-year numerical experiments starting from 1851 was conducted using a global climate model coupled with an aerosol-cloud-radiation model to investigate the response of the Asian summer monsoon to variations in the secondary organic aerosol (SOA) flux induced by two different estimations of biogenic volatile organic compound (BVOC) emissions. One estimation was obtained from a pre-existing archive and the other was generated by a next-generation model (the Model of Emissions of Gases and Aerosols from Nature, MEGAN). The use of MEGAN resulted in an overall increase of the SOA production through a higher rate of gasto-particle conversion of BVOCs. Consequently, the atmospheric loading of organic carbon (OC) increased due to the contribution of SOA to OC aerosol. The increase of atmospheric OC aerosols was prominent in particular in the Indian subcontinent and Indochina Peninsula (IP) during the pre- and early-monsoon periods because the terrestrial biosphere is the major source of BVOC emissions and the atmospheric aerosol concentration diminishes rapidly with the arrival of monsoon rainfall. As the number of atmospheric OC particles increased, the number concentrations of cloud droplets increased, but their size decreased. These changes represent a combination of aerosol-cloud interactions that were favorable to rainfall suppression. However, the modeled precipitation was slightly enhanced in May over the oceans that surround the Indian subcontinent and IP. Further analysis revealed that a compensating updraft in the surrounding oceans was induced by the thermally-driven downdraft in the IP, which was a result of surface cooling associated with direct OC aerosol radiative forcing, and was able to surpass the aerosolcloud interactions. The co-existence of oceanic ascending motion with the maximum convective available potential energy was also found to be crucial for rainfall formation. Although the model produced statistically significant rainfall changes with locally organized patterns, the suggested pathways should be considered guardedly because in the simulation results, 1) the BVOC-induced aerosol direct effect was marginal; 2) cloud-aerosol interactions were modeldependent; and 3) Asian summer monsoons were biased to a nonnegligible extent.     

Star formation triggered by supernova explosions in young galaxies

We study the evolution of supernova remnants in a low-metallicity medium   Z /Z_⊙= 10⁻⁴ to 10⁻²  in the early universe, using one-dimensional hydrodynamics with non-equilibrium chemistry. Once a post-shock layer is able to cool radiatively, a dense shell forms behind the shock. If this shell becomes gravitationally unstable and fragments into pieces, next-generation stars are expected to form from these fragments. To explore the possibility of this triggered star formation, we apply a linear perturbation analysis of an expanding shell to our results and constrain the parameter range of ambient density, explosion energy and metallicity where fragmentation of the shell occurs. For the explosion energy of  10⁵¹ erg (10⁵² erg)  , the shell fragmentation occurs for ambient densities higher than  ≳10² cm⁻³ (10 cm⁻³  ), respectively. This condition depends little on the metallicity in the ranges we examined. We find that the mode of star formation triggered occurs only in massive  (≳10⁸ M_⊙)  haloes.     

Foldback reflectors near methane hydrate bottom‐simulating reflectors: Indicators of gas distribution from 3D seismic images in the eastern Nankai Trough

Understanding of fluid behavior and gas distribution in the shallow subsurface are important considerations in gas hydrate formation and the global carbon cycle. Estimation of gas distribution based on reflection seismic surveys, however, is difficult because the boundary of a gas‐bearing zone is indistinct and not systematically defined. This study reports distinctive features related to gas‐hydrate distribution and possible fluid migration in high‐resolution 3D seismic‐reflection data from sediments of the eastern Nankai Trough. These features, here termed foldback reflectors (FBRs), descend in accordion shaped reflectors near the edges of bottom‐simulating reflectors (BSRs). FBRs generally correspond to lateral boundaries between two seismic facies, a ‘dimmed’ facies with relatively low amplitude and subdued high‐frequency components beneath the BSR and the contrasting facies around the BSR. The dimmed facies corresponds to areas of anomalously low velocity consistent with a small amount of free gas. FBR is mostly developed in well‐stratified formations in uplifted regions. Dip directions of the FBR appear to be restricted by orientation of the host formations. Edges of the FBR often correspond to high‐amplitude layers. Such occurrences of FBR suggest that regional uplift and layer‐parallel fluid migration are related to the formation of FBR as well as BSR.     

Climatic physical snowpack properties for large-scale modeling examined by observations and a physical model

Here we have conducted an integral study using site observations and a model with detailed snow dynamics, to examine the capability of the model for deriving a simple relationship between the density and thermal conductivity of the snowpack within different climatic zones used in large-scale climate modeling. Snow and meteorological observations were conducted at multiple sites in different climatic regions (two in Interior Alaska, two in Japan). A series of thermal conductivity measurements in snow pit observations done in Alaska provided useful information for constructing the relationship. The one-dimensional snow dynamics model, SNOWPACK, simulated the evolution of the snowpack and compared observations between all sites. Overall, model simulations tended to underestimate the density and overestimate the thermal conductivity, and failed to foster the relationship evident in the observations from the current and previous research. The causes for the deficiency were analyzed and discussed, regarding a low density of the new snow layer and a slow compaction rate. Our working relationships were compared to the equations derived by previous investigators. Discrepancy from the regression for the melting season observations in Alaska was found in common.     

The countermeasures against coastal hazards in Japan

The total length of Japanese coastline is about 34,500 km. Most of the coastal lowlands are now so intensively used that a large proportion of the Japanese coast has been artificially modified. Therefore, coastal countermeasures against several kinds of natural and man-induced hazards are the serious problems for Japanese. When we go down to the beaches, there are many kinds of armor blocks, breakwaters and dykes for the shore protection. Among the Japanese coastline, 15,900 km needs some protection of which 9,400 km was protected with some kinds of artificial structures in 1992. It is probable that no nation has used armor blocks so extensively as Japan.     

Fouling community dynamics in Lynnhaven Bay,Virginia

Acrylic panels were employed to follow fouling community development at Channel Points, Lynnhaven Bay, Virginia, from April 1977 through March 1979. Panel treatments were devised to study the effects of predator exclusion and panel area on community development. An additional panel was introduced in April 1978 to test the repeatability of initial community development. Early development was similar both years with hydroids, barnacles, and the polychaetePolydora ligni the initial macroinvertebrate colonists. Early community development in both years was dominated by the solitary tunicateMolgula manhattensis. Molgula was able to invade and rapidly overgrow all resident speices except largeBalanus eburneus. Panel area was shown to affect the community structure and development. The community on caged panels (predator exclusion) was significantly different from controls, although the differences could not be unequivocally attributed to predation.     

Pressure dependence of electrical conductivity of (Mg,Fe)SiO<Subscript>3</Subscript> ilmenite

The electrical conductivity of (Mg_0.93Fe_0.07)SiO₃ ilmenite was measured at temperatures of 500–1,200 K and pressures of 25–35 GPa in a Kawai-type multi-anvil apparatus equipped with sintered diamond anvils. In order to verify the reliability of this study, the electrical conductivity of (Mg_0.93Fe_0.07)SiO₃ perovskite was also measured at temperatures of 500–1,400 K and pressures of 30–35 GPa. The pressure calibration was carried out using in situ X-ray diffraction of MgO as pressure marker. The oxidation conditions of the samples were controlled by the Fe disk. The activation energy at zero pressure and activation volume for ilmenite are 0.82(6) eV and −1.5(2) cm³/mol, respectively. Those for perovskite were 0.5(1) eV and −0.4(4) cm³/mol, respectively, which are in agreement with the experimental results reported previously. It is concluded that ilmenite conductivity has a large pressure dependence in the investigated P–T range.