Array programming provides a powerful, compact and expressive syntax for accessing, manipulating and operating on data in vectors, matrices and higher-dimensional arrays. NumPy is the primary array programming library for the Python language. It has an essential role in research analysis pipelines in fields as diverse as physics, chemistry, astronomy, geoscience, biology, psychology, materials science, engineering, finance and economics. For example, in astronomy, NumPy was an important part of the software stack used in the discovery of gravitational waves1 and in the first imaging of a black hole2. Here we review how a few fundamental array concepts lead to a simple and powerful programming paradigm for organizing, exploring and analysing scientific data. NumPy is the foundation upon which the scientific Python ecosystem is constructed. It is so pervasive that several projects, targeting audiences with specialized needs, have developed their own NumPy-like interfaces and array objects. Owing to its central position in the ecosystem, NumPy increasingly acts as an interoperability layer between such array computation libraries and, together with its application programming interface (API), provides a flexible framework to support the next decade of scientific and industrial analysis.

Abstract:Rapidly depleting unconfined aquifers are the primary source of water for irrigation on the North China Plain. Yet, despite its critical importance, groundwater recharge to the Plain remains an enigma. We introduce a one-dimensional soil-water-balance model to estimate precipitation-and irrigation-generated areal recharge from commonly available crop and soil characteristics and climate data. To limit input data needs and to simplify calculations, the model assumes that water flows vertically downward under a unit gradient; infiltration and evapotranspiration are separate, sequential processes; evapotranspiration is allocated to evaporation and transpiration as a function of leaf-area index and is limited by soil-moisture content; and evaporation and transpiration are distributed through the soil profile as exponential functions of soil and root depth, respectively. For calibration, model-calculated water contents of 11 soildepth intervals from 0 to 200 cm were compared with measured water contents of loam soil at four sites in Luancheng County, Hebei Province, over 3 years (1998)(1999)(2000)(2001). Each 50-m 2 site was identically cropped with winter wheat and summer maize, but received a different irrigation treatment. Average root mean-squared error between measured and model-calculated water content of the top 180 cm was 4Ð2 cm, or 9Ð3% of average total water content. In addition, model-calculated evapotranspiration compared well with that measured by a large-scale lysimeter. To test the model, 12 additional sites were simulated successfully. Model results demonstrate that drainage from the soil profile is not a constant fraction of precipitation and irrigation inputs, but rather the fraction increases as the inputs increase. Because this drainage recharges the underlying aquifer, improving irrigation efficiency by reducing seepage will not reverse water-table declines.

No abstract

Abstract:Because the traditional Soil Conservation Service curve-number (SCS-CN) approach continues to be used ubiquitously in water quality models, new application methods are needed that are consistent with variable source area (VSA) hydrological processes in the landscape. We developed and tested a distributed approach for applying the traditional SCS-CN equation to watersheds where VSA hydrology is a dominant process. Predicting the location of source areas is important for watershed planning because restricting potentially polluting activities from runoff source areas is fundamental to controlling non-point-source pollution. The method presented here used the traditional SCS-CN approach to predict runoff volume and spatial extent of saturated areas and a topographic index, like that used in TOPMODEL, to distribute runoff source areas through watersheds. The resulting distributed CN-VSA method was applied to two subwatersheds of the Delaware basin in the Catskill Mountains region of New York State and one watershed in south-eastern Australia to produce runoff-probability maps. Observed saturated area locations in the watersheds agreed with the distributed CN-VSA method. Results showed good agreement with those obtained from the previously validated soil moisture routing (SMR) model. When compared with the traditional SCS-CN method, the distributed CN-VSA method predicted a similar total volume of runoff, but vastly different locations of runoff generation. Thus, the distributed CN-VSA approach provides a physically based method that is simple enough to be incorporated into water quality models, and other tools that currently use the traditional SCS-CN method, while still adhering to the principles of VSA hydrology.

[1] Effective control of nonpoint source contaminants in runoff from urbanized watersheds requires knowledge about the locations of runoff source areas under a variety of conditions. Physical monitoring of spatially variant processes, such as runoff production form variable source areas, is time-consuming and expensive. Thus modeling the processes may provide a valuable cost-effective alternative. In this paper we adapt and validate a variable source model for an urban watershed to predict areas of the landscape prone to elevated soil moisture levels and saturation excess runoff. We modified the soil moisture distribution and routing (SMDR) model to simulate hydrologic processes in an urban upstate New York watershed by considering the impact of impervious surfaces, hydraulic control structures (detention basins), and land use on the water balance. In our model, infiltration excess runoff from impervious surfaces infiltrates in the surrounding soils. Simulated and observed streamflow agreed well, and more importantly, the distribution of soil moisture levels and overland flow throughout the watershed were well predicted. Removing urban features from the model resulted in substantially lower peak stormflows than observed, and the base and interflows increased accordingly. Both modeled and measured distributed results indicated that the more urbanized areas of the study site had both higher soil moistures and runoff losses due to shallower soils, greater upslope contributing areas, and a larger area of impervious surfaces generating runoff. The model results helped identify processes that describe how urbanization impacts integrated and distributed hydrology, which provides useful information for targeted water quality management.