Optimization problems can be found in many areas of scienceand technology. Not only the global optimum, but also a(large) number of near-optima are often of interest. Thisgives rise to what are referred to as multimodal optimizationproblems. In most cases, the number and quality of the optimaare unknown and assumptions cannot be made about theobjective functions. In this paper, we focus on continuous,unconstrained optimization in moderately high-dimensionalcontinuous spaces (d ≤ 10).We present a scalablealgorithm with virtually no parameters, which performs wellfor general objective functions (non-convex, discontinuous).It is based on two well-established algorithms (CMA-ES,deterministic crowding). Novel elements of the algorithminclude the detection of seed points for local searches andcollision avoidance, both based on nearest neighbors, and astrategy for semi-sequential optimization to realize scalability.The performance of the proposed algorithm is numericallyevaluated using the CEC2013 niching benchmark suite for1 − 20 dimensional functions, and a 9 dimensional real-worldproblem from constraint optimization in climate research.The algorithm performs well in relation to the CEC2013benchmarks and only falls short on higher dimensional andstrongly inisotropic problems. In the case of the climate-related problem, the algorithm is able to find a high number(> 150) of optima of relevance to climate research. Theproposed algorithm does not require special configuration forthe optimization problems considered in this paper, i.e., itshows good black-box behavior.