Faculty and Staff

This paper studies the effect of targeted observations on state and parameter estimates determined with Kalman filter data assimilation (DA) techniques. We first provide an analytical result demonstrating that targeting observations within the Kalman filter for a linear model can significantly reduce state estimation error as opposed to fixed or randomly located observations. We next conduct observing system simulation experiments for a chaotic model of meteorological interest, where we demonstrate that the local ensemble transform Kalman filter (LETKF) with targeted observations based on largest ensemble variance is skillful in providing more accurate state estimates than the LETKF with randomly located observations. Additionally, we find that a hybrid ensemble Kalman filter parameter estimation method accurately updates model parameters within the targeted observation context to further improve state estimation.

The effects of urbanization on ozone levels have been widely investigated over cities primarily located in temperate and/or humid regions. In this study, nested WRF-Chem simulations with a finest grid resolution of 1 km are conducted to investigate ozone concentrations O₃ due to urbanization within cities in arid/semi-arid environments. First, a method based on a shape preserving Monotonic Cubic Interpolation (MCI) is developed and used to downscale anthropogenic emissions from the 4 km resolution 2005 National Emissions Inventory (NEI05) to the finest model resolution of 1 km. Using the rapidly expanding Phoenix metropolitan region as the area of focus, we demonstrate the proposed MCI method achieves ozone simulation results with appreciably improved correspondence to observations relative to the default interpolation method of the WRF-Chem system. Next, two additional sets of experiments are conducted, with the recommended MCI approach, to examine impacts of urbanization on ozone production: (1) the urban land cover is included (i.e., urbanization experiments) and, (2) the urban land cover is replaced with the region's native shrubland. Impacts due to the presence of the built environment on O₃ are highly heterogeneous across the metropolitan area. Increased near surface O₃ due to urbanization of 10–20 ppb is predominantly a nighttime phenomenon while simulated impacts during daytime are negligible. Urbanization narrows the daily O₃ range (by virtue of increasing nighttime minima), an impact largely due to the region's urban heat island. Our results demonstrate the importance of the MCI method for accurate representation of the diurnal profile of ozone, and highlight its utility for high-resolution air quality simulations for urban areas.

An interesting occurrence of a Rossby wave breaking event observed during the VORCORE experiment is presented and explained. Twenty-seven balloons were launched inside the Antarctic polar vortex. Almost all of these balloons evolved in the stratosphere around 500K within the vortex, except the one launched on 28 October 2005. In this case, the balloon was caught within a tongue of high potential vorticity (PV), and was ejected from the polar vortex. The evolution of this event is studied for the period between 19 and 25 November 2005. It is found that at the beginning of this period, the polar vortex experienced distortions due to the presence of Rossby waves. Then, these waves break and a tongue of high PV develops. On 25 November, the tongue became separated from the vortex and the balloon was ejected into the surf zone. Lagrangian simulations demonstrate that the air masses surrounding the balloon after its ejection were originating from the vortex edge. The wave breaking and the development of the tongue are confined within a region where a planetary Quasi-Stationary Wave 1 (QSW1) induces wind speeds with weaker values. The QSW1 causes asymmetry in the wind speed and the horizontal PV gradient along the edge of the polar vortex, resulting in a localized jet. Rossby waves with smaller scales propagating on top of this jet amplify as they enter the jet exit region and then break. The role of the QSW1 on the formation of the weak flow conditions that caused the non-linear wave breaking observed near the vortex edge is confirmed by three-dimensional numerical simulations using forcing with and without the contribution of the QSW1.

Physical mechanisms of incongruency between observations and Weather Research and Forecasting (WRF) Model predictions are examined. Limitations of evaluation are constrained by (i) parameterizations of model physics, (ii) parameterizations of input data, (iii) model resolution, and (iv) flux observation resolution. Observations from a new 22.1-m flux tower situated within a residential neighborhood in Phoenix, Arizona, are utilized to evaluate the ability of the urbanized WRF to resolve finescale surface energy balance (SEB) when using the urban classes derived from the 30-m-resolution National Land Cover Database. Modeled SEB response to a large seasonal variation of net radiation forcing was tested during synoptically quiescent periods of high pressure in winter 2011 and premonsoon summer 2012. Results are presented from simulations employing five nested domains down to 333-m horizontal resolution. A comparative analysis of model cases testing parameterization of physical processes was done using four configurations of urban parameterization for the bulk urban scheme versus three representations with the Urban Canopy Model (UCM) scheme, and also for two types of planetary boundary layer parameterization: the local Mellor–Yamada–Janjić scheme and the nonlocal Yonsei University scheme. Diurnal variation in SEB constituent fluxes is examined in relation to surface-layer stability and modeled diagnostic variables. Improvement is found when adapting UCM for Phoenix with reduced errors in the SEB components. Finer model resolution is seen to have insignificant (<1 standard deviation) influence on mean absolute percent difference of 30-min diurnal mean SEB terms.

Forecasts of noise pollution from a highway line segment noise source are obtained from a sound propagation model utilizing effective sound speed profiles derived from a Numerical Weather Prediction (NWP) limited area forecast with 1 km horizontal resolution and near-ground vertical resolution finer than 20 m. Methods for temporal along with horizontal and vertical spatial nesting are demonstrated within the NWP model for maintaining forecast feasibility. It is shown that vertical nesting can improve the prediction of finer structures in near-ground temperature and velocity profiles, such as morning temperature inversions and low level jet-like features. Accurate representation of these features is shown to be important for modeling sound refraction phenomena and for enabling accurate noise assessment. Comparisons are made using the parabolic equation model for predictions with profiles derived from NWP simulations and from field experiment observations during mornings on November 7 and 8, 2006 in Phoenix, Arizona. The challenges faced in simulating accurate meteorological profiles at high resolution for sound propagation applications are highlighted and areas for possible improvement are discussed.

Explicit solutions of the inhomogeneous paraxial wave equation in a linear and quadratic approximation are applied to wave fields with invariant features, such as oscillating laser beams in a parabolic waveguide and spiral light beams in varying media. A similar effect of superfocusing of particle beams in a thin monocrystal film, harmonic oscillations of cold trapped atoms, and motion in magnetic field are also mentioned.

In this paper, we study oscillating solutions of the 1D-quintic nonlinear Schrödinger equation with the help of Wigner's quasiprobability distribution in quantum phase space. An "absolute squeezing property," namely a periodic in time total localization of wave packets at some finite spatial points without violation of the Heisenberg uncertainty principle, is analyzed in this nonlinear model.

Oppia Quadrangle Av-10 (288–360°E, ±22°) is a junction of key geologic features that preserve a rough history of Asteroid (4) Vesta and serves as a case study of using geologic mapping to define a relative geologic timescale. Clear filter images, stereo-derived topography, slope maps, and multispectral color-ratio images from the Framing Camera on NASA’s Dawn spacecraft served as basemaps to create a geologic map and investigate the spatial and temporal relationships of the local stratigraphy. Geologic mapping reveals the oldest map unit within Av-10 is the cratered highlands terrain which possibly represents original crustal material on Vesta that was then excavated by one or more impacts to form the basin Feralia Planitia. Saturnalia Fossae and Divalia Fossae ridge and trough terrains intersect the wall of Feralia Planitia indicating that this impact basin is older than both the Veneneia and Rheasilvia impact structures, representing Pre-Veneneian crustal material. Two of the youngest geologic features in Av-10 are Lepida (∼45 km diameter) and Oppia (∼40 km diameter) impact craters that formed on the northern and southern wall of Feralia Planitia and each cross-cuts a trough terrain. The ejecta blanket of Oppia is mapped as ‘dark mantle’ material because it appears dark orange in the Framing Camera ‘Clementine-type’ color-ratio image and has a diffuse, gradational contact distributed to the south across the rim of Rheasilvia. Mapping of surface material that appears light orange in color in the Framing Camera ‘Clementine-type’ color-ratio image as ‘light mantle material’ supports previous interpretations of an impact ejecta origin. Some light mantle deposits are easily traced to nearby source craters, but other deposits may represent distal ejecta deposits (emplaced >5 crater radii away) in a microgravity environment.

We produced two 1:250,000 scale geologic maps of the adjacent quadrangles Av-6 Gegania and Av-7 Lucaria, located in the equatorial region of (4) Vesta (0–144°E, 22°S to 22°N). The mapping is based on clear and color filter images of the Framing Camera (FC) onboard the Dawn spacecraft, which has captured the entire illuminated surface of Vesta with high spatial resolution (up to ∼20 m/pixel), and on a digital terrain model derived from FC imagery. Besides the geologic mapping itself, a secondary purpose of this work is to investigate one of the most prominent morphological features on Vesta, namely the aggregation of several giant equatorial troughs termed the Divalia Fossae, most probably formed during the Rheasilvia impact near Vesta’s south pole. The up to 465 km long and 22 km wide troughs show height differences of up to 5 km between adjacent troughs and ridges. Another imprint of the Rheasilvia impact is the >350 km long and ∼250 km wide swath of ejecta crossing quadrangle Av-6 Gegania. This lobe shows a distinct appearance in FC color ratios and a high albedo in FC images, indicating a mineralogical similarity to material typically found within the Rheasilvia basin, in particular composed of diogenite-rich howardites. Almost the entire northern half of the mapping area shows the oldest surface, being dominated by upper crustal basaltic material. To the south, increasingly younger formations related to the Rheasilvia impact occur, either indicated by the troughs formed by Rheasilvia or by the Rheasilvia ejecta itself. Only medium sized impact craters with diameters less than 22 km occur within the two mapped quadrangles. Some of the craters exhibit ejecta blankets and/or distinctly dark or bright ejecta material in ejecta rays outside and exposures within the crater, and mass-wasting deposits down crater slopes, forming the youngest surfaces.

We used Dawn spacecraft data to identify and delineate geological units and landforms in the Marcia quadrangle of Vesta as a means to assess the role of the large, relatively young impact craters Marcia (∼63 km diam.) and Calpurnia (∼53 km diam.) and their surrounding ejecta field on the local geology. We also investigated a local topographic high with a dark-rayed crater named Aricia Tholus, and the impact crater Octavia that is surrounded by a distinctive diffuse mantle. Crater counts and stratigraphic relations suggest that Marcia is the youngest large crater on Vesta, in which a putative impact melt on the crater floor ranges in age between ∼40 and 60 Ma (depending upon choice of chronology system), and Marcia’s ejecta blanket ranges in age between ∼120 and 390 Ma (depending upon choice of chronology system).

We interpret the geologic units in and around Marcia crater to mark a major vestan time-stratigraphic event, and that the Marcia Formation is one of the geologically youngest formations on Vesta. Marcia crater reveals pristine bright and dark material in its walls and smooth and pitted terrains on its floor. The smooth unit we interpret as evidence of flow of impact melts and (for the pitted terrain) release of volatiles during or after the impact process. The distinctive dark ejecta surrounding craters Marcia and Calpurnia is enriched in OH- or H-bearing phases and has a variable morphology, suggestive of a complex mixture of impact ejecta and impact melts including dark materials possibly derived from carbonaceous chondrite-rich material. Aricia Tholus, which was originally interpreted as a putative vestan volcanic edifice based on lower resolution observations, appears to be a fragment of an ancient impact basin rim topped by a dark-rayed impact crater. Octavia crater has a cratering model formation age of ∼280–990 Ma based on counts of its ejecta field (depending upon choice of chronology system), and its ejecta field is the second oldest unit in this quadrangle. The relatively young craters and their related ejecta materials in this quadrangle are in stark contrast to the surrounding heavily cratered units that are related to the billion years old or older Rheasilvia and Veneneia impact basins and Vesta’s ancient crust preserved on Vestalia Terra.