The Fragmentation in a Protoplanetary Disk and the Properties of the Resultant Self-Gravitating Clumps

We study the fragmentation properties in the protoplanetary disk and properties of the resultant self-gravitating clumps using our newly constructed disk model. Our disk model includes the mass inflall term from a molecular cloud core and the photoevaporation winds effect. We adopt the conventional fragmentation criterion to judge whether a protoplanetary disk can fragment. In this work, we follow our previous work to investigate the properties of the resultant self-gravitating clumps. In our calculation, the initial masses of the resultant self-gravitating clumps lie in the range of tens of M_J to more than one hundred of M_J, where M_J is the Jupiter mass. These initial masses can seemingly account for the masses of extrasolar planets in magnitude. We also calculate the subsequent gas accretion of clumps in 1.27 × 10⁴ yr after the formation of self-gravitating clumps. We find that the subsequent gas accretion of self-gravitating clumps is very efficient, and the clump masses grow to hundreds of M_J and the physical radii R_c of clumps increase to about 10 AU. Additionally, we also calculate the orbital migration of clumps. We find that most clumps have short migration timescale to be accreted onto the protostar, and only a small fraction of clumps have long migration timescale (>10⁶ yr) to successfully become gas giant planets. These results are consistent with previous studies.     

Ensemble clustering for efficient robust optimization of naturally fractured reservoirs

In the development of naturally fractured reservoirs (NFRs), the existence of natural fractures induces severe fingering and breakthrough. To manage the flooding process and improve the ultimate recovery, we propose a numerical workflow to generate optimal production schedules for smart wells, in which the inflow control valve (ICV) settings can be controlled individually. To properly consider the uncertainty introduced by randomly distributed natural fractures, the robust optimization would require a large ensemble size and it would be computationally demanding. In this work, a hierarchical clustering method is proposed to select representative models for the robust optimization in order to avoid redundant simulation runs and improve the efficiency of the robust optimization. By reducing the full ensemble of models into a small subset ensemble, the efficiency of the robust optimization algorithm is significantly improved. The robust optimization is performed using the StoSAG scheme to find the optimal well controls that maximize the net-present-value (NPV) of the NFR’s development. Due to the discrete property of a natural fracture field, traditional feature extraction methods such as model-parameter-based clustering may not be directly applicable. Therefore, two different kinds of clustering-based optimization methods, a state-based (e.g., s _w profiles) clustering and a response-based (e.g., production rates) clustering, are proposed and compared. The computational results show that the robust clustering optimization could increase the computational efficiency significantly without sacrificing much expected NPV of the robust optimization. Moreover, the performance of different clustering algorithms varies widely in correspondence to different selections of clustering features. By properly extracting model features, the clustered subset could adequately represent the uncertainty of the full ensemble.     

Identifying influence areas with connectivity analysis – application to the local perturbation of heterogeneity distribution for history matching

Numerical representations of a target reservoir can help to assess the potential of different development plans. To be as predictive as possible, these representations or models must reproduce the data (static, dynamic) collected on the field. However, constraining reservoir models to dynamic data – the history-matching process – can be very time consuming. Many uncertain parameters need to be taken into account, such as the spatial distribution of petrophysical properties. This distribution is mostly unknown and usually represented by millions of values populating the reservoir grid. Dedicated parameterization techniques make it possible to investigate many spatial distributions from a small number of parameters. The efficiency of the matching process can be improved from the perturbation of specific regions of the reservoir. Distinct approaches can be considered to define such regions. For instance, one can refer to streamlines. The leading idea is to identify areas that influence the production behavior where the data are poorly reproduced. Here, we propose alternative methods based on connectivity analysis to easily provide approximate influence areas for any fluid-flow simulation. The reservoir is viewed as a set of nodes connected by weighted links that characterize the distance between two nodes. The path between nodes (or grid blocks) with the lowest cumulative weight yields an approximate flow path used to define influence areas. The potential of the approach is demonstrated on the basis of 2D synthetic cases for the joint integration of production and 4D saturation data, considering several formulations for the weights attributed to the links.     

Hydrogeology of a montane headwater groundwater system downgradient of a coal-mine waste rock dump: Elk Valley,British Columbia,Canada

Steelmaking-coal waste rock placed in mountain catchments in the Elk Valley, British Columbia, Canada, drain constituents of interest (CIs) to surface water downgradient of the waste rock dumps. The role of groundwater in transporting CIs in the headwaters of mountain catchments is not well understood. This study characterizes the physical hydrogeology of a portion of a 10-km² headwater catchment (West Line Creek) downgradient of a 2.7-km² waste rock dump placed over a natural headwater valley-bottom groundwater system. The study site was instrumented with 13 monitoring wells. Drill core samples were collected to determine subsurface lithology and geotechnical properties. The groundwater system was characterized using field testing and water-level monitoring. The valley-bottom sediments were composed of unconsolidated glacial and meltwater successions (<64 m thick) deposited as a series of cut and fill structures overlying shale bedrock. An unconfined basal alluvial aquifer located above fractured bedrock was identified as the primary conduit for groundwater flow toward Line Creek (650 m from the toe of the dump). Discharge through the basal alluvial aquifer was estimated using the geometric mean hydraulic conductivity (±1 standard deviation). These calculations suggest groundwater discharge could account for approximately 15% (ranging from 2 to 60%) of the total water discharged from the watershed. The residence time from the base of the waste rock dump to Line Creek was estimated at <3 years. The groundwater system was defined as a snowmelt (i.e., nival) regime dominated by direct recharge (percolation of precipitation) across the catchment.     

Evaluating Tidal Wetland Restoration Performance Using National Estuarine Research Reserve System Reference Sites and the Restoration Performance Index (RPI)

Evaluations of tidal wetland restoration efforts suffer from a lack of appropriate reference sites and standardized methods among projects. To help address these issues, the National Estuarine Research Reserve System (NERRS) and the NOAA Restoration Center engaged in a partnership to monitor ecological responses and evaluate 17 tidal wetland restoration projects associated with five reserves. The goals of this study were to (1) determine the level of restoration achieved at each project using the restoration performance index (RPI), which compares change in parameters over time between reference and restoration sites, (2) compare hydrologic and excavation restoration projects using the RPI, (3) identify key indicator parameters for assessing restoration effectiveness, and (4) evaluate the value of the NERRS as reference sites for local restoration projects. We found that the RPI, modified for this study, was an effective tool for evaluating relative differences in restoration performance; most projects achieved an intermediate level of restoration from 2008 to 2010, and two sites became very similar to their paired reference sites, indicating that the restoration efforts were highly effective. There were no differences in RPI scores between hydrologic and excavation restoration project types. Two abiotic parameters (marsh platform elevation and groundwater level) were significantly correlated with vegetation community structure and thus can potentially influence restoration performance. Our results highlight the value of the NERRS as reference sites for assessing tidal wetland restoration projects and provide improved guidance for scientists and restoration practitioners by highlighting the RPI as a trajectory analysis tool and identifying key monitoring parameters.     

Tidal Habitats Support Large Numbers of Invasive Blue Catfish in a Chesapeake Bay Subestuary

The introduction of a non-native freshwater fish, blue catfish Ictalurus furcatus, in tributaries of Chesapeake Bay resulted in the establishment of fisheries and in the expansion of the population into brackish habitats. Blue catfish are an invasive species in the Chesapeake Bay region, and efforts are underway to limit their impacts on native communities. Key characteristics of the population (population size, survival rates) are unknown, but such knowledge is useful in understanding the impact of blue catfish in estuarine systems. We estimated population size and survival rates of blue catfish in tidal habitats of the James River subestuary. We tagged 34,252 blue catfish during July–August 2012 and 2013; information from live recaptures (n = 1177) and dead recoveries (n = 279) were used to estimate annual survival rates and population size using Barker’s Model in a Robust Design and allowing for heterogeneity in detection probabilities. The blue catfish population in the 12-km study area was estimated to be 1.6 million fish in 2013 (95% confidence interval [CI] adjusted for overdispersion: 926,307–2,914,208 fish). Annual apparent survival rate estimates were low: 0.16 (95% CI 0.10–0.24) in 2012–2013 and 0.44 (95% CI 0.31–0.58) in 2013–2014 and represent losses from the population through mortality, permanent emigration, or both. The tagged fish included individuals that were large enough to exhibit piscivory and represented size classes that are likely to colonize estuarine habitats. The large population size that we estimated was unexpected for a freshwater fish in tidal habitats and highlights the need to effectively manage such species.     

Do Macroalgal Mats Impact Microphytobenthos on Mudflats? Evidence from a Meta-Analysis,Comparative Survey,and Large-Scale Manipulation

Eutrophication and species introductions have resulted in increased macroalgal biomass in coastal ecosystems around the globe. Macroalgal mats may compete with microphytobenthos (MPB) for light and nutrients and, due to their position in the canopy, have a negative impact on MPB biomass. We tested this effect by conducting a meta-analysis of prior experiments, as well as a comparative survey, and a macroalgal-removal manipulation in the coastal lagoons of the Virginia Coastal Reserve (VCR) on the eastern shore of Virginia (USA). In all cases, MPB biomass was estimated using benthic chlorophyll as a proxy. While prior individual studies documented impacts of macroalgae, when effect sizes were averaged across studies, there was no consistent effect of macroalgal biomass on MPB biomass. In the VCR, a non-native red macroalga, Gracilaria vermiculophylla, dominates intertidal mats and attains high biomasses at some sites. Nevertheless, MPB biomass was unrelated to macroalgal mass based on a survey of mudflats. Further, when macroalgae were experimentally manipulated on a mudflat using a before and after impact design, there was no change in MPB. Based on the meta-analysis, survey, and manipulation we conducted, macroalgal mats do not have a generalizable effect on MPB, interactions seem context-dependent, and in the VCR, the effects on MPB appear neutral. This finding is important given the significance of MPB in supporting food webs and other estuarine ecosystem functions, as well as the increasing frequency and intensity of macroalgal blooms.     

Great South Bay After Sandy: Changes in Circulation and Flushing due to New Inlet

The coastal ocean model FVCOM is applied to quantify the changes in circulation, flushing, and exposure time in Great South Bay, New York, after Superstorm Sandy breached the barrier island in 2012. Since then, the lagoon system is connected to the Atlantic via five instead of four inlets. The model simulations are run on two high-resolution unstructured grids, one for the pre-breach configuration, one including the new inlet, with tidal-only forcing, and summer and winter forcing conditions. Despite its small cross-sectional size, the breach has a relatively large net inflow that leads to a strengthening of the along-bay through-flow in Great South Bay (GSB); the tidally driven volume transport in central GSB quadrupled. The seasonal forcing scenarios show that the southwesterly sea breeze in summer slows down the tidally driven flow, while the forcing conditions in winter are highly variable, and the circulation is dependent on wind direction and offshore sea level. Changes in flushing and exposure time associated with the modified transport patterns are evaluated using a Eulerian passive tracer technique. Results show that the new inlet produced a significant decrease in flushing time (approximately 35% reduction under summer wind conditions and 20% reduction under winter wind conditions). Maps of exposure time reflect the local changes in circulation and flushing.     

Weathering characterization for landslides modeling in granitoid rock masses of the Capo Vaticano promontory (Calabria,Italy)

This work focuses on developing multidisciplinary researches concerning weathering profiles related to landscape evolution of the Capo Vaticano promontory on the Calabria Tyrrhenian side (southern Italy). In this area, the tectonic uplift, occurred at least since Pleistocene, together with the Mediterranean climatic conditions, is the main cause of deep weathering and denudation processes. The latter occurred on the outcropping rocks of the crystalline-metamorphic basement, made up of weathered granitoids, in turn belonging to the Monte Poro granitoid complex (intermediate to felsic plutonic rocks covered by Cenozoic sedimentary successions). Field observations coupled to borehole explorations, geophysical surveys, and minero-petrographical analyses allowed the characterization of the granitoid outcrops typical of the studied area in terms of kind and degree of slope instability. This characterization was based on suitable correlations verified between several factors as weathering degree, elastic properties of rocks, and discontinuity features. Weathering profiles are mainly composed by rock masses varying from completely weathered rock with corestones of highly weathered rock (classes IV–V) to slightly weathered rocks (class II). The weathered rocks are involved in several landslide typologies such as debris flow (frequency 48.5%), translational slide (frequency 33.3%), and minor rock fall and rotational slide (frequency 9%). The achieved data allowed the establishment of a general correlation between weathering degree and type of slope instability. Debris flow-type instabilities are predominant on the steeper slopes, involving very poor rock masses ascribed to the shallowest portions of the weathering class IV. Translational slides are less widespread than the previous ones and often involve a mixture of soil and highly weathered rocks. Rotational slides are more frequently close to the top of the slopes, where the thicknesses of more weathered rocks increase, and involve mainly rock masses belonging to the weathering classes IV and V. Rock falls mostly occur on the vertical escarpments of the road cuts and are controlled by the characteristics of the main discontinuities. The assessment of rock mass rating and slope mass rating, based on the application of the discontinuity data, allowed respectively an evaluation of the quality of rock masses and of the susceptibility of rock slopes to failure. The comparison between the last one and the real stability conditions along the cut slopes shows a good correspondence. Finally, the geological strength index system was also applied for the estimation of rock mass properties. The achieved results give a worthy support for a better understanding of the relationship between the distribution of landslides and the geological features related to different weathering degrees. Therefore, they can provide a reliable tool to evaluate the potential stability conditions of the rock slopes in the studied area and a general reference framework for the study of weathering processes in other regions with similar geological features.