MedFlux: Investigations of particle flux in the Twilight Zone

The MedFlux project was devised to determine and model relationships between organic matter and mineral ballasts of sinking particulate matter in the ocean. Specifically we investigated the ballast ratio hypothesis, tested various commonly used sampling and modeling techniques, and developed new technologies that would allow better characterization of particle biogeochemistry. Here we describe the rationale for the project, the biogeochemical provenance of the DYFAMED site, the international support structure, and highlights from the papers published here. Additional MedFlux papers can be accessed at the MedFlux web site (http://msrc.sunysb.edu/MedFlux/).     

Some statistical characteristics of large deepwater waves around the Korean Peninsula

The relationship between significant wave height and period, the variability of significant wave period, the spectral peak enhancement factor, and the directional spreading parameter of large deepwater waves around the Korean Peninsula have been investigated using various sources of wave measurement and hindcasting data. For very large waves comparable to design waves, it is recommended to use the average value of the empirical formulas proposed by Shore Protection Manual in 1977 and by Goda in 2003 for the relationship between significant wave height and period. The standard deviation of significant wave periods non-dimensionalized with respect to the mean value for a certain significant wave height varies between 0.04 and 0.21 with a typical value of 0.1 depending upon different regions and different ranges of significant wave heights. The probability density function of the peak enhancement factor is expressed as a lognormal distribution, with its mean value of 2.14, which is somewhat smaller than the value in the North Sea. For relatively large waves, the probability density function of the directional spreading parameter at peak frequency is also expressed as a lognormal distribution.     

Benthic habitat quality change as measured by macroinfauna community in a tidal flat on the west coast of Korea

Yellow Sea tidal flats are internationally recognised for their contribution to biological diversity and yet are under enormous pressure from reclamation, pollution and overexploitation. The benthic macroinfauna community is the dominant community on these tidal flats and a reliable indicator of benthic environmental changes. We surveyed the current benthic macroinfauna community of the Ganghwa Southern Tidal Flat, the largest remaining Korean mud flat in the Yellow Sea, in order to examine changes in the environmental situation of this benthic ecosystem. The results show a significant decline in species diversity from the last survey made in 2003, and a shift in species composition with appearances of polychaetes indicative of pollution and physical disturbances and other opportunistic species becoming dominant in both density and biomass. The benthic community shift observed during the two study periods may be associated with increased nutrient pollution as well as increased physical disturbances in this area. However, we recognise the limitations of the data both in frequency and scope but believe the significant changes to the composition of the benthic fauna are sufficient to warrant concern. Observations are required to examine the extent to which these human activities induce benthic community shift in this tidal flat.     

Optimal design of two-dimensional wings in ground effect using multi-objective genetic algorithm

The shape optimization of the 2-dimensional wing in ground effect (WIG) has been performed by the integration of CFD (computational fluid dynamics) and MOGA (multi-objective genetic algorithm). Because of the trade-off between the aerodynamic forces and the height stability, it is difficult to satisfy the design requirements of efficiency and stability at the same time. In this study, the lift coefficient, the lift-drag ratio and the static height stability are chosen as the objective functions to obtain the optimal wing profiles of a WIG craft. An NACA0015 airfoil is used for the baseline model; the aerodynamic characteristics of the base model are compared with that of the optimal solutions. The profile of the airfoil is constructed by four Bezier curves with fourteen control points resulting in the eighteen coordinates, which are adopted as the design variables. The optimal solutions of the multi-objective optimization are not unique but a set of the non-dominated optima: the Pareto frontiers or a Pareto set. As the results of the multi-objective optimization, the forty Pareto optima, which include high-lift, high-efficiency, and more stable airfoils on the edge of the 3-dimensional objective space, are obtained at thirty evolutions of the generation.     

Designation processes for marine protected areas in the coastal wetlands of South Korea

The history of Korean tidal flat management and the process for designating Coastal Wetland Protected Areas (CWPAs) are described. Korean coastal wetlands have a long history of intensive use through reclamation for agricultural and industrial uses in the 20th century. Recently, the management policy is shifting away from intensive use towards the conservation of wetlands. This shift is caused by increasing public awareness of the value of wetlands and strong institutional support from the government. Since the Wetlands Conservation Act was passed in 1999, a total of twelve CWPAs have been designated through both top-down and bottom-up processes. Three designation paths are classified based on the relevant drivers, namely government-driven designations (seven CWPAs), local community driven designations (three CWPAs), and conflict resolution (trade-offs) driven designation (two CWPAs). The lessons learned from the designation of Korean CWPAs is that diversification of designation process could facilitate voluntary participation of local stakeholders and thereby enhance the chance of successful implementation of wise use strategy of tidal flats.     

Application of the weighted residual method to diffraction by 2-D canyons of arbitrary shape: I. Incident SH waves

A solution for the two-dimensional scattering and diffraction of plane SH waves by canyons of arbitrary shape in an elastic half space is presented. The wave field for arbitrary geometry in this paper is computed numerically by the method of weighted residues (moment method). The wave displacement field computed by the present residual method for the case of a semi-circular canyon was shown to agree analytically and numerically with that computed by the exact closed form series solution. The same observations about ground amplifications, their dependence on frequencies and orientations of the incident waves, can be stated here for canyons of arbitrary shape as previously made for circular canyons.     

Coastal environmental assessment and management by ecological simulation in Yeoja Bay, Korea

An eco-hydrodynamic model was used to estimate the carrying capacity of pollutant loads and response of water quality to environmental change in Yeoja Bay, Korea. An energy-system model also was used to simulate the fluctuation in nutrients and organic matter in the bordering wetland. Most water quality factors showed a pulsed pattern, and the concentrations of nutrients and organic matter of seawater increased when input loads of nutrients increased due to freshwater discharge. The well-developed tidal zones and wetlands in the northern area of the bay were highly sensitive to input loads. Residence times of water, chemical oxygen demand (COD), and dissolved inorganic nitrogen (DIN) within the bay were estimated to be about 16 days, 43.2 days, and 50.2 days, respectively. Water quality reacted more sensitively to the effects of nitrogen and phosphorus input than to COD. A plan to reduce the present levels of COD and dissolved inorganic phosphorus (DIP) by 20–30% and DIN by at least 50% in pollutant loads is needed for satisfying the target water quality criteria. The natural removal rate of nutrients in wetlands by reeds was assessed to be approximately 10%.     

Origin of the hydrocarbon component of carbonaceous chondrites: the star-meteorite connection

We report the results of an experiment that produced a residue which closely matches the hydrocarbon component of the Murchison carbonaceous chondrite. This experiment suggests that the parent material of the meteoritic component originated as polycyclic aromatic hydrocarbon species in carbon stars during their later stages of evolution. The experiments also indicate that the pathway from those formation sites to eventual incorporation into the meteorite parent body involved hydrogenation in a plasma in the solar nebula or in H II regions prior to the solar nebula. This model is consistent with what is known about the meteoritic hydrocarbon component including deuterium abundance, the observation of cosmic infrared emission bands best attributed to polycyclic aromatic hydrocarbon molecules, and the inherent stability of these molecules that allows their formation in stars and subsequent survival in the interstellar medium.     

Effects of Geometries and Substructures of ICMEs on Geomagnetic Storms

To better understand geomagnetic storm generations by ICMEs, we consider the effect of substructures (magnetic cloud, MC, and sheath) and geometries (impact location of flux-rope at the Earth) of the ICMEs. We apply the toroidal magnetic flux-rope model to 59 CDAW CME–ICME pairs to identify their substructures and geometries, and select 20 MC-associated and five sheath-associated storm events. We investigate the relationship between the storm strength indicated by minimum Dst index \((\mathrm{Dst}_{\mathrm{min}})\) and solar wind conditions related to a southward magnetic field. We find that all slopes of linear regression lines for sheath-storm events are steeper (\({\geq}\,1.4\)) than those of the MC-storm events in the relationship between \(\mathrm{Dst}_{\mathrm{min}}\) and solar wind conditions, implying that the efficiency of sheath for the process of geomagnetic storm generations is higher than that of MC. These results suggest that different general solar wind conditions (sheaths have a higher density, dynamic and thermal pressures with a higher fluctuation of the parameters and higher magnetic fields than MCs) have different impact on storm generation. Regarding the geometric encounter of ICMEs, 100% (2/2) of major storms (\(\mathrm{Dst}_{\mathrm{min}} \leq -100~\mbox{nT}\)) occur in the regions at negative \(P_{Y}\) (relative position of the Earth trajectory from the ICME axis in the \(Y\) component of the GSE coordinate) when the eastern flanks of ICMEs encounter the Earth. We find similar statistical trends in solar wind conditions, suggesting that the dependence of geomagnetic storms on 3D ICME–Earth impact geometries is caused by asymmetric distributions of the geoeffective solar wind conditions. For western flank events, 80% (4/5) of the major storms occur in positive \(P_{Y}\) regions, while intense geoeffective solar wind conditions are not located in the positive \(P_{Y}\). These results suggest that the strength of geomagnetic storms depends on ICME–Earth impact geometries as they determine the solar wind conditions at Earth.     

Lie group variational integrators for the full body problem in orbital mechanics

Equations of motion, referred to as full body models, are developed to describe the dynamics of rigid bodies acting under their mutual gravitational potential. Continuous equations of motion and discrete equations of motion are derived using Hamilton’s principle. These equations are expressed in an inertial frame and in relative coordinates. The discrete equations of motion, referred to as a Lie group variational integrator, provide a geometrically exact and numerically efficient computational method for simulating full body dynamics in orbital mechanics; they are symplectic and momentum preserving, and they exhibit good energy behavior for exponentially long time periods. They are also efficient in only requiring a single evaluation of the gravity forces and moments per time step. The Lie group variational integrator also preserves the group structure without the use of local charts, reprojection, or constraints. Computational results are given for the dynamics of two rigid dumbbell bodies acting under their mutual gravity; these computational results demonstrate the superiority of the Lie group variational integrator compared with integrators that are not symplectic or do not preserve the Lie group structure.