Dynamic p-enrichment schemes for multicomponent reactive flows

We present a family of p-enrichment schemes. These schemes may be separated into two basic classes: the first, called fixed tolerance schemes, rely on setting global scalar tolerances on the local regularity of the solution, and the second, called dioristic schemes, rely on time-evolving bounds on the local variation in the solution. Each class of p-enrichment scheme is further divided into two basic types. The first type (the Type I schemes) enrich along lines of maximal variation, striving to enhance stable solutions in “areas of highest interest.” The second type (the Type II schemes) enrich along lines of maximal regularity in order to maximize the stability of the enrichment process. Each of these schemes are tested on three model systems. The first is an academic exact system where basic analysis is easily performed. Then we discuss a pair of application model problems arising in coastal hydrology. The first being a contaminant transport model, which addresses a declinature problem for a contaminant plume with respect to a bay inlet setting. And the second, a multicomponent chemically reactive flow model of estuary eutrophication arising in the Gulf of Mexico.     

The properties of convective storms in central Mexico: A radar and lightning approach

Radar data from Cerro Catedral (a peak close to Mexico City) were used to investigate the properties of convective storms over central Mexico, a region with complex orography. The spatial distribution shows a preference for storms to form and move to the west of radar, over a narrow band of high terrain. However, the storms with the higher volumes and echo-top heights tend to be located southwestward over lower terrain. Each radar feature was matched with the number of cloud-to-ground (CG) lightning produced inside it, as retrieved from the World Wide Lightning Location Network dataset. The storms in which lightning was detected, with an average of more than six lightning bolts, clearly outperform in size and intensity the group of storms in which lightning was not detected, and tend to lie over lower terrain. The sample of over 98 000 identified cells was divided into four elevation groups to look for elevation trends in the mean properties, as reported for other Mexican regions. While the number of storms per unit area increases with terrain height, the average values for properties related to both size (area, volume, echo-top height) and intensity (maximum reflectivity, number of CG bolts, height of maximum reflectivity, maximum height of 30 dBZ echo) decrease. These results could be related to the possible shallower warm-cloud depths over the higher elevations. The diurnal cycles of convection and lightning north of the radar show a nearly typical continental regime of precipitation in that zone, with maxima at 18:00 LT in both variables. However, south of the radar, a maximum in lightning activity occurs during late night and early morning, which is linked with the deeper nocturnal convection over the lower terrain in that zone.     

ENSEMBLE PREDICTION EXPERIMENTS OF TRACKS OF TROPICAL CYCLONES BY USING MULTIPLE CUMULUS PARAMETERIZATION SCHEMES

Ensemble prediction experiments of the tracks of eight tropical cyclones occurring between 2004 -2006 over the western Pacific have been performed by using MM5 with five cumulus parameterizations chemes. The results show that the predictions of the tracks of the tropical cyclones are sensitive to the selection of cumulus parameterization schemes. Each scheme has its own advantage and disadvantage, and the predications without cumulus parameterization schemes are not the worst, sometimes even better than the others. And all of the three ensemble methods improve the predictions of the tracks significantly, among which the ensemble method without parameterization schemes, the Grell, Betts-Miller and Kain-Fritsch schemes are the best.     

A Simulating Study on Resolvable-Scale Microphysical Parameterization in a Mesoscale Model

The Penn State/ NCAR Mesoscale Model (MM5) is used to simulate the precipitation event that oc-curred during 1-2 May 1994 to the south of the Yangtze River. In five experiments the Kain-Fritsch scheme is made use of for the subgrid-scale convective precipitation, but five different resolvable-scale microphysical parameterization schemes are employed. They are the simple super-saturation removal scheme, the warm rain scheme of Hsie et al. (1984), the simple ice scheme of Dudhia (1989), the complex mixed-phase scheme developed by Reisner et al. (1993). and the GSFC microphysical scheme with graupel. Our interest is how the various resolvable-scale schemes affect the domain-averaged precipitation, the pre-cipitation distribution, the sea level pressure, the cloud water and the cloud ice.Through a series of experiments about a warm sector rainfall case, results show that although the dif-ferent resolvable-scale scheme is used, the differences of the precipitation characteristics among all five runs are not very obvious. However, the precipitation is over-predicted and the strong mesoscale low is produced by the simple super-saturation removal scheme. The warm rain scheme with the inclusion of condensation and evaporation under-predicts the precipitation and allows the cloud water to reach the 300 hPa level The scheme of the addition of graupel increases the resolvable-scale precipitation by about 20%–30%. The inclusion of supercooled liquid water in the grid-scale scheme does not affect significantly the results.     

Impacts of Two Ice Parameterization Schemes on the Cloud Microphysical Processes and Precipitation of a Severe Storm in Northern China

A severe storm that occurred over Beijing in northern China on 23 June 2011 was simulated with two different ice crystal parameterization schemes(the DeMott scheme and Meyers scheme) by using the Regional Atmospheric Modeling System. Compared with the DeMott scheme, the simulation results with the Meyers scheme have the following characteristics:(1) Updrafts are stronger and more numerous;(2) The cloud is better organized and contains a greater peak of ice-phase hydrometeor mixing ratios;(3) Cloud water and hail mixing ratios increase while graupel mixing ratios decrease;(4) The surface precipitation is initially greater. However, at the end of the simulation, less precipitation is produced. In short, the differences between the two schemes are not obvious, but the De Mott scheme has a relatively more reasonable result.     

总体最小二乘方法在三维坐标转换中的应用

工程测量中经常需要实现不同坐标系下成果的相互转换,而高精度的转换参数是完成这一工作的基础。获取基准转换参数的实质就是利用公共点在两套坐标系下的坐标,根据相似变换原理建立误差方程求解。传统的最小二乘（LS）相似变换法只考虑了公共点在一套坐标系下的误差,与实际情况不符。基于此,探讨了坐标参数化的平面基准转换方法,解决了考虑公共点在两套坐标系下坐标都含有误差时高斯-马尔科夫（Gauss-Markov,G-M）模型不成立的问题,以相似变换原理为基础,采取通用的最小二乘方法解算基准转换参数。     

Cases-97: Late-Morning Warming And Moistening Of The Convective Boundary Layer Over The Walnut River Watershed

Aircraft, radiosonde, surface-flux, and boundary-layer windprofiler data from the Cooperative Atmosphere Surface Exchange Study's 1997 field project, CASES-97, are combined with synoptic data to study the evolution of the vertically-averaged mixed-layerpotential temperature []and mixing-ratio [Q] onthree nearly-cloudless days from 1000 CST to 1200CST (local noon is approximately 1230 CST). This was achieved through examination of the terms in the time-tendency (`budget')equations for []and [Q]. We estimate three of the terms –local time rate of change, vertical flux divergence, andhorizontal advection. For the [Q]-budget, vertical flux divergence usually dominates, buthorizontal advection is significant on one of the three days. The [Q]-budget balances for two of the three days to within the large experimental error. For the -budget,vertical flux divergence accounts for most of the morningwarming, with horizontal advection of secondary importance.The residual in the -budget has the same sign for all three days, indicating that not all the heating is accounted for. We can balance the []-budgets to within experimental error on two of the three days by correcting the vertical-flux divergence for apparent low biases in the flux measurements of one of the aircraft and in the surface fluxes, and accounting for direct heating of the mixed layer by radiative flux divergence allowing for the effects of carbonaceous aerosols. The [];-budget with these corrections also balances on the third day if horizontal gradients from synoptic maps are used to estimate the horizontal advection. However, the corrected budget for this day does not balance if the horizontal gradient in the advection term is estimated using CASES-97aircraft and radiosondes; we suggest that persistent mesoscale circulations led to an overestimate of the horizontal gradient andhence horizontal advection.     

Performance of shape indices and classification schemes for characterising perceptual shape complexity of building footprints in GIS

Shape characterisation is important in many fields dealing with spatial data. For this purpose, numerous shape analysis and recognition methods with different degrees of complexity have so far been developed. Among them, relatively simple indices are widely used in spatial applications, but their performance has not been investigated sufficiently, particularly for building footprints (BFs). Therefore, this article focuses on BF shape characterisation with shape indices and classification schemes in a GIS environment. This study consists of four phases. In the first phase, the criteria for BF shape complexity were identified, and accordingly, benchmark data was constructed by human experts in three shape complexity categories. In the second phase, 18 shape indices were selected from the literature and automatically computed in GIS. The performance of these indices was then statistically assessed with histograms, correlation matrix and boxplots, and consequently four indices were found to be appropriate for further investigation. In the third phase, two new indices (Equivalent Rectangular index and Roughness index) were proposed with the objective to measure some BF shape characteristics more efficiently. The proposed indices also were found to be appropriate with the same statistical assessment procedures. In the final phase, BF shape complexity categories were created with the pairs of six appropriate indices and four choropleth mapping classification schemes (equal intervals, natural break, standard deviation, and custom) in GIS. The performance of the index–scheme pairs was assessed against the benchmark data. The findings demonstrated that both new indices and two of the selected indices (Convexity and Rectangularity) delivered higher performance. The custom classification scheme was found more ideal to reveal absolute shape complexity with the index value ranges derived from the boxplots while the other classification schemes were more appropriate to reveal relative shape complexity.     

Application of Dynamic Subgrid-scale Models for Large-eddy Simulation of the Daytime Convective Boundary Layer over Heterogeneous Surfaces

The sensitivity of large-eddy simulation (LES) to the representation of subgrid-scale (SGS) processes is explored for the case of the convective boundary layer (CBL) developing over surfaces with varying degrees of spatial heterogeneity. Three representations of SGS processes are explored: the traditional constant Smagorinsky–Lilly model and two other dynamic models with Lagrangian averaging approaches to calculate the Smagorinsky coefficient (C _S ) and SGS Prandtl number (Pr). With initial data based roughly on the observed meteorology, simulations of daytime CBL growth are performed over surfaces with characteristics (i.e. fluxes and roughness) ranging from homogeneous, to striped heterogeneity, to a realistic representation of heterogeneity as derived from a recent field study. In both idealized tests and the realistic case, SGS sensitivities are mostly manifest near the surface and entrainment zone. However, unlike simulations over complex domains or under neutral or stable conditions, these differences for the CBL simulation, where large eddies dominate, are not significant enough to distinguish the performance of the different SGS models, irrespective of surface heterogeneity.     

A sensitivity of squall-line intensity to environmental static stability under various shear and moisture conditions

Squall lines develop in various climate regions having diverse environmental profiles of wind shear, moisture, and temperature. In order to explore the sensitivity of squall lines to these environmental profiles, we have performed an extensive set of numerical simulations under various shear and moisture conditions in midlatitude-continental and tropical–oceanic temperature environments. From the results of the sensitivity simulations and the analyses of the environmental parameters, it is found that the static stability in a convectively unstable layer is of primary importance in determining the strength of squall lines. Under temperature environments having the same static stability, convective available potential energy (CAPE) and precipitable water content (PWC) well describe the squall-line intensity. Vertical shear also plays an important role in determining the squall-line structure as well as the intensity through the interaction with surface cold pool. The combination of the static stability in a convectively unstable layer, CAPE, and PWC should be examined in diagnosing the intensity of squall lines that develop with an optimal shear for their environment.