SEARCH · Search NASA

Results for “Convective mixing layer”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 73 records · Page 4

All sources NASA NTRS DOE OSTI

NASA NTRS2016

Small Scale Motions Observed by Aircraft in the Tropical Tropopause Layer - Convective and Non-Convective Environments

Indirect evidence indicates a role for vertical mixing in the Tropical Tropopause Layer (TTL). In particular, detailed model studies suggest that such vertical mixing may be required to explain the value of the water vapor minimum in the TTL. There have been previous observations during the STEP Tropical aircraft campaign (1987) of bursts of high frequency activity associated with convectively generated gravity waves in the tropical western Pacific. Higher frequency, higher quality measurements from NASA high altitude aircraft (ER-2, WB-57, and Global Hawk) have been made available in the last 20 years. These include measurements of vertical velocity and other meteorological parameters. Most recently, during the ATTREX Global Hawk aircraft mission (Airborne Tropical TRopopause EXperiment), there have been extensive measurements at all altitudes of the TTL in both convective (winter western Pacific) and less convective (winter eastern Pacific) regions. This presentation represents an initial analysis of high frequency small scale (a few km max) meteorological measurements from the ATTREX dataset. We obtain some basic information about the distribution and character of high frequency activity in vertical velocity in the TTL. In particular, we focus on relating the high frequency activity to nearby tropical convection and to vertical shears associated with gravity and inertia-gravity waves.

convection↗

NASA NTRSDate not supplied

Understanding and Utilizing PBL Height Data from Multiple Observing Systems in the GEOS System

The accuracy of PBL height simulation is a key issue in many applications including forecasting near surface meteorology and air quality, however, it is a very challenging problem due to the lack of not only comprehensive, global Planetary Boundary Layer (PBL) observations but also a strategy and infrastructure to utilize PBL height data from a variety of sensors. Following the designation of PBL as an incubation class observable in the 2017 Decadal Survey, the PBL Incubation Study Team Report [14] made clear that “a future global PBL observing system requires modeling and data assimilation as essential components.” There is an urgent need for global modeling development in order to utilize Program of Record (POR) observations, assess their impacts, and identify gaps to be filled by future PBL missions. Our overall objective is to develop PBL data assimilation capabilities in the NASA Global Earth Observing System (GEOS), focusing on PBL height from multiple observing systems, to support the assessment and use of future PBL observations. The NASA GEOS system is composed of the GEOS global atmospheric general circulation model (AGCM) and the atmospheric data assimilation system (ADAS). The PBL parameterizations include the “Lock” K-profile scheme driven by surface and cloud-top buoyancy fluxes ([4]), and the “Louis” local scheme for stable conditions based on the Richardson number ([5]). Above the mixed layer defined by the Lock surface plume, shallow cumulus convection is represented by the mass flux scheme of [9]. Additional parameterizations are summarized in [1]. The ADAS employs the hybrid 4D Ensemble- Variational (EnVar) configuration ([15]), with the ensemble providing flow-dependent background error covariance information. The resultant analysis increments are fed back to the forecast model through the 4D incremental analysis update (IAU) approach ([11]). In this study, PBL height data are being or have been generated from radiosondes, GNSS RO, satellite (CATS, CALIPSO and ICESat-2) and ground-based (MPLNET) lidars, and wind profiler. Investigations have been conducted to specify quality marks for PBL height retrievals for the data assimilation purpose. These PBL height data have different strengths and weaknesses ([2], [3], [6], [7], [8], [10]), and the satellite PBL height data provide better global coverage and complement in-situ PBL height data. Radiosondes offer high accuracy and in situ measurement of temperature and humidity profiles, but with poor spatio-temporal sampling. The in-situ observing systems like MPLNET and wind profiler provide long history of PBL height records at each station. The GNSS RO based PBL height is retrieved based on the sharp gradients in refractivity profile that represent the fine vertical structure of temperature and moisture changes above the PBL. However, not all RO refractivity profiles reach the surface depending on location and regime, and RO refractivity retrievals can be negatively biased below 2km. The PBL height data from satellite lidars provide high resolution along track PBL height retrievals, but over land they are affected by previous day convective PBL aerosol and strongly associated with mixing layer and retrievals cannot be made below thick, attenuating clouds. A successful assimilation of PBL height data requires a thorough understanding of the observing method and the retrieval algorithm for each observing system in order to use the PBL height data from multiple observing systems properly. Due to the sensitivity of PBL height data to the observing method and choice of algorithm, it is important to use a model definition appropriate for each observation type to compute differences between PBL height data and model PBL height (OmFs). The GEOS model currently includes two PBL height definitions suitable for direct comparison with observed PBL height, and additional definitions are being added in this study. Evaluation of different model PBL height definitions is underway. Meanwhile, efforts have been made in the GEOS data assimilation system to develop PBL height data assimilation capability. PBL height data can be assimilated using two different approaches. The traditional approach is to construct an observation operator and its tangent linear and adjoint, which link control variables to PBL height data from each observing system. This observation operator can be very complicated, e.g., the lidar-based PBL height observation operator includes the backscatter lidar forward observation operator, the algorithm to derive PBL height from attenuated total backscatter, interpolation, and calculations handling the mismatch between observed and model scales. The other approach is to augment PBL height to the control variable vector, and it is adopted in this study. The latter approach was also used in previous studies, e.g., the assimilation of PBL height data from radiosonde and aircraft in the Real Time Mesoscale Analysis (RTMA) system for a dispersion modelling study ([13]); the PBL height assimilation study using lidar PBL height data at Greensburg, Kansas for a field campaign ([12]). The PBL height assimilation from multiple observing systems in this study allows us to take advantage of the diverse PBL height data that provide much better global coverage collectively under different meteorological conditions and with different temporal and spatial scales. As all the PBL heights are tightly coupled with the PBL thermodynamic variables, the strong correlations, which are provided by the 4D ensemble forecast, enable PBL height data from various sources to interact and combine coherently and provide additional information for PBL temperature and moisture fields. The results of comparisons among PBL height data from different sources and the evaluation of the model PBL height definitions with the PBL height data will be presented, and the PBL height data synergy strategies and preliminary results will also be discussed at the conference.

Y. Zhu↗

NASA NTRS1992

Integration of satellite and surface data using a radiative-convective oceanic boundary-layer model

A mixing-line boundary-layer model is used to retrieve cloud-top height from satellite-derived cloud-top temperatures, using 700-hPa National Meteorological Center (NMC) analyses and the Comprehensive Ocean and Atmosphere Data Set (COADS) surface data as supporting datasets. Results are compared with the fixed-lapse-rate method of retrieving boundary-layer depth from sea surface temperatures (SST) and cloud-top temperatures. A radiative-convective equilibrium boundary-layer model is used to retrieve boundary-layer structure given SST and surface wind, satellite cloud-top temperatures and cloud fraction, and the 700-hPa NMC thermodynamic analyses. Good agreement is found between the COADS data and the model solutions for low-level temperature and moisture. This suggests that equilibrium boundary-layer models may be of use over remote oceans in the retrieval of boundary-layer structure.

Betts, Alan K.↗

DOE OSTI2020

Observations of Disequilibrium CO Chemistry in the Coldest Brown Dwarfs

Cold brown dwarfs are excellent analogs of widely separated, gas giant exoplanets, and provide insight into the potential atmospheric chemistry and physics we may encounter in objects to be discovered by future direct imaging surveys. We present a low-resolution, R ∼ 300, M-band spectroscopic sequence of seven brown dwarfs with effective temperatures between 750 and 250 K along with Jupiter. These spectra reveal disequilibrium abundances of carbon monoxide (CO) produced by atmospheric quenching. We use the eddy diffusion coefficient (K {sub zz}) to estimate the strength of vertical mixing in each object. The K {sub zz} values of cooler gaseous objects are close to their theoretical maximum, and warmer objects show weaker mixing, likely due to less efficient convective mixing in primarily radiative layers. The CO-derived K {sub zz} values imply that disequilibrium phosphine (PH{sub 3}) should be easily observable in all of the brown dwarfs, but none as yet show any evidence for PH{sub 3} absorption. We find that ammonia is relatively insensitive to atmospheric quenching at these effective temperatures. We are able to improve the fit to WISE 0855's M-band spectrum by including both CO and water clouds in the atmospheric model.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS2002

Simulation of TRMM Microwave Imager Brightness Temperature using Precipitation Radar Reflectivity for Convective and Stratiform Rain Areas over Land

Rain is highly variable in space and time. In order to measure rainfall over global land with satellites, we need observations with very high spatial resolution and frequency in time. On board the Tropical Rainfall Measuring Mission (TRMM) satellite, the Precipitation Radar (PR) and Microwave Imager (TMI) are flown together for the purpose of estimating rain rate. The basic method to estimate rain from PR has been developed over the past several decades. On the other hand, the TMI method of rain estimation is still in the state development, particularly over land. The objective of this technical memorandum is to develop a theoretical framework that helps relate the observations made by these two instruments. The principle result of this study is that in order to match the PR observations with the TMI observations in convective rain areas, a mixed layer of graupel and supercooled water drops above the freezing level is needed. On the other hand, to match these observations in the stratiform region, a layer of snowflakes with appropriate densities above the freezing level, and a melting layer below the freezing level, are needed. This understanding can lead to a robust rainfall estimation technique from the microwave radiometer observations.

Prabhakara, C.↗

NASA NTRS2011

LES of Temporally Evolving Mixing Layers by an Eighth-Order Filter Scheme

An eighth-order filter method for a wide range of compressible flow speeds (H.C. Yee and B. Sjogreen, Proceedings of ICOSAHOM09, June 22-26, 2009, Trondheim, Norway) are employed for large eddy simulations (LES) of temporally evolving mixing layers (TML) for different convective Mach numbers (Mc) and Reynolds numbers. The high order filter method is designed for accurate and efficient simulations of shock-free compressible turbulence, turbulence with shocklets and turbulence with strong shocks with minimum tuning of scheme parameters. The value of Mc considered is for the TML range from the quasi-incompressible regime to the highly compressible supersonic regime. The three main characteristics of compressible TML (the self similarity property, compressibility effects and the presence of large-scale structure with shocklets for high Mc) are considered for the LES study. The LES results using the same scheme parameters for all studied cases agree well with experimental results of Barone et al. (2006), and published direct numerical simulations (DNS) work of Rogers & Moser (1994) and Pantano & Sarkar (2002).

Hadjadj, A↗

NASA NTRS2012

Comparative Study of Three High Order Schemes for LES of Temporally Evolving Mixing Layers

Three high order shock-capturing schemes are compared for large eddy simulations (LES) of temporally evolving mixing layers (TML) for different convective Mach numbers (Mc) ranging from the quasi-incompressible regime to highly compressible supersonic regime. The considered high order schemes are fifth-order WENO (WENO5), seventh-order WENO (WENO7) and the associated eighth-order central spatial base scheme with the dissipative portion of WENO7 as a nonlinear post-processing filter step (WENO7fi). This high order nonlinear filter method (H.C. Yee and B. Sjogreen, Proceedings of ICOSAHOM09, June 22-26, 2009, Trondheim, Norway) is designed for accurate and efficient simulations of shock-free compressible turbulence, turbulence with shocklets and turbulence with strong shocks with minimum tuning of scheme parameters. The LES results by WENO7fi using the same scheme parameter agree well with experimental results of Barone et al. (2006), and published direct numerical simulations (DNS) work of Rogers & Moser (1994) and Pantano & Sarkar (2002), whereas results by WENO5 and WENO7 compare poorly with experimental data and DNS computations.

Yee, Helen M. C.↗

NASA NTRS2019

The Role of Oceanic Processes in the Initiation of Boreal Winter Intraseasonal Oscillations over the Indian Ocean

Observational analyses and a hierarchy of ocean general circulation model (OGCM) experiments were performed to understand the influence of oceanic processes on the warm sea surface temperature anomalies (SSTAs) prior to the convection initiation of boreal winter intraseasonal oscillations (ISOs), including the Madden-Julian Oscillation (MJO), in the equatorial Indian Ocean. We found 39 strong ISOs that passed over the Indian Ocean Warm Pool region during the November-April season of the 2001-2012 period. 17/39 ISO events initiated in the Seychelles-Chagos Thermocline Ridge (SCTR) before propagating eastward; the remaining events initiated in the southern Arabian Sea (6) or Warm Pool (16) regions. The SCTR event set was notable in that it contained more global-scale MJOs (71-76%), as defined by the RMM and OMI indices, than the WP events (25-44%). Additionally, ~24% (44%) of the SCTR (Warm Pool) events were preceded by strong oceanic process-induced SSTAs of similar magnitude to those of shortwave radiative and turbulent heat fluxes. The Arabian Sea events, however, were not associated with statistically significant SSTA signals prior to convection. Based on a mixed layer heat budget analysis, entrainment and upwelling reduction were the dominant oceanic processes contributing to the warming, in contrast with boreal summer, when horizontal advection dominated. We examined several case studies, including primary MJO events, where oceanic Rossby waves were associated with the entrainment and upwelling reduction. Two simple atmospheric boundary layer convergence models revealed that the SSTAs contributed at least half of the total convergence and suggested that the ocean dynamical effect was responsible for the majority of SSTA-forced convergence for those case studies. These results underscore the need for climate prediction models to accurately represent the ocean structure and processes to include the effects of oceanic predictors.

West, B. Jason↗

NASA NTRS1993

Turbulence modeling for high speed compressible flows

The following grant objectives were delineated in the proposal to NASA: to offer course work in computational fluid dynamics (CFD) and related areas to enable mechanical engineering students at North Carolina A&T State University (N.C. A&TSU) to pursue M.S. studies in CFD, and to enable students and faculty to engage in research in high speed compressible flows. Since no CFD-related activity existed at N.C. A&TSU before the start of the NASA grant period, training of students in the CFD area and initiation of research in high speed compressible flows were proposed as the key aspects of the project. To that end, graduate level courses in CFD, boundary layer theory, and fluid dynamics were offered. This effort included initiating a CFD course for graduate students. Also, research work was performed on studying compressibility effects in high speed flows. Specifically, a modified compressible dissipation model, which included a fourth order turbulent Mach number term, was incorporated into the SPARK code and verified for the air-air mixing layer case. The results obtained for this case were compared with a wide variety of experimental data to discern the trends in the mixing layer growth rates with varying convective Mach numbers. Comparison of the predictions of the study with the results of several analytical models was also carried out. The details of the research study are described in the publication entitled 'Compressibility Effects in Modeling Turbulent High Speed Mixing Layers,' which is attached to this report.

Chandra, Suresh↗

NASA NTRS1984

Lidar observations of mixed layer dynamics - Tests of parameterized entrainment models of mixed layer growth rate

Ground based lidar measurements of the atmospheric mixed layer depth, the entrainment zone depth and the wind speed and wind direction were used to test various parameterized entrainment models of mixed layer growth rate. Six case studies under clear air convective conditions over flat terrain in central Illinois are presented. It is shown that surface heating alone accounts for a major portion of the rise of the mixed layer on all days. A new set of entrainment model constants was determined which optimized height predictions for the dataset. Under convective conditions, the shape of the mixed layer height prediction curves closely resembled the observed shapes. Under conditions when significant wind shear was present, the shape of the height prediction curve departed from the data suggesting deficiencies in the parameterization of shear production. Development of small cumulus clouds on top of the layer is shown to affect mixed layer depths in the afternoon growth phase.

Boers, R.↗

DOE OSTI2021

Idealized Large-Eddy Simulations of Stratocumulus Advecting over Cold Water. Part I: Boundary Layer Decoupling

Abstract We explore the decoupling physics of a stratocumulus-topped boundary layer (STBL) moving over cooler water, a situation mimicking warm-air advection (WADV). We simulate an initially well-mixed STBL over a doubly periodic domain with the sea surface temperature decreasing linearly over time using the System for Atmospheric Modeling large-eddy model. Due to the surface cooling, the STBL becomes increasingly stably stratified, manifested as a near-surface temperature inversion topped by a well-mixed cloud-containing layer. Unlike the stably stratified STBL in cold-air advection (CADV) that is characterized by cumulus coupling, the stratocumulus deck in the WADV is unambiguously decoupled from the sea surface, manifested as weakly negative buoyancy flux throughout the subcloud layer. Without the influxes of buoyancy from the surface, the convective circulation in the well-mixed cloud-containing layer is driven by cloud-top radiative cooling. In such a regime, the downdrafts propel the circulation, in contrast to that in CADV regime for which the cumulus updrafts play a more determinant role. Such a contrast in convection regime explains the difference in many aspects of the STBLs including the entrainment rate, cloud homogeneity, vertical exchanges of heat and moisture, and lifetime of the stratocumulus deck, with the last being subject to a more thorough investigation in Part II. Finally, we investigate under what conditions a secondary stratus near the surface (or fog) can form in the WADV. We found that weaker subsidence favors the formation of fog whereas a more rapid surface cooling rate does not. Significant Statement The low-lying blanketlike clouds, called stratocumulus (Sc), reflect much incoming sunlight, substantially modulating Earth’s temperature. While much is known about how the Sc evolves when it moves over warmer water, few studies examine the opposite situation of Sc moving over colder water. We used a high-resolution numerical model to simulate such a case. When moving over cold water, the Sc becomes unambiguously decoupled from the water surface, distinctive from its warm counterpart in which the Sc interacts with the water surface via intermittent cauliflower-like clouds called cumulus clouds. Such decoupling influences many aspects of the Sc–sea surface system, which combine to alter the ability of the Sc to reflect sunlight, thereby influencing the climate. This work laid the foundation for future work that quantifies the contribution of such a decoupled Sc regime to Earth’s radiative budget and climate change.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS1997

Vertical Transport Processes for Inert and Scavenged Species: TRACE-A Measurements

The TRACE-A mission of the NASA DC-8 aircraft made a large-scale survey of the tropical and subtropical atmosphere in September and October of 1992. Both In-situ measurements of CO (G. Sachsen NASA Langley) and aerosol size (J. Browell group, NASA Langley) provide excellent data sets with which to constrain vertical transport by planetary boundary layer mixing and deep-cloud cumulus convection. Lidar profiles of aerosol-induced scattering and ozone (also by Bremen) are somewhat require more subtle interpretation as tracers, but the vertical information on layering largely compensates for these complexities. The reason this DC-8 dataset is so useful is that very large areas of biomass burning over Africa and South America provide surface sources of appropriate sizes with which to characterize vertical and horizontal motions; the major limitation of our source description is that biomass burning patterns move considerably every few days, and daily burning inventories are a matter of concurrent, intensive research. We use the Penn State / NCAR MM5 model in an assimilation mode on the synoptic and intercontinental scale, and assess the success it shows in vertical transport descriptions. We find that the general level of emissions suggested by the climatological approach (Will. Has, U. of Montana) appears to be approximately correct, possibly a bit low, for this October, 1992, time period. Vertical transport in planetary boundary layer mixing to 5.5 kin was observed and reproduced in our simulations. Furthermore we find evidence that Blackader "transilient" or matrix-transport scheme is needed, but may require some adaptation in our tracer model: CO seems to exhibit very high values at the top of the planetary boundary layer, a process that stretches the eddy-diffusion parameterization. We will report on progress in improving the deep convective transport of carbon monoxide: the Grail scheme as we used it at 100 kin resolution did not transport enough material to the upper troposphere. We expect to be able to attribute this to either parameterization reasons (inadequacy of this parameterization at the large 100km scale) or other reasons. Nevertheless, the qualitative nature of deep transport by clouds shows up well in the simulations. As for scavengable species, the simulations predict tens of micrograms per standard cubic meter of smoke aerosol in the boundary layer. In a straightforward illustration of our simple bulk-mass scavenging parameterization, to one or two micrograms per standard cubic meter of smoke aerosol in the free troposphere just above the source regions: very high concentrations for the free troposphere. We expect to report on comparisons of these predictions to a variety of observations.

Chatfield, Robert B.↗

DOE OSTI2026

The Spatial and Temporal Variability of the Clear Convective Boundary Layer at the ARM SGP Supersite

The convective boundary layer (CBL), also known as the mixing layer, constitutes the critical lower segment of the atmosphere that significantly influences daily human activities. The growing demand for precise weather forecasts is driven by the requirements of agriculture, transportation, and routine societal functions. Here, to enhance understanding of the CBL, this study investigates the spatiotemporal variability in the CBL and its controlling factors using four-year Doppler lidar, surface flux, and profiling measurements at five ARM Southern Great Plains sites within a 100 km radius. This investigation utilizes data collected exclusively under clear-sky conditions or scattered low-cloud conditions. Results reveal significant spatial differences in CBL evolutions. Daily mixing layer heights (MLHs) vary up to 1 km (30% of the mean) in late afternoon. There is a clear east–west contrast: western sites (C1, E32, E37) exhibit higher summer MLH (1.9–2.1 km) and vertical velocity variances (1.0–1.2 m 2 s −2 ) than eastern sites (1.6–1.8 km), reversing in winter. Temporally, the MLH peaks at 70% of the sunrise–sunset interval, the lagging heat flux (HF) peaks at 50%; and the seasonal MLH maxima lag the HF by approximately one month, influenced by nighttime PBL (planetary boundary layer) properties. The HF and lower tropospheric stability are the main factors of influence for the CBL, but site-specific dependencies highlight the critical roles of local factors, underscoring the need for including them in CBL modeling.

ARM SGP site↗

NASA NTRS1978

Stack plume characterization and model assessment with lidar data

A mobile lidar system was used to obtain plume rise and dispersion measurements for validating plume models. Predictions of rise and spread from the two-thirds law were found to agree well with observations. These predictions were used in combination with Lamb's results for neutrally buoyant particle diffusion in the mixed layer to estimate dispersion when convective turbulence became important.

Weil, J. C.↗

NASA NTRS1990

Interface shapes during vertical Bridgman growth of (Pb, Sn)Te crystals

Melt-solid interfaces obtained during vertical Bridgman growth of (Pb, Sn)Te crystals were investigated with a quenching technique. The shapes of these interfaces, revealed by etching longitudinally cut sections, were correlated with the composition variations determined by EMPA. These experiments demonstrated that the interface shape can be changed from concave to convex by moving its location from the edge of the cold zone into the hot zone. The metallography and microsegregation near the melt-solid interface were analyzed in detail. A sharp change in composition above the interface indicated the existence of a diffusion boundary layer 40-90 microns thick. This small diffusion boundary layer is consistent with strong convective mixing in the (Pb, Sn)Te melt.

Huang, YU↗

NASA NTRS1993

Distribution of tropical tropospheric water vapor

Utilizing a conceptual model for tropical convection and observational data for water vapor, the maintenance of the vertical distribution of the tropical tropospheric water vapor is discussed. While deep convection induces large-scale subsidence that constrains the turbulent downgradient mixing to within the convective boundary layer and effectively dries the troposphere through downward advection, it also pumps hydrometeors into the upper troposphere, whose subsequent evaporation appears to be the major source of moisture for the large-scale subsiding motion. The development of upper-level clouds and precipitation from these clouds may also act to dry the outflow, thus explaining the low relative humidity near the tropopause. A one-dimensional model is developed to simulate the mean vertical structure of water vapor in the tropical troposphere. It is also shown that the horizontal variation of water vapor in the tropical troposphere above the trade-wind boundary layer can be explained by the variation of a moisture source that is proportional to the amount of upper-level clouds. Implications for the nature of water vapor feedback in global warming are discussed.

Sun, De-Zheng↗

NASA NTRSDate not supplied

Dynamic Adjustment of Model Parameters Using Assimilated Boundary Layer Height

The NASA Global Modeling and Assimilation Office is exploring the assimilation of boundary layer height (PBLH) derived from a variety of observations, including radiosondes, GNSS-RO, space- and ground-based lidar, and radar wind profilers. Previous work using the Goddard Earth Observing System (GEOS) has shown that within the lower troposphere, the representation of thermodynamic structure is strongly affected by model parameterizations. In assimilation experiments, parameterizations can compensate for analysis thermodynamic increments and reduce the observation information carried forward in time. To mitigate this effect, the present study explores a dynamic parameter adjustment based on the PBLH analysis field. The PBLH analysis increment is used to adjust model parameters that influence the depth of parameterized boundary layer mixing in both stable and convective regimes. Replay (nudging) experiments and cycled data assimilation experiments are used to evaluate the approach. The parameter adjustment is shown to reduce model compensation of the analysis thermodynamic tendencies and bring model parameterized mixing depths into closer agreement with observed PBLH.

Nathan Arnold↗

NASA NTRS2000

The Northeast Monsoon's Impact on Mixing, Phytoplankton Biomass and Nutrient Cycling in the Arabian Sea

In the northern Arabian Sea, atmospheric conditions during the Northeast (winter) Monsoon lead to deep convective mixing. Due to the proximity of the permanent pyncnocline to the sea surface, this mixing does not penetrate below 125 m. However, a strong nitracline is also present and the deep convection results in significant nitrate flux into the surface waters. This leads to nitrate concentrations over the upper 100 m that exceed 4 micrometers toward the end of the Monsoon. During the 1994/1995 US JGOFS/Arabian Sea expedition, the mean areal gross primary production over two successive Northeast Monsoons was determined to be 1.35gC/sq m/d. Thus, despite the deep penetrative convection, high rates of primary productivity were maintained. An interdisciplinary model was developed to elucidate the biogeochemical processes involved in supporting the elevated productivity. This model consists of a 1-D mixed-layer model coupled to a set of equations that tracked phytoplankton growth and the concentration of the two major nutrients (nitrate and ammonium). Zooplankton grazing was parameterized by rate constant determined by shipboard experiments. Model boundary conditions consist of meteorological time-series measured from the surface buoy that was part of the ONR Arabian Sea Experiment's central mooring. Our numerical experiments show that elevated surface evaporation, and the associated salinization of the mixed layer, strongly contributes to the frequency and penetration depth of the observed convective mixing. Cooler surface temperatures, increased nitrate entrainment, reduced water column stratification, and lower near-surface chlorophyll a concentrations all result from this enhanced mixing. The model also captured a dependence on regenerated nitrogen observed in nutrient uptake experiments performed during the Northeast Monsoon. Our numerical experiments also indicate that variability in mean pycnocline depth causes up to a 25% reduction in areal chlorophyll a concentration. We hypothesize that such shifts in pycnocline depth may contribute to the interannual variations in primary production and surface chlorophyll a concentration that have been previously observed in this region.

Wiggert, J. D.↗