Pollutant Plume Dispersion in the Atmospheric Boundary Layer over Idealized Urban Roughness

The Gaussian model of plume dispersion is commonly used for pollutant concentration estimates. However, its major parameters, dispersion coefficients, barely account for terrain configuration and surface roughness. Large-scale roughness elements (e.g. buildings in urban areas) can substantially modify the ground features together with the pollutant transport in the atmospheric boundary layer over urban roughness (also known as the urban boundary layer, UBL). This study is thus conceived to investigate how urban roughness affects the flow structure and vertical dispersion coefficient in the UBL. Large-eddy simulation (LES) is carried out to examine the plume dispersion from a ground-level pollutant (area) source over idealized street canyons for cross flows in neutral stratification. A range of building-height-to-street-width (aspect) ratios, covering the regimes of skimming flow, wake interference, and isolated roughness, is employed to control the surface roughness. Apart from the widely used aerodynamic resistance or roughness function, the friction factor is another suitable parameter that measures the drag imposed by urban roughness quantitatively. Previous results from laboratory experiments and mathematical modelling also support the aforementioned approach for both two- and three-dimensional roughness elements. Comparing the UBL plume behaviour, the LES results show that the pollutant dispersion strongly depends on the friction factor. Empirical studies reveal that the vertical dispersion coefficient increases with increasing friction factor in the skimming flow regime (lower resistance) but is more uniform in the regimes of wake interference and isolated roughness (higher resistance). Hence, it is proposed that the friction factor and flow regimes could be adopted concurrently for pollutant concentration estimate in the UBL over urban street canyons of different roughness.     

Future trend of extreme value distributions of wintertime surface air temperatures over Korea and the associated physical changes

Daily winter temperatures in Korea have been analyzed via CSEOF analysis. Then, each PC time series was detrended and was fitted to an AR (autoregressive) model. Based on the identified AR model, an artificial time series of arbitrary length can be generated by using an arbitrary white-noise time series. In this way, one hundred new sets of PC time series were generated over the period of 1973–2058. Then, the trend for each PC time series was added back to the artificial PC time series extending the trend until 2058. Ultimately, artificial daily winter temperatures in Korea have been constructed by using the artificial PC time series and the original loading vectors derived from the observational data. The 100 new data sets have been investigated in order to understand the winter temperature variability 50 years into the future. Regression analysis in CSEOF space shows that temperature increase in Korea is associated with increased 850-hPa air temperature over most of the Asian domain (97°-153°E × 22°-73°N) and increased 850-hPa geopotential height in the southern part of the domain. As a result, southerly and southeasterly wind anomalies develop carrying positive temperature anomalies northward and northwestward. Both the 200-hPa air temperature and geopotential height changes indicate that there will be fairly significant northward shift of the jet stream in future. The standard deviation of the 200-hPa potential vorticity increases implying that shortwave trough and henceforth baroclinic instability will increase in future. Finally, GEV (Generalized Extreme Value) distribution and GPD (Generalized Pareto distribution) distribution have been compared between the observational records and the future records of the same length. The extreme value distributions based on the synthetic datasets show that warm extreme events will be more extreme in future and cold extreme events, on the other hand, will be less extreme. This study provides an estimate of future temperatures based on the observational data and serves as an independent baseline solution for comparisons with numerical model solutions.     

Changes in future precipitation over South Korea using a global high-resolution climate model

We used daily precipitation data from a global high-resolution climate scenario to analyze the features of future precipitation including extreme and heavy rainfall. The scenario shows that the model reproduces the daily precipitation over South Korea well. The projections show an increase in annual precipitation of approximately 18% in the late 21st century, with the highest increase (38%) occurring in winter. The number of days with daily precipitation of less than 5 mm decreases, but that of daily precipitation of more than 5 mm increases slightly in the latter part of the 21st century. The peak of precipitation days shifted from July to August. The number of days with relatively small amounts of precipitation (10 and 30 mm d^?1) increases most substantially in the winter season, but that for large amounts of precipitation (50, 80, 100, and 130 mm d^?1) increases most in the summer season. Events with heavy precipitation rates of 100 and 130 mm d^?1 are expected to occur in the winter season in the late 21st century, although no such events occurred during the winter season in the reference period.     

Recent changes of Mixed Layer Depth in the East/Japan Sea: 1994–2007

Recent changes of Mixed Layer Depth (MLD) in the East/Japan Sea (EJS) is investigated by using ocean reanalysis data. There is a significant shoaling of MLD in the southern EJS for 1994-2007, showing the change rate of - 1.24 m yr^?1. The present study demonstrates that the wind stress curl mostly attributes to the MLD changes. The long-term variation of wind stress curl with a positive trend in the southern EJS is generally well consistent with the MLD variation. In addition, upward displacement estimated from the wind stress curl reaches 14.6 m during 14 years, which is comparable to the shoaling of MLD (i.e., 17.4 m). Thus, the wind stress curl with positive trend would induce the increase of upward Ekman pumping which in turn lead to the shoaling of MLD. The change of temperature structure in the subsurface layer, that is shrinking and upward displacement of thermal layer from the reanalysis and observational data, also supports this conclusion.     

The effect of optimization and the nesting domain on carbon flux analyses in Asia using a carbon tracking system based on the ensemble Kalman filter

To estimate the surface carbon flux in Asia and investigate the effect of the nesting domain on carbon flux analyses in Asia, two experiments with different nesting domains were conducted using the CarbonTracker developed by the National Oceanic and Atmospheric Administration. CarbonTracker is an inverse modeling system that uses an ensemble Kalman filter (EnKF) to estimate surface carbon fluxes from surface CO₂ observations. One experiment was conducted with a nesting domain centered in Asia and the other with a nesting domain centered in North America. Both experiments analyzed the surface carbon fluxes in Asia from 2001 to 2006. The results showed that prior surface carbon fluxes were underestimated in Asia compared with the optimized fluxes. The optimized biosphere fluxes of the two experiments exhibited roughly similar spatial patterns but different magnitudes. Weekly cumulative optimized fluxes showed more diverse patterns than the prior fluxes, indicating that more detailed flux analyses were conducted during the optimization. The nesting domain in Asia produced a detailed estimate of the surface carbon fluxes in Asia and exhibited better agreement with the CO₂ observations. Finally, the simulated background atmospheric CO₂ concentrations in the experiment with the nesting domain in Asia were more consistent with the observed CO₂ concentrations than those in the experiment with the nesting domain in North America. The results of this study suggest that surface carbon fluxes in Asia can be estimated more accurately using an EnKF when the nesting domain is centered in Asian regions.     

Tectonometamorphic evolution of the Chuncheon amphibolite, central Gyeonggi massif, South Korea

Abstract The Chuncheon amphibolite, part of the Gubongsan Group which overlies the Yongduri gneiss complex, is interlayered with calc-silicate rock, marble, quartzite, biotite schist and quartzofeldspathic gneiss in the central Gyeonggi massif, South Korea. Metamorphic pressures and temperatures estimated from the amphibolite are 5.5–10.6 kbar and 615–714°C. These P—T conditions are close to those defined by the reaction curve between kyanite and sillimanite, and suggest medium-pressure-type metamorphism of the Chuncheon amphibolite. For two metapelites intercalated with the amphibolite, temperatures are estimated to be 607–699° C, consistent with those obtained from the amphibolite. On the other hand, pressures estimated from these metapelites are significantly different, 4–6 kbar and 9–13 kbar, when rim and core compositions of garnet are, respectively, used. These P—T estimates obtained from the amphibolite and metapelite suggest a nearly isothermal decompression of 3–7 kbar during denudation. Rapid decompression is likely on the basis of the results of mineral chemistry, phase equilibria and geothermobarometer. Moreover, in conjunction with the occurrence of kyanite in the adjacent Gyeonggi gneiss complex, P—T estimates of the Chuncheon amphibolite and metapelite suggest a clockwise P—T—t path. This evolutionary path may be related to the amalgamation of continents during the late Proterozoic event which corresponds to the Jinningian orogeny in the Qinling belt of China.     

Transition from the Zeolite to Prehnite-Pumpellyite Facies in the Karmutsen Metabasites, Vancouver Island, British Columbia

The upper Triassic Karmutsen metabasites from northeast VancouverIsland, B.C., are thermally metamorphosed by the intrusion ofthe Coast Range Batholith. The amygdaloidal metabasites developedin the outer portion of the contact aureole show a progressivemetamorphism from zeolite to prehnite-pumpellyite facies. Thesize of an equilibrium domain is extremely small for these metabasites,and the individual amygdule assemblages are assumed to be inequilibrium. Two major calcite-free assemblages (+chlorite+quartz)are characteristic: (i) laumontite+pumpellyite+epidote in thezeolite facies and (ii) prehnite+pumpellyite+epidote in theprehnite-pumpellyite facies. The assemblages and compositionsof Ca-Al silicates are chemographically and theoretically interpretedon the basis of the predicted P-T grid for the model basalticsystem, CaO-MgO-A1₂O₃-Fe₂O₃-SiO₂-H₂O. The results indicate:(1) local equilibrium has been approached in mineral assemblagesand compositions; (2) the X_Fe³⁺ values in the coexisting Ca-Alsilicates decrease from epidote, through pumpellyite to prehnite;(3) with increasing metamorphic grade, the Fe³⁺ contents ofepidotes in reaction assemblages decrease in the zeolite facies,then increase in the prehnite-pumpellyite facies rocks. Suchvariations in the assemblages and mineral compositions are controlledby a sequence of continuous and discontinuous reactions, andallow delineation of T-X_Fe³⁺ relations at constant pressure.The transition from the zeolite to prehnite-pumpellyite faciesof the Karmutsen metabasites is defined by a discontinuous reaction:0·18 laumontite+pumpellyite+0·15 quartz = 1·31prehnite+ 0·78 epidote+0·2 chlorite+ 1·72H₂O, where the X_Fe³⁺ values of prehnite, pumpellyite and epidoteare 0·03, 0·10 and 0·18, respectively.These values together with available thermodynamic data andour preliminary experimental data are used to calculate theP-T condition for the discontinuous reaction as P = 1·1±0·5 kb and T = 190±30°C. The effectsof pressure on the upper stability of the zeolite facies assemblagesare discussed utilizing T-X_Fe³⁺ diagrams. The stability of thelaumontite-bearing assemblages for the zeolite facies metamorphismof basaltic rocks may be defined by either continuous or discontinuousreactions depending on the imposed metamorphic field gradient.Hence, the zeolite and prehnite-pumpellyite facies transitionboundary is multivariant.     

A wind tunnel study of flows over idealised urban surfaces with roughness sublayer corrections

Variation of soil moisture during active and weak phases of summer monsoon JJAS (June, July, August, and September) is very important for sustenance of the crop and subsequent crop yield. As in situ observations of soil moisture are few or not available, researchers use data derived from remote sensing satellites or global reanalysis. This study documents the intercomparison of soil moisture from remotely sensed and reanalyses during dry spells within monsoon seasons in central India and central Myanmar. Soil moisture data from the European Space Agency (ESA)—Climate Change Initiative (CCI) has been treated as observed data and was compared against soil moisture data from the ECMWF reanalysis-Interim (ERA-I) and the climate forecast system reanalysis (CFSR) for the period of 2002–2011. The ESA soil moisture correlates rather well with observed gridded rainfall. The ESA data indicates that soil moisture increases over India from west to east and from north to south during monsoon season. The ERA-I overestimates the soil moisture over India, while the CFSR soil moisture agrees well with the remotely sensed observation (ESA). Over Myanmar, both the reanalysis overestimate soil moisture values and the ERA-I soil moisture does not show much variability from year to year. Day-to-day variations of soil moisture in central India and central Myanmar during weak monsoon conditions indicate that, because of the rainfall deficiency, the observed (ESA) and the CFSR soil moisture values are reduced up to 0.1 m³/m³ compared to climatological values of more than 0.35 m³/m³. This reduction is not seen in the ERA-I data. Therefore, soil moisture from the CFSR is closer to the ESA observed soil moisture than that from the ERA-I during weak phases of monsoon in the study region.

     

The diel vertical migration of the sound-scattering layer in the Yellow Sea Bottom Cold Water of the southeastern Yellow Sea：focus on its relationship with a temperature structure

Using the hydroacoustic method with a 200 kHz scientific echo sounding system,the diel vertical migration(DVM) of the sound-scattering layer(SSL) in the Yellow Sea Bottom Cold Water(YSBCW) of the southeastern Yellow Sea was studied in April(spring) and August(summer) of 2010 and 2011.For each survey,13–27 hours of acoustic data were continuously collected at a stationary station.The acoustic volume scattering strength(Sv) data were analyzed with temperature profile data.In the spring of both 2010 and 2011,the SSL clearly showed the vertical migration throughout the entire water column,moving from the surface layer at night to near the bottom during the day.Conductivity,temperature,and depth data indicated that the entire water column was well mixed with low temperature of about 8 C.However,the SSL showed different patterns in the summers of 2010 and 2011.In the summer of 2010(28 C at the surface),the SSL migrated to near the bottom during the day,but there were two SSLs above and below the thermocline at depth of 10–30 m at night.In the summer of 2011(20 C at the surface),the SSL extended throughout the entire water column at night,possibly owing to an abrupt change in sea weather conditions caused by the passage of a Typhoon Muifa over the study area.It was concluded that the DVM patterns in summer in the YSBCW area may be greatly influenced by a strengthened or weakened thermocline.     

Statistical multisite simulations of summertime precipitation over South Korea and its future change based on observational data

In this study, a weather generator for summer (May 19 – September 15) precipitation over South Korea is developed. Precipitation data for 33 years (1979–2011) observed at 57 stations of Korea Meteorological Administration (KMA) are used to develop a new weather generator. Using the cyclostationary empirical orthogonal function (CSEOF) technique, the observed precipitation data is described as a linear combination of deterministic evolution patterns and corresponding stochastic amplitude (principal component) time series. An autoregressive-moving average (ARMA) model is used to generate one hundred sets of synthetic amplitude time series for the period of 1979–2061 (83 years) with similar statistical properties of the original amplitude time series. Based on these synthetic time series and the annually repeating evolution patterns, one hundred sets of synthetic summer precipitation were generated. Statistical characteristics of the synthetic datasets are examined in comparison with those of the KMA observational record for the period of the observational record. Characteristic changes of synthetic precipitations for a future period are also examined. The seasonal cycle in the synthetic precipitation is reproduced faithfully with typical bimodal peaks of summer precipitation. The spatial correlation patterns of the synthetic precipitation are fairly similar to that of the observational data. The frequency-intensity relationship of the synthetic precipitation also looks similar to that of the observational data. In the future period, precipitation amount increases except in the precipitation range of (0,10) mm day^?1 with nearly no change in the frequency of no-rain days; frequency increase is particularly conspicuous in the range of (100,500) mm day^?1.