Air parcel motions during a simulated sudden stratospheric warming
Explore the source record for details and available documents.
SEARCH · Search NASA
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Explore the source record for details and available documents.
Midstratospheric trajectories for February and March 1979 are calculated using geopotential analyses derived from limb infrared monitor of the stratosphere data. These trajectories are compared with the corresponding results using stratospheric sounding unit data. The trajectories are quasi-isentropic in that a radiation scheme is used to simply cross-isentrope flow. The results show that in disturbed conditions, quantitative agreement the trajectories, that is, within 25 great circle degrees (GCD) (one GCD about 110 km) may be valid for only 3 or 4 days, whereas during quiescent periods, quantitative agreement may last up to 10 days. By comparing trajectories calculated with different data some insight can be gained as to errors due to vertical resolution and horizontal resolution (due to infrequent sampling) in the analyzed geopotential height fields. For the disturbed trajectories described in this paper the horizontal resolution of the data was more important than vertical resolution; however, for the quiescent trajectories, which could be calculated accurately for a longer duration because of the absence of appreciable transients, the vertical resolution of the data was found to be more important than the horizontal resolution. It is speculated that these characteristics are also applicable to trajectories calculated during disturbed and quiescent periods in general. A review of some recently published trajectories shows that the qualitative conclusions of such works remains unaffected when the calculations are repeated using different data.
An overview of Airborne Arctic Stratospheric Expedition 2 is given. Effects of Pinatubo aerosol on stratospheric ozone at mid-latitudes, in situ measurements of ClO and ClO/HCl ratio, balloon-borne measurements of ClO, NO, and O3 in a volcanic cloud, and new observations of the NO(y)/N2O correlation in the lower stratosphere are discussed. Among other topics addressed are the following: in situ tracer correlations of methane, nitrous oxide, and ozone as observed aboard the DC-8, in situ measurements of changes in stratospheric aerosol and the N2O-aerosol relationship inside and outside of the polar vortex, measurements of halogenated organic compounds near the tropical tropopause, and airborne brightness measurements of the polar winter troposphere.
During January and February 1987, as part of the Stratosphere-Troposphere Exchange Project, the NASA ER-2 made 11 flights from Darwin, Australia to investigate dehydration mechanisms in the vicinity of the tropical tropopause. After the monsoon onset in the second week of January, steady easterly flow of 15-25 ms (exp -1) was established in the upper troposphere and lower stratosphere over northern Australia and adjacent seas. Penetrating into this regime were elements of the monsoon convection such as overshooting convective turrets and extensive anvils including cyclone cloud shields. In cases of the latter, the resulting flow obstructions tended to produce mesoscale gravity waves. In several instances the ER- 2 meteorological and trace constituent measurements provide a detailed description of the structure of these gravity waves. Among these was STEP Flight 6, 22-23 January. It is of particular interest to STEP because of the close proximity of ice-laden and dehydrated air on the same isentropic surfaces. Convective events inject large amounts of ice into the upper troposphere and lower stratosphere which may not be completely removed by local precipitation processes. In the present instance, a gravity wave for removed from the source region appears to induce relativity rapid upward motion in the ice-laden air and subsequent dessication. Potential mechanisms for such a localized removal process are under investigation.
Explore the source record for details and available documents.
Initial droplet spectra produced upon activation impact the ensuing chain of microphysical processes andtherefore play a crucial role in cloud evolution. This work re-examines dependencies of newly formed clouddroplet size distribution (CDSD) characteristics on environmental and aerosol properties via parcel model simulationsthat serve as the basis for a multi-moment bulk microphysics droplet activation scheme suitable for acloud-resolving model (CRM). It is found that applying a fixed size threshold to define activated droplets versusemploying physical considerations can lead to erroneous activation and overly broad CDSDs for high aerosolconcentration and weak updraft conditions. Aerosol distributions characterized by larger median sizes and/orincreased solubility can result in greater activated droplet numbers, whereas impacts of these parameters onCDSD spectral width depend on both aerosol number concentration and updraft velocity. An expansion of theactivation scheme to include CDSD spectral width is proposed to aid efforts to extend high-order momentprediction to cloud droplet categories in CRMs as well as better represent variability in the activation process onthe cloud scale.simulations to investigate the regime dependence of the relative dispersion(d)1 of newly activated CDSDs, where d is the ratio of dropletradius standard deviation (σ) to the mean radius (r ). C16 demonstratedthat increasing Na resulted in increasing (decreasing) d values via reducedcondensational narrowing (spectral broadening) rates within theAL (UL) regime, with d values peaking in the TR regime. Their findingssuggest a similar regime dependence for d as R09 noted for Nc and helpexplain reportedly conflicting relationships between Na and CDSDspectral characteristics (cf. Hudson and Noble, 2014; Liu et al., 2014),although the applicability of these results within bulk microphysicalschemes was not addressed.Simulating aerosol-cloud interactions with CRMs employing bulkmicrophysics requires that the scheme minimally predict two CDSDparameters, namely mass and number concentrations, and represent thedroplet activation process. Various activation schemes aim to determineNc from aerosol and environmental properties and include analyticalexpressions (e.g., Abdul-Razzak et al., 1998; Morrison et al., 2005) aswell as lookup tables (LUTs) based on detailed parcel model calculations(e.g., Saleeby and Cotton, 2004, hereafter SC04; Segal and Khain,2006; Thompson and Eidhammer, 2014). Expressions to diagnose CDSDspectral width from Nc (Grabowski, 1998; Liu et al., 2006; Morrison andGrabowski, 2007) or cloud water content (Geoffroy et al., 2010) havealso been developed, although more robust methods to obtain CDSDspectral width upon activation are presently lacking. This latter point isrelevant for triple-moment (3 M) bulk microphysics that aim to predictdistribution spectral width alongside number and mass concentrations(e.g., Loftus et al., 2014; Milbrandt and Yau, 2005).The current work extends the findings of C16 to the current LUTbasedaerosol activation scheme used in the Regional AtmosphericModeling System (RAMS) (Cotton et al., 2003; SC04; Saleeby and vanden Heever, 2013, hereafter SvdH13) and additionally examinesaerosol size and solubility impacts on newly activated CDSD properties.Because early cloud development processes such as condensationalgrowth, evaporation, and droplet self-collection depend on and impactCDSD spectral width (Hudson and Yum, 1997; Seifert and Beheng 2001;Lu and Seinfeld, 2006; Igel and van den Heever, 2017), an expansion ofthe activation LUTs to include CDSD spectral width is proposed as apreliminary step for extending 3M prediction to CDSDs in CRMs forimproved simulations of aerosol-cloud interactions.2. MethodologyThe current RAMS two-moment microphysics module determinesthe fractional number of aerosol particles that activate to cloud dropletsfrom five-dimensional LUTs based on model predicted air temperature(T), w, Na, and the geometric median radius (rg) and soluble fraction (ε)of the aerosol size distribution (SvdH13). These LUTs are created offlineusing a one-dimensional Lagrangian adiabatic parcel model (Feingoldand Heymsfield, 1992; Heymsfield and Sabin, 1989; SC04) to simulateexplicit droplet activation and initial CDSD growth for a range of ambientatmospheric conditions [T, w] and binned lognormal aerosol sizedistributions given by= ⎡⎣ ⎢− ⎤⎦ ⎥N r Nr π σr rσ( )2 lnexp[ln( / )]2(ln )aggg22 (1)where r is the dry aerosol particle bin radius and σg is the geometricstandard deviation of the distribution. As the parcel model simulationsfocus on the activation process, other processes such as coalescence,sedimentation, and mixing are not considered. Details of the parcelmodel can be found in SC04 and SvdH13, and only a brief description isprovided here. At the onset of parcel model calculations, the initiallydry aerosol particles in all bins first deliquesce and reach theirequilibrium diameters in a sub-saturated environment based on theKöhler equation for solution droplets. The parcel is then lifted at a fixedupward velocity w, and particle growth by vapor diffusion, along withconcurrent changes in the ambient environment, are iteratively computedusing the Variable-coefficient Ordinary Differential Equation(VODE) solver (Brown et al., 1989). The time resolution of these calculationsis determined within the VODE solver, and the frequency atwhich the solver is called is controlled by a longer model time stepbased on fixed upward parcel displacement increments (Δz) at thespecified w (Δt=Δz/w). Model calculations proceed until the parcelreaches a height 50m beyond the level of maximum saturation ratio(Smax) or total parcel displacement exceeds 2 km. Upon model termination,Smax and the fractional number of aerosols (factv) resulting innewly formed cloud droplets, defined as particles having diameters of atleast 2 μm, are cataloged in the LUTs according to the specified T, w, Na,rg, and ε parameter values.A critical point regarding the creation of these LUTs is the use of afixed minimum diameter (Dmin) to define cloud droplets in the parcelmodel, which can produce erroneous CDSD characteristics, particularlywithin the UL regime. For aerosol distributions with large rg valuesunder low SS conditions, for example, deliquesced aerosols within thelarge tail of the distribution can exceed 2 μm in diameter yet remainunactivated as ‘haze’ particles (Levin and Cotton, 2009; McFigganset al., 2006). For this study, aerosol particles activate to cloud dropletsbased on the critical diameter Dcrit as a function of parcel supersaturationratio (Sr) as in R09:D = σ MS RTρ83 ln( ) critsol wr w (2)where σsol is the surface tension of a solution droplet, Mw and ρw are themolar mass and density of liquid water, respectively, and R is theuniversal gas constant. Additionally, at relatively large w values withinthe AL regime, Nc stabilizes shortly after reaching supersaturation.However, parcel ascent and condensational growth continue beyondthe level of Smax, potentially causing additional narrowing of the CDSD.In the current work, model calculations terminate upon reaching Smaxas changes in Nc are negligible with continued ascent (Peng et al., 2007;R09).Parcel model simulations are performed to examine the sensitivitiesof CDSD characteristics to w, Na, rg, and ε, with the ranges for theseparameters listed in Table 1. Aerosols are assumed to be a mix of solubleand insoluble material of equal density, specified by ε, where fullysoluble aerosols correspond to ammonium sulfate with hygroscopicityparameter κ=0.61 (Petters and Kreidenweis, 2007). FollowingSvdH13, aerosol geometric standard deviation is fixed at σg=1.8, andaerosol distributions (Eq. 1) are partitioned into 100 logarithmicallyspacedbins spanning a size range specific to each rg value. For all simulations,Δz=1 m, and initial values of relative humidity, air temperatureand pressure are set to RH=0.99, T=10 °C and p=900 hPa,respectively.
The present study deals with several issues involving the improvement of physical submodels and the computational efficiency in modeling dense fuel sprays. To improve the computational efficiency, a parcel PDF approach is implemented which can account for turbulent dispersion within each computational parcel. The advantage of a parcel PDF tracking method is to reduce the number of computational parcels representing the spray dynamics as well as to obtain grid-independent solutions for two-phase flows. To account for the dense spray effects, an existing drop collision and coalescence model, two breakup models, and a Reitz's wave instability model were used. These models were incorporated into a state-of-the-art multiphase all-speed transient flow solution procedure. Comparative performance for two breakup models as well as the turbulence modulation effects are also studied. Validation cases include the nonevaporating and evaporating solid-cone dense sprays. The predictions show a reasonably good agreement with available experimental results in terms of spray penetration, drop sizes, gas and drop mean velocities, and gas and drop rms velocities. The numerical results indicate that the present parcel PDF model has the capability of accurately representing drop dispersion in dense sprays with manageable number of computational parcels.
In this paper, we have looked for enhancements of the O/N2 ratio in data measured by the Dynamics Explorer 2 (DE 2) satellite in the middle latitudes of the winter hemisphere, based on a prediction that was made by the National Center for Atmospheric Research thermosphere/tonosphere general circulation model (NCAR-TIGCM) that such increases occur. The NCAR-TIGCM predicts that these enhancements should be seen throughout the low latitude region and in many middle latitude locations, but that the enhancements in O/N2 are particularly strong in the middle-latitude, evening-to-midnight sector of the winter hemisphere. When this prediction was used to look for these effects in DE 2 NACS (neutral atmosphere composition spectrometer) data, large enhancements in the O/N2 ratio (approx. 50 to 90%) were seen. These enhancements were observed during the main phase of a storm that occurred on November 24, 1982, and were seen in the same region of the winter hemisphere predicted by the NCAR-TIGCM. They are partially the result of the depletion of N2 and, as electron loss is dependent on dissociative recombination at F(sub 2) altitudes, they have implications for electron densities in this area. Parcel trajectories, which have been followed through the NCAR-TIGCM history file for this event, show that large O/N2 enhancements occur in this limited region in the winter hemisphere for two reasons. First, these parcels of air are decelerated by the antisunward edge of the ion convection pattern; individual parcels converge and subsidence occurs. Thus molecular-nitrogen-poor air is brought from higher to lower heights. Because neutral parcels that are found a little poleward of the equatorial edge of the eveningside convection pattern are swept inward toward the center of the auroral oval, the enhancements occur only in a very limited range of latitudes. Second, nitrogen-poor air is transported from regions close to the magnetic pole in the winter hemisphere. During geomagnetic storms, enhanced meridional winds are driven by the increased pressure-gradient force that is associated with intensified Joule heating in the auroral oval. These pressure-driven winds decrease rapidly on the dayside beyond the auroral oval where the parcels originate, limiting the region into which the parcels can be transported. Thus these two processes drive values of O/N2 in a limited region of the winter hemisphere, and reinforce only in the evening sector, causing large changes in this region.
The existence of the multi-year HALOE CH4 data set, together with some comparisons of forward with back trajectory calculations which we have carried out, has motivated us to reexamine the question of polar vortex descent. Three-dimensional diabatic trajectory calculations have been carried out for the seven month fall to spring period in both the northern hemisphere (NH) and southern hemisphere (SH) polar stratosphere for the years 1992-1999. These computations are compared to fixed descent computations where the parcels were fixed at their latitude-longitude locations and allowed to descend without circulating. The forward trajectory computed descent is always less than the fixed descent due to horizontal parcel motions and variations in heating rates with latitude and longitude. Although the forward calculations estimate the maximum amount of descent that can occur, they do not necessarily indicate the actual origin of springtime vortex air. This is because more equator-ward air can be entrained within the vortex during its formation. To examine the origin of the springtime vortex air, the trajectory model was run backward for seven months from spring to fall. The back trajectories show a complex distribution of parcels in which one population originates in the upper stratosphere and mesosphere and experiences considerable descent in the polar regions, while the remaining parcels originate at lower altitudes of the middle and lower stratosphere and are mixed into the polar regions during vortex formation without experiencing as much vertical transport. The amount of descent experienced by the first population shows little variability from year to year, while the computed descent and mixing of the remaining parcels show considerable interannual variability due to the varying polar meteorology. Because of this complex parcel distribution it is not meaningful to speak of a net amount of descent experienced over the entire winter period. Since the back trajectories indicate that much of the air can come from lower altitudes than would be implied by the forward calculations, using a comparison between pre-winter and post-winter tracer profiles to estimate the amount of descent over this period will give erroneous descent amounts. In order to evaluate the computed descent, spring methane amounts were computed by mapping HALOE fall observations onto the final latitude-altitude locations of the back trajectories. These locations indicate the origin of the spring vortex air. The agreement between the computed means and the spring HALOE means is generally within 0.1-0.2 ppmv in the NH and 0.1-0.4 ppmv in the SH.
We evaluate evidence for stratospheric mass transport into, and mass remaining in, the troposphere during an intense midlatitude cyclone. Mesoscale forecast model analysis fields from the Mesoscale Analysis and Prediction System were matched with total ozone observations from the Total Ozone Measurement Spectrometer. Combined with parcel back trajectory calculations, the analyses imply that two mechanisms contributed to the mass exchange: (1) An area of dynamically induced exchange was observed on the cyclone's southern edge. Parcels originally in the stratosphere crossed the jet core and were diluted through turbulent mixing with tropospheric air; (2) Diabetic effects reduced parcel potential vorticity (PC) for trajectories traversing precipitation regions, creating a 'PV hole' signature in the center of the cyclone. Air with characteristics of ozone and water vapor found in the lower stratosphere remained in the troposphere. The strength of the latter process may be unusual. Combined with other research, these results suggest that precipitation-induced diabetic effects can significantly modify (either decreasing or increasing) parcel potential vorticity, depending on parcel trajectory configuration with respect to maximum heating regions and jet core. The diabetic heating effect on stratosphere-troposphere exchange (STE) is more important to tropopause erosion than to altering parcel trajectories. In addition, these results underline the importance of using not only PC but also chemical constituents for diagnoses of STE.
Phoenix, Arizona is the hottest large city in the United States with an average summer daytime temperature of 106°F. Temperatures in Phoenix continue to climb due to increasing global greenhouse gas concentrations and regional urbanization. The impacts of high temperatures, including heat-related illnesses and deaths, are disproportionately concentrated in low-income neighborhoods often characterized by little tree canopy, lack of green space, and insufficient access to shade. The City of Phoenix’s Office of Heat Response and Mitigation and Arizona State University’s Urban Climate Research Center, partnered with NASA DEVELOP to identify residential neighborhoods and parcels within qualified census tracts (QCTs) to be prioritized for tree planting initiatives using funding from the American Rescue Plan Act (ARPA). This project conducted analyses using NASA Earth observations, socioeconomic data from the 2019 American Community Survey, and local tree canopy and mobility data. For Earth observations, daytime land surface temperature, vegetation, and land cover were obtained from the Landsat 8 Thermal Infrared Sensor (TIRS) and Operational Land Imager (OLI). The project team incorporated these data into a heat vulnerability index (HVI) with an emphasis on tree canopy and social vulnerability to rank block groups within QCTs and focused on the resulting top 25 block group HVI scores. These top 25 most vulnerable block groups were then processed through a parcel analysis to determine the feasibility of residential tree planting based on building footprints on each parcel. Within the top 25 block groups, 2,411 parcels were analyzed, and 3,133 existing trees were identified averaging 1.3 trees per parcel with 90% of parcels having 3 trees or less. Homes with 2 trees or fewer were considered high priority for future planting efforts. Based on the City’s goals for increased tree canopy the project team determined that <10,000 additional trees would need to be planted within the most vulnerable 25 block groups. The project findings have helped initiate community engagement efforts and have contributed to the approval of tree planting funds in Phoenix.
The central role of sea ice within the Arctic climate system is clear, as is the importance of sea ice properties in the Arctic Ocean’s evolution and response to climate change. Less clear is how the history of the sea ice floe (in winter and spring) precondition the summer thickness evolution. Accurate accounting for the collective influences of processes that affect the summer sea ice thickness evolution is necessary for determining sea ice survivability and the rate of future sea ice loss. In this study, we use a satellite-derived sea ice parcel database co-located with the MOSAiC floe to explore the factors that precondition the summer sea ice thickness evolution. First. we compare the >500 collocated satellite-tracked sea ice parcels that intersected the MOSAiC drift track with observations to evaluate the satellite-derived surface energy budget and surface properties, finding agreement with correlations >0.9. Approximately 60% of the satellite tracked parcels survived the summer melt season. Secondly, we explore the spread in sea ice thickness and survivability across the satellite-tracked parcels, finding contributions from spring/summer surface albedo evolution, surface temperature, clouds, and parcel sea ice concentration. The satellite tracked sea ice parcel analysis provides insights into the spatial coherence of the factors that influence sea ice thickness evolution.
The state of Arctic sea ice influences aspects of many global systems including ecosystems, economies, geopolitics, and climate. And, it is in severe decline. Numerous studies have considered the factors that influence Arctic sea ice decline using monthly mean or gridded data. However, these methods struggle to cleanly separate the influences on sea ice melt because the sea ice cover is always moving; thus, different sea ice floes and sea ice regimes constitute the monthly average sea ice state within the Eulerian frame. This study uses a LaGrangian sea ice parcel tracking satellite database to investigate the factors that influence sea ice parcel survivability. The concept of sea ice survivability is defined as the likelihood that a sea ice parcel will last through the summer melt season. Survivability can be stratified in many ways to analyze how it changes as a function of sea ice regime, region, and sea ice characteristics. How is the survivability of first year (FY) and multi-year (MY) ice classes influenced by sea ice motion? This is the central question addressed in this presentation. Specifically, this analysis determines sea ice survivability as a function of distance traveled by sea ice parcels and compares the survivability statistics for sea ice parcels that stay in their starting region versus those that leave their starting region. The preliminary results indicate that the influence of sea ice movement on survivability depends strongly on which region the sea ice resides in at the beginning of the growth season. Lastly, we find substantial inter-annual variability in the regional distribution of sea ice survivability that seems to be strongly linked to sea ice parcel drift between regions of the Arctic Ocean. Moreover, we use these results to discuss the potential contribution of changes in sea ice motion on the observed slowdown in Arctic sea ice extent loss.
Simulations of Antarctic ozone photochemistry performed using a photochemical model integrated along air parcel trajectories are described. This type of model has a major advantage at high latitudes of being able to simulate correctly the complex interaction between photolysis and temperature fields, which, because of the polar night cannot be represented accurately in a zonally averaged framework. Isentropic air parcel trajectories were computed using Meteorological Office global model analyses and forecast fields from positions along the ER-2 flight paths during the Airborne Antarctic Ozone Experiment in Austral Spring 1987. A photochemical model is integrated along these trajectories using the aircraft observations to initialize constituent concentrations. The model includes additional reactions of the ClO dimer and also bromine reactions, which are thought to play a significant role in Antarctica. The model also includes heterogeneous reactions which are invoked when the air parcel passes through a polar stratospheric cloud (PSC). The existence of a PSC is determined throughout the course of the model integration from the parcel temperature and the saturated vapour pressure of water over an assumed H2O/HNO3 mixture. The air parcel temperature is used to determine the saturated vapor pressure of HNO3 over the same mixture. Mixing ratios which exceed saturation result in condensation of the excess in the model and hence lead to a reduction of the amount of gas phase NO2 available for chemical reaction.
We present here the chemical composition of outflow from the Asian continent to the atmosphere over the western Pacific basin during the Pacific Exploratory Mission-West (PEM-West B) in February-March 1994. Comprehensive measurements of important tropospheric trace gases and aerosol particulate matter were performed from the NASA DC-8 airborne laboratory. Backward 5 day isentropic trajectories were used to partition the outflow from two major source regions- continental north (greater than 20 deg N) and continental south (less than 20 deg N). Air parcels that had not passed over continental areas for the previous 5 days were classified as originating from an aged marine source. The trajectories and the chemistry together indicated that there was extensive rapid outflow of air parcels at altitudes below 5 km, while aged marine air was rarely encountered and only at less than 20 deg N latitude. The outflow at low altitudes had enhancements in common industrial solvent vapors such as C2Cl4, CH3CCl3, and C6H6, intermixed with the combustion emission products C2H2, C2H6, CO, and NO. The mixing ratios of all species were up to tenfold greater in outflow from the continental north compared to the continental south source region, with Pb-210 concentrations reaching 38 fCi (10(exp -15) curies) per standard cubic meter. In the upper troposphere we again observed significant enhancements in combustion-derived species in the 8-10 km altitude range, but water-soluble trace gases and aerosol species were depleted. These observations suggest that ground level emissions were lofted to the upper troposphere by wet convective systems which stripped water-soluble components from these air parcels. There were good correlations between C2H2 and CO and C2H6 (r(sup 2) = 0.70 - 0.97) in these air parcels and much weaker ones between C2H2 and H2O2 or CH3OOH (r(sup 2) = 0.50). These correlations were the strongest in the continental north outflow where combustion inputs appeared to be recent (1 - 2 days old). Ozone and PAN showed general correlation in these same air parcels but not with the combustion products. It thus appears that several source inputs were intermixed in these upper tropospheric air masses, with possible contributions from European or Middle Eastern source regions. In aged marine air mixing ratios of 03 (approximately equals 20 parts per billion by volume) and PAN (less than or equal to 10 parts per trillion by volume) were nearly identical at less than 2 km and 10 - 12 km altitudes due to extensive convective uplifting of marine boundary layer air over the equatorial Pacific even in wintertime. Comparison of the Pacific Exploratory Mission-West A and PEM-West B data sets shows significantly larger mixing ratios of SO2 and H2O2 during PEM-West A. Emissions from eruption of Mount Pinatubo are a likely cause for the former, while suppressed photochemical activity in winter was probably responsible for the latter. This comparison also highlighted the twofold enhancement in C2H2, C2H6, and C3H8 in the continental north outflow during /PEM-West B. Although this could be due to reduced OH oxidation rates of these species in wintertime, we argue that increased source emissions are primarily responsible.
We present preliminary results from a systematic study using simultaneous data from three spacecraft, Wind, IMP 8 (Interplanetary Monitoring Platform) and Geotail to examine interplanetary length scales and their implications on predictability for magnetic field parcels in the typical solar wind. Time periods were selected when the plane formed by the three spacecraft included the GSE (Ground Support Equipment) x-direction so that if the parcel fronts were strictly planar, the two adjacent spacecraft pairs would determine the same phase front angles. After correcting for the motion of the Earth relative to the interplanetary medium and deviations in the solar wind flow from radial, we used differences in the measured front angle between the two spacecraft pairs to determine structure radius of curvature. Results indicate that the typical radius of curvature for these IMF parcels is of the order of 100 R (Sub E). This implies that there are two important IMF (Interplanetary Magnetic Field) scale lengths relevant to predictability: (1) the well-established scale length over which correlations observed by two spacecraft decay along a given IMF parcel, of the order of a few tens of Earth radii and (2) the scale length over which two spacecraft are unlikely to even observe the same parcel because of its curvature, of the order of a hundred Earth radii.
We present evidence for stratospheric mass transport into, and remaining in, the troposphere in an intense midlatitude cyclone. Mesoscale forecast model analysis fields from the Mesoscale Analysis and Prediction System (MAPS) were compared with total ozone observations from the Total Ozone Measurement Spectrometer (TOMS). Coupled with parcel back-trajectory calculations, the analyses suggest two mechanisms contributed to the mass exchange: (1) A region of dynamical ly-induced exchange occurred on the cyclone's southern edge. Parcels originally in the stratosphere crossed the jet core and experienced dilution by turbulent mixing with tropospheric air. (2) Diabatic effects reduced parcel potential vorticity (PV) for trajectories traversing precipitation regions, resulting in a "PV-hole" signature in the cyclone center. Air with lower-stratospheric values of ozone and water vapor was left in the troposphere. The strength of the latter process may be atypical. These results, combined with other research, suggest that precipitation-induced diabatic effects can significantly modify, (either decreasing or increasing) parcel potential vorticity, depending on parcel trajectory configuration with respect to jet core and maximum heating regions. In addition, these results underscore the importance of using not only PV but also chemical constituents for diagnoses of stratosphere-troposphere exchange (STE).
This project was focused on modeling the processes by which increasing demand for developed land uses, brought about by changes in the regional economy and the socio-demographics of the region, are translated into a changing spatial pattern of land use. Our study focused on a portion of the Chesapeake Bay Watershed where the spatial patterns of sprawl represent a set of conditions generally prevalent in much of the U.S. Working in the region permitted us access to (i) a time-series of multi-scale and multi-temporal (including historical) satellite imagery and (ii) an established network of collaborating partners and agencies willing to share resources and to utilize developed techniques and model results. In addition, a unique parcel-level tax assessment database and linked parcel boundary maps exists for two counties in the Maryland portion of this region that made it possible to establish a historical cross-section time-series database of parcel level development decisions. Scenario analyses of future land use dynamics provided critical quantitative insight into the impact of alternative land management and policy decisions. These also have been specifically aimed at addressing growth control policies aimed at curbing exurban (sprawl) development. Our initial technical approach included three components: (i) spatial econometric modeling of the development decision, (ii) remote sensing of suburban change and residential land use density, including comparisons of past change from Landsat analyses and more traditional sources, and (iii) linkages between the two through variable initialization and supplementation of parcel level data. To these we added a fourth component, (iv) cellular automata modeling of urbanization, which proved to be a valuable addition to the project. This project has generated both remote sensing and spatially explicit socio-economic data to estimate and calibrate the parameters for two different types of land use change models and has undertaken analyses of these models. One (the CA model) is driven largely by observations on past patterns of land use change, while the other (the EC model) is driven by mechanisms of the land use change decision at the parcel level. Our project may be the first serious attempt at developing both types of models for the same area, using as much common data as possible. We have identified the strengths and weaknesses of the two approaches and plan to continue to revise each model in the light of new data and new lessons learned through continued collaboration. Questions, approaches, findings, publication and presentation lists concerning the research are also presented.